Changeset 3145 for trunk/LMDZ.MARS/deftank/pem

Timestamp:

Dec 1, 2023, 1:48:25 PM (13 months ago)

Author:

jbclement

Message:

PEM:

• Addition of a script in LMDZ.MARS/deftank/pem/ to launch a chained simulation of 1D PCM runs which follow, year by year, the orbital parameters (obliquity, eccentricity, Ls perihelion) given in a specified file.
• Small changes to other files of the deftank directory (check and cosmetic).

JBC

Location:

trunk/LMDZ.MARS/deftank/pem

Files:

• concat_diagpem.sh (modified) (1 diff)
• launch_orb_1Dchained.sh (added)
• launch_pem.sh (modified) (3 diffs)
• modify_startfi_orbit.sh (modified) (3 diffs)

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

@@ -19 +19 @@
 
 # Checking if everything is ok
+if [ ! -d "$directory" ]; then
+    echo "Error: directory \"$diags\" not found!"
+    exit 1
+fi
 if [ ${#files[@]} -lt 2 ]; then
-  echo "Error: there are not enough files for concatenation!"
-  exit 1
+    echo "Error: there are not enough files for concatenation!"
+    exit 1
 fi
 if [ -f "$output_file" ]; then

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

@@ -38 +38 @@
 exeReshape=reshape_XIOS_output_64x48x29_phymars_seq.e
 
-#------------------------------ Parameters -----------------------------
-myear=686.9725      # Number of Earth days in Martian year
-eyear=365.256363004 # Number of days in Earth year
-
 #------ Check if files/directories necessary for the script exist ------
 dir=`pwd`
 machine=`hostname`
 address=`whoami`
+if [ ! -f "$exePEM" ]; then
+    echo "Error: file \"$exePEM\" does not exist in $dir!"
+    exit 1
+fi
+if [ ! -f "$exeReshape" ]; then
+    echo "Error: file \"$exeReshape\" does not exist in $dir!"
+    exit 1
+fi
 if [ ! -f "exePCM.sh" ]; then
     echo "Error: file \"exePCM.sh\" does not exist in $dir!"
@@ -92 +96 @@
 address=`whoami`
 echo "This is a chained simulation for PEM and PCM runs in $dir on $machine."
+myear=686.9725      # Number of Earth days in Martian year
+eyear=365.256363004 # Number of days in Earth year
 convert_years=$(echo "$myear/$eyear" | bc -l)
 convert_years=$(printf "%.4f" $convert_years) # Rounding to the 4th decimal to respect the precision of Martian year
@@ -110 +116 @@
 if [ -f "start.nc" ]; then
     cp start.nc starts/
+elif [ -f "star1D.nc" ]; then
+    cp star1D.txt starts/
 fi
 

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

@@ -15 +15 @@
 
 # Calculate modified values
+if [ ! -f "$name_file" ]; then
+    echo "Error: file \"$name_file\" not found!"
+    exit 1
+fi
 year_day=669.
 halfaxe=227.94
@@ -26 +30 @@
 zx0=$(echo "-2.*a($tan*sqrt((1. - $new_ecc)/(1. + $new_ecc)))" | bc -l)
 if [ $(echo "$zx0 <= 0." | bc -l) -eq 1 ]; then
-  zx0=$(echo "$zx0 + 2.*$pi" | bc -l)
+    zx0=$(echo "$zx0 + 2.*$pi" | bc -l)
 fi
 peri_day=$(echo "$year_day*(1. - ($zx0 - $new_ecc*s($zx0))/(2.*$pi))" | bc -l)
@@ -36 +40 @@
 new_inisol=$(echo "$peri_day + $xref*$year_day/(2.*$pi)" | bc -l)
 if [ $(echo "$new_inisol < 0." | bc -l) -eq 1 ]; then
-  new_inisol=$(echo "$new_inisol + $year_day" | bc -l)
+    new_inisol=$(echo "$new_inisol + $year_day" | bc -l)
 fi
 if [ $(echo "$new_inisol >= $year_day" | bc -l) -eq 1 ]; then
-  new_inisol=$(echo "$new_inisol - $year_day" | bc -l)
+    new_inisol=$(echo "$new_inisol - $year_day" | bc -l)
 fi