Skip to content

Add RUC LSM support #297

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
grantfirl merged 14 commits into from
Feb 11, 2022
Merged

Add RUC LSM support #297

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
14 commits
Select commit Hold shift + click to select a range
f830d3c
save work for RUC LSM
grantfirl Oct 6, 2021
3ee0eb3
Merge branch 'main' into RUC_LSM_support
grantfirl Oct 14, 2021
ea38552
more RUC LSM support work
grantfirl Oct 19, 2021
6ac4fd7
update gabls3 case data files
Oct 19, 2021
1009219
add GABLS3 case info for RUC LSM (before tuning surface fluxes)
grantfirl Oct 22, 2021
0994509
retune GABLS3 NOAH and RUC LSM cases
Oct 22, 2021
9220b29
update *typedefs to handle lsm_cold_start in CCPP
grantfirl Oct 25, 2021
6a61f7e
set Sfcprop%zorl instead of zorlw for specified surface flux case
grantfirl Oct 28, 2021
ac5de24
add RUC LSM soil temperature/moisture to output routines; fix example…
grantfirl Oct 29, 2021
e29e258
fix a few bugs in RUC LSM initialization caught in debug mode
grantfirl Oct 30, 2021
80a1b65
retune gabls3_ruc soil moisture with proper initialization
Oct 30, 2021
e4e2e75
Merge branch 'main' into RUC_LSM_support
grantfirl Feb 2, 2022
d3c3ae1
retune gabls3_ruc case
Feb 2, 2022
522cfa6
update ccpp/physics submodule pointer
grantfirl Feb 9, 2022
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
2 changes: 1 addition & 1 deletion ccpp/physics
Open in desktop
Submodule physics updated 8 files
+3 −3 physics/GFS_phys_time_vary.fv3.F90
+3 −3 physics/GFS_phys_time_vary.fv3.meta
+3 −3 physics/GFS_phys_time_vary.scm.F90
+3 −3 physics/GFS_phys_time_vary.scm.meta
+59 −52 physics/gwdps.f
+4 −4 physics/module_sf_ruclsm.F90
+18 −18 physics/sfc_drv_ruc.F90
+6 −6 physics/sfc_drv_ruc.meta
Binary file modified scm/data/processed_case_input/gabls3.nc
View file
Open in desktop
Binary file not shown.
Binary file modified scm/data/processed_case_input/gabls3_noahmp.nc
View file
Open in desktop
Binary file not shown.
Binary file added scm/data/processed_case_input/gabls3_ruc.nc
View file
Open in desktop
Binary file not shown.
29 changes: 29 additions & 0 deletions scm/etc/case_config/gabls3_ruc.nml
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,29 @@
$case_config
model_name = 'FV3',
n_columns = 1,
case_name = 'gabls3_ruc',
dt = 600.0,
time_scheme = 1,
runtime = 86400,
n_itt_out = 1,
n_itt_diag = 6,
n_levels = 64,
output_file = 'output',
case_data_dir = '../data/processed_case_input',
vert_coord_data_dir = '../data/vert_coord_data',
thermo_forcing_type = 2,
mom_forcing_type = 2,
relax_time = 7200.0,
sfc_flux_spec = .false.,
lsm_ics = .true.,
do_spinup = .true.,
spinup_timesteps = 12,
sfc_roughness_length_cm = 15.0,
sfc_type = 1,
reference_profile_choice = 2,
year = 2006,
month = 7,
day = 1,
hour = 12,
column_area = 1.45E8,
$end
2 changes: 1 addition & 1 deletion scm/etc/scm_qsub_example.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -21,7 +21,7 @@
USER = os.getenv('USER')
STRINGS = [USER, 'ucar.edu']
MY_EMAIL = '@'.join(STRINGS)
PROC = Popen('qsub', shell=True, stdin=PIPE, stdout=PIPE, stderr=PIPE, close_fds=True)
PROC = Popen('qsub -V', shell=True, stdin=PIPE, stdout=PIPE, stderr=PIPE, close_fds=True)

### Begin User-editable section ###
JOB_NAME = "test_job"
Expand Down
18 changes: 13 additions & 5 deletions scm/etc/scripts/GABLS3_LSM.ncl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -431,7 +431,11 @@ begin
tg3@units = "K"
tg3@long_name = "deep soil temperature"

zorlw = 15.0 ;(cm) from case specs
zorl = 15.0 ;(cm) from case specs
zorl@long_name = "composite surface roughness length"
zorl@units = "cm"

zorlw = zorl ;(cm) from case specs
zorlw@long_name = "surface roughness length over ocean"
zorlw@units = "cm"

Expand Down Expand Up @@ -482,7 +486,7 @@ begin

;derive friction velocity from known values
wind1 = sqrt(u_specified(1)^2 + v_specified(1)^2)
uustar = von_K*wind1/log(u_h(1)/(0.01*zorlw))
uustar = von_K*wind1/log(u_h(1)/(0.01*zorl))
uustar@units = "m s-1"
uustar@long_name = "friction velocity"

Expand Down Expand Up @@ -533,19 +537,19 @@ begin
tsfcl@long_name = "surface skin temperature over land"
tsfcl@units = "K"

zorll = zorlw ;(cm) from case specs
zorll = zorl ;(cm) from case specs
zorll@long_name = "surface roughness length over land"
zorll@units = "cm"

zorli = zorlw ;(cm) from case specs
zorli = zorl ;(cm) from case specs
zorli@long_name = "surface roughness length over ice"
zorli@units = "cm"

stc = (/ 292.55, 289.9, 285.35, 283.15 /)
stc@units = "K"
stc@long_name = "initial profile of soil temperature"

smc = (/ 0.197, 0.197, 0.197, 0.197 /)
smc = (/ 0.203, 0.203, 0.203, 0.203 /)
smc@units = "m3 m-3"
smc@long_name = "initial profile of soil moisture"

Expand Down Expand Up @@ -636,6 +640,10 @@ begin
filevarattdef(g1,"tg3",tg3)
g1->tg3 = tg3

filevardef(g1,"zorl",typeof(zorl),"ncl_scalar")
filevarattdef(g1,"zorl",zorl)
g1->zorl = zorl

filevardef(g1,"zorlw",typeof(zorlw),"ncl_scalar")
filevarattdef(g1,"zorlw",zorlw)
g1->zorlw = zorlw
Expand Down
16 changes: 12 additions & 4 deletions scm/etc/scripts/GABLS3_LSM_NoahMP.ncl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -431,7 +431,11 @@ begin
tg3@units = "K"
tg3@long_name = "deep soil temperature"

zorlw = 15.0 ;(cm) from case specs
zorl = 15.0 ;(cm) from case specs
zorl@long_name = "composite surface roughness length"
zorl@units = "cm"

zorlw = zorl ;(cm) from case specs
zorlw@long_name = "surface roughness length over ocean"
zorlw@units = "cm"

Expand Down Expand Up @@ -482,7 +486,7 @@ begin

;derive friction velocity from known values
wind1 = sqrt(u_specified(1)^2 + v_specified(1)^2)
uustar = von_K*wind1/log(u_h(1)/(0.01*zorlw))
uustar = von_K*wind1/log(u_h(1)/(0.01*zorl))
uustar@units = "m s-1"
uustar@long_name = "friction velocity"

Expand Down Expand Up @@ -533,11 +537,11 @@ begin
tsfcl@long_name = "surface skin temperature over land"
tsfcl@units = "K"

zorll = zorlw ;(cm) from case specs
zorll = zorl ;(cm) from case specs
zorll@long_name = "surface roughness length over land"
zorll@units = "cm"

zorli = zorlw ;(cm) from case specs
zorli = zorl ;(cm) from case specs
zorli@long_name = "surface roughness length over ice"
zorli@units = "cm"

Expand Down Expand Up @@ -636,6 +640,10 @@ begin
filevarattdef(g1,"tg3",tg3)
g1->tg3 = tg3

filevardef(g1,"zorl",typeof(zorl),"ncl_scalar")
filevarattdef(g1,"zorl",zorl)
g1->zorl = zorl

filevardef(g1,"zorlw",typeof(zorlw),"ncl_scalar")
filevarattdef(g1,"zorlw",zorlw)
g1->zorlw = zorlw
Expand Down
Loading