Changeset 3349 for trunk/LMDZ.COMMON/libf

Timestamp:

May 31, 2024, 10:46:55 AM (7 months ago)

Author:

jbclement

Message:

PEM:
New script called "launchPEM.sh" to launch a chained simulation of PCM/PEM runs:

• It is far more efficient since it submits a cycle of PCM/PEM dependent jobs waiting in the SLURM queue. Plus, it does not involve a endless running master script;
• It is more robust and user-friendly since it holds lots of warning/error messages handling different situations. It can take arguments to adapt to the needs of the user. It also creates automotically a script called "kill_launchPEM.sh" to cancel all the current dependent jobs.

There are 3 other files to make it work:

• "lib_launchPEM.sh" which is the library containing the bash functions used by the script;
• "jobPCM.slurm" which is the SLURM job script for the PCM used as a template by the script;
• "jobPEM.slurm" which is the SLURM job script for the PEM used as a template by the script.

More info in the "README" of the deftank.
JBC

Location:

trunk/LMDZ.COMMON/libf/evolution

Files:

• changelog.txt (modified) (1 diff)
• conf_pem.F90 (modified) (2 diffs)
• deftank/README (modified) (2 diffs)
• deftank/exePCM.sh (deleted)
• deftank/jobPCM.slurm (added)
• deftank/jobPEM.slurm (added)
• deftank/launchPEM.sh (added)
• deftank/launch_pem.sh (deleted)
• deftank/lib_launchPEM.sh (added)
• info_PEM_mod.F90 (modified) (3 diffs)

trunk/LMDZ.COMMON/libf/evolution/changelog.txt

@@ -327 +327 @@
 == 27/05/2024 == LL
 Fixing bug in update_soilthermal properties: in case of an ice table, the dry regolith was wrongly initialized.
+
+== 31/05/2024 == JBC
+New script called "launchPEM.sh" to launch a chained simulation of PCM/PEM runs:
+  - It is far more efficient since it submits a cycle of PCM/PEM dependent jobs waiting in the SLURM queue. Plus, it does not involve a endless running master script;
+  - It is more robust and user-friendly since it holds lots of warning/error messages handling different situations. It can take arguments to adapt to the needs of the user. It also creates automotically a script called "kill_launchPEM.sh" to cancel all the current dependent jobs.
+There are 3 other files to make it work:
+  - "lib_launchPEM.sh" which is the library containing the bash functions used by the script;
+  - "jobPCM.slurm" which is the SLURM job script for the PCM used as a template by the script;
+  - "jobPEM.slurm" which is the SLURM job script for the PEM used as a template by the script.
+More info in the "README" of the deftank.

trunk/LMDZ.COMMON/libf/evolution/conf_pem.F90

@@ -30 +30 @@
 use ice_table_mod,         only: icetable_equilibrium, icetable_dynamic
 use layering_mod,          only: layering_algo
+use info_PEM_mod,          only: iPCM, iPEM, nPCM, nPCM_ini
 
 implicit none
@@ -48 +49 @@
     stop
 else
-    read(100,*) i_myear, n_myears, convert_years
+    read(100,*) i_myear, n_myears, convert_years, iPCM, iPEM, nPCM, nPCM_ini
 endif
 close(100)

trunk/LMDZ.COMMON/libf/evolution/deftank/README

@@ -1 +1 @@
 This folder contains all the files needed to run the Planetary Evolution Model (PEM).
 
-# launch_pem.sh:
+# launchPEM.sh:
   Bash script file to launch the chained simulation of PEM and PCM runs.
   The user has to specify:
       (i)   n_mars_years, n_earth_years -> the number of Mars/Earth years to be simulated in total;
-      (ii)  nPCM -> the number of PCM calls between PEM runs (usually nPCM=2);
-      (iii) exePEM -> the name of executable file for the PEM run;
-      (iv)  exeReshape -> the name of executable file for the data reshaping task with XIOS.
+      (ii)  nPCM_ini -> the number of initial PCM runs (at least 2);
+      (iii) nPCM -> the number of PCM runs between each PEM run (usually 2).
+  The script can take an argument:
+      - If there is no argument, then the script initiates a PEM simulation from scratch.
+      - If the argument is 're', then the script relaunches an existing PEM simulation. It will ask for parameters to know the starting point that you want to. 
+
+# liblaunchPEM.sh:
+  Library of bash functions used by the PEM launching script "launchPEM.sh".
+
+# jobPCM.slurm:
+  Bash script file to run the PCM with a SLURM job. The name of the PCM executable file should be adapted. The header corresponds to the ADASTRA supercomputer and should be changed for other machines.
   The path to source the arch file should be adapted to the machine.
 
-# exePCM.sh:
-  Bash script file to run the PCM with a SLURM job. The name of the PCM executable file should be adapted. The header corresponds to the ADASTRA supercomputer and should be changed for other machines.
+# jobPEM.slurm:
+  Bash script file to run the PEM with a SLURM job. The name of the PEM executable file and Reshaping executable file should be adapted. The header corresponds to the ADASTRA supercomputer and should be changed for other machines.
+  The path to source the arch file should be adapted to the machine.
 
 # run_PEM.def
@@ -23 +32 @@
   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's data in "obl_ecc_lsp.asc" to start a chained simulation at a specific date.
 
+# modify_startfi_var.sh:
+  Bash script file to modify the value of a variable in a file "startfi.nc".
+
+# inipem_orbit.sh:
+  Bash script file to modify the orbital parameters of a file "startfi.nc" according to the date set in the file "run_PEM.def".
+
 # concat_diagpem.sh:
   Bash script file to concatenate along the variable 'Time' all the "diagpem.nc" files of the PEM into one NetCDF file. 'Time' is re-indexed with the numbering of Martian years simulated by the PEM run.
+
+# output_layering.py:
+  Python script file to output the stratification data from the "startpem.nc" files.
 
 Note:

trunk/LMDZ.COMMON/libf/evolution/info_PEM_mod.F90

@@ -2 +2 @@
 
 implicit none
+
+integer :: iPCM, iPEM, nPCM, nPCM_ini ! Data about the chained simulation of PCM/PEM runs
 
 !=======================================================================
@@ -28 +30 @@
 logical       :: ok
 integer       :: cstat
-character(10) :: ich1, ich2, fch
+character(10) :: ich1, ich2, ich3, ich4, ich5, ich6, fch
 
 !----- Code
@@ -36 +38 @@
     write(ich2,'(i0)') n_myear
     write(fch,'(f0.4)') convert_years ! 4 digits to the right of the decimal point to respect the precision of Martian year in "launch_pem.sh"
-    call execute_command_line('sed -i "1s/.*/'//trim(ich1)//' '//trim(ich2)//' '//trim(fch)//'/" info_PEM.txt',cmdstat = cstat)
+    write(ich3,'(i0)') iPCM
+    write(ich4,'(i0)') iPEM + 1
+    write(ich5,'(i0)') nPCM
+    write(ich6,'(i0)') nPCM_ini
+    call execute_command_line('sed -i "1s/.*/'//trim(ich1)//' '//trim(ich2)//' '//trim(fch)//' '//trim(ich3)//' '//trim(ich4)//' '//trim(ich5)//' '//trim(ich6)//'/" info_PEM.txt',cmdstat = cstat)
     if (cstat > 0) then
         error stop 'info_PEM: command execution failed!'