Adding diurnal cycle simulations

scripts/modifysam/prmtemplates/KGW/D1282km300V64.prm

@@ -0,0 +1,110 @@
+&PARAMETERS
+
+caseid ='ExploreWTGSpace-D1282km300V64-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     = .false.,
+doseasons       = .false.,
+dosolarconstant = .true.,
+zenith_angle    = 0.0,
+solar_constant  = 1361.,
+
+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.
+
+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
+
+/

scripts/modifysam/prmtemplates/TDG/D1282km300V64.prm

@@ -0,0 +1,111 @@
+&PARAMETERS
+
+caseid ='ExploreWTGSpace-D1282km300V64-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     = .false.,
+doseasons       = .false.,
+dosolarconstant = .true.,
+zenith_angle    = 0.0,
+solar_constant  = 1361.,
+
+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
+
+/