Changeset 128 in lmdz_wrf

Timestamp:

Jul 28, 2014, 5:51:20 PM (10 years ago)

Author:

lfita

Message:

Adding more conditions to the criteria

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.  &
-            (pplay(ig,klcl(ig)+1)-pplay(ig,klcl(ig))) .ne. 0.) then
+        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)))