Changeset 79 in lmdz_wrf for branches

Timestamp:

Jul 25, 2014, 10:06:36 AM (10 years ago)

Author:

lfita

Message:

Adding checks due to NaNs? in Alp_bl!

File:

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

branches/LMDZ_WRFmeas/WRFV3/lmdz/thermcell_main.F90

@@ -621 +621 @@
       lvarname = 'zgla'
       CALL check_var3D(lfname, lvarname, fm, ngrid, nlay, largest, .FALSE.)
-      lvarname = 'fm'
-      CALL check_var3D(lfname, lvarname, fm, ngrid, nlay, largest, .FALSE.)
 
       CALL thermcell_flux2(ngrid,nlay,ptimestep,masse, &
@@ -1000 +998 @@
       lvarname = 'pdtadj'
       CALL check_var3D(lfname, lvarname, pdtadj, klon, klev, largest, .FALSE.)
-
+      lvarname = 'alp_bl_det'
+      CALL check_var(lfname, lvarname, alp_bl_det,klon, largest, .FALSE.)
+      lvarname = 'rhobarz0'
+      CALL check_var(lfname, lvarname, rhobarz0,klon, largest, .FALSE.)
+      lvarname = 'fraca0'
+      CALL check_var(lfname, lvarname, fraca0,klon, largest, .FALSE.)
+      lvarname = 'w_conv'
+      CALL check_var(lfname, lvarname, w_conv,klon, largest, .FALSE.)
+      lvarname = 'w0'
+      CALL check_var(lfname, lvarname, w0,klon, largest, .FALSE.)
+      lvarname = 'alp_bl_fluct_m'
+      CALL check_var(lfname, lvarname, alp_bl_fluct_m,klon, largest, .FALSE.)
+      lvarname = 'alp_bl_fluct_tke'
+      CALL check_var(lfname, lvarname, alp_bl_fluct_tke,klon, largest, .FALSE.)
+      lvarname = 'therm_tke_max0'
+      CALL check_var(lfname, lvarname, therm_tke_max0,klon, largest, .FALSE.)
+      lvarname = 'env_tke_max0'
+      CALL check_var(lfname, lvarname, env_tke_max0,klon, largest, .FALSE.)
+      lvarname = 'pbl_tke_max0'
+      CALL check_var(lfname, lvarname, pbl_tke_max0,klon, largest, .FALSE.)
+      lvarname = 'alp_bl_conv'
+      CALL check_var(lfname, lvarname, alp_bl_conv,klon, largest, .FALSE.)
+      lvarname = 'alp_bl_stat'
+      CALL check_var(lfname, lvarname, alp_bl_stat,klon, largest, .FALSE.)
 !------------Closure------------------
 
@@ -1027 +1048 @@
 
   ENDIF ! (iflag_clos_bl.ge.1)
+      lfname='after closure'
+      lvarname = 'pt'
+      CALL check_var3D(lfname, lvarname, pt, klon, klev, largest, .FALSE.)
+      lvarname = 'pdtadj'
+      CALL check_var3D(lfname, lvarname, pdtadj, klon, klev, largest, .FALSE.)
+      lvarname = 'alp_bl_det'
+      CALL check_var(lfname, lvarname, alp_bl_det,klon, largest, .FALSE.)
+      lvarname = 'rhobarz0'
+      CALL check_var(lfname, lvarname, rhobarz0,klon, largest, .FALSE.)
+      lvarname = 'fraca0'
+      CALL check_var(lfname, lvarname, fraca0,klon, largest, .FALSE.)
+      lvarname = 'w_conv'
+      CALL check_var(lfname, lvarname, w_conv,klon, largest, .FALSE.)
+      lvarname = 'w0'
+      CALL check_var(lfname, lvarname, w0,klon, largest, .FALSE.)
+      lvarname = 'alp_bl_fluct_m'
+      CALL check_var(lfname, lvarname, alp_bl_fluct_m,klon, largest, .FALSE.)
+      lvarname = 'alp_bl_fluct_tke'
+      CALL check_var(lfname, lvarname, alp_bl_fluct_tke,klon, largest, .FALSE.)
+      lvarname = 'therm_tke_max0'
+      CALL check_var(lfname, lvarname, therm_tke_max0,klon, largest, .FALSE.)
+      lvarname = 'env_tke_max0'
+      CALL check_var(lfname, lvarname, env_tke_max0,klon, largest, .FALSE.)
+      lvarname = 'pbl_tke_max0'
+      CALL check_var(lfname, lvarname, pbl_tke_max0,klon, largest, .FALSE.)
+      lvarname = 'alp_bl_conv'
+      CALL check_var(lfname, lvarname, alp_bl_conv,klon, largest, .FALSE.)
+      lvarname = 'alp_bl_stat'
+      CALL check_var(lfname, lvarname, alp_bl_stat,klon, largest, .FALSE.)
 
 !!! fin nrlmd le 10/04/2012
@@ -1063 +1113 @@
          enddo
       enddo
+      lfname='after calculation of Al[p/e]_bl'
+      lvarname = 'fm'
+      CALL check_var3D(lfname, lvarname, fm, klon, klev, largest, .FALSE.)
+      lvarname = 'rhobarz'
+      CALL check_var3D(lfname, lvarname, rhobarz, klon, klev, largest, .FALSE.)
+      lvarname = 'Alp_bl'
+      CALL check_var(lfname, lvarname, Alp_bl,klon, largest, .FALSE.)
+      lvarname = 'Ale_bl'
+      CALL check_var(lfname, lvarname, Ale_bl,klon, largest, .FALSE.)
 
 ! assez bizarre car, si on est dans la couche d'alim et que alim_star et