Skip to content

Commit

Permalink
Update, split Kuang experiments into time-dependent and non-time-depe…
Browse files Browse the repository at this point in the history
…ndent experiments
  • Loading branch information
natgeo-wong committed Nov 2, 2024
1 parent dd1570f commit 2a3d946
Show file tree
Hide file tree
Showing 5 changed files with 223 additions and 0 deletions.
111 changes: 111 additions & 0 deletions scripts/modifysam/prmtemplates/TDG/P1282km300V64.prm
Original file line number Diff line number Diff line change
@@ -0,0 +1,111 @@
&PARAMETERS

caseid ='ExploreWTGSpace-P1282km300V64-member[xx]',

dx = 2000.,
dy = 2000.,
dt = 10.,

doensemble = .true.,
nensemble = [en],
nrestart = 0,
nrestart_skip = 239,
perturb_type = 0,

LES_S = .false.,
OCEAN = .true.,

dosgs = .true.,
dodamping = .true.,
doupperbound = .true.,
docloud = .true.,
doprecip = .true.,
dolongwave = .true.,
doshortwave = .true.,
dosurface = .true.,
dodynamicocean = .false.,
dolargescale = .true.,
doradforcing = .false.,
doradhomo = .false.,
dosfcforcing = .false.,
docoriolis = .false.,
donudging_uv = .true.,
donudging_tq = .false.,

tauls = 3600.,

SFC_FLX_FXD = .false.,
SFC_TAU_FXD = .false.,

longitude0 = 0.,
latitude0 = 0.,
day0 = 0.,
tabs_s = 300.,
Szero = 0,

depth_slab_ocean = 1.,

nrad = 12,

! note: insolation = solar_constant*cos(zenith_angle)
! below values are equivalent to insolation = 423.4773 W/m2
! coszrs = 0.732650

doperpetual = .true.,
doseasons = .false.,
dosolarconstant = .true.,
zenith_angle = 0.0,
solar_constant = 431.06,

doisccp = .false.,
domodis = .false.,
domisr = .false.,

nstop = 2160000,
nprint = 8640,
nstat = 360, ! Calculate statistics every hour!
nstatfrq = 60, ! Collect 60 samples per statistical calculation (every 1 minute)

doSAMconditionals = .false.,
dosatupdnconditionals = .false.,

nsave2D = 120,
nsave2Dstart = 9999999,
nsave2Dend = 999999999,
save2Davg = .true.,

nsave3D = 120,
nsave3Dstart = 9999999,
nsave3Dend = 999999999,
rad3Dout = .true.,

/

&KUANG_PARAMS

doradtendency = .false.,
dobulksfcflx = .false.,
bulksfcflx_u = 5.

dowtg_blossey_etal_JAMES2009 = .false.
dowtg_kuang_JAS2008 = .true.
dowtg_raymondzeng_QJRMS2005 = .false.
dowtg_hermanraymond_JAMES2014 = .false.
dowtg_decompdgw = .false.
dowtg_decomptgr = .false.
dowtg_timedependence = .true.

wtgscale_time = 0.1

am_wtg = [am] ! momentum damping rate in day^-1, default= 1.,
am_wtg_exp = 0., ! exponent of p/p0 in momentum damping rate.
lambda_wtg = 650.e3, ! quarter wavelength in m. default = 650.e3 m

tau_wtg = [tau] ! temperature relaxation in hours, default= 1.,
dthetadz_min = 1.e-3
dowtgLBL = .false.

wtgscale_vertmodenum = 2 ! number of vertical modes
wtgscale_vertmodescl = 1., 1. ! wtg strength scaling for vertical modes

/
111 changes: 111 additions & 0 deletions scripts/modifysam/prmtemplates/TDG/T1282km300V64.prm
Original file line number Diff line number Diff line change
@@ -0,0 +1,111 @@
&PARAMETERS

caseid ='ExploreWTGSpace-T1282km300V64-member[xx]',

dx = 2000.,
dy = 2000.,
dt = 10.,

doensemble = .true.,
nensemble = [en],
nrestart = 0,
nrestart_skip = 239,
perturb_type = 0,

LES_S = .false.,
OCEAN = .true.,

dosgs = .true.,
dodamping = .true.,
doupperbound = .true.,
docloud = .true.,
doprecip = .true.,
dolongwave = .false.,
doshortwave = .false.,
dosurface = .true.,
dodynamicocean = .false.,
dolargescale = .true.,
doradforcing = .false.,
doradhomo = .false.,
dosfcforcing = .false.,
docoriolis = .false.,
donudging_uv = .true.,
donudging_tq = .false.,

tauls = 3600.,

SFC_FLX_FXD = .false.,
SFC_TAU_FXD = .false.,

longitude0 = 0.,
latitude0 = 0.,
day0 = 0.,
tabs_s = 300.,
Szero = 0,

depth_slab_ocean = 1.,

nrad = 12,

! note: insolation = solar_constant*cos(zenith_angle)
! below values are equivalent to insolation = 423.4773 W/m2
! coszrs = 0.732650

doperpetual = .true.,
doseasons = .false.,
dosolarconstant = .true.,
zenith_angle = 0.0,
solar_constant = 431.06,

doisccp = .false.,
domodis = .false.,
domisr = .false.,

nstop = 2160000,
nprint = 8640,
nstat = 360, ! Calculate statistics every hour!
nstatfrq = 60, ! Collect 60 samples per statistical calculation (every 1 minute)

doSAMconditionals = .false.,
dosatupdnconditionals = .false.,

nsave2D = 120,
nsave2Dstart = 9999999,
nsave2Dend = 999999999,
save2Davg = .true.,

nsave3D = 120,
nsave3Dstart = 9999999,
nsave3Dend = 999999999,
rad3Dout = .true.,

/

&KUANG_PARAMS

doradtendency = .true.,
dobulksfcflx = .false.,
bulksfcflx_u = 5.

dowtg_blossey_etal_JAMES2009 = .false.
dowtg_kuang_JAS2008 = .true.
dowtg_raymondzeng_QJRMS2005 = .false.
dowtg_hermanraymond_JAMES2014 = .false.
dowtg_decompdgw = .false.
dowtg_decomptgr = .false.
dowtg_timedependence = .true.

wtgscale_time = 0.1

am_wtg = [am] ! momentum damping rate in day^-1, default= 1.,
am_wtg_exp = 0., ! exponent of p/p0 in momentum damping rate.
lambda_wtg = 650.e3, ! quarter wavelength in m. default = 650.e3 m

tau_wtg = [tau] ! temperature relaxation in hours, default= 1.,
dthetadz_min = 1.e-3
dowtgLBL = .false.

wtgscale_vertmodenum = 2 ! number of vertical modes
wtgscale_vertmodescl = 1., 1. ! wtg strength scaling for vertical modes

/
Original file line number Diff line number Diff line change
Expand Up @@ -93,6 +93,7 @@ dowtg_raymondzeng_QJRMS2005 = .false.
dowtg_hermanraymond_JAMES2014 = .false.
dowtg_decompdgw = .false.
dowtg_decomptgr = .false.
dowtg_timedependence = .true.

wtgscale_time = 0.1

Expand Down

0 comments on commit 2a3d946

Please sign in to comment.