Changeset 104 in lmdz_wrf for branches/LMDZ_WRFmeas

Timestamp:

Jul 25, 2014, 4:23:01 PM (11 years ago)

Author:

lfita

Message:

Adding artesanal chck for NaNs?

File:

• branches/LMDZ_WRFmeas/WRFV3/lmdz/thermcell_flux2.F90 (modified) (2 diffs)

branches/LMDZ_WRFmeas/WRFV3/lmdz/thermcell_flux2.F90

@@ -248 +248 @@
              zzz=fm(ig,l)*rhobarz(ig,l+1)*zw2(ig,l+1)  &
      &                          /(rhobarz(ig,l)*zw2(ig,l))
+             IF (zzz /= zzz .OR. zw2(ig,l) == 0. ) 
+               PRINT *,'  Lluis ', ig,',',l,' wrong zzz: ',zzz,' zw2: ',zw2(ig,l),' fm: ',fm(ig,l)
+             END IF
              if (fm(ig,l+1).gt.zzz) then
                 detr(ig,l)=detr(ig,l)+fm(ig,l+1)-zzz
@@ -261 +264 @@
       lvarname = 'fm'
       CALL check_var(lfname, lvarname, fm(:,l), ngrid, klev+1, largest, .FALSE.)
-      lvarname = '1/zw2'
-      CALL check_var(lfname, lvarname, 1./zw2(:,l), ngrid, klev+1, largest, .FALSE.)
-      lvarname = '1/rhobarz'
-      CALL check_var(lfname, lvarname, 1./rhobarz(:,l), ngrid, klev+1, largest, .FALSE.)
+!      lvarname = '1/zw2'
+!      CALL check_var(lfname, lvarname, 1./zw2(:,l), ngrid, klev+1, largest, .FALSE.)
+!      lvarname = '1/rhobarz'
+!      CALL check_var(lfname, lvarname, 1./rhobarz(:,l), ngrid, klev+1, largest, .FALSE.)
       lvarname = 'entr'
       CALL check_var(lfname, lvarname, entr(:,l), ngrid, klev, largest, .FALSE.)