Changeset 3766 for trunk/LMDZ.MARS/deftank/ciclad/run_month1

Timestamp:

May 15, 2025, 6:26:24 PM (6 weeks ago)

Author:

jbclement

Message:

Mars PCM:
Cleaning and improvement of robustness for "run0" and "run_month1" scripts with file checks and clearer errors.
JBC

File:

• trunk/LMDZ.MARS/deftank/ciclad/run_month1 (modified) (2 diffs)

trunk/LMDZ.MARS/deftank/ciclad/run_month1

@@ -8 +8 @@
 #PBS -l walltime=20:00:00
 
-## Script to run chained simulations on Ciclad
+## Script to run chained simulations
 ## (uses script "run0" and reference file "run.def.ref")
 ## Set values of "num_now" and "num_end" in the script below
@@ -16 +16 @@
 cd $PBS_O_WORKDIR
 
-# source the environment (you might need to adapt this path to your case)
+# A few parameters that might need be changed depending on your setup:
+# Path to the arch.env to source
 source ../trunk/LMDZ.COMMON/arch.env
+#
+########################################################################
+set -exv
+ls -al
+trap 'echo -e "Error at line $LINENO!"' ERR
 
-date
-
-## set starting month and ending month below:
+# Set starting month and ending month below:
 num_now=1
 num_end=12
-(( num_previous = $num_now - 1 ))
+num_previous=$(( num_now - 1 ))
 
 echo "$num_previous" > num_run
 # next month number
-(( num_next = $num_now + 1 ))
+num_next=$(( num_now + 1 ))
 # true (i.e. modulo 12) month number
-(( true_num = $num_now % 12 ))
+true_num=$(( num_now % 12 ))
 
+# Check if required files exist
 \rm -f  error; touch error
 
+if [ ! -f run.def.ref ]; then
+  echo "Error: file \"run.def.ref\" not found in current directory!" > error
+  exit 1
+fi
+
+if [ ! -x run0 ]; then
+  echo "Error: file \"run0\" not found or not executable in current directory!" > error
+  exit 1
+fi
+
+if [ ! -f run_month$num_now ]; then
+  echo "Error: file \"run_month${num_now}\" not found in current directory!" > error
+  exit 1
+fi
+
+# Run model depending on current month
 case $true_num in
-  1 ) sed s/9999/61/ run.def.ref > run.def ; run0 >> error ;;    #1
-  2 ) sed s/9999/66/ run.def.ref > run.def ; run0 >> error ;;    #2
-  3 ) sed s/9999/66/ run.def.ref > run.def ; run0 >> error ;;    #3
-  4 ) sed s/9999/65/ run.def.ref > run.def ; run0 >> error ;;    #4
-  5 ) sed s/9999/60/ run.def.ref > run.def ; run0 >> error ;;    #5
-  6 ) sed s/9999/54/ run.def.ref > run.def ; run0 >> error ;;    #6
-  7 ) sed s/9999/50/ run.def.ref > run.def ; run0 >> error ;;    #7
-  8 ) sed s/9999/46/ run.def.ref > run.def ; run0 >> error ;;    #8
-  9 ) sed s/9999/47/ run.def.ref > run.def ; run0 >> error ;;    #9
- 10 ) sed s/9999/47/ run.def.ref > run.def ; run0 >> error ;;    #10
- 11 ) sed s/9999/51/ run.def.ref > run.def ; run0 >> error ;;    #11 
-  0 ) sed s/9999/56/ run.def.ref > run.def ; run0 >> error ;;    #12
-  * ) echo "error" ;;
+  1 ) sed s/9999/61/ run.def.ref > run.def ; ./run0 >> error ;;    #1
+  2 ) sed s/9999/66/ run.def.ref > run.def ; ./run0 >> error ;;    #2
+  3 ) sed s/9999/66/ run.def.ref > run.def ; ./run0 >> error ;;    #3
+  4 ) sed s/9999/65/ run.def.ref > run.def ; ./run0 >> error ;;    #4
+  5 ) sed s/9999/60/ run.def.ref > run.def ; ./run0 >> error ;;    #5
+  6 ) sed s/9999/54/ run.def.ref > run.def ; ./run0 >> error ;;    #6
+  7 ) sed s/9999/50/ run.def.ref > run.def ; ./run0 >> error ;;    #7
+  8 ) sed s/9999/46/ run.def.ref > run.def ; ./run0 >> error ;;    #8
+  9 ) sed s/9999/47/ run.def.ref > run.def ; ./run0 >> error ;;    #9
+ 10 ) sed s/9999/47/ run.def.ref > run.def ; ./run0 >> error ;;    #10
+ 11 ) sed s/9999/51/ run.def.ref > run.def ; ./run0 >> error ;;    #11
+  0 ) sed s/9999/56/ run.def.ref > run.def ; ./run0 >> error ;;    #12
+  * ) echo "Error: Invalid value of true_num ($true_num)" ; exit 1 ;;
 esac
 
-# launch job for next month
-if (( $num_next <= $num_end )) ; then
+# Launch job for next month
+if (( num_next <= num_end )) ; then
   cp -f run_month$num_now tmp
   sed -e "s@run_month${num_now}@run_month${num_next}@" \