Microphysics tuning PR

FV3/atmos_cubed_sphere/model/fv_cmp.F90

@@ -57,7 +57,7 @@ module fv_cmp_mod
     use constants_mod, only: rvgas, rdgas, grav, hlv, hlf, cp_air
     use fv_mp_mod, only: is_master
     use fv_arrays_mod, only: r_grid
-    use gfdl_cloud_microphys_mod, only: ql_gen, qi_gen, qi0_max, ql_mlt, ql0_max, qi_lim, qs_mlt
+    use gfdl_cloud_microphys_mod, only: ql_gen, qi0_max, ql_mlt, ql0_max, qi_lim, qs_mlt
     use gfdl_cloud_microphys_mod, only: icloud_f, sat_adj0, t_sub, cld_min
     use gfdl_cloud_microphys_mod, only: tau_r2g, tau_smlt, tau_i2s, tau_v2l, tau_l2v, tau_imlt, tau_l2r
     use gfdl_cloud_microphys_mod, only: rad_rain, rad_snow, rad_graupel, dw_ocean, dw_land, tintqs
@@ -109,6 +109,7 @@ module fv_cmp_mod
 
     real :: d0_vap !< the same as dc_vap, except that cp_vap can be cp_vap or cv_vap
     real :: lv00 !< the same as lv0, except that cp_vap can be cp_vap or cv_vap
+    real :: qi_gen ! TODO this is no longer provided in Linjiongs fork.
 
     real, allocatable :: table (:), table2 (:), tablew (:), des2 (:), desw (:)
 

FV3/gfsphysics/physics/gfdl_cloud_microphys.F90

@@ -304,6 +304,11 @@ module gfdl_cloud_microphys_mod
     real :: resmin = 150.0, resmax = 10000.0
     real :: regmin = 300.0, regmax = 10000.0
 
+    ! Linjiong-Noah custom tunining parameters
+    ! TODO give these better names
+    real :: rs_factor = 0.2
+    real :: rh_factor = 10.0
+
     ! -----------------------------------------------------------------------
     ! namelist
     ! -----------------------------------------------------------------------
@@ -320,7 +325,7 @@ module gfdl_cloud_microphys_mod
         rad_snow, rad_graupel, rad_rain, cld_min, use_ppm, mono_prof,         &
         do_sedi_heat, sedi_transport, do_sedi_w, de_ice, icloud_f, irain_f,   &
         mp_print, reiflag, rewmin, rewmax, reimin, reimax, rermin, rermax,    &
-        resmin, resmax, regmin, regmax, tintqs
+        resmin, resmax, regmin, regmax, tintqs, rs_factor, rh_factor
 
     public                                                                    &
         mp_time, t_min, t_sub, tau_r2g, tau_smlt, tau_g2r, dw_land, dw_ocean, &
@@ -334,7 +339,7 @@ module gfdl_cloud_microphys_mod
         rad_snow, rad_graupel, rad_rain, cld_min, use_ppm, mono_prof,         &
         do_sedi_heat, sedi_transport, do_sedi_w, de_ice, icloud_f, irain_f,   &
         mp_print, reiflag, rewmin, rewmax, reimin, reimax, rermin, rermax,    &
-        resmin, resmax, regmin, regmax, tintqs 
+        resmin, resmax, regmin, regmax, tintqs, rs_factor, rh_factor
 
 contains
 
@@ -2083,7 +2088,7 @@ subroutine subgrid_z_proc (ktop, kbot, p1, den, denfac, dts, rh_adj, tz, qv, &
             ! factor = min (1., fac_l2v * sqrt (max (0., ql (k)) / 1.e-5) * 10. * dq0 / qsw)
             ! factor = fac_l2v
             ! factor = 1
-            factor = min (1., fac_l2v * (10. * dq0 / qsw)) ! the rh dependent factor = 1 at 90%
+            factor = min (1., fac_l2v * (rh_factor * dq0 / qsw)) ! the rh dependent factor = 1 at 90%
             evap = min (ql (k), factor * dq0 / (1. + tcp3 (k) * dwsdt))
         else ! condensate all excess vapor into cloud water
             ! -----------------------------------------------------------------------
@@ -2181,7 +2186,7 @@ subroutine subgrid_z_proc (ktop, kbot, p1, den, denfac, dts, rh_adj, tz, qv, &
                 qi_crt = qi_gen * min (qi_lim, 0.1 * tmp) / den (k)
                 sink = min (sink, max (qi_crt - qi (k), pidep), tmp / tcpk (k))
             else ! ice -- > vapor
-                pidep = pidep * min (1., dim (tz (k), t_sub) * 0.2)
+                pidep = pidep * min (1., dim (tz (k), t_sub) * rs_factor)
                 sink = max (pidep, sink, - qi (k))
             endif
             qv (k) = qv (k) - sink