The run icosa.def Input File

1. 1. mpi_threading_mode=serialized
2. mpi_threading_mode=funneled

1. halo_i=0
2. halo_j=1

1. -------------------------------------------------------------------------------------
2. --------------------------------------- Mesh ----------------------------------------
3. -------------------------------------------------------------------------------------

1. Number of subdivision on a main triangle (nbp) : integer (default=40)

nbp = 160

1. nbp = 40

1. nbp 20 40 80 160
2. T-edge length (km) 500 250 120 60

1. 1. sub splitting of main rhombus : integer (default=1)
2. nsplit_i=1
3. nsplit_j=1
4. omp_level_size=1
   2. There must be less MPIxOpenMP processes than the 10 x nsplit_i x nsplit_j tiles
   3. typically for pure MPI runs, let nproc = 10 x nsplit_i x nsplit_j
   4. it is better to have nbp/split >~ 10
   5. 2 procs
5. nsplit_i = 1
6. nsplit_j = 1
   1. 40 procs
7. nsplit_i = 2
8. nsplit_j = 2
   1. 1. 1. 40 noeuds de 24 processeurs = 960 procs
9. nsplit_i = 8
10. nsplit_j = 12
    1. 1. 1. 30 noeuds de 24 processeurs = 720 procs
11. nsplit_i = 8
12. nsplit_j = 9
    1. 1. 1. 20 noeuds de 24 processeurs = 480 procs

nsplit_i = 8 nsplit_j = 6

1. Number of vertical layer (llm) : integer (default=19)

llm = 64

1. disvert : vertical discretisation : string (default='std') : std, ncar, ncar30l

disvert=read_apbp

1. optim_it : mesh optimisation : number of iteration : integer (default=0)
2. optim_it = 1000

optim_it = 0

1. -------------------------------------------------------------------------------------
2. --------------------------------------- Time ----------------------------------------
3. -------------------------------------------------------------------------------------

1. 1. scheme type : string ( default='runge_kutta') euler, leapfrog_matsuno, runge_kutta )
2. scheme = runge_kutta leapfrog_matsuno

1. 1. matsuno period : integer ( default=5)
2. matsuno_period = 5

1. number of dynamical time step per day

day_step = 80

1. advection called every itau_adv time steps : integer (default 2)
2. standard : umax=100m/s vs c=340m/s (ratio 1:3)
3. in JW06 umax=35m/s vs c=340m/s (ratio 1:10)

itau_adv=3

1. number of days to run (substitute run_length) !!!!

ndays = 5

1. 1. activate IO (default = true)
2. enable_io = false

1. 1. output with XIOS (only if compiled with XIOS): true/false (default true)
2. xios_output=false

1. output field period (only when not using XIOS) : integer (default none)
2. write_period=7200
3. write_period=14400
4. write_period=118.9125

write_period = 446.75

1. itau_write_etat0=380520

1. -------------------------------------------------------------------------------------
2. -------------------------------------- Misc -----------------------------------------
3. -------------------------------------------------------------------------------------

1. number of tracer (nqtot) : integer (default 1)

nqtot=2

1. pression value where output is interpolated : real (default=0, no output)

out_pression_level=85000

1. etat0 : initial state : string (default=jablonowsky06) :
2. jablonowsky06, academic, ncar

etat0=start_file etat0_start_file_colocated=false

1. 1. -- to cross the 2y limit
2. etat0_start_iteration=0
3. run_length=38052

1. etat0=isothermal
2. etat0_isothermal_temp=175

1. start file name (default is start.nc)
2. restart file name (default is restart.nc)

start_file_name=start_icosa restart_file_name=restart_icosa

1. 1. 2. to run a one-day simulation
      3. using a ASCII profile
      4. to create start files
2. etat0=temperature_profile
3. temperature_profile_file=temp_profile.txt

1. -------------------------------------------------------------------------------------
2. ----------------------------------- Dynamics ----------------------------------------
3. -------------------------------------------------------------------------------------

1. caldyn : computation type for gcm equation : string (default=gcm) : gcm, adv

caldyn=gcm

1. caldyn_conserv : string (default=energy) : energy,enstrophy

caldyn_conserv=energy

1. caldyn_exner : scheme for computing Exner function : string (default=direct) : direct,lmdz

caldyn_exner=direct

1. caldyn_hydrostat : scheme for computing geopotential : string (default=direct) : direct,lmdz

caldyn_hydrostat=direct

1. guided_type : string (default=none) : none, ncar

guided_type=none

1. Sponge layer
   1. 1. iflag_sponge=0 for no sponge
      2. iflag_sponge=1 for sponge over 4 topmost layers
      3. iflag_sponge=2 for sponge from top to ~1% of top layer pressure
      4. tau_sponge --> damping frequency at last layer
      5. e-5 medium / e-4 strong yet reasonable / e-3 very strong
      6. mode_sponge=1 for u,v --> 0
      7. mode_sponge=2 for u,v --> zonal mean (NOT IMPLEMENTED)
      8. mode_sponge=3 for u,v,h --> zonal mean (NOT IMPLEMENTED)
2. iflag_sponge = 1

iflag_sponge = 0

1. tau_sponge = 1.e-4
2. mode_sponge = 1

1. -------------------------------------------------------------------------------------
2. ---------------------------------- Dissipation --------------------------------------
3. -------------------------------------------------------------------------------------

1. dissipation time graddiv : real (default=5000)

tau_graddiv = 10000

1. number of iteration for graddiv : integer (default=1)

nitergdiv = 2

1. dissipation time nxgradrot (default=5000)

tau_gradrot = 10000

1. number of iteration for nxgradrot : integer (default=1)

nitergrot=2

1. dissipation time divgrad (theta) (default=5000)

tau_divgrad= 10000

1. number of iteration for divgrad : integer (default=1)

niterdivgrad=2

1. Rayleigh friction : string (default=none) :
2. none, dcmip2_schaer_noshear, dcmip2_schaer_shear, giant_liu_schneider

rayleigh_friction_type = giant_liu_schneider rayleigh_limlat = 16. rayleigh_friction_tau = 8640000.

1. -------------------------------------------------------------------------------------
2. ------------------------------------- Physics ---------------------------------------
3. -------------------------------------------------------------------------------------

1. (itau_physics=160)*(dt=111.6875) --> half a jupiter day as physical timestep
2. ---- for some reason (day consistency, physical timestep must be less a day)
3. ---- change also in run.def
4. itau_physics = 160

itau_physics = 40

1. kind of physics : string : none, dcmip (default=none)

physics=phys_external

1. 1. testcase physics for dcmip : INTEGER : 1, 2 (default=1)
2. dcmip_physics=1

1. ---------------------- parameters for NCAR test cases -------------------------

1. NCAR advection test, initial tracer : string ( default='cos_bell')
2. const, slotted_cyl, cos_bell, dbl_cos_bell_q1, dbl_cos_bell_q2, complement, hadley

ncar_adv_shape=cos_bell

1. NCAR advection test, wind field : string (default='deform') : solid, deform, hadley

ncar_adv_wind=solid

1. ncar_T0 : reference temperature for NCAR test cases : real (default=300)
2. also used by disvert if disvert=ncar

ncar_T0=300

1. ncar_p0 : reference pressure for NCAR test cases : real (default=1e5)
2. also used by disvert if disvert=ncar

ncar_p0=1e5

1. ncar_disvert_c : exponent for B(eta) : integer (default=1)
2. used by disvert if disvert=ncar

ncar_disvert_c=1

1. dcmip 4 testcase : integer (default=1) : 1, 2

dcmip4_testcase=1

1. -------------------------------------------------------------------------------------
2. -------------------------------------------------------------------------------------