Run.def for Pluto-DYNAMICO - Revision history

Emillour at 16:00, 11 April 2025

2025-04-11T16:00:09Z

Aurelien.falco at 14:22, 26 November 2024

2024-11-26T14:22:42Z

Aurelien.falco: run.def without a startfi

2024-11-26T14:10:34Z

run.def without a startfi

New page

#---------------- 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] (default=std)
# disvert = std
disvert = plugin
hybrid = .false.

# Mesh optimisation : number of iterations : integer (default=0)
optim_it = 1000

# Sub splitting of main rhombus : integer (default=1)
nsplit_i = 10
nsplit_j = 10

#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 = 12000

# Number of tracers : integer (default=1)
nqtot = 7

#---------------- Time and output ----------------

# Time style : [none|dcmip] (default=dcmip)
time_style = none

# Run length in s : real (default=??)
run_length = 120
# Alternative to specifying "run_length", specify number of days to run : ndays
# ndays=1

# Interval in s between two outputs : integer (default=??)
write_period = 120

#---------------- Planet ----------------

INCLUDEDEF=pluto_const.def

#---------------- Physical parameters ----------------

# Initial state :
# [jablonowsky06|academic|dcmip[1-4]|heldsz|dcmip2_schaer_noshear] (default=jablonowsky06)
etat0 = isothermal
etat0_isothermal_temp=50
etat0_ps_white_noise=0.01
etat0_theta_rhodz_white_noise=0.01

# Physics package : [none|held_suarez|dcmip] (default=none)
physics = phys_external

startphy_file = false

# 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

INCLUDEDEF=callphys.def

← Older revision		Revision as of 14:22, 26 November 2024
Line 1:		Line 1:
		+	Example of a run.def to run without a startfi input file:
		+
		+	<pre>
	#---------------- Mesh ----------------		#---------------- Mesh ----------------

Line 87:		Line 90:

	INCLUDEDEF=callphys.def		INCLUDEDEF=callphys.def
		+	</pre>

@@ Line 1: / Line 1: @@
-Example of a run.def to run without a startfi input file:
+An example of a run_icosa.def to use for a Pluto-DYNAMICO run is given in the LMDZ.PLUTO/deftank/dynamico subdirectory (see e.g. https://trac.lmd.jussieu.fr/Planeto/browser/trunk/LMDZ.PLUTO/deftank/dynamico/run_icosa.def ).
+It should look something like:
 <pre>
 #---------------- Mesh ----------------
 # Number of subdivisions on a main triangle : integer (default=40)
-nbp = 40
+nbp = 20
 # Number of vertical layers : integer (default=19)
-llm = 19
+llm = 27
 # Vertical grid : [std|ncar|ncarl30] (default=std)
@@ Line 19: / Line 19: @@
 # Sub splitting of main rhombus : integer (default=1)
-nsplit_i = 10
+# NB: total number of computational subdomains is 10*nsplit_i*nsplit_j
-nsplit_j = 10
+nsplit_i = 1
 #number of openmp task on vertical level
@@ Line 33: / Line 34: @@
 # dt = 720
 # Alternative to specifying "dt", specify number of steps per day : day_step
-day_step = 12000
+day_step = 2400
 # Number of tracers : integer (default=1)
@@ Line 44: / Line 45: @@
 # Run length in s : real (default=??)
-run_length = 120
+# run_length = 100
 # Alternative to specifying "run_length", specify number of days to run : ndays
-# ndays=1
+ndays=1
 # Interval in s between two outputs : integer (default=??)
-write_period = 120
+write_period = 1000
@@ Line 61: / Line 62: @@
 # Initial state :
 #   [jablonowsky06|academic|dcmip[1-4]|heldsz|dcmip2_schaer_noshear] (default=jablonowsky06)
-etat0 = isothermal
+#etat0=isothermal
-etat0_isothermal_temp=50
+#etat0_isothermal_temp=80
-etat0_ps_white_noise=0.01
+etat0 = start_file
-etat0_theta_rhodz_white_noise=0.01
+etat0_start_file_colocated=true
 # Physics package : [none|held_suarez|dcmip] (default=none)
 physics = phys_external
-startphy_file = false
+# Call physics every itau_physics dynamical steps
 # Dissipation time for grad(div) : real (default=5000)
@@ Line 75: / Line 94: @@
 # Exponent of grad(div) disspation : integer (default=1)
-nitergdiv = 2
+nitergdiv = 1
 # Dissipation time for curl(curl) : real (default=5000)
@@ Line 89: / Line 108: @@
 niterdivgrad = 2
 </pre>