Changeset 3088 for trunk/LMDZ.MARS/deftank

Timestamp:

Oct 16, 2023, 10:31:19 AM (14 months ago)

Author:

jbclement

Message:

PEM:

• Addition of the bash script "concat_diagpem.sh" in deftank to concatenate all the "diagpem.nc" outputs into one file;
• Addition of "writediagpem_mod.F90" inspired by the subroutine "writediagfi.F" in phymars/;
• Renaming of PEM outputs ('startfi_PEM' -> 'startpem' and 'diagfi' -> 'diagpem');
• Gathering of all "diagpem.nc" outputs at the same place at the end of the main PEM loop;
• Update of "launch_pem.sh", "README" and "modify_startfi_orbit.sh" in the deftank.

JBC

Location:

trunk/LMDZ.MARS/deftank/pem

Files:

• README (modified) (1 diff)
• concat_diagpem.sh (added)
• launch_pem.sh (modified) (1 diff)
• modify_startfi_orbit.sh (modified) (2 diffs)

trunk/LMDZ.MARS/deftank/pem/README

@@ -23 +23 @@
   Bash script file to modify orbital parameters in the file "startfi.nc" (obliquity, eccentricity, Ls perihelion and possibly initial Ls). They should be consistent with Laskar data in "obl_ecc_lsp.asc" to start the chained simulation.
 
+# concat_diagpem.sh:
+  Bash script file to concatenate along the variable 'Time' all the "diagfi.nc" files of the PEM into one NetCDF file. 'Time' is reindexed with the numbering of Martian years simulated by the PEM.
+
 Note:
   To run the chained simulation, be careful to have:

trunk/LMDZ.MARS/deftank/pem/launch_pem.sh

@@ -179 +179 @@
     # Copy data files and prepare the next run
     mv out_runPEM${iPEM} out_PEM/run${iPEM}
-    mv diagfi.nc diags/diagfi_PEM${iPEM}.nc
-    cp restartfi_PEM.nc starts/startfi_PEM${iPEM}.nc
-    mv restartfi_PEM.nc startfi_PEM.nc
+    mv diagpem.nc diags/diagpem_${iPEM}.nc
+    cp restartpem.nc starts/startpem_${iPEM}.nc
+    mv restartpem.nc startpem.nc
     cp restartfi_evol.nc starts/startfi_postPEM${iPEM}.nc
     mv restartfi_evol.nc startfi.nc

trunk/LMDZ.MARS/deftank/pem/modify_startfi_orbit.sh

@@ -1 +1 @@
 #!/bin/bash
+######################################################################
+### Script to modify the orbital parameters of a file "startfi.nc" ###
+######################################################################
 
-# Bash script to modify the orbital parameters of a file "startfi.nc"
+#Name of the file to be modified
 name_file="startfi.nc"
 
+# New values for the orbital parameters
 new_obl=25.18941689
 new_ecc=0.09340902
@@ -49 +53 @@
          -s "controle(16)=$peri_day" \
          -s "controle(2)=$new_inisol" \
-         $name_file $name_file.temp
-
-# Rename the temporary file back to the original filename
-mv $name_file.temp $name_file
+         $name_file $name_file
 
 echo "In \"$name_file\":"