Skip to content

Removed SNET. #15

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
climbfuji merged 1 commit into from
May 25, 2021
Merged

Removed SNET. #15

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
30 changes: 6 additions & 24 deletions physics/sfc_drv_ruc.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -277,7 +277,6 @@ end subroutine lsm_ruc_finalize
! sigmaf - real, areal fractional cover of green vegetation im !
! dlwflx - real, total sky sfc downward lw flux ( w/m**2 ) im !
! dswflx - real, total sky sfc downward sw flux ( w/m**2 ) im !
! snet - real, total sky sfc netsw flx into ground(w/m**2) im !
! delt - real, time interval (second) 1 !
! tg3 - real, deep soil temperature (k) im !
! cm - real, surface exchange coeff for momentum (m/s) im !
Expand Down Expand Up @@ -337,7 +336,7 @@ subroutine lsm_ruc_run & ! inputs
& imp_physics, imp_physics_gfdl, imp_physics_thompson, &
& do_mynnsfclay, lsoil_ruc, lsoil, rdlai, zs, &
& t1, q1, qc, soiltyp, vegtype, sigmaf, laixy, &
& dlwflx, dswsfc, snet, tg3, coszen, land, icy, lake, &
& dlwflx, dswsfc, tg3, coszen, land, icy, lake, &
& rainnc, rainc, ice, snow, graupel, &
& prsl1, zf, wind, shdmin, shdmax, &
& srflag, sfalb_lnd_bck, snoalb, &
Expand Down Expand Up @@ -389,7 +388,7 @@ subroutine lsm_ruc_run & ! inputs
integer, intent(in) :: imp_physics, imp_physics_gfdl, imp_physics_thompson

real (kind=kind_phys), dimension(:), intent(in) :: &
& t1, sigmaf, laixy, dlwflx, dswsfc, snet, tg3, &
& t1, sigmaf, laixy, dlwflx, dswsfc, tg3, &
& coszen, prsl1, wind, shdmin, shdmax, &
& sfalb_lnd_bck, snoalb, zf, qc, q1, &
! for land
Expand Down Expand Up @@ -938,17 +937,8 @@ subroutine lsm_ruc_run & ! inputs

snoalb1d_lnd(i,j) = snoalb(i)
albbck_lnd(i,j) = albbcksol(i) !sfalb_lnd_bck(i)
! alb_lnd takes into account snow on the ground
!if (kdt == 1) then
! if (dswsfc(i) > 0.) then
! alb_lnd(i,j) = max(0.01, 1. - snet(i)/dswsfc(i))
! else
! alb_lnd(i,j) = albbck_lnd(i,j) * (1.-sncovr_lnd(i,j)) + snoalb(i) * sncovr_lnd(i,j)
! endif
!else
alb_lnd(i,j) = albbck_lnd(i,j) * (1.-sncovr_lnd(i,j)) + snoalb(i) * sncovr_lnd(i,j) ! sfalb_lnd(i)
!endif
solnet_lnd(i,j) = snet(i) !dswsfc(i)*(1.-alb_lnd(i,j)) !..net sw rad flx (dn-up) at sfc in w/m2
solnet_lnd(i,j) = dswsfc(i)*(1.-alb_lnd(i,j)) !..net sw rad flx (dn-up) at sfc in w/m2

cmc(i,j) = canopy(i) ! [mm]
soilt_lnd(i,j) = tsurf_lnd(i) ! clu_q2m_iter
Expand Down Expand Up @@ -1239,21 +1229,13 @@ subroutine lsm_ruc_run & ! inputs
!-- alb_ice* is computed in setalb called from rrtmg_sw_pre.
snoalb1d_ice(i,j) = 0.75 !alb_ice_snow(i) !0.75 is RAP value for max snow alb on ice
albbck_ice(i,j) = 0.55 !alb_ice_snowfree(i) !0.55 is RAP value for ice alb
if (kdt == 1) then
if (dswsfc(i) > 0.) then
alb_ice(i,j) = max(0.01, 1. - snet(i)/dswsfc(i))
else
alb_ice(i,j) = albbck_ice(i,j) * (1.-sncovr_ice(i,j)) + snoalb1d_ice(i,j) * sncovr_ice(i,j)
endif
else
alb_ice(i,j) = sfalb_ice(i)
endif
solnet_ice(i,j) = snet(i) !dswsfc(i)*(1.-alb_ice(i,j))
alb_ice(i,j) = sfalb_ice(i)
solnet_ice(i,j) = dswsfc(i)*(1.-alb_ice(i,j))
qvg_ice(i,j) = sfcqv_ice(i)
qsfc_ice(i,j) = sfcqv_ice(i)/(1.+sfcqv_ice(i))
qsg_ice(i,j) = rslf(prsl1(i),tsurf_ice(i))
qcg_ice(i,j) = sfcqc_ice(i)
semis_bck(i,j) = 0.99
semis_bck(i,j) = 0.99
if (kdt == 1) then
sfcems_ice(i,j) = semisbase(i) * (1.-sncovr_ice(i,j)) + 0.99 * sncovr_ice(i,j)
else
Expand Down
9 changes: 0 additions & 9 deletions physics/sfc_drv_ruc.meta
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -810,15 +810,6 @@
kind = kind_phys
intent = in
optional = F
[snet]
standard_name = surface_net_downwelling_shortwave_flux
long_name = surface net downwelling shortwave flux at current time
units = W m-2
dimensions = (horizontal_loop_extent)
type = real
kind = kind_phys
intent = in
optional = F
[tg3]
standard_name = deep_soil_temperature
long_name = deep soil temperature
Expand Down