Changeset 116 in lmdz_wrf

Timestamp:

Jul 28, 2014, 10:51:37 AM (10 years ago)

Author:

lfita

Message:

Adding new checks. NaN values on 'w_con'

Location:

branches/LMDZ_WRFmeas/WRFV3/lmdz

Files:

• thermcell_flux2.F90 (modified) (1 diff)
• thermcell_main.F90 (modified) (5 diffs)

branches/LMDZ_WRFmeas/WRFV3/lmdz/thermcell_flux2.F90

@@ -481 +481 @@
       CALL check_var3D(lfname, lvarname, fm, ngrid, klev+1, largest, .FALSE.)
       lvarname = 'rhobarz'
-      CALL check_var3D(lfname, lvarname, rhobarz, ngrid, klev+1, largest, .FALSE.)
+      CALL check_var3D(lfname, lvarname, rhobarz, ngrid, klev, largest, .FALSE.)
       lvarname = 'zw2'
       CALL check_var3D(lfname, lvarname, zw2, ngrid, klev+1, largest, .FALSE.)

branches/LMDZ_WRFmeas/WRFV3/lmdz/thermcell_main.F90

@@ -429 +429 @@
       lvarname = 'zw2'
       CALL check_var3D(lfname, lvarname, zw2, ngrid, nlay+1, largest, .FALSE.)
-      lvarname = 'alim_star'
-      CALL check_var3D(lfname, lvarname, alim_star, ngrid, nlay, largest, .FALSE.)
 
       if (prt_level.ge.1) print*,'thermcell_main apres initialisation'
@@ -902 +900 @@
       enddo
     enddo
+      lfname='thermcell_main calculation of LCL'
+      lvarname = 'inerp'
+      CALL check_var3D(lfname, lvarname, interp, ngrid, nlay, largest, .FALSE.)
+      lvarname = 'klcl'
+      CALL check_var(lfname, lvarname, REAL(klcl), ngrid, largest, .FALSE.)
+      lvarname = 'pcon'
+      CALL check_var(lfname, lvarname, pcon, ngrid, largest, .FALSE.)
+      lvarname = 'pplay'
+      CALL check_var3D(lfname, lvarname, pplay, ngrid, nlay, largest, .FALSE.)
 
 !------------Hauteur des thermiques
@@ -1058 +1065 @@
       lvarname = 'pdtadj'
       CALL check_var3D(lfname, lvarname, pdtadj, ngrid, nlay, largest, .FALSE.)
-      lvarname = 'alp_bl_det'
-      CALL check_var(lfname, lvarname, alp_bl_det,ngrid, largest, .FALSE.)
       lvarname = 'rhobarz0'
       CALL check_var(lfname, lvarname, rhobarz0,ngrid, largest, .FALSE.)
@@ -1066 +1071 @@
       lvarname = 'w_conv'
       CALL check_var(lfname, lvarname, w_conv,ngrid, largest, .FALSE.)
+      lvarname = 'interp'
+      CALL check_var(lfname, lvarname, interp,ngrid, largest, .FALSE.)
       lvarname = 'w0'
       CALL check_var(lfname, lvarname, w0,ngrid, largest, .FALSE.)
-      lvarname = 'alp_bl_fluct_m'
-      CALL check_var(lfname, lvarname, alp_bl_fluct_m,ngrid, largest, .FALSE.)
-      lvarname = 'alp_bl_fluct_tke'
-      CALL check_var(lfname, lvarname, alp_bl_fluct_tke,ngrid, largest, .FALSE.)
       lvarname = 'therm_tke_max0'
       CALL check_var(lfname, lvarname, therm_tke_max0,ngrid, largest, .FALSE.)
@@ -1078 +1081 @@
       lvarname = 'pbl_tke_max0'
       CALL check_var(lfname, lvarname, pbl_tke_max0,ngrid, largest, .FALSE.)
-      lvarname = 'alp_bl_conv'
-      CALL check_var(lfname, lvarname, alp_bl_conv,ngrid, largest, .FALSE.)
-      lvarname = 'alp_bl_stat'
-      CALL check_var(lfname, lvarname, alp_bl_stat,ngrid, largest, .FALSE.)
 !------------Closure------------------