
&SGS_TKE
dosmagor = .false.
/


 &MICRO_DRIZZLE

 Nc0 = 100.
 sigmag = 1.2

/

 &MICRO_M2005

 Nc0 = 100.,
 doicemicro = .false.

/

 &MICRO_M2005_PA

 Nc0 = 100.,
 dofix_pgam = .true.
 pgam_fixed = 10.3 ! Geoffroy et al (2010, doi:10.5194/acp-10-4835-2010)

 doicemicro = .false.
 dototalwater = .true.

 dopredictNc = .true.
 dospecifyaerosol = .true.
 doprogaerosol = .true.
 aerinitmode = 0
 doevapnc = .true.
 docloudedgeactivation = .true.
 douse_reffc = .true.
 massfactor = 1.

 do_m2011_scavenge = .true.
 nscv = 3

 aer_rm1 = 0.127
 aer_sig1 = 1.5
 aer_n1 = 33.  !below inversion, #/mg

 aer_n2 = 100. !above inversion, #/mg
 aer_rm2 = 0.102
 aer_sig2 = 1.5

 aerfluxmode = 2
 sflux_nacc_coef = 44440000.
 sflux_racc = 0.13

 doshipv2 = .true.
 shipv2_radius = 100.
 shipv2_aer_n = 1.e16 ! total aerosol input in #/s
 shipv2_aer_r = 100.e-9
 shipv2_time0 = 196.75004, !-99999.
 shipv2_timef = 196.78695, !-99999.

 do_output_micro_process_rates = .true.

 do_aerosol_chunk_budgets = .true.

/

 &MICRO_THOM

 Nc0 = 100.,
 doicemicro = .false.

/

 &PARAMETERS

 caseid ='3072x128x288_100m_ship_effrad',

 nrestart = 2,

 caseid_restart = 'spinup_3072x128x288_100m_nudge',
 !caseid_restart = '3072x128x288_100m_ship_control',
 case_restart = 'SCT_jeyun',
 nrestart_skip = 3, 
 restart_sep = .true.,

 CEM = .true.,

 OCEAN = .true.,

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

 ! apply temperature nudging starting 500m above inversion
 dovariable_tauz = .true.
 variable_tauz_offset_above_inversion = 500.
 variable_tauz_minimum_height = 1200.
 tautqls = 3600.,

 tauls = 7200.,

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

 latitude0 = 34.0,
 longitude0 = -25.,

 ug = 0.,
 vg = -4.,

 dx = 	100.,
 dy = 	100.,

 dt = 3.0,

 nxco2 = 1.

 doperpetual = .false.
 dosolarconstant = .true.
 solar_constant = 593.6336
 zenith_angle = 38.7340
 nrad = 5,
 day0= 196.,

nstop 	= 30000,
nprint  = 300,
nstat 	= 300,
nstatfrq= 60,

 nsave2D	= 600,
 nsave2Dstart	= 0,
 nsave2Dend	= 99997200,
 save2Dbin      = .true.
 save2Davg = .false.

 nsave3D	= 1200,
 nsave3Dstart	= 0,
 nsave3Dend	= 999999999,
 save3Dbin      = .true.

 doSAMconditionals = .true.
 dosatupdnconditionals = .true.
 doShipTrackconditionals = .true.

 doisccp = .true.
 domodis = .true.
 domisr = .true.

 doSCAMiopdata = .true.
 dozero_out_day0 = .false.
 iopfile = 'composite_ref_CTL_Aerosol_BL33_FT100_D200um_q4_WithUgVg_rot.nc'
 use_scam_reference_sounding_for_wtg=.false.
 use_scam_reference_sounding=.true.,

 perturb_type = 6

 ncycle_min = 2,
 ncycle_max = 40,
 cfl_safety_factor = 1.5

 do_chunked_energy_budgets = .true.
 nchunk_x_gl = 192
 nchunk_y_gl = 16
 nsaveMSE = 1200

/ 

 &UWOPTIONS

 nelapsemin = 2800,

 donudging_transient = .true.
 transient_nudging_start = 194.75 ! day
 transient_nudging_end = 196.75 ! day
 transient_nudging_ramp = 0.0625 ! day
 tau_transient_nudging = 600. ! second

 dowtg_blossey_etal_JAMES2009 = .false.
! am_wtg = 1.,
! am_wtg_exp = 1.,
! lambda_wtg = 200.e3, ! This is stronger than default to strenthen "inversion nudging"
! wtg_end_time = 196.75 ! when SST starts increasing (day)
! wtg_end_ramp = 0.0625 ! time over which WTG falls off (day)

! do_smooth_reference_density = .true.
! lz_smooth = 600.

/
