Changeset 110 in lmdz_wrf

Timestamp:

Jul 25, 2014, 5:46:02 PM (10 years ago)

Author:

lfita

Message:

Searching for sources of NaN alim_star

File:

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

branches/LMDZ_WRFmeas/WRFV3/lmdz/thermcell_plume.F90

@@ -174 +174 @@
                alim_star(ig,l)=alim_star(ig,l)/alim_star_tot(ig)
             endif
+            IF (zlev(ig,l+1) < 0. .OR. alim_star_tot(ig) == 0.) THEN
+              PRINT *,'Lluis ',ig,',',l,' negative zlev: ',zlev(ig,l+1),' null: ',alim_start_tot(ig)
+            END IF
          enddo
       enddo
       alim_star_tot(:)=1.
 
+      lfname='thermcell_plume after computing alim_star'
+      lvarname = 'alim_star'
+      CALL check_var3D(lfname, lvarname, alim_star, ngrid, klev, largest, .FALSE.)
+      lvarname = 'ztv'
+      CALL check_var3D(lfname, lvarname, ztv, ngrid, klev, largest, .FALSE.)
+      lvarname = 'zlev'
+      CALL check_var3D(lfname, lvarname, zlev, ngrid, klev+1, largest, .FALSE.)
 
 !------------------------------------------------------------------------------