fix ibas_prf=3 initialization when zf(k) is exactly 11000 (1st value …

…in the lapse rate table) - fixes #41

src/modstartup.f90

@@ -1192,7 +1192,7 @@ subroutine baseprofs
             tb(k)=tsurf+lapserate(1)*(zf(k)-zsurf)
           else
             j=1
-            do while(zf(k)>zmat(j))
+            do while(zf(k)>=zmat(j))
               j=j+1
             end do
             tb(k)=tmat(j-1)+lapserate(j)*(zf(k)-zmat(j-1))