Changeset 184 in lmdz_wrf for trunk/WRFV3/dyn_em

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...

File:

• trunk/WRFV3/dyn_em/solve_em.F (modified) (3 diffs)

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)