-
Notifications
You must be signed in to change notification settings - Fork 249
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.
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
Add a WAM regression test with the MULTI_GASES and MOLECULAR_DIFFUSION options on #714
Changes from all commits
521e530
8a33360
651fe72
8e3a118
6434f2b
89bf150
da602ea
6fd63b7
241c032
a9de025
a9646e5
f7d84bb
bcc17d5
37c58ec
0c3e334
a4e01c5
777724b
8061851
415fa47
61cca3c
41b2b9b
5a83e59
7ec07e3
fa1f1dc
4ba420c
1cefc23
2f6cc5d
1f5283b
247cfd3
99fb3b8
ed04c62
91a9aa7
e556523
3d7dc97
d26b1f7
7712255
50a5b2d
5af402b
3f4c6d6
05833df
bd62e72
208a045
614e5bc
b1b67b0
ae671ac
be2358d
061b4bf
d465a10
File filter
Filter by extension
Conversations
Jump to
Diff view
Diff view
There are no files selected for viewing
+4 −0 | CMakeLists.txt | |
+1 −1 | atmos_cubed_sphere | |
+3 −0 | ccpp/CMakeLists.txt | |
+21 −10 | ccpp/data/GFS_typedefs.F90 | |
+1 −1 | ccpp/physics | |
+7 −6 | ccpp/suites/suite_FV3_GFS_v16_fv3wam.xml | |
+4 −0 | io/FV3GFS_io.F90 | |
+4 −0 | io/post_gfs.F90 |
Large diffs are not rendered by default.
Large diffs are not rendered by default.
Large diffs are not rendered by default.
Large diffs are not rendered by default.
Large diffs are not rendered by default.
Large diffs are not rendered by default.
Large diffs are not rendered by default.
Large diffs are not rendered by default.
Original file line number | Diff line number | Diff line change |
---|---|---|
|
@@ -33,15 +33,15 @@ | |
ntiles = 6, | ||
npz = @[NPZ] | ||
grid_type = -1 | ||
ncnst = 5 | ||
ncnst = 11 | ||
make_nh = @[MAKE_NH] | ||
fv_debug = .false. | ||
|
||
molecular_diffusion=.true. | ||
range_warn = .false. | ||
reset_eta = .false. | ||
n_sponge = -1 | ||
nudge_qv = .true. | ||
tau = 0.015 | ||
tau = -5.0e1 | ||
rf_cutoff = 2.e2 | ||
d2_bg_k1 = 0.20 | ||
d2_bg_k2 = 0.005 | ||
|
@@ -57,7 +57,7 @@ | |
p_fac = 0.1 | ||
k_split = 2 | ||
n_split = 6 | ||
nwat = 2 | ||
nwat = 6 | ||
na_init = @[NA_INIT] | ||
d_ext = 0.0 | ||
dnats = 0 | ||
|
@@ -83,7 +83,8 @@ | |
hord_dp = 6 | ||
hord_tr = 8 | ||
adjust_dry_mass = .false. | ||
consv_te = 1. | ||
consv_te = 0. | ||
do_sat_adj = .false. | ||
consv_am = .false. | ||
fill = .true. | ||
dwind_2d = .false. | ||
|
@@ -124,7 +125,7 @@ | |
nst_anl = .true. | ||
use_ufo = .true. | ||
pre_rad = .false. | ||
imp_physics = 99 | ||
imp_physics = 11 | ||
pdfcld = .false. | ||
fhswr = 3600. | ||
fhlwr = 3600. | ||
|
@@ -142,16 +143,16 @@ | |
cal_pre = .true. | ||
redrag = .true. | ||
dspheat = .true. | ||
hybedmf = @[HYBEDMF] | ||
satmedmf = @[SATMEDMF] | ||
hybedmf = .false. | ||
satmedmf = .true. | ||
lheatstrg = @[LHEATSTRG] | ||
lseaspray = @[LSEASPRAY] | ||
random_clds = .true. | ||
trans_trac = .false. | ||
cnvcld = .true. | ||
imfshalcnv = 2 | ||
imfdeepcnv = 2 | ||
isatmedmf = @[ISATMEDMF] | ||
isatmedmf = 1 | ||
cdmbgwd = 3.5,0.25 | ||
prslrd0 = 0. | ||
ivegsrc = 1 | ||
|
@@ -170,8 +171,11 @@ | |
iopt_tbot = 2 | ||
iopt_stc = 1 | ||
debug = .false. | ||
levr = 90 | ||
There was a problem hiding this comment. Choose a reason for hiding this commentThe reason will be displayed to describe this comment to others. Learn more. Good to see that the code runs fine with |
||
h2o_phys = @[H2O_PHYS] | ||
nstf_name = @[NSTF_NAME] | ||
oz_phys = .false. | ||
oz_phys_2015 = .true. | ||
cplflx = @[CPLFLX] | ||
cplwav = @[CPLWAV] | ||
iau_delthrs = 6 | ||
|
@@ -183,6 +187,56 @@ | |
lndp_type = @[LNDP_TYPE] | ||
n_var_lndp = @[N_VAR_LNDP] | ||
/ | ||
&gfdl_cloud_microphysics_nml | ||
sedi_transport = .true. | ||
do_sedi_heat = .false. | ||
rad_snow = .true. | ||
rad_graupel = .true. | ||
rad_rain = .true. | ||
const_vi = .false. | ||
const_vs = .false. | ||
const_vg = .false. | ||
const_vr = .false. | ||
vi_max = 1. | ||
vs_max = 2. | ||
vg_max = 12. | ||
vr_max = 12. | ||
qi_lim = 1. | ||
prog_ccn = .false. | ||
do_qa = .true. | ||
fast_sat_adj = .true. | ||
tau_l2v = 225. | ||
tau_v2l = 150. | ||
tau_g2v = 900. | ||
rthresh = 10.e-6 ! This is a key parameter for cloud water | ||
dw_land = 0.16 | ||
dw_ocean = 0.10 | ||
ql_gen = 1.0e-3 | ||
ql_mlt = 1.0e-3 | ||
qi0_crt = 8.0E-5 | ||
qs0_crt = 1.0e-3 | ||
tau_i2s = 1000. | ||
c_psaci = 0.05 | ||
c_pgacs = 0.01 | ||
rh_inc = 0.30 | ||
rh_inr = 0.30 | ||
rh_ins = 0.30 | ||
ccn_l = 300. | ||
ccn_o = 100. | ||
c_paut = 0.5 | ||
c_cracw = 0.8 | ||
use_ppm = .false. | ||
use_ccn = .true. | ||
mono_prof = .true. | ||
z_slope_liq = .true. | ||
z_slope_ice = .true. | ||
de_ice = .false. | ||
fix_negative = .true. | ||
icloud_f = 1 | ||
mp_time = 150. | ||
reiflag = 2 | ||
/ | ||
|
||
|
||
&interpolator_nml | ||
interp_method = 'conserve_great_circle' | ||
|
@@ -245,8 +299,8 @@ | |
SPPT=@[SPPT], | ||
SPPT_TAU=21600, | ||
SPPT_LSCALE=500000, | ||
SPPT_LOGIT=.TRUE., | ||
SPPT_SFCLIMIT=.TRUE., | ||
SPPT_LOGIT=.true., | ||
SPPT_SFCLIMIT=.true., | ||
ISEED_SHUM=1, | ||
ISEED_SKEB=2, | ||
ISEED_SPPT=3, | ||
|
@@ -277,6 +331,16 @@ | |
launch_level = 25 | ||
/ | ||
&multi_gases_nml | ||
ri = 287.050, 461.50, 0.0, 173.2247, 0.0, 0.0 | ||
cpi = 1004.600, 1846.00, 0.0, 820.2391, 0.0, 0.0 | ||
ri= 295.3892, 461.50, 0., 0., 0., 0., 0., 173.2247, 519.674, 259.8370 ,0.,0. | ||
cpi=1031.1083, 1846.00, 0., 0., 0., 0., 0., 820.2391, 1299.185, 918.0969, 0.,0. | ||
/ | ||
&molecular_diffusion_nml | ||
tau_visc = 1.0 | ||
tau_cond = 1.0 | ||
tau_diff = 1.0 | ||
md_layers = 70 | ||
md_tadj_layers = 0 | ||
md_consv_te = 1.0 | ||
md_wait_hr = 0.0 | ||
/ | ||
|
This file was deleted.
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -0,0 +1,43 @@ | ||
############################################################################### | ||
# | ||
# Global control multigases test | ||
# | ||
############################################################################### | ||
|
||
export TEST_DESCR="Compare global control multigases results with previous trunk version" | ||
|
||
export CNTL_DIR=control_wam | ||
|
||
export LIST_FILES="sfcf024.nc \ | ||
atmf024.nc" | ||
|
||
export_fv3 | ||
export NPZ=149 | ||
export NPZP=150 | ||
export DT_ATMOS="225" | ||
export SYEAR="2021" | ||
export SMONTH="03" | ||
export SDAY="22" | ||
export SHOUR="06" | ||
export FHMAX=24 | ||
export WARM_START=.true. | ||
export NA_INIT=0 | ||
export EXTERNAL_IC=.false. | ||
export NGGPS_IC=.false. | ||
export MAKE_NH=.false. | ||
export MOUNTAIN=.true. | ||
export FDIAG=3 | ||
export NSTF_NAME=0,0,1,0,5 | ||
export INPUT_NML=wam.nml.IN | ||
export FV3_RUN=wam_run.IN | ||
export CCPP_SUITE=FV3_GFS_v16_fv3wam | ||
export FIELD_TABLE=field_table_wam | ||
export IAER=5111 | ||
export QUILTING=.true. | ||
export OUTPUT_GRID="gaussian_grid" | ||
export OZ_PHYS_OLD=.false. | ||
export OZ_PHYS_NEW=.true. | ||
export H2O_PHYS=.true. | ||
export NFHOUT=6 | ||
export NFHMAX=24 | ||
export NFHOUT_HF=6 |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
What does a negative
tau
do to the dycore?There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
tau<0 uses the GFS type RF damping, instead of the GFDL type.