Changeset 3432 for trunk/LMDZ.COMMON/libf/evolution/deftank

Timestamp:

Sep 20, 2024, 5:18:49 PM (11 months ago)

Author:

jbclement

Message:

PEM:
Small corrections for the launching script (run numbers must be integer) and the PEM stopping criteria algorithm related to r3430.
JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/deftank/lib_launchPEM.sh (modified) (5 diffs)

trunk/LMDZ.COMMON/libf/evolution/deftank/lib_launchPEM.sh

@@ -173 +173 @@
         echo "Run PCM $iPCM: call $ii/$2..."
         if [ $1 -eq 1 ]; then # 1D model
-            sed -i "s/^k=[0-9]\+$/k=$(echo "3 - $nPCM_ini" | bc -l)/" PCMrun.job
+            sed -i "s/^k=[0-9]\+$/k=$(echo "3 - $nPCM_ini" | bc)/" PCMrun.job
             ./PCMrun.job
             if [ $? -ne 0 ]; then
@@ -181 +181 @@
             cp PCMrun.job PCMrun${iPCM}.job
             sed -i -E "s/($name_job[^0-9]*[0-9]*[^0-9]*)[0-9]+$/\1${iPCM}/" PCMrun${iPCM}.job
-            sed -i "s/^k=[0-9]\+$/k=$(echo "3 - $nPCM_ini" | bc -l)/" PCMrun${iPCM}.job
+            sed -i "s/^k=[0-9]\+$/k=$(echo "3 - $nPCM_ini" | bc)/" PCMrun${iPCM}.job
             jobID=$(eval "$submit_job PCMrun${iPCM}.job")
             # Create a file to cancel the dependent jobs of the cycle
@@ -198 +198 @@
             echo "Run PCM $iPCM: call $i/$2..."
             if [ $1 -eq 1 ]; then # 1D model
-                sed -i "s/^k=[0-9]\+$/k=$(echo "$i + 2 - $nPCM_ini" | bc -l)/" PCMrun.job
+                sed -i "s/^k=[0-9]\+$/k=$(echo "$i + 2 - $nPCM_ini" | bc)/" PCMrun.job
                 ./PCMrun.job
                 if [ $? -ne 0 ]; then
@@ -206 +206 @@
                 cp PCMrun.job PCMrun${iPCM}.job
                 sed -i -E "s/($name_job[^0-9]*[0-9]*[^0-9]*)[0-9]+$/\1${iPCM}/" PCMrun${iPCM}.job
-                sed -i "s/^k=[0-9]\+$/k=$(echo "$i + 2 - $nPCM_ini" | bc -l)/" PCMrun${iPCM}.job
+                sed -i "s/^k=[0-9]\+$/k=$(echo "$i + 2 - $nPCM_ini" | bc)/" PCMrun${iPCM}.job
                 jobID=$(eval "$submit_dependjob=afterok:${jobID} PCMrun${iPCM}.job")
                 echo $kill_job $jobID >> kill_launchPEM.sh
@@ -363 +363 @@
 # arg1: model dimension
 relaunchPEM() {
-    iPEM=$(($irelaunch + 1))
-    iPCM=$(($nPCM_ini + $nPCM*($irelaunch - 1) + 1))
+    iPEM=$(echo "$irelaunch + 1" | bc)
+    iPCM=$(echo "$nPCM_ini + $nPCM*($irelaunch - 1) + 1" | bc)
     i_myear=$(awk "NR==$(($iPEM + 1)) {print \$1}" "info_PEM.txt")
     sed -i "1s/.*/$i_myear $n_myear $convert_years $iPCM $iPEM $nPCM $nPCM_ini/" info_PEM.txt