Run.def for Venus-DYNAMICO
From Planets
Revision as of 17:37, 18 October 2024 by Aurelien.falco (talk | contribs) (Created page with "Example of run.def for Venus-DYNAMICO (type Held&Suarez): <pre> #---------------- Mesh -...")
Example of run.def for Venus-DYNAMICO (type Held&Suarez):
#---------------- Mesh ---------------- # Number of subdivisions on a main triangle : integer (default=40) nbp = 40 # Number of vertical layers : integer (default=19) llm = 19 # Vertical grid : [std|ncar|ncarl30;strato;strato_custom;ncar;dcmip31;dcmip200;read_apbp;plugin] (default=std) disvert = std # Mesh optimisation : number of iterations : integer (default=0) optim_it = 1000 # Sub splitting of main rhombus : integer (default=1) nsplit_i = 1 nsplit_j = 1 #number of openmp task on vertical level omp_level_size=1 #---------------- Numerics ---------------- # Advection called every itau_adv time steps : integer (default=2) itau_adv = 1 # Time step in s : real (default=480) # dt = 720 # Alternative to specifying "dt", specify number of steps per day : day_step day_step = 240000 # Number of tracers : integer (default=1) nqtot = 1 #---------------- Time and output ---------------- # Time style : [none|dcmip] (default=dcmip) time_style = none # Run length in s : real (default=??) # run_length = 1036800 # Alternative to specifying "run_length", specify number of days to run : ndays ndays=1 # Interval in s between two outputs : integer (default=??) write_period = 314928 #---------------- Planet ---------------- INCLUDEDEF=venus_const.def #---------------- Physical parameters ---------------- # Initial state : # [jablonowsky06|academic|dcmip[1-4]|heldsz|dcmip2_schaer_noshear] (default=jablonowsky06) etat0 = venus # Si on souhaite repartir de la "fin" d'une précédente simulation #etat0=start_file #start_file_name=start # Physics package : [none|held_suarez|dcmip] (default=none) physics = Lebonnois2012 # Dissipation time for grad(div) : real (default=5000) tau_graddiv = 18000 # Exponent of grad(div) disspation : integer (default=1) nitergdiv = 2 # Dissipation time for curl(curl) : real (default=5000) tau_gradrot = 18000 # Exponent of curl(curl) disspation : integer (default=1) nitergrot = 2 # Dissipation time for div(grad) : real (default=5000) tau_divgrad = 18000 # Exponent of div(grad) disspation : integer (default=1) niterdivgrad = 2