Changeset 127 in lmdz_wrf

Timestamp:

Jul 28, 2014, 5:42:58 PM (10 years ago)

Author:

lfita

Message:

Adding condition to avoid pintep == NaN

File:

• branches/LMDZ_WRFmeas/WRFV3/lmdz/thermcell_main.F90 (modified) (3 diffs)

branches/LMDZ_WRFmeas/WRFV3/lmdz/thermcell_main.F90

@@ -914 +914 @@
       do ig=1,ngrid
         if (ok_lcl(ig)) then 
-          if ((pplay(ig,l) .ge. pcon(ig)) .and. (pplay(ig,l+1) .le. pcon(ig))) then
+!          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
           klcl(ig)=l
           interp(ig)=(pcon(ig)-pplay(ig,klcl(ig)))/(pplay(ig,klcl(ig)+1)-pplay(ig,klcl(ig)))
@@ -922 +925 @@
     enddo
       lfname='thermcell_main calculation of LCL'
-      lvarname = 'inerp'
+      lvarname = 'interp'
       CALL check_var3D(lfname, lvarname, interp, ngrid, nlay, largest, .FALSE.)
       lvarname = 'klcl'
@@ -930 +933 @@
       lvarname = 'pplay'
       CALL check_var3D(lfname, lvarname, pplay, ngrid, nlay, largest, .FALSE.)
+      lvarname = '1/pplay'
+      CALL check_var3D(lfname, lvarname, 1./pplay, ngrid, nlay, largest, .FALSE.)
+
 
 !------------Hauteur des thermiques