Skip to content

Commit

Permalink
Merge pull request #4 from tanyasmirnova/PR_moorthi_tanya_30sept21
Browse files Browse the repository at this point in the history 
Addition to Moorthi's PR
  • Loading branch information
@SMoorthi-emc
SMoorthi-emc authored Sep 30, 2021
2 parents d6f8ad4 + 19e675a commit 6c709c9
Show file tree
Hide file tree
Showing 2 changed files with 7 additions and 65 deletions.
50 changes: 5 additions & 45 deletions physics/sfc_drv_ruc.F90
View file
Original file line number Diff line number Diff line change
Expand Up @@ -353,8 +353,6 @@ subroutine lsm_ruc_run & ! inputs
! --- constants
& con_cp, con_rd, con_rv, con_g, con_pi, con_hvap, &
& con_fvirt, &
! for water
& ch_wat, tskin_wat, &
! --- in/outs for ice and land
& semisbase, semis_lnd, semis_ice, sfalb_lnd, sfalb_ice, &
& sncovr1_lnd, weasd_lnd, snwdph_lnd, tskin_lnd, &
Expand All @@ -370,8 +368,7 @@ subroutine lsm_ruc_run & ! inputs
& albdvis_lnd, albdnir_lnd, albivis_lnd, albinir_lnd, &
! for ice
& sfcqc_ice, sfcqv_ice, &
& tsurf_ice, tsnow_ice, z0rl_ice, &
! & tice, tsurf_ice, tsnow_ice, z0rl_ice, &
& tsurf_ice, tsnow_ice, z0rl_ice, &
& qsurf_ice, gflux_ice, evap_ice, ep1d_ice, hflx_ice, &
& cm_ice, ch_ice, snowfallac_ice, &
& albdvis_ice, albdnir_ice, albivis_ice, albinir_ice, &
Expand All @@ -388,7 +385,6 @@ subroutine lsm_ruc_run & ! inputs
! --- constant parameters:
real(kind=kind_phys), parameter :: rhoh2o = 1000.0
real(kind=kind_phys), parameter :: stbolt = 5.670400e-8
real(kind=kind_phys), parameter :: con_tice = 271.2

! --- input:
integer, intent(in) :: me, master
Expand All @@ -403,7 +399,7 @@ subroutine lsm_ruc_run & ! inputs
! for land
& cm_lnd, ch_lnd, &
! for water
& ch_wat, tskin_wat, oceanfrac, &
& oceanfrac, &
! for ice
& cm_ice, ch_ice

Expand Down Expand Up @@ -557,7 +553,7 @@ subroutine lsm_ruc_run & ! inputs

do i = 1, im ! i - horizontal loop
flag_ice(i) = .false.
if (icy(i) .and. .not. flag_cice(i)) then
if (icy(i) .and. .not. flag_cice(i)) then ! flag_cice(i)=.true. when coupled to CICE
! - uncoupled ice model
if (oceanfrac(i) > zero) then
cimin(i) = min_seaice
Expand All @@ -569,8 +565,8 @@ subroutine lsm_ruc_run & ! inputs
flag_ice(i) = .true.
endif
endif
! - 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.
! - Ice points for uncoupled model
! - Exclude ice on the lakes if the lake model is turned on: lake(i)=.true.
flag_ice_uncoupled(i) = (flag_ice(i) .and. .not. lake(i))
!> - Set flag for land and ice points.
!- 10may19 - ice points are turned off.
Expand Down Expand Up @@ -1388,42 +1384,6 @@ subroutine lsm_ruc_run & ! inputs
enddo ! j
enddo ! i

!-- Take care of fractional sea ice for uncoupled run with frac_grid=.false.
!-- When frac_grid=.true. GFS_surface_composite will take care of this.
do i = 1, im ! i - horizontal loop
if ( flag_iter(i) .and. flag(i) ) then
! Do this only when the fractional grid is not turned on!
! Compute composite for a fractional sea ice: fice(i) < 1.
! This is needed for the 2-way coupling
! in the upcoupled case (when sfc_cice is not used).
if(.not. frac_grid) then
if( flag_ice_uncoupled(i) .and. fice(i) < 1.) then
!write (0,*)'Fractional sea ice at i', i, fice(i)
fwat = 1.0 - fice(i)
! Check if ice fraction is below the minimum value: 15% in GFS
! physics.
if (fice(i) < cimin(i)) then ! cimin - minimal ice fraction
write (0,*)'warning: ice fraction is low:', fice(i)
fice(i) = cimin(i)
fwat = 1.0 - cimin(i)
write (0,*)'fix ice fraction: reset it to:', fice(i), tskin_wat(i)
endif

! Compute the composite of ice and open water for 2-way coupling in the
! uncoupled sea-ice model. Use ice variables for the composite.
tsurf_ice(i) = tsurf_ice(i) * fice(i) + min(con_tice,tskin_wat(i)) * fwat
chh_ice(i) = chh_ice(i) * fice(i) + ch_wat(i) * wind(i) * rho(i) * fwat
hfxw = ch_wat(i) * wind(i) * (min(con_tice,tskin_wat(i)) - t1(i))
hflx_ice(i) = hflx_ice(i) * fice(i) + hfxw * fwat
qsw = rslf(prsl1(i),min(con_tice,tskin_wat(i)))
evapw = ch_wat(i) * wind(i) * (qsw - q0(i))
evap_ice(i) = evap_ice(i) * fice(i) + evapw * fwat
qsurf_ice(i) = q1(i) + evap_ice(i) * rho(i) / chh_ice(i)
endif ! flag_ice_uncoupled(i) .and. fice(i) < 1.
endif ! flag_iter, icy, not frac_grid
endif
enddo ! i

!> - Restore land-related prognostic fields for guess run.
do j = 1, 1
do i = 1, im
Expand Down
22 changes: 2 additions & 20 deletions physics/sfc_drv_ruc.meta
View file
Original file line number Diff line number Diff line change
Expand Up @@ -192,8 +192,8 @@
intent = in
optional = F
[tsfc_ice]
standard_name = surface_skin_temperature_over_ice
long_name = surface skin temperature over ice
standard_name = sea_ice_temperature
long_name = sea ice surface skin temperature
units = K
dimensions = (horizontal_dimension)
type = real
Expand Down Expand Up @@ -1102,24 +1102,6 @@
kind = kind_phys
intent = in
optional = F
[ch_wat]
standard_name = surface_drag_coefficient_for_heat_and_moisture_in_air_over_water
long_name = surface exchange coeff heat surface exchange coeff heat & moisture over ocean moisture over water
units = none
dimensions = (horizontal_loop_extent)
type = real
kind = kind_phys
intent = in
optional = F
[tskin_wat]
standard_name = surface_skin_temperature_over_water
long_name = surface skin temperature over water
units = K
dimensions = (horizontal_loop_extent)
type = real
kind = kind_phys
intent = in
optional = F
[semisbase]
standard_name = baseline_surface_longwave_emissivity
long_name = baseline surface lw emissivity in fraction
Expand Down

0 comments on commit 6c709c9

Please sign in to comment.