Skip to content
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

Jump to bottom

Add option to derive the override for the net surface shortwave radiative flux #253

Draft
wants to merge 3 commits into
base: master
Choose a base branch
from
Draft

Add option to derive the override for the net surface shortwave radiative flux #253

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
8c749c1
Add option to derive net shortwave radiative flux in physics driver
spencerkclark Dec 17, 2021
89c4b30
Update name of namelist parameter
spencerkclark Dec 17, 2021
1ef20c8
Be sure to compute derived net shortwave after adjsfcdsw and adjsfcns…
spencerkclark Dec 19, 2021
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
8 changes: 8 additions & 0 deletions FV3/gfsphysics/GFS_layer/GFS_physics_driver.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1457,6 +1457,14 @@ subroutine GFS_physics_driver &
!---------------------------------------------------------------------
endif
!
if (Model%override_surface_radiative_fluxes .and. Model%derive_net_surface_shortwave_radiative_flux) then
where (adjsfcdsw .gt. adjsfcnsw)
Statein%adjsfcnsw_override = (adjsfcnsw / adjsfcdsw) * Statein%adjsfcdsw_override
elsewhere
Statein%adjsfcnsw_override = 0.0
endwhere
endif

if (Model%lsidea) then !idea jw
dtdt(:,:) = zero
endif
Expand Down
10 changes: 8 additions & 2 deletions FV3/gfsphysics/GFS_layer/GFS_typedefs.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1084,6 +1084,7 @@ module GFS_typedefs
logical :: iau_filter_increments
real(kind=kind_phys) :: sst_perturbation ! Sea surface temperature perturbation to climatology or nudging SST (default 0.0 K)
logical :: override_surface_radiative_fluxes ! Whether to use Statein to override the surface radiative fluxes
logical :: derive_net_surface_shortwave_radiative_flux ! Whether to compute the net shortwave radiative flux using override downward shortwave flux and the RRTMG effective albedo
logical :: use_climatological_sst ! Whether to allow the Python wrapper to override the sea surface temperature
logical :: emulate_zc_microphysics ! Use an emulator in place of ZC microphysics
logical :: save_zc_microphysics ! Save ZC microphysics state
Expand Down Expand Up @@ -3162,6 +3163,7 @@ subroutine control_initialize (Model, nlunit, fn_nml, me, master, &

real(kind=kind_phys) :: sst_perturbation = 0.0 ! Sea surface temperature perturbation [K]
logical :: override_surface_radiative_fluxes = .false.
logical :: derive_net_surface_shortwave_radiative_flux = .false.
logical :: use_climatological_sst = .true.
logical :: emulate_zc_microphysics = .false.
logical :: save_zc_microphysics = .false.
Expand Down Expand Up @@ -3256,8 +3258,10 @@ subroutine control_initialize (Model, nlunit, fn_nml, me, master, &
!--- aerosol scavenging factors ('name:value' string array)
fscav_aero, &
sst_perturbation, &
override_surface_radiative_fluxes, use_climatological_sst, &
emulate_zc_microphysics, save_zc_microphysics
override_surface_radiative_fluxes, &
derive_net_surface_shortwave_radiative_flux, &
use_climatological_sst, emulate_zc_microphysics, &
save_zc_microphysics

!--- other parameters
integer :: nctp = 0 !< number of cloud types in CS scheme
Expand Down Expand Up @@ -3726,6 +3730,7 @@ subroutine control_initialize (Model, nlunit, fn_nml, me, master, &

Model%sst_perturbation = sst_perturbation
Model%override_surface_radiative_fluxes = override_surface_radiative_fluxes
Model%derive_net_surface_shortwave_radiative_flux = derive_net_surface_shortwave_radiative_flux
Model%use_climatological_sst = use_climatological_sst

!--- emulation parameters
Expand Down Expand Up @@ -4544,6 +4549,7 @@ subroutine control_print(Model)
print *, ' isot : ', Model%isot
print *, ' sst_perturbation : ', Model%sst_perturbation
print *, ' override_surface_radiative_fluxes: ', Model%override_surface_radiative_fluxes
print *, ' derive_net_surface_shortwave_radiative_flux: ', Model%derive_net_surface_shortwave_radiative_flux
print *, ' use_climatological_sst: ', Model%use_climatological_sst
if (Model%lsm == Model%lsm_noahmp) then
print *, ' Noah MP LSM is used, the options are'
Expand Down