From 56a8f01b494433f6c24d7b0cb99549d3bcf29123 Mon Sep 17 00:00:00 2001
From: Nathanael Wong <natgeo.wong@outlook.com>
Date: Sat, 23 Mar 2024 16:09:39 -0400
Subject: [PATCH] Removing 2Doutput, adding new experiments

---
 .../prmtemplates/DGW_P1282km300V64.prm        |   7 +-
 .../prmtemplates/DGW_P1282km300V64_FSF.prm    | 111 ++++++++++++++++++
 .../prmtemplates/DGW_T0642km300V64.prm        |   7 +-
 .../prmtemplates/DGW_T1282km300V64.prm        |   7 +-
 .../prmtemplates/DGW_T1282km300V64_FSF.prm    | 111 ++++++++++++++++++
 .../prmtemplates/SPC_P1282km300V64_FSF.prm    | 111 ++++++++++++++++++
 .../prmtemplates/SPC_T1282km300V64_FSF.prm    | 111 ++++++++++++++++++
 7 files changed, 450 insertions(+), 15 deletions(-)
 create mode 100644 scripts/modifysam/prmtemplates/DGW_P1282km300V64_FSF.prm
 create mode 100644 scripts/modifysam/prmtemplates/DGW_T1282km300V64_FSF.prm
 create mode 100644 scripts/modifysam/prmtemplates/SPC_P1282km300V64_FSF.prm
 create mode 100644 scripts/modifysam/prmtemplates/SPC_T1282km300V64_FSF.prm

diff --git a/scripts/modifysam/prmtemplates/DGW_P1282km300V64.prm b/scripts/modifysam/prmtemplates/DGW_P1282km300V64.prm
index b8f3b82..adca3b3 100644
--- a/scripts/modifysam/prmtemplates/DGW_P1282km300V64.prm
+++ b/scripts/modifysam/prmtemplates/DGW_P1282km300V64.prm
@@ -1,9 +1,6 @@
 &PARAMETERS
 
-caseid ='ExploreWTGSpace-P1282km300V64-member[xx]-OUT2D',
-
-case_restart = 'DGW',
-caseid_restart = 'ExploreWTGSpace-P1282km300V64-member[xx]',
+caseid ='ExploreWTGSpace-P1282km300V64-member[xx]',
 
 dx = 2000.,
 dy = 2000.,
@@ -11,7 +8,7 @@ dt = 10.,
 
 doensemble    = .true.,
 nensemble     = [en],
-nrestart      = 2,
+nrestart      = 0,
 nrestart_skip = 239,
 perturb_type  = 0,
 
diff --git a/scripts/modifysam/prmtemplates/DGW_P1282km300V64_FSF.prm b/scripts/modifysam/prmtemplates/DGW_P1282km300V64_FSF.prm
new file mode 100644
index 0000000..ccc9bdd
--- /dev/null
+++ b/scripts/modifysam/prmtemplates/DGW_P1282km300V64_FSF.prm
@@ -0,0 +1,111 @@
+&PARAMETERS
+
+caseid ='ExploreWTGSpace-P1282km300V64_FSF-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 = .true.,
+SFC_TAU_FXD = .false.,
+fluxt0 = 10.,
+fluxq0 = 90.,
+
+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    = 43200,
+nprint   = 8640,
+nstat    = 360,     ! Calculate statistics every hour!
+nstatfrq = 60,     ! Collect 60 samples per statistical calculation (every 1 minute)
+
+doSAMconditionals     = .false.,
+dosatupdnconditionals = .false.,
+
+nsave2D      = 90,
+nsave2Dstart = 0,
+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  = .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
+
+/
\ No newline at end of file
diff --git a/scripts/modifysam/prmtemplates/DGW_T0642km300V64.prm b/scripts/modifysam/prmtemplates/DGW_T0642km300V64.prm
index d3ccc96..759df12 100644
--- a/scripts/modifysam/prmtemplates/DGW_T0642km300V64.prm
+++ b/scripts/modifysam/prmtemplates/DGW_T0642km300V64.prm
@@ -1,9 +1,6 @@
 &PARAMETERS
 
-caseid ='ExploreWTGSpace-T0642km300V64-member[xx]-OUT2D',
-
-case_restart = 'DGW',
-caseid_restart = 'ExploreWTGSpace-T0642km300V64-member[xx]',
+caseid ='ExploreWTGSpace-T0642km300V64-member[xx]',
 
 dx = 2000.,
 dy = 2000.,
@@ -11,7 +8,7 @@ dt = 10.,
 
 doensemble    = .true.,
 nensemble     = [en],
-nrestart      = 2,
+nrestart      = 0,
 nrestart_skip = 239,
 perturb_type  = 0,
 
diff --git a/scripts/modifysam/prmtemplates/DGW_T1282km300V64.prm b/scripts/modifysam/prmtemplates/DGW_T1282km300V64.prm
index fe1611a..e181974 100644
--- a/scripts/modifysam/prmtemplates/DGW_T1282km300V64.prm
+++ b/scripts/modifysam/prmtemplates/DGW_T1282km300V64.prm
@@ -1,9 +1,6 @@
 &PARAMETERS
 
-caseid ='ExploreWTGSpace-T1282km300V64-member[xx]-OUT2D',
-
-case_restart = 'DGW',
-caseid_restart = 'ExploreWTGSpace-T1282km300V64-member[xx]',
+caseid ='ExploreWTGSpace-T1282km300V64-member[xx]',
 
 dx = 2000.,
 dy = 2000.,
@@ -11,7 +8,7 @@ dt = 10.,
 
 doensemble    = .true.,
 nensemble     = [en],
-nrestart      = 2,
+nrestart      = 0,
 nrestart_skip = 239,
 perturb_type  = 0,
 
diff --git a/scripts/modifysam/prmtemplates/DGW_T1282km300V64_FSF.prm b/scripts/modifysam/prmtemplates/DGW_T1282km300V64_FSF.prm
new file mode 100644
index 0000000..a68e6d9
--- /dev/null
+++ b/scripts/modifysam/prmtemplates/DGW_T1282km300V64_FSF.prm
@@ -0,0 +1,111 @@
+&PARAMETERS
+
+caseid ='ExploreWTGSpace-T1282km300V64_FSF-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 = .true.,
+SFC_TAU_FXD = .false.,
+fluxt0 = 20.,
+fluxq0 = 120.,
+
+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    = 43200,
+nprint   = 8640,
+nstat    = 360,     ! Calculate statistics every hour!
+nstatfrq = 60,     ! Collect 60 samples per statistical calculation (every 1 minute)
+
+doSAMconditionals     = .false.,
+dosatupdnconditionals = .false.,
+
+nsave2D      = 90,
+nsave2Dstart = 0,
+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  = .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
+
+/
\ No newline at end of file
diff --git a/scripts/modifysam/prmtemplates/SPC_P1282km300V64_FSF.prm b/scripts/modifysam/prmtemplates/SPC_P1282km300V64_FSF.prm
new file mode 100644
index 0000000..546dde7
--- /dev/null
+++ b/scripts/modifysam/prmtemplates/SPC_P1282km300V64_FSF.prm
@@ -0,0 +1,111 @@
+&PARAMETERS
+
+caseid ='ExploreWTGSpace-P1282km300V64_FSF-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 = .true.,
+SFC_TAU_FXD = .false.,
+fluxt0 = 10.,
+fluxq0 = 90.,
+
+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_raymondzeng_QJRMS2005   = .false.
+dowtg_hermanraymond_JAMES2014 = .true.
+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             = .true.
+
+wtgscale_vertmodenum = 2      ! number of vertical modes
+wtgscale_vertmodescl = 1., 1. ! wtg strength scaling for vertical modes
+
+/
\ No newline at end of file
diff --git a/scripts/modifysam/prmtemplates/SPC_T1282km300V64_FSF.prm b/scripts/modifysam/prmtemplates/SPC_T1282km300V64_FSF.prm
new file mode 100644
index 0000000..904d4d0
--- /dev/null
+++ b/scripts/modifysam/prmtemplates/SPC_T1282km300V64_FSF.prm
@@ -0,0 +1,111 @@
+&PARAMETERS
+
+caseid ='ExploreWTGSpace-T1282km300V64_FSF-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 = .true.,
+SFC_TAU_FXD = .false.,
+fluxt0 = 20.,
+fluxq0 = 120.,
+
+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_raymondzeng_QJRMS2005   = .false.
+dowtg_hermanraymond_JAMES2014 = .true.
+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             = .true.
+
+wtgscale_vertmodenum = 2      ! number of vertical modes
+wtgscale_vertmodescl = 1., 1. ! wtg strength scaling for vertical modes
+
+/
\ No newline at end of file