New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
trusting_func.sh in branches/2015/dev_r5092_CNRS_SETTE/NEMOGCM/TRUST – NEMO

source: branches/2015/dev_r5092_CNRS_SETTE/NEMOGCM/TRUST/trusting_func.sh @ 5672

Last change on this file since 5672 was 5672, checked in by nicolasmartin, 9 years ago

dev_r5092_CNRS_SETTE Trusting functions: bugfix & rebuild all restart files

  • Property eol-style set to native
  • Property svn:executable set to *
  • Property svn:keywords set to Author Date Id Rev URL
File size: 10.1 KB
Line 
1#!/bin/bash
2
3
4## Messenger filenames
5FILE_DATE=mesg_01_date_$PATTERNAME.txt    ; FILE_STAT=mesg_02_status_$PATTERNAME.txt
6FILE_RESU=mesg_03_result_$PATTERNAME.txt  ; FILE_NEMO=mesg_04_nemogcm_$PATTERNAME.txt
7FILE_XIOS=mesg_05_xios_$PATTERNAME.txt    ; FILE_COMP=mesg_06_compiler_$PATTERNAME.txt
8FILE_MPIN=mesg_07_mpi_$PATTERNAME.txt     ; FILE_NCDF=mesg_08_netcdf_$PATTERNAME.txt
9FILE_INPT=mesg_09_inputs_$PATTERNAME.txt  ; FILE_CPUT=mesg_10_cputime_$PATTERNAME.txt
10FILE_COMM=mesg_11_comments_$PATTERNAME.txt;
11
12## Trusting summary & archive timestamped filenames
13FILE_TRUS=trusting_${DATE}_$PATTERNAME.txt; FILE_ARCH=trusting_${DATE}_$PATTERNAME.tar.gz
14
15print_step() { printf "\nStep.....\n$1\n"; }
16
17init_files() {
18    echo 'Date'          >  ${FILE_DATE}; echo 'Status'           >  ${FILE_STAT}
19    echo 'Result'        >  ${FILE_RESU}; echo 'NEMOGCM rev.'     >  ${FILE_NEMO}
20    echo 'XIOS rev.'     >  ${FILE_XIOS}; echo 'Fortran compiler' >  ${FILE_COMP}
21    echo 'MPI libs'      >  ${FILE_MPIN}; echo 'NetCDF libs'      >  ${FILE_NCDF}
22    echo 'Input files'   >  ${FILE_INPT}; echo 'Real CPU time'    >  ${FILE_CPUT}
23    echo 'Comments'      >  ${FILE_COMM}
24
25    # 'Failed' status with 'Unknown error' by default
26    echo $ST             >> ${FILE_STAT}; echo 'Unknown error'    >> ${FILE_RESU}
27}
28
29get_date() {
30    local dat=$( date -ud ${DATE} +"%F %R %Z" )
31
32    echo $dat >> ${FILE_DATE}
33}
34
35get_nemo_rev() {
36    local rev_loc svn_cmd
37    local rev=0 svn_opt='status'
38
39    ## If -v|--version option has been set, modify svn command
40    if   [ $( echo ${NEMO_VERS} | grep  "HEAD\|up\|update"                     ) ]; then
41   svn_cmd='svn update -r HEAD'
42    elif [ $( echo ${NEMO_VERS} | grep -o '{[0-9]\{4\}-[0-9]\{2\}-[0-9]\{2\}}' ) ]; then
43   svn_cmd='svn update -r '$( echo ${NEMO_VERS} | grep -o '{[0-9]\{4\}-[0-9]\{2\}-[0-9]\{2\}}' )
44    elif [ $( echo ${NEMO_VERS} | grep -o '[0-9]*'                             ) ]; then
45   svn_cmd='svn update -r '$( echo ${NEMO_VERS} | grep -o '[0-9]*'                             )
46    else
47   svn_cmd='svn '${svn_opt}
48    fi
49
50    for dir in ${NEMO_ARCH} ${NEMO_CONF} ${NEMO_ENGI} \
51               ${NEMO_EIOI} ${NEMO_EFCM}              \
52          ${NEMO_TCMP} ${NEMO_TRBD}              \
53          ${DIR_XIOS}                             ; do
54
55   ## For time being, just get revision number from XIOS
56   if [ "$dir" == "${DIR_XIOS}" ]; then
57       rev_loc=$( svn info $dir | awk '(NR == 9) {print $NF}' )
58       echo 'XIOS '${rev_loc} >> model.log
59       echo "<a href=\"https://forge.ipsl.jussieu.fr/ioserver/changeset/${rev_loc}\" target=\"_blank\">${rev_loc}</a>" \
60      >> ${FILE_XIOS}
61       continue
62   fi
63
64   echo $dir
65   ${svn_cmd} $dir
66   rev_loc=$( svn info $dir | awk '(NR == 9) {print $NF}' )
67   (( ${rev_loc} >= $rev )) && rev=${rev_loc}
68    done
69
70    echo 'NEMOGCM '$rev   >> model.log
71    echo "<a href=\"https://forge.ipsl.jussieu.fr/nemo/changeset/$rev\" target=\"_blank\">$rev</a>" \
72   >> ${FILE_NEMO}
73}
74
75get_soft_rel() {
76    local soft_rel
77
78    for soft in $CDO $COMPILER $MPI $NETCDF; do
79   soft_rel=''
80
81        if (( ${I_MODULE} == 1 )); then
82       soft_rel=$( echo $LOADEDMODULES | sed "s/.*$soft\/\([0-9.a-z_]*\).*/\1/i" )
83   else
84       soft_rel=$( echo $PATH          | sed "s/.*$soft\([0-9.a-z_]*\).*/\1/i"   )
85   fi
86
87   [ "$soft" == "$COMPILER" ] && soft_rel=$( $soft --version | grep -m1 -oe '\<[0-9. ]*\>' )
88   ## Cleaning characters string to display proper soft name
89   soft=$( echo $soft | sed 's#\\##g; s#[/-]$##' )
90   echo $soft ${soft_rel} >> model.log
91    done
92
93    sed -n 4p model.log >> ${FILE_COMP}; sed -n 5p model.log >> ${FILE_MPIN}
94    sed -n 6p model.log >> ${FILE_NCDF}
95}
96
97get_inputs() {
98    local cmd_iol="tar -tvf ${NEMO_FORC}/${NEMO_TARF}" cmd_iof="tar -vxf ${NEMO_FORC}/${NEMO_TARF}"
99
100    [ -z "${NEMO_TARF}" ] && cmd_iol="ls ${NEMO_FORC}/*" && cmd_iof="\cp ${NEMO_FORC}/* ."
101    ${cmd_iol} > inputs_list.txt && ${cmd_iof} > /dev/null
102}
103
104diff_inputs() {
105    local dif
106    local files_list='' mesg='Same' 
107
108    for file in inputs_list.txt namelist_* *.xml; do
109   dif=''
110   if [ -e ${REFE_DIR}/$file ]; then dif=$( diff -q $file ${REFE_DIR}/$file ); else dif=0; fi
111
112   if [ -n "$dif" ]; then
113       mesg='Different'
114       (( "$dif" != '0' )) && echo $dif && files_list+=$file' '
115   fi
116
117    done
118
119    [ "$mesg" == 'Same' ] && echo $mesg
120    echo $mesg >> ${FILE_INPT}
121
122    [ -n "${files_list}" ] && echo 'Inputs  : '${files_list}'differ<br>' \
123   >> temp_${FILE_COMM}
124}
125
126job_pending() {
127    local outline=$( printf "%100s" ) time_elapsed=0 time_increment=30
128
129    sleep ${time_increment}
130
131    while [[ $( ${JOB_LIST} | grep ${JOB_ID} ) && ${time_elapsed} -lt ${TIME_LIMI} ]]; do
132   printf "\n%s\n" ${outline// /#} >> computation.log
133   ${JOB_INFO} ${JOB_ID}           >> computation.log
134   sleep ${time_increment}
135   time_elapsed=$(( ${time_elapsed} + ${time_increment} ))
136    done
137
138    sleep ${time_increment}
139
140    if (( ${time_elapsed} == ${TIME_LIMI} )); then
141   ${JOB_DELE} ${JOB_ID} &> /dev/null
142   let TIME_LIMIT/=3600
143   export ${TIME_LIMIT}
144   get_out 6
145    fi
146}
147
148diff_results() {
149    local files_list=''
150
151    for file in ocean.output *.stat; do
152   [ ! -e ${REFE_DIR}/$file ] && export ST='FAILED' && get_out 7
153   diff -q $file ${REFE_DIR}/$file
154   (( $? != 0 )) && export ST='FAILED' && files_list+=$file' '
155    done
156
157    [ -n "${files_list}" ] && echo 'Results : '${files_list}'differ<br>' \
158   >> temp_${FILE_COMM}
159}
160
161diff_restart() {
162    local base_name dif dift file list_comp list_tmsp nb_dom
163    local files_list='' 
164
165    [ ! -e  ${REFE_DIR}/time.step ] && export ST='FAILED' && get_out 7
166    export TIME_STEP=$( cat ${REFE_DIR}/time.step | tr -d [:space:] )
167    echo 'Last time step of benchmark run: '${TIME_STEP}
168
169    ## Find all restart files to rebuild
170    if [ $( find -regex ".*_restart.*[0-9]\.nc" -print -quit ) ]; then
171   base_name=$( find -regex ".*_restart.*[0-9]\.nc"                      \
172                | sed "s#^\./\(.*\)_[0-9]*_restart.*#\1#"       | sort -u )
173   list_comp=$( find -regex ".*_restart.*[0-9]\.nc"                      \
174                | sed "s#^.*\(restart[a-z_]*\)_[0-9].*\.nc#\1#" | sort -u )
175   list_tmsp=$( find -regex ".*_restart.*[0-9]\.nc"                      \
176                | sed "s#^.*\([0-9]\{8\}\)_restart.*#\1#"       | sort -u )
177
178   for tmsp in ${list_tmsp}; do
179       for comp in ${list_comp}; do
180      file=${base_name}_${tmsp}_${comp}
181      nb_dom=$( find -name "${file}_[0-9]*.nc" | wc -l | awk '{ print $1 }' )
182
183      if   (( ${nb_dom} >  1 )); then
184          ${NEMO_TRBD}/rebuild_nemo -t $NPROC $file ${nb_dom} > /dev/null
185          [ $? -eq 0 ] && rm -f ${file}_[0-9]*.nc             > /dev/null
186      elif (( ${nb_dom} == 0 )); then
187          export ST='FAILED' && get_out 8
188      fi
189
190      ## Compare restart files at same time step
191      if (( $tmsp == ${TIME_STEP} )); then
192          if [ -e ${REFE_DIR}/$file.nc ]; then
193                   ## UNIX `cmp` not suitable (timestamp in .nc file)
194         dif=$( $CDO diffn $file.nc ${REFE_DIR}/$file.nc 2> /dev/null      \
195                | awk '/records/ {print $0}' | sed '2 s/^/,/' | tr -d '\n'  )
196
197         if [[ -n "${dif}" && $( echo $dif | awk '{ print $1 }') -ne 0 ]]; then
198             export ST='FAILED' && files_list+=$file' '
199             let dift+=$( echo $dif | awk '{print $1}' )
200             echo $file'.nc: '$dif
201         else
202             echo 'identical'
203         fi
204
205          else
206         export ST='FAILED' && get_out 7
207          fi
208
209      else
210          continue
211      fi
212
213       done
214   done
215
216   [ -n "${files_list}" ] && echo 'Restarts: '${files_list}$dift' record(s) differ<br>' \
217       >> temp_${FILE_COMM}
218    else
219   export ST='FAILED'
220    fi
221
222}
223
224get_cpu_time() {
225    local real_cpu_time=$( eval ${JOB_TIME} )
226
227    echo ${real_cpu_time} | tee -a ${FILE_CPUT}
228}
229
230comments() {
231    local state=$1
232
233    if [ "$state" == 'E R R O R' ]; then
234   [ -e time.step ] && export TIME_STEP=$( cat time.step | tr -d ' ' )
235   [ -e ocean.output ] && line=$( grep -m1 -A4 "$state"   ocean.output | tr -d '\n' )
236    else
237   [ -e ocean.output ] && line=$( grep -m1 -A2 "^ $state" ocean.output | tr -d '\n' )
238    fi
239
240    echo   $line
241    printf "$line<br>" >> temp_${FILE_COMM}
242}
243
244mesg_make() {
245    ## Format comments for web
246    [ -e temp_${FILE_COMM} ] && cat temp_${FILE_COMM} | tr -d '\n' | sed 's/<br>$//' \
247   >> ${FILE_COMM}
248
249    paste -d ';' mesg_*.txt | tee ${FILE_TRUS}
250}
251
252mesg_publish() {
253    local cmd
254
255    if [ $PUBLISH -eq 1 ]; then
256
257   if [ -f ${REFE_DIR}/trusting_$PATTERNAME.txt ]; then cmd='tail -1'; else cmd='cat'; fi
258
259   $cmd ${FILE_TRUS} >> ${REFE_DIR}/trusting_$PATTERNAME.txt
260
261        ## Send mail only when FAILED
262   if [[ ! -z "$EMAIL" && "$ST" == 'FAILED' ]]; then
263       cat << END_MAIL > trusting.mail
264XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
265
266Dear $USER,
267
268The trusting sequence for ${REFE_CONF} has failed.
269Directory: ${TEST_DIR}
270Archive created: ${FILE_ARCH} in ${REFE_DIR}
271
272XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
273
274END_MAIL
275       [ -e ${TEST_DIR}/${FILE_TRUS} ] && cat ${TEST_DIR}/${FILE_TRUS}                  \
276      >> trusting.mail
277       mail -s "[NEMO Trusting ${NEMO_HPCC} ${REFE_CONF} ${NEMO_BRAN}] $ST $ERR" $EMAIL \
278      <  trusting.mail
279   fi
280
281    fi
282}
283
284get_out() {
285    ERR=$1
286
287    cd ${TEST_DIR}
288
289    if [ "$ST" == 'FAILED' ]; then
290
291        ## Error identification
292   case $ERR in
293       ## Compilation
294       1) ERR='XIOS compilation failed'          ;; 2) ERR='NEMO compilation failed'            ;;
295       ## Submission
296       3) ERR='Missing input files'              ;; 4) ERR='Job submission error'               ;;
297       ## Running
298       5) ERR='Crashed at time step '${TIME_STEP};; 6) ERR='Exceeded time limit '${TIME_LIMI}'h';;
299       ## Results
300       7) ERR='Missing previous outputs'         ;; 8) ERR='New outputs differ'                 ;;
301       ## Undefined
302       *) ERR='Unknown error'                    ;;
303   esac
304
305   export ERR
306    fi
307
308    sed -i "2 s/.*/$ST/"  ${FILE_STAT} && sed -i "2 s/.*/$ERR/" ${FILE_RESU}
309
310    ## Save tested configuration if trusting failed
311    if [[ "$ST" == 'FAILED' && $PUBLISH -eq 1 && $DEBUG -eq 0 ]]; then
312   echo 'Creating archive '${FILE_ARCH}' under '${REFE_DIR}
313   tar -czf ${REFE_DIR}/${FILE_ARCH} * ${TEST_DIR}/${NEMO_CONF}/${TEST_CONF}/cpp_*  \
314                                        ${TEST_DIR}/${NEMO_CONF}/${TEST_CONF}/MY_SRC \
315                                       ${TEST_DIR}/${NEMO_CONF}/${TEST_CONF}/WORK
316    fi
317
318    mesg_make
319    mesg_publish
320
321    exit 0
322}
Note: See TracBrowser for help on using the repository browser.