Changeset 3899 for LMDZ5

Timestamp:

May 14, 2021, 7:12:13 PM (4 years ago)

Author:

lguez

Message:

Old work on cold pool death.

Location:

LMDZ5/branches/Cold_pool_death

Files:

• bld.cfg (modified) (2 diffs)
• libf/phylmd/Cold_pool_death (added)
• libf/phylmd/Cold_pool_death/cold_pool_death_condition_m.F90 (added)
• libf/phylmd/Cold_pool_death/cold_pool_death_m.F90 (added)
• libf/phylmd/Cold_pool_death/derivs_m.F90 (added)
• libf/phylmd/Cold_pool_death/derivs_param.F90 (added)
• libf/phylmd/calwake.F90 (modified) (3 diffs)
• libf/phylmd/phys_output_ctrlout_mod.F90 (modified) (1 diff)
• libf/phylmd/phys_output_var_mod.F90 (modified) (3 diffs)
• libf/phylmd/phys_output_write_mod.F90 (modified) (3 diffs)
• libf/phylmd/physiq_mod.F90 (modified) (1 diff)
• libf/phylmd/wake.F90 (modified) (6 diffs)

LMDZ5/branches/Cold_pool_death/bld.cfg

@@ -37 +37 @@
 src::cosp    %COSP
 src::ext_src %EXT_SRC
+src::cold_pool_death %SRC_PATH/%PHYS/Cold_pool_death
 
 bld::lib            lmdz
@@ -98 +99 @@
 bld::excl_dep        use::ifile_attr
 bld::excl_dep        use::ixml_tree
+bld::excl_dep        use::numer_rec_95
 
 # Don't generate interface files

LMDZ5/branches/Cold_pool_death/libf/phylmd/calwake.F90

@@ -14 +14 @@
     wake_omg, wake_dp_deltomg, &
     wake_spread, wake_cstar, wake_d_deltat_gw, &
-    wake_ddeltat, wake_ddeltaq, wake_ds, wake_ddens)
+    wake_ddeltat, wake_ddeltaq, wake_ds, wake_ddens, itap)
   ! **************************************************************
   ! *
@@ -68 +68 @@
   REAL, DIMENSION(klon, klev),   INTENT (OUT)        :: wake_ddeltat, wake_ddeltaq
   REAL, DIMENSION(klon),         INTENT (OUT)        :: wake_ds, wake_ddens
+  integer, intent(in):: itap ! index of call to physiq
 
 
@@ -211 +212 @@
     dtke, dqke, omg, dp_deltomg, spread, cstar, &
     d_deltat_gw, &
-    d_deltatw, d_deltaqw, d_sigmaw, d_wdens)                     ! tendencies
+    d_deltatw, d_deltaqw, d_sigmaw, d_wdens, itap)                     ! tendencies
 
 !

LMDZ5/branches/Cold_pool_death/libf/phylmd/phys_output_ctrlout_mod.F90

@@ -1476 +1476 @@
   TYPE(ctrl_out), SAVE :: o_wake_omg = ctrl_out((/ 4, 5, 10, 10, 10, 10, 11, 11, 11, 11/), &
     'wake_omg', 'wake_omg', ' ', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_wake_lifetime = ctrl_out((/ 4, 5, 10, 10, 10, 10, &
+       11, 11, 11, 11/), 'wake_lifetime', 'lifetime of wake', 's', &
+       (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_wake_final_radius = ctrl_out((/ 4, 5, 10, 10, 10, &
+       10, 11, 11, 11, 11/), 'wake_final_radius', 'final radius of wake', 'm', &
+       (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_wake_final_height = ctrl_out((/ 4, 5, 10, 10, 10, &
+       10, 11, 11, 11, 11/), 'wake_final_height', 'final height of wake', 'm', &
+       (/ ('', i=1, 10) /))
   TYPE(ctrl_out), SAVE :: o_wdtrainA = ctrl_out((/ 4, 5, 10,  4, 10, 10, 11, 11, 11, 11 /), &
     'wdtrainA', 'precipitation from AA', '-', (/ ('', i=1, 10) /))

LMDZ5/branches/Cold_pool_death/libf/phylmd/phys_output_var_mod.F90

@@ -114 +114 @@
  !$OMP THREADPRIVATE(sens_prec_liq_o, sens_prec_sol_o,lat_prec_liq_o,lat_prec_sol_o)
 
+  REAL, allocatable:: wake_lifetime(:) ! (klon) in s
+  !$OMP THREADPRIVATE(wake_lifetime)
+
+  REAL, allocatable:: wake_final_radius(:), wake_final_height(:) ! (klon) in m
+  !$OMP THREADPRIVATE(wake_final_radius, wake_final_height)
+
 CONTAINS
 
@@ -123 +129 @@
 
     include "clesphys.h"
+
+    include "thermcell.h"
+    ! (for iflag_wake)
 
     !------------------------------------------------
@@ -163 +172 @@
     IF (ok_gwd_rando) allocate(zustr_gwd_rando(klon), zvstr_gwd_rando(klon))
 
+    if (iflag_wake>=1) allocate(wake_lifetime(klon), wake_final_radius(klon), &
+         wake_final_height(klon))
+
   END SUBROUTINE phys_output_var_init
 

LMDZ5/branches/Cold_pool_death/libf/phylmd/phys_output_write_mod.F90

@@ -182 +182 @@
          o_map_emis_Anv, o_map_pcld_Anv, o_map_tcld_Anv, &
          o_map_ntot, o_map_hc,o_map_hist,o_map_Cb,o_map_ThCi,o_map_Anv, &
-         o_alt_tropo
+         o_alt_tropo, o_wake_lifetime, o_wake_final_radius, o_wake_final_height
 
 #ifdef CPP_StratAer
@@ -317 +317 @@
          map_emis_ThCi,map_pcld_ThCi,map_tcld_ThCi,&
          map_emis_Anv,map_pcld_Anv,map_tcld_Anv, &
-         alt_tropo
+         alt_tropo, wake_lifetime, wake_final_radius, wake_final_height  
 
     USE ocean_slab_mod, ONLY: nslay, tslab, slab_bils, slab_bilg, tice, &
@@ -1057 +1057 @@
              CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)
              CALL histwrite_phy(o_dqwak2d, zx_tmp_fi2d)
+             CALL histwrite_phy(o_wake_lifetime, wake_lifetime)
+             CALL histwrite_phy(o_wake_final_radius, wake_final_radius)
+             CALL histwrite_phy(o_wake_final_height, wake_final_height)
           ENDIF ! iflag_wake>=1
           CALL histwrite_phy(o_ftd, ftd)

LMDZ5/branches/Cold_pool_death/libf/phylmd/physiq_mod.F90

@@ -2772 +2772 @@
                wake_omg, wake_dp_deltomg,  &
                wake_spread, wake_Cstar, d_deltat_wk_gw,  &
-               d_deltat_wk, d_deltaq_wk, d_s_wk, d_dens_wk)
+               d_deltat_wk, d_deltaq_wk, d_s_wk, d_dens_wk, itap)
           !
           !jyg    Reinitialize itapwk when wakes have been called

LMDZ5/branches/Cold_pool_death/libf/phylmd/wake.F90

@@ -11 +11 @@
                 dtke, dqke, omg, dp_deltomg, spread, cstar, &
                 d_deltat_gw, &
-                d_deltatw2, d_deltaqw2, d_sigmaw2, d_wdens2)                ! tendencies
+                d_deltatw2, d_deltaqw2, d_sigmaw2, d_wdens2, itap)                ! tendencies
 
 
@@ -26 +26 @@
   USE dimphy
   use mod_phys_lmdz_para
+  use cold_pool_death_m, only: cold_pool_death
   USE print_control_mod, ONLY: prt_level
+  use write_field_phy, only: WriteField_phy
+  use phys_local_var_mod, only: zxtsol, q_seri, sens_w
+  use phys_output_var_mod, only: wake_lifetime, wake_final_radius, &
+       wake_final_height
+
   IMPLICIT NONE
   ! ============================================================================
@@ -158 +164 @@
   REAL, DIMENSION (klon, klev),     INTENT(OUT)         :: d_deltatw2, d_deltaqw2
   REAL, DIMENSION (klon),           INTENT(OUT)         :: d_sigmaw2, d_wdens2
+  integer, intent(in):: itap ! index of call to physiq
 
   ! Variables internes
@@ -272 +279 @@
 
   REAL, DIMENSION(klon)                                 :: sigmaw_in   ! pour les prints
+
+  real:: z_cp0  = 1000. ! initial height of cold pool, in m
+  !$omp THREADPRIVATE(z_cp0)
+
+  real:: r0  = 2000. ! initial radius of cold pool, in m
+  !$omp THREADPRIVATE(r0)
+
+  real:: h0  = 1000. ! initial height of environmental boundary layer, in m
+  !$omp THREADPRIVATE(h0)
+
+  real:: eps_cp = 0.06 ! entrainment coefficient at the border of the cold pool
+  !$omp THREADPRIVATE(eps_cp)
+
+  real:: beta = 0.2 ! entrainment at the top of the mixed layer
+  !$omp THREADPRIVATE(beta)
+
+  logical:: av_thvu_constant = .true.
+  !$omp THREADPRIVATE(av_thvu_constant)
+
+  real:: k_star = 0.68
+  !$omp THREADPRIVATE(k_star)
+
+  real rate_uns(klon) ! mass rate in one unsaturated downdraft, in kg s-1
+  real theta_vs(klon) ! virtual potential temperature at the surface, in K
+
+  real theta_v_uns(klon)
+  ! virtual potential temperature of the unsaturated downdraft, in K
+    
+  integer return_code(klon)
 
   ! -------------------------------------------------------------------------
@@ -348 +384 @@
   WRITE(*,*) 'flag_wk_check_trgl=', flag_wk_check_trgl
 
+  CALL getin_p('z_cp0', z_cp0)
+  CALL getin_p('r0', r0)
+  CALL getin_p('h0', h0)
+  CALL getin_p('eps_cp', eps_cp)
+  CALL getin_p('beta', beta)
+  CALL getin_p('av_thvu_constant', av_thvu_constant)
+  CALL getin_p('k_star', k_star)
   first=.false.
  endif
@@ -1970 +2013 @@
 !>jyg
 
-
+  do i = 1, klon
+     if (ktopw(i) /= 0) then
+        rate_uns(i) = - amdwn(i, ktopw(i)) / wdens(i)
+     else
+        rate_uns(i) = -1.
+     end if
+  end do
+
+  theta_vs = zxtsol * (1e5 / ph(:, 1))**rkappa * (1 + epsim1 * q_seri(:, 1))
+
+  forall (i = 1:klon, ktopw(i) /= 0 .and. theta_vs(i) >= av_thvu(i)) &
+       theta_v_uns(i) = the(i, ktopw(i)) * (1. + epsim1 * qe(i, ktopw(i)))
+  
+  call WriteField_phy("thu", thu, klev)
+  call WriteField_phy("qu", qu, klev)
+  call WriteField_phy("av_thvu", av_thvu, 1)
+  call WriteField_phy("ktopw", real(ktopw), 1)
+  call WriteField_phy("wdens", wdens, 1)
+  call WriteField_phy("av_rho", av_rho, 1)
+  call WriteField_phy("rate_uns", rate_uns, 1)
+  call WriteField_phy("theta_vs", theta_vs, 1)
+  call WriteField_phy("the", the, klev)
+  call WriteField_phy("qe", qe, klev)
+
+  do i = 1, klon
+     call cold_pool_death(z_cp0, r0, h0, eps_cp, beta, av_thvu_constant, &
+          k_star, av_thvu(i), rate_uns(i), theta_vs(i), theta_v_uns(i), &
+          sens_w(i) / (av_rho(i) * rcpd * (theta_vs(i) - theta_v_uns(i))), &
+          av_rho(i), return_code(i), wake_lifetime(i), wake_final_radius(i), &
+          wake_final_height(i))
+  end do
+  if (any(return_code /= 0)) print *, &
+       "cold_pool_death: return_code /= 0 somewhere, itap  = ", itap
 
   RETURN