diff --git a/physics/sfc_drv_ruc.F90 b/physics/sfc_drv_ruc.F90 index 985bee414..ec3a3c62b 100644 --- a/physics/sfc_drv_ruc.F90 +++ b/physics/sfc_drv_ruc.F90 @@ -16,8 +16,6 @@ module lsm_ruc public :: lsm_ruc_init, lsm_ruc_run, lsm_ruc_finalize - real(kind=kind_phys), parameter :: zero = 0.0d0, one = 1.0d0, epsln = 1.0d-10 - contains !> This subroutine calls set_soilveg_ruc() to specify vegetation and soil parameters for @@ -260,7 +258,7 @@ subroutine lsm_ruc_run & ! inputs & srflag, snoalb, isot, ivegsrc, fice, smcwlt2, smcref2, & ! --- constants & con_cp, con_rd, con_rv, con_g, con_pi, con_hvap, & - & con_fvirt, & + & con_fvirt, min_seaice, con_tice, rhowater, & ! for water & ch_wat, tskin_wat, & ! --- in/outs for ice and land @@ -291,11 +289,7 @@ subroutine lsm_ruc_run & ! inputs implicit none ! --- constant parameters: - real(kind=kind_phys), parameter :: rhoh2o = 1000.0 real(kind=kind_phys), parameter :: stbolt = 5.670400e-8 - real(kind=kind_phys), parameter :: cimin = 0.15 !--- in GFS - !real(kind=kind_phys), parameter :: cimin = 0.02 !--- minimum ice concentration, 0.15 in GFS - real(kind=kind_phys), parameter :: con_tice = 271.2 ! --- input: integer, intent(in) :: me, master @@ -317,7 +311,8 @@ subroutine lsm_ruc_run & ! inputs real (kind=kind_phys), intent(in) :: delt real (kind=kind_phys), intent(in) :: con_cp, con_rv, con_g, & con_pi, con_rd, & - con_hvap, con_fvirt + con_hvap, con_fvirt, & + min_seaice, con_tice, rhowater logical, dimension(im), intent(in) :: flag_iter, flag_guess logical, dimension(im), intent(in) :: land, icy, lake @@ -419,7 +414,7 @@ subroutine lsm_ruc_run & ! inputs & qsfc_ice, qsg_ice, qvg_ice, qcg_ice, soilt1_ice - real (kind=kind_phys) :: xice_threshold + real (kind=kind_phys) :: xice_threshold, cimin, rhoh2o real (kind=kind_phys) :: fwat, qsw, evapw, hfxw character(len=256) :: llanduse !< Land-use dataset. Valid values are : @@ -454,6 +449,9 @@ subroutine lsm_ruc_run & ! inputs chklowq = 1. + cimin = min_seaice + rhoh2o = rhowater + do i = 1, im ! i - horizontal loop ! - Set flag for ice points for uncoupled model (islmsk(i) == 4 when coupled to CICE) ! - Exclude ice on the lakes if the lake model is turned on. diff --git a/physics/sfc_drv_ruc.meta b/physics/sfc_drv_ruc.meta index ba61ed899..0017c5763 100644 --- a/physics/sfc_drv_ruc.meta +++ b/physics/sfc_drv_ruc.meta @@ -798,6 +798,33 @@ kind = kind_phys intent = in optional = F +[min_seaice] + standard_name = sea_ice_minimum + long_name = minimum sea ice value + units = frac + dimensions = () + type = real + kind = kind_phys + intent = in + optional = F +[con_tice] + standard_name = freezing_point_temperature_of_seawater + long_name = freezing point temperature of seawater + units = K + dimensions = () + type = real + kind = kind_phys + intent = in + optional = F +[rhowater] + standard_name = liquid_water_density + long_name = density of liquid water + units = kg m-3 + dimensions = () + type = real + kind = kind_phys + intent = in + optional = F [ch_wat] standard_name = surface_drag_coefficient_for_heat_and_moisture_in_air_over_ocean long_name = surface exchange coeff heat & moisture over ocean