Changeset 3720

Timestamp:

Apr 11, 2025, 5:26:37 PM (2 months ago)

Author:

emillour

Message:

Pluto PCM:
Some update/cleanup of the "reference" files in deftank/dynamico
and of reference constants at the interface.
EM

Location:

trunk

Files:

• ICOSA_LMDZ/src/phypluto/interface_icosa_lmdz.f90 (modified) (1 diff)
• LMDZ.PLUTO/changelog.txt (modified) (1 diff)
• LMDZ.PLUTO/deftank/dynamico/README (added)
• LMDZ.PLUTO/deftank/dynamico/iodef.xml (modified) (2 diffs)
• LMDZ.PLUTO/deftank/dynamico/nudging_dynamico.xml (deleted)
• LMDZ.PLUTO/deftank/dynamico/pluto_const.def (modified) (1 diff)
• LMDZ.PLUTO/deftank/dynamico/run.def.newstart (deleted)
• LMDZ.PLUTO/deftank/dynamico/run_icosa.def (modified) (7 diffs)
• LMDZ.PLUTO/deftank/dynamico/traceur.def (added)

trunk/ICOSA_LMDZ/src/phypluto/interface_icosa_lmdz.f90

@@ -362 +362 @@
     ENDIF
 
-    day_length=86400
+    day_length=551854.0625 ! Pluto day length (s)
     CALL getin('day_length',day_length)
-    year_length = 31557600 ! 365.25 * 86400
+    ! Pluto year length (s) 14178.3037109375 Pluto days
+    year_length = 7824354502.23969
     CALL getin('year_length',year_length)
     ndays=nint(itaumax*(dt/day_length))! number of days to run

trunk/LMDZ.PLUTO/changelog.txt

@@ -1866 +1866 @@
 Follow-up of reindexing of columns for restartfi files: In 1D one should also
 initialize the "ind_cell_glo" array.
+
+== 11/04/2025 == EM
+Some update/cleanup of the "reference" files in deftank/dynamico
+and of reference constants at the interface.

trunk/LMDZ.PLUTO/deftank/dynamico/iodef.xml

@@ -27 +27 @@
   <!-- <context id="icosagcm"        src="./tracers_dynamico.xml" /> -->
 
-  <!-- For running with LMDZ, following line is added by dynamicolmdz.driver while using libIGCM:   -->
+  <!-- For running with physics:   -->
   <context id="LMDZ" src="./context_pcm_physics.xml" />
 
@@ -33 +33 @@
 </simulation>
 
-<!-- This file was generated by : python iodef.xml.py -->
-<!-- For help : python iodef.xml.py -h -->
-

trunk/LMDZ.PLUTO/deftank/dynamico/pluto_const.def

@@ -33 +33 @@
 # Reference pressure : real (default=1.1)
 preff = 1.1
+
+# Reference "pa" for hybrid vertical coordinate (Pa)
+pa = 0.2
+
+# Atmospheric Scale height (m) # should match the one in z2sig.def !!
+scale_height=18000.

trunk/LMDZ.PLUTO/deftank/dynamico/run_icosa.def

@@ -2 +2 @@
 
 # Number of subdivisions on a main triangle : integer (default=40)
-nbp = 40
+nbp = 20
 
 # Number of vertical layers : integer (default=19)
@@ -16 +16 @@
 
 # Sub splitting of main rhombus : integer (default=1)
-nsplit_i = 10
-nsplit_j = 10
+# NB: total number of computational subdomains is 10*nsplit_i*nsplit_j
+nsplit_i = 1
+nsplit_j = 1
 
 #number of openmp task on vertical level
@@ -30 +31 @@
 # dt = 720
 # Alternative to specifying "dt", specify number of steps per day : day_step
-day_step = 8000
+day_step = 2400
 
 # Number of tracers : integer (default=1)
@@ -58 +59 @@
 # Initial state :
 #   [jablonowsky06|academic|dcmip[1-4]|heldsz|dcmip2_schaer_noshear] (default=jablonowsky06)
+#etat0=isothermal
+#etat0_isothermal_temp=80
 etat0 = start_file
 etat0_start_file_colocated=true
+# start file name (default: start)
+start_file_name = start_icosa
+
+# restart file name (default: restart)
+restart_file_name = restart_icosa
+
+# optional perturbations to add to initial state:
+#etat0_ps_white_noise=0.01
+#etat0_theta_rhodz_white_noise=0.01
+#etat0_u_white_noise=0.01
 
 # Physics package : [none|held_suarez|dcmip] (default=none)
@@ -65 +78 @@
 # physics = none
 physics = phys_external
+
+# Call physics every itau_physics dynamical steps
+itau_physics=5
+# try to be consistant with LMDZ 
+iphysiq=5
 
 startphy_file = true
@@ -73 +91 @@
 
 # Exponent of grad(div) disspation : integer (default=1)
-nitergdiv = 2
+nitergdiv = 1
 
 # Dissipation time for curl(curl) : real (default=5000)
@@ -86 +104 @@
 # Exponent of div(grad) disspation : integer (default=1)
 niterdivgrad = 2
-