Changeset 2239 for trunk/LMDZ.TITAN/libf/phytitan/bilinearbig.F90

Timestamp:

Feb 6, 2020, 4:49:54 PM (5 years ago)

Author:

jvatant

Message:

Minor correction in bilinearbig, fixing a bug
when reaching exactly the T range limit
causing nasty overflow error.
--JVO

File:

• trunk/LMDZ.TITAN/libf/phytitan/bilinearbig.F90 (modified) (3 diffs)

trunk/LMDZ.TITAN/libf/phytitan/bilinearbig.F90

@@ -55 +55 @@
 
 !     ... and for y within the temperature range
-      if ((y.lt.y_arr(1)).or.(y.gt.y_arr(nY))) then
-         print*,y_arr(1),y_arr(nY)
+      if ((y.le.y_arr(1)).or.(y.ge.y_arr(nY))) then
+         !print*,y_arr(1),y_arr(nY)
          !write(*,*) 'Warning from bilinearbig routine:'
          !write(*,*) 'Outside continuum temperature range!'
-         if(y.lt.y_arr(1))then
+         if(y.le.y_arr(1))then
             y=y_arr(1)+0.01
             b=1
@@ -65 +65 @@
             y2=y_arr(b+1)
          endif
-         if(y.gt.y_arr(nY))then
+         if(y.ge.y_arr(nY))then
             y=y_arr(nY)-0.01
             b=nY-1
@@ -74 +74 @@
         j=1
         y2=y_arr(j)
-        do while ( y2 .le. y )
+        do while ( y2 .lt. y )
           y1=y2
           j=j+1