Changeset 129 in lmdz_wrf

Timestamp:

Jul 28, 2014, 6:04:06 PM (10 years ago)

Author:

lfita

Message:

Adding particular checking

File:

• branches/LMDZ_WRFmeas/WRFV3/lmdz/thermcell_main.F90 (modified) (1 diff)

branches/LMDZ_WRFmeas/WRFV3/lmdz/thermcell_main.F90

@@ -916 +916 @@
 !          if ((pplay(ig,l) .ge. pcon(ig)) .and. (pplay(ig,l+1) .le. pcon(ig))) then
 ! L. Fita, LMD July 2014. Adding avoiding divisons by zero...
-        if ((pplay(ig,l) .ge. pcon(ig)) .and. (pplay(ig,l+1) .le. pcon(ig)) .and.    &
-          (klcl(ig) .gt. 0) .and. (klcl(ig) + 1 .le. nlay-1) .and.                   &
-          (pplay(ig,klcl(ig)+1)-pplay(ig,klcl(ig)) .ne. 0.) ) then
-          klcl(ig)=l
-          interp(ig)=(pcon(ig)-pplay(ig,klcl(ig)))/(pplay(ig,klcl(ig)+1)-pplay(ig,klcl(ig)))
+          if ((pplay(ig,l).ge.pcon(ig)) .and. (pplay(ig,l+1).le.pcon(ig)) .and.      &
+            (klcl(ig).gt.0).and.(klcl(ig)+1.le.nlay-1) .and. (klcl(ig)+1.gt.0) then
+!            (pplay(ig,klcl(ig)+1)-pplay(ig,klcl(ig)).ne.0.) ) then
+            klcl(ig)=l
+            interp(ig)=(pcon(ig)-pplay(ig,klcl(ig)))/(pplay(ig,klcl(ig)+1)-pplay(ig,klcl(ig)))
+            IF interp(ig) /= interp(ig) THEN
+              PRINT *,'  Lluis wrong interp= ',interp(ig),' at ', ig,' klcl(ig): ',  &
+                klcl(ig),' klcl(ig)+1: ', klcl(ig)+1
+            END IF
           endif
         endif