Changeset 184 in lmdz_wrf for trunk

Timestamp:

Oct 15, 2014, 3:54:12 PM (10 years ago)

Author:

lfita

Message:

Adding 'qc' on mp==0 in Registry/Registry?.EM.LMDZ. With this now we have also

'QCLOUD' even with mp=0. This is the way to keep liquid water for the next
time-step and have large scale precipitation!

Adding some checks around moisture values and tendencies...

Location:

trunk/WRFV3

Files:

• Registry/Registry.EM.LMDZ (modified) (2 diffs)
• Registry/Registry.backup (modified) (1 diff)
• dyn_em/solve_em.F (modified) (3 diffs)
• phys/module_lmdz_phys.F (modified) (1 diff)

trunk/WRFV3/Registry/Registry.EM.LMDZ

@@ -86 +86 @@
 state    real   v_gc           igj      dyn_em      1        YZ    i1  "VV"     "y-wind component"    "m s-1"
 state    real   t_gc           igj      dyn_em      1        Z     i1  "TT"     "temperature"         "K"
-state    real   rh_gc          igj      dyn_em      1        Z     i1  "RH"    "relative humidity"   "%"
+state    real   rh_gc          igj      dyn_em      1        Z     i01 "RH"    "relative humidity"   "%"
 state    real   ght_gc         igj      dyn_em      1        Z     i1  "GHT"   "geopotential height" "m"
 state    real   p_gc           igj      dyn_em      1        Z     i1  "PRES"   "pressure"            "Pa"
@@ -1961 +1961 @@
 
 #package   passivec1     chem_opt==0                  -             
-package   passiveqv     mp_physics==0                -             moist:qv
+# L. Fita, LMD October 2014. We require also qc in order to properly have liquid
+#   water which is required in LMDZ phsics
+#package   passiveqv     mp_physics==0                -             moist:qv
+package   passiveqv     mp_physics==0                -             moist:qv,qc
 package   kesslerscheme mp_physics==1                -             moist:qv,qc,qr
 package   linscheme     mp_physics==2                -             moist:qv,qc,qr,qi,qs,qg

trunk/WRFV3/Registry/Registry.backup

@@ -1962 +1962 @@
 
 #package   passivec1     chem_opt==0                  -             
-package   passiveqv     mp_physics==0                -             moist:qv
+# L. Fita, LMD October 2014. We require also qc in order to properly have liquid
+#   water which is required in LMDZ phsics
+#package   passiveqv     mp_physics==0                -             moist:qv
+package   passiveqv     mp_physics==0                -             moist:qv,qc
 package   kesslerscheme mp_physics==1                -             moist:qv,qc,qr
 package   linscheme     mp_physics==2                -             moist:qv,qc,qr,qi,qs,qg

trunk/WRFV3/dyn_em/solve_em.F

@@ -2494 +2494 @@
 !</DESCRIPTION>
 
+       PRINT *,'  Lluis PARAM_FIRST_SCALAR: ', PARAM_FIRST_SCALAR,    &
+         ' p_qv: ',p_qv, ' p_qc: ',p_qc,' U(moist): ', UBOUND(moist)
+       PRINT *,'   moist: ',moist(config_flags%i_check_point, 2,   &
+               config_flags%j_check_point,:)
 
        moist_scalar_advance: IF (num_3d_m >= PARAM_FIRST_SCALAR )  THEN
@@ -2510 +2514 @@
                CALL wrf_debug ( 200 , ' call rk_scalar_tend' )
                tenddec = .false.
+           PRINT *,'  Lluis before rk_update_scalar moist im=',im,         &
+             ' moist_old moist_tend moist________'
+           DO iz = kms, kme
+             PRINT *,moist_old(config_flags%i_check_point, iz,             &
+               config_flags%j_check_point,im),                             &
+               moist_tend(config_flags%i_check_point, iz,                  &
+               config_flags%j_check_point,im),                             &
+               moist(config_flags%i_check_point, iz,                       &
+               config_flags%j_check_point,im)
+           END DO
 
 BENCH_START(rk_scalar_tend_tim)
@@ -2601 +2615 @@
                                kts=k_start    , kte=k_end                              )
 BENCH_END(update_scal_tim)
+           PRINT *,'  Lluis after rk_update_scalar moist im=',im,         &
+             ' moist_old moist_tend moist________'
+           DO iz = kms, kme
+             PRINT *,moist_old(config_flags%i_check_point, iz,             &
+               config_flags%j_check_point,im),                             &
+               moist_tend(config_flags%i_check_point, iz,                  &
+               config_flags%j_check_point,im),                             &
+               moist(config_flags%i_check_point, iz,                       &
+               config_flags%j_check_point,im)
+           END DO
 
 BENCH_START(flow_depbdy_tim)

trunk/WRFV3/phys/module_lmdz_phys.F

@@ -377 +377 @@
       ' wrf_xtime: ',wrf_xtime,' wrf_rst: ', wrf_isrestart,' wrf_lwbdy: ',           &
       wrf_islowbdyin
+    PRINT *,'  Number of water species from WRF: ',wnum3dm
 
     phykeys=0.