Ignore:
Timestamp:
May 15, 2025, 6:26:24 PM (3 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

Location:
trunk/LMDZ.MARS/deftank/irene
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.MARS/deftank/irene/run0

    r2776 r3766  
    11#!/bin/bash
    2 ###########################################################################
     2########################################################################
    33# Script to perform several chained LMD Mars GCM simulations
    44# SET HERE the maximum total number of simulations
    … …  
    99gcm=gcm_64x48x73_phymars_para.e
    1010#
    11 ###########################################################################
     11########################################################################
    1212set -x
    1313
    1414echo "---------------------------------------------------------"
    15 echo "starting run0"
     15echo "Starting run0..."
    1616
    1717dir=`pwd`
    … …  
    1919address=`whoami`
    2020
    21 # Look for file "num_run" which should contain
     21# Check if GCM executable exists and is executable
     22if [ ! -x $gcm ]; then
     23  echo "Error: file \"$gcm\" not found or not executable in $dir!"
     24  exit 1
     25fi
     26
     27# Look for file "num_run" which should contain
    2228# the value of the previously computed season
    2329# (defaults to 0 if file "num_run" does not exist)
    2430if [[ -r num_run ]] ; then
    25   echo "found file num_run"
     31  echo "Found file \"num_run\""
    2632  numold=`cat num_run`
    2733else
    2834  numold=0
    2935fi
    30 echo "numold is set to" ${numold}
     36echo "numold is set to" $numold
    3137
    32 
    33 # Set value of current season
    34 (( numnew = ${numold} + 1 ))
    35 echo "numnew is set to" ${numnew}
     38# Set value of current season
     39numnew=$((numold + 1 ))
     40echo "numnew is set to" $numnew
    3641
    3742# Look for initialization data files (exit if none found)
    … …  
    4045   \cp -f startfi${numold}.nc startfi.nc
    4146else
    42    if (( ${numold} == 99999 )) ; then
    43     echo "No run because previous run crashed ! (99999 in num_run)"
    44     exit
     47   if (( numold == 99999 )) ; then
     48     echo "Error: no run because previous run crashed! (99999 in \"num_run\")"
     49     exit 1
    4550   else
    46    echo "Where is file start"${numold}".nc??"
    47    exit
     51     echo "Error: missing input files \"start${numold}.nc\" or \"startfi${numold}.nc\" in $dir!"
     52     exit 1
    4853   fi
    4954fi
    … …  
    5156
    5257# Run GCM
    53 ccc_mprun -l $gcm > lrun${numnew} 2>&1 
     58ccc_mprun -l $gcm > lrun${numnew} 2>&1
    5459
    5560
    5661# Check if run ended normaly and copy datafiles
    5762if [[ ( -r restartfi.nc  &&  -r restart.nc ) ]] ; then
    58   echo "Run seems to have ended normaly"
     63  echo "Run seems to have ended normally."
    5964  \mv -f restartfi.nc startfi${numnew}.nc
    6065  \mv -f restart.nc start${numnew}.nc
    6166else
     67  echo "Error: Run crashed or incomplete output!"
    6268  if [[ -r num_run ]] ; then
    6369    \mv -f num_run num_run.crash
    6470  else
    65     echo "No file num_run to build num_run.crash from !!"
     71    echo "No file num_run to build \"num_run.crash\" from!"
    6672    # Impose a default value of 0 for num_run
    6773    echo 0 > num_run.crash
    6874  fi
    69  echo 99999 > num_run
     75  echo 99999 > num_run
    7076############## To receive an Email message if the run crashes ########
    71 mail -s "crash run GCM" $address <<ENDMAIL
    72 The run on $machine in $dir has just crashed.
     77mail -s "Crash in GCM run" $address <<ENDMAIL
     78The GCM run on $machine in $dir has just crashed.
     79Check the output logs for more information.
    7380ENDMAIL
    74 ############################################""
    75  exit
     81#############################################
     82  exit 1
    7683fi
    7784
    … …  
    93100# Prepare things for upcoming runs by writing
    94101# value of computed season in file num_run
    95 echo ${numnew} > num_run
     102echo $numnew > num_run
    96103
    97104# If we are over nummax : stop
    98105if (( $numnew + 1 > $nummax )) ; then
    99    exit
     106  exit 0
    100107else
    101    \cp -f run0 exe_mars
    102    ./exe_mars
    103 fi 
     108  \cp -f run0 exe_mars
     109  ./exe_mars
     110fi

  • trunk/LMDZ.MARS/deftank/irene/run_month1

    r2816 r3766  
    2323#MSUB -x
    2424
     25## Script to run chained simulations
     26## (uses script "run0" and reference file "run.def.ref")
     27## Set values of "num_now" and "num_end" in the script below
     28## to set initial month # and final month # of the simulation
     29
    2530# A few parameters that might need be changed depending on your setup:
    2631# Path to the arch.env to source
    2732source ../trunk/LMDZ.COMMON/arch.env
    28 # Number of threads to use (must be the same as "#MSUB -c" above) 
     33# Number of threads to use (must be the same as "#MSUB -c" above)
    2934export OMP_NUM_THREADS=5
    3035export OMP_STACKSIZE=400M
     36#
     37########################################################################
     38set -exv
     39ls -al
     40trap 'echo -e "Error at line $LINENO!"' ERR
    3141
    32 set -xv
    33 ls -al
    34 
    35 ## set starting month and ending month below:
     42# Set starting month and ending month below:
    3643num_now=1
    3744num_end=12
    38 (( num_previous = $num_now - 1 ))
     45num_previous=$(( num_now - 1 ))
    3946
    4047echo "$num_previous" > num_run
    4148# next month number
    42 (( num_next = $num_now + 1 ))
     49num_next=$(( num_now + 1 ))
    4350# true (i.e. modulo 12) month number
    44 (( true_num = $num_now % 12 ))
     51true_num=$(( num_now % 12 ))
    4552
     53# Check if required files exist
    4654\rm -f  error; touch error
    4755
     56if [ ! -f run.def.ref ]; then
     57  echo "Error: file \"run.def.ref\" not found in current directory!" > error
     58  exit 1
     59fi
     60
     61if [ ! -x run0 ]; then
     62  echo "Error: file \"run0\" not found or not executable in current directory!" > error
     63  exit 1
     64fi
     65
     66if [ ! -f run_month$num_now ]; then
     67  echo "Error: file \"run_month${num_now}\" not found in current directory!" > error
     68  exit 1
     69fi
     70
     71# Run model depending on current month
    4872case $true_num in
    49   1 ) sed s/9999/61/ run.def.ref > run.def ; run0 >> error ;;    #1
    50   2 ) sed s/9999/66/ run.def.ref > run.def ; run0 >> error ;;    #2
    51   3 ) sed s/9999/66/ run.def.ref > run.def ; run0 >> error ;;    #3
    52   4 ) sed s/9999/65/ run.def.ref > run.def ; run0 >> error ;;    #4
    53   5 ) sed s/9999/60/ run.def.ref > run.def ; run0 >> error ;;    #5
    54   6 ) sed s/9999/54/ run.def.ref > run.def ; run0 >> error ;;    #6
    55   7 ) sed s/9999/50/ run.def.ref > run.def ; run0 >> error ;;    #7
    56   8 ) sed s/9999/46/ run.def.ref > run.def ; run0 >> error ;;    #8
    57   9 ) sed s/9999/47/ run.def.ref > run.def ; run0 >> error ;;    #9
    58  10 ) sed s/9999/47/ run.def.ref > run.def ; run0 >> error ;;    #10
    59  11 ) sed s/9999/51/ run.def.ref > run.def ; run0 >> error ;;    #11
    60   0 ) sed s/9999/56/ run.def.ref > run.def ; run0 >> error ;;    #12
    61   * ) echo "error" ;;
     73  1 ) sed s/9999/61/ run.def.ref > run.def ; ./run0 >> error ;;    #1
     74  2 ) sed s/9999/66/ run.def.ref > run.def ; ./run0 >> error ;;    #2
     75  3 ) sed s/9999/66/ run.def.ref > run.def ; ./run0 >> error ;;    #3
     76  4 ) sed s/9999/65/ run.def.ref > run.def ; ./run0 >> error ;;    #4
     77  5 ) sed s/9999/60/ run.def.ref > run.def ; ./run0 >> error ;;    #5
     78  6 ) sed s/9999/54/ run.def.ref > run.def ; ./run0 >> error ;;    #6
     79  7 ) sed s/9999/50/ run.def.ref > run.def ; ./run0 >> error ;;    #7
     80  8 ) sed s/9999/46/ run.def.ref > run.def ; ./run0 >> error ;;    #8
     81  9 ) sed s/9999/47/ run.def.ref > run.def ; ./run0 >> error ;;    #9
     82 10 ) sed s/9999/47/ run.def.ref > run.def ; ./run0 >> error ;;    #10
     83 11 ) sed s/9999/51/ run.def.ref > run.def ; ./run0 >> error ;;    #11
     84  0 ) sed s/9999/56/ run.def.ref > run.def ; ./run0 >> error ;;    #12
     85  * ) echo "Error: Invalid value of true_num ($true_num)" ; exit 1 ;;
    6286esac
    6387
    64 # launch job for next month
    65 if (( $num_next <= $num_end )) ; then
     88# Launch job for next month
     89if (( num_next <= num_end )) ; then
    6690  cp -f run_month$num_now tmp
    6791  sed -e "s@run_month${num_now}@run_month${num_next}@" \
    6892      -e "s@num_now=${num_now}@num_now=${num_next}@" tmp > run_month$num_next
    6993  rm tmp
    70  
    7194  ccc_msub run_month$num_next
    7295fi
Note: See TracChangeset for help on using the changeset viewer.