Updates of interstitial *pre for MM5 SFCLAY_REV

modified: gfs_mmm_sfclayrev_pre.F90 modified: gfs_mmm_sfclayrev_pre.meta
bluefinweiwei · Mar 7, 2024 · 738ccd7 · 738ccd7
1 parent a598993
commit 738ccd7
Show file tree

Hide file tree

Showing 2 changed files with 56 additions and 23 deletions.
diff --git a/physics/gfs_mmm_sfclayrev_pre.F90 b/physics/gfs_mmm_sfclayrev_pre.F90
@@ -17,15 +17,16 @@ module gfs_mmm_sfclayrev_pre
   !
   ! #########################################################################################
   subroutine gfs_mmm_sfclayrev_pre_run(do_mmm_sfclayrev, nCol, zorl, znt, wetness, mavail, pbl_regime, fm, fh, cpm, slimask, xland, &
-             wet, dry, icy, tsfc_wat, tsfc_lnd, tsfc_ice, qsfc_wat, qsfc_lnd, qsfc_ice, tsfc, qsfc, &
-             ep1, t, q, ps, con_rd, con_cp, hfx_wat, hfx_lnd, hfx_ice, hfx, qfx_wat, qfx_lnd, qfx_ice,qfx, scm_force_flux, &
-             qgh, gz1oz0, water_depth, shalwater_depth, svp1, svp2, svp3, svpt0, p1000mb, its, ite, errmsg, errflg)
+             wet, dry, icy, tsfc_wat, tsfc_lnd, tsfc_ice, qsfc_wat, qsfc_lnd, qsfc_ice, tsfc, qsfc, shum1d, rv1d, &
+             ep1, t, q, ps, con_rd, con_cp, hfx_wat, hfx_lnd, hfx_ice, hfx, qfx_wat, qfx_lnd, qfx_ice,qfx, scm_force_flux, isfflx, &
+             qgh, gz1oz0, water_depth, shalwater_depth, con_psat, svp1, svp2, svp3, p1000mb, its, ite, errmsg, errflg)
 
     ! Input
     logical,                       intent(in)           :: &
         do_mmm_sfclayrev  ! Flag for MMM SFCLAYREV scheme
     logical,                       intent(out)          :: &
-        scm_force_flux    ! Flag for not computing surface fluxes but prescribing
+        scm_force_flux, & ! Flag for not computing surface fluxes but prescribing them from a SCM
+        isfflx            ! Flag for computing surface heat and moisture fluxes in surface layer scheme  
     logical,         dimension(:), intent(in)           :: &
         wet,            & ! Flag indicating presence of some ocean or lake surface area fraction
         dry,            & ! Flag indicating presence of some land surface area fraction
@@ -39,7 +40,8 @@ subroutine gfs_mmm_sfclayrev_pre_run(do_mmm_sfclayrev, nCol, zorl, znt, wetness,
     real(kind_phys),               intent(in)           :: &
         ep1,            & ! (rv/rd) - 1
         con_rd,         & ! rd
-        con_cp            ! cp
+        con_cp,         & ! cp
+        con_psat          ! con_psat pres at H2O 3pt (Pa)
     real(kind_phys), dimension(:), intent(in)           :: &
         zorl,           & ! Surface roughness-length (cm)
         tsfc_wat,       & ! Surface skin temperature over water (K)
@@ -54,6 +56,7 @@ subroutine gfs_mmm_sfclayrev_pre_run(do_mmm_sfclayrev, nCol, zorl, znt, wetness,
         qfx_wat,        & ! Kinematic surface upward latent heat flux over water (kg kg-1 m s-1)
         qfx_lnd,        & ! Kinematic surface upward latent heat flux over land (kg kg-1 m s-1)
         qfx_ice,        & ! Kinematic surface upward latent heat flux over ice (kg kg-1 m s-1)
+        shum1d,         & ! Water vapor specific humidity at lowest model layer (kg kg-1)
         ps                ! Surface pressure (Pa)
     real(kind_phys), dimension(:), intent(in), optional :: &
         wetness           ! Normalized soil wetness (frac)  
@@ -63,10 +66,11 @@ subroutine gfs_mmm_sfclayrev_pre_run(do_mmm_sfclayrev, nCol, zorl, znt, wetness,
         q                 ! Model layer mean tracer concentration (kg kg-1)
 
     ! Outputs
-    real(kind=kind_phys),           intent(out) :: shalwater_depth, svp1, svp2, svp3, svpt0, p1000mb
+    real(kind=kind_phys),          intent(out) :: shalwater_depth, svp1, svp2, svp3, p1000mb
     real(kind_phys), dimension(:), intent(out) :: &
         tsfc,           & ! Surface skin temperature (K)
         qsfc,           & ! Surface specific humidity (kg kg-1)
+        rv1d,           & ! Water vapor mixing ratio at the lowest model level (kg kg-1)
         znt,            & ! Surface roughness-length (m)
         mavail,         & ! Surface moisture availability (frac; [0-1])
         pbl_regime,     & ! PBL regime categories
@@ -109,26 +113,32 @@ subroutine gfs_mmm_sfclayrev_pre_run(do_mmm_sfclayrev, nCol, zorl, znt, wetness,
     ! Pass in flag for computing or prescribing surface fluxes
     if (do_mmm_sfclayrev) then
         scm_force_flux = .false.
+        isfflx = .true.
     endif
 
-    ! Constants from module_sf_mynn.F90
     shalwater_depth=1.0
-    svp1 = 0.6112
+
+    ! Constants for calculating saturation vapor pressure of water (Bolton 1980)
+    svp1 = con_psat * 0.001 ! convert unit from Pa to kPa
     svp2 = 17.67
     svp3 = 29.65
-    svpt0 = 273.15
     p1000mb = 100000.
 
     ! Convert roughness length (cm -> m)
     do iCol = 1,nCol
         znt(iCol) = zorl(iCol) * 0.01
     end do
     ! Calculate surface available moisture (fraction [0-1])
-    mavail(:) = 1.0     ! WL2023 - should mavail be separated for dry, icy and wet?
-    if (present(wetness)) then
+    mavail(:) = 1.0 ! WL2023: lower values may be used - TBD    
+    if (present(wetness)) then ! only when RUC LSM used
         mavail(:) = wetness(:)
     end if
 
+    ! Convertion between specific humidity of water vapor to water vapor mixing ratio
+    do iCol=1,nCol
+       rv1d(iCol) = shum1d(iCol) / (1.0 - shum1d(iCol))
+    enddo
+
     ! Compute land/sea mask convention from (0-sea/1-land/2-ice) ---> (1-land/2-sea)
     do iCol=1,nCol
         if(slimask(iCol).eq.0)then
@@ -139,7 +149,7 @@ subroutine gfs_mmm_sfclayrev_pre_run(do_mmm_sfclayrev, nCol, zorl, znt, wetness,
     enddo
 
     do iCol = 1,nCol
-       tvcon = (1. + ep1 * q(iCol,1,1))
+       tvcon = (1.0 + ep1 * q(iCol,1,1))
        rho   = ps(iCol) / (con_rd * t(iCol,1) * tvcon)
         if (dry(iCol)) then
             hfx(iCol)  = rho * con_cp * hfx_lnd(iCol)

diff --git a/physics/gfs_mmm_sfclayrev_pre.meta b/physics/gfs_mmm_sfclayrev_pre.meta
@@ -185,6 +185,22 @@
   type = real
   kind = kind_phys
   intent = out
+[shum1d]
+  standard_name = specific_humidity_at_surface_adjacent_layer
+  long_name = water vapor specific humidity at surface adjacent layer
+  units = kg kg-1
+  dimensions = (horizontal_loop_extent)
+  type = real
+  kind = kind_phys
+  intent = in
+[rv1d]
+  standard_name = water_vapor_mixing_ratio_at_surface_adjacent_layer
+  long_name =  water vapor mixing ratio at surface adjacent layer
+  units = kg kg-1
+  dimensions = (horizontal_loop_extent)
+  type = real
+  kind = kind_phys
+  intent = out
 [ep1]
   standard_name = ratio_of_vapor_to_dry_air_gas_constants_minus_one
   long_name = rv/rd - 1 (rv = ideal gas constant for water vapor)
@@ -304,6 +320,13 @@
   dimensions = ()
   type = logical
   intent = out
+[isfflx]
+  standard_name = flag_to_compute_surface_heat_and_moisture_fluxes
+  long_name = flag to compute surface heat and moisture fluxes
+  units = flag
+  dimensions = ()
+  type = logical
+  intent = out
 [qgh]
   standard_name = saturation_mixing_ratio_at_lowest_model_level
   long_name = saturation water vapor mixing ratio at lowest model level
@@ -336,10 +359,18 @@
   type = real
   kind = kind_phys
   intent = out
+[con_psat]
+  standard_name = saturation_pressure_at_triple_point_of_water
+  long_name = saturation pressure at triple point of water
+  units = Pa
+  dimensions = ()
+  type = real
+  kind = kind_phys
+  intent = in
 [svp1]
-  standard_name = saturation_vapor_pressure_constant_1
-  long_name = constant for saturation vapor pressure calculation
-  units = k Pa
+  standard_name = saturation_pressure_at_triple_point_of_water_in_kPa
+  long_name = saturation pressure at triple point of water in kPa
+  units = kPa
   dimensions = ()
   type = real
   kind = kind_phys
@@ -360,14 +391,6 @@
   type = real
   kind = kind_phys
   intent = out
-[svpt0]
-  standard_name = saturation_vapor_pressure_constant_0
-  long_name = constant for saturation vapor pressure calculation
-  units = K
-  dimensions = ()
-  type = real
-  kind = kind_phys
-  intent = out
 [p1000mb]
   standard_name = pressure_value_at_1000_hPa
   long_name = pressure value at 1000 hPa