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

Timestamp:

May 15, 2025, 6:26:24 PM (9 months 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/irene/run_month1 (modified) (1 diff)

trunk/LMDZ.MARS/deftank/irene/run_month1

@@ -23 +23 @@
 #MSUB -x
 
+## 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
+## to set initial month # and final month # of the simulation
+
 # 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
-# Number of threads to use (must be the same as "#MSUB -c" above) 
+# Number of threads to use (must be the same as "#MSUB -c" above)
 export OMP_NUM_THREADS=5
 export OMP_STACKSIZE=400M
+#
+########################################################################
+set -exv
+ls -al
+trap 'echo -e "Error at line $LINENO!"' ERR
 
-set -xv
-ls -al
-
-## 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}@" \
       -e "s@num_now=${num_now}@num_now=${num_next}@" tmp > run_month$num_next
   rm tmp
-  
   ccc_msub run_month$num_next
 fi