source: trunk/libIGCM/libIGCM_ensemble/libIGCM_ensemble.ksh @ 737

Last change on this file since 737 was 737, checked in by sdipsl, 12 years ago
  • cosmetics
  • Property svn:keywords set to Revision Author Date
File size: 18.2 KB
Line 
1#!/bin/ksh
2
3#**************************************************************
4# Author: Sonia Labetoulle
5# Contact: Sonia.Labetoulle__at__locean-ipsl.upmc.fr
6# $Revision::                                          $ Revision of last commit
7# $Author::                                            $ Author of last commit
8# $Date::                                              $ Date of last commit
9# IPSL (2012)
10#  This software is governed by the CeCILL licence see libIGCM/libIGCM_CeCILL.LIC
11#
12#**************************************************************
13
14function IGCM_ensemble_CastInit
15{
16  IGCM_debug_PushStack "IGCM_ensemble_CastInit"
17
18  set -A Alpha      A B C D E F G H I J K L M N O P Q R S T U V W X Y Z
19  set -A AlphaMonth a b c d e f g h i j k l
20
21  IGCM_sys_Mkdir ${RUN_DIR}
22
23  IGCM_sys_Cp ${SUBMIT_DIR}/config.card   ${RUN_DIR}
24  IGCM_sys_Cp ${SUBMIT_DIR}/ensemble.card ${RUN_DIR}
25  IGCM_sys_Cp ${SUBMIT_DIR}/Job_*         ${RUN_DIR}
26  IGCM_sys_Cp ${SUBMIT_DIR}/run.card.init ${RUN_DIR}
27
28  if [ -f  ${SUBMIT_DIR}/CreatedDir.txt ] ; then
29    IGCM_sys_Cp ${SUBMIT_DIR}/CreatedDir.txt ${RUN_DIR}
30  fi
31  if [ -f  ${SUBMIT_DIR}/Qsub.txt ] ; then
32    IGCM_sys_Cp ${SUBMIT_DIR}/Qsub.txt ${RUN_DIR}
33  fi
34
35  IGCM_card_DefineVariableFromOption ${FileName} Ens_PERTURB actif
36  IGCM_card_DefineVariableFromOption ${FileName} Ens_PERTURB NAME
37  IGCM_card_DefineVariableFromOption ${FileName} Ens_PERTURB DEBUT_INIT
38  IGCM_card_DefineVariableFromOption ${FileName} Ens_PERTURB FIN_INIT
39  IGCM_card_DefineVariableFromOption ${FileName} Ens_PERTURB PERIODICITE
40  IGCM_card_DefineArrayFromOption    ${FileName} Ens_PERTURB NONPERIODIC
41  IGCM_card_DefineVariableFromOption ${FileName} Ens_PERTURB DUREE
42  IGCM_card_DefineArrayFromOption    ${FileName} Ens_PERTURB DUREE_NONPERIODIC
43  IGCM_card_DefineVariableFromOption ${FileName} Ens_PERTURB MEMBER
44  IGCM_card_DefineArrayFromOption    ${FileName} Ens_PERTURB PERTURB_BIN
45  IGCM_card_DefineVariableFromOption ${FileName} Ens_PERTURB INITFROM
46  IGCM_card_DefineVariableFromOption ${FileName} Ens_PERTURB INITPATH
47  IGCM_card_DefineVariableFromOption config.card UserChoices TagName
48  IGCM_card_DefineVariableFromOption config.card UserChoices CalendarType
49  IGCM_card_DefineArrayFromSection   config.card ListOfComponents
50
51  echo
52  IGCM_debug_Print 1 "ACTIVE            = ${ensemble_Ens_PERTURB_actif}"
53  IGCM_debug_Print 1 "NAME              = ${ensemble_Ens_PERTURB_NAME}"
54  IGCM_debug_Print 1 "DEBUT_INIT        = ${ensemble_Ens_PERTURB_DEBUT_INIT}"
55  IGCM_debug_Print 1 "FIN_INIT          = ${ensemble_Ens_PERTURB_FIN_INIT}"
56  IGCM_debug_Print 1 "PERIODICITE       = ${ensemble_Ens_PERTURB_PERIODICITE}"
57  IGCM_debug_Print 1 "NONPERIODIC       = ${ensemble_Ens_PERTURB_NONPERIODIC[*]}"
58  IGCM_debug_Print 1 "DUREE             = ${ensemble_Ens_PERTURB_DUREE}"
59  IGCM_debug_Print 1 "DUREE_NONPERIODIC = ${ensemble_Ens_PERTURB_DUREE_NONPERIODIC[*]}"
60  IGCM_debug_Print 1 "MEMBER            = ${ensemble_Ens_PERTURB_MEMBER}"
61  IGCM_debug_Print 1 "PERTURB_BIN       = ${ensemble_Ens_PERTURB_PERTURB_BIN[*]}"
62  IGCM_debug_Print 1 "INITFROM          = ${ensemble_Ens_PERTURB_INITFROM}"
63  IGCM_debug_Print 1 "INITPATH          = ${ensemble_Ens_PERTURB_INITPATH}"
64  IGCM_debug_Print 1 "TagName           = ${config_UserChoices_TagName}"
65  IGCM_debug_Print 1 "CalendarType      = ${config_UserChoices_CalendarType}"
66  IGCM_debug_Print 1 "ListOfComponents  = ${config_ListOfComponents[*]}"
67
68  PerturbExe=${ensemble_Ens_PERTURB_PERTURB_BIN[0]}
69  PerturbComp=${ensemble_Ens_PERTURB_PERTURB_BIN[1]}
70  PerturbFile=${ensemble_Ens_PERTURB_PERTURB_BIN[2]}
71  PerturbVar=${ensemble_Ens_PERTURB_PERTURB_BIN[3]}
72  PerturbAmp=${ensemble_Ens_PERTURB_PERTURB_BIN[4]}
73
74  IGCM_debug_Print 1 "PerturbExe  = ${PerturbExe}"
75  IGCM_debug_Print 1 "PerturbFile = ${PerturbFile}"
76  IGCM_debug_Print 1 "PerturbComp = ${PerturbComp}"
77  IGCM_debug_Print 1 "PerturbVar  = ${PerturbVar}"
78  IGCM_debug_Print 1 "PerturbAmp  = ${PerturbAmp}"
79
80  # A few checks Period case:
81
82  # ... Check PERIODICITE ...
83  case ${ensemble_Ens_PERTURB_PERIODICITE} in
84  _0_)
85    IGCM_debug_Print 1 "periodic start not active"
86    CastPeriodicStart=false
87    return
88    ;;
89  *[Yy]|*[Mm])
90    CastPeriodicStart=true
91    IGCM_debug_Print 1 "Periodic length : ${ensemble_Ens_PERTURB_PERIODICITE}" ;;
92  *)
93    IGCM_debug_Exit "IGCM_ensemble_CastInit ${ensemble_Ens_PERTURB_PERIODICITE} : invalid PERIODICITE"
94    IGCM_debug_Exit "Choose a value in *Y or *M"
95    IGCM_debug_Verif_Exit ;;
96  esac
97# ... Check DUREE ...
98  case ${ensemble_Ens_PERTURB_DUREE} in
99  *[Yy]|*[Mm])
100    IGCM_debug_Print 1 "Periodic duration : ${ensemble_Ens_PERTURB_DUREE}" ;;
101  *)
102    IGCM_debug_Exit "IGCM_ensemble_CastInit ${ensemble_Ens_PERTURB_DUREE} invalid DUREE"
103    IGCM_debug_Exit "Choose a value in choose in *Y or *M"
104    IGCM_debug_Verif_Exit ;;
105  esac
106
107  # A few checks for the Non-Periodic case:
108  DateNum=0
109  while [ $DateNum -lt ${#ensemble_Ens_PERTURB_NONPERIODIC[*]} ] ; do
110
111    # - Check DUREE_NONPERIODIC
112    case ${ensemble_Ens_PERTURB_DUREE_NONPERIODIC[$DateNum]} in
113    _0_)
114      IGCM_debug_Print 1 "non-periodic start not active"
115      CastNonPeriodicStart=false
116      return
117      ;;
118    *[Yy]|*[Mm])
119      IGCM_debug_Print 1 "Non-periodic duration : ${ensemble_Ens_PERTURB_DUREE_NONPERIODIC[$DateNum]}"
120      CastNonPeriodicStart=true
121      ;;
122    *)
123      IGCM_debug_Exit "IGCM_ensemble_CastInit ${ensemble_Ens_PERTURB_DUREE_NONPERIODIC[$DateNum]} : invalid DUREE"
124      IGCM_debug_Exit "choose in *Y or *M"
125      IGCM_debug_Verif_Exit ;;
126    esac
127    (( DateNum = DateNum + 1 ))
128  done
129
130  IGCM_debug_PopStack "IGCM_ensemble_CastInit"
131}
132
133function IGCM_ensemble_CastPeriodicStarts
134{
135  IGCM_debug_PushStack "IGCM_ensemble_CastPeriodicStarts"
136
137  [ ${CastPeriodicStart} = false ] && return
138
139  echo
140  IGCM_debug_Print 1 "Manage periodic starts"
141
142#.. Manage periodic starts ..
143#   ======================
144
145# ... Loop over DateBegin ...
146  eval DateBegin=\${ensemble_Ens_PERTURB_DEBUT_INIT}
147
148  while [ ${DateBegin} -le ${ensemble_Ens_PERTURB_FIN_INIT} ] ; do
149    IGCM_date_GetYearMonth ${DateBegin} year month
150
151  # - Determine number of day(s) in PERIODICITE
152    PeriodLengthInDays=$( IGCM_date_DaysInCurrentPeriod ${DateBegin} ${ensemble_Ens_PERTURB_PERIODICITE} )
153
154  # - Determine number of day(s) in DUREE
155    DureeLengthInDays=$(( $( IGCM_date_DaysInCurrentPeriod ${DateBegin} ${ensemble_Ens_PERTURB_DUREE} ) - 1 ))
156
157  # - Determine DateEnd
158    (( DateEnd = $( IGCM_date_AddDaysToGregorianDate ${DateBegin} ${DureeLengthInDays} ) ))
159
160  # - Build directory name
161    IGCM_ensemble_CastDirectoryName ${ensemble_Ens_PERTURB_NAME} ${ensemble_Ens_PERTURB_PERIODICITE} $year $month $StartDir
162
163  # - Determine RestartDate
164    (( Offset = -1 ))
165    (( RestartDate = $( IGCM_date_AddDaysToGregorianDate ${DateBegin} ${Offset} ) ))
166
167    IGCM_debug_Print 2 "${DateBegin} => ${DateEnd} : ${StartDir}"
168    echo "${DateBegin} ${DateEnd} ${StartDir}" >> CreatedDir.txt
169
170  # - Create directory for current DateBegin
171    if [ ! -d  ${StartDir} ] ; then
172      IGCM_sys_Mkdir ${SUBMIT_DIR}/${StartDir}
173      IGCM_sys_Cd ${SUBMIT_DIR}/${StartDir}
174      ln -s ../../.resol .
175      ln -s ../../.libmpi .
176      IGCM_sys_Cd ${RUN_DIR}
177    fi
178
179  # - Create directory in which to store new restart files
180    RestartDir=${STORAGE}/IGCM_IN/${config_UserChoices_TagName}/${StartDir}
181    IGCM_sys_MkdirArchive ${RestartDir}
182
183  # - Loop over members
184    i=0
185    while [ $i -lt ${ensemble_Ens_PERTURB_MEMBER} ] ; do
186      MemberDir="${StartDir}${Alpha[$i]}"
187      echo
188      IGCM_debug_Print 3 "${MemberDir}"
189
190      JobName="Job_${MemberDir}"
191
192    # * Create directory if it doesn't exist and copy/link files
193      if [ ! -d  ${SUBMIT_DIR}/${StartDir}/${MemberDir} ] ; then
194        IGCM_sys_Mkdir ${SUBMIT_DIR}/${StartDir}/${MemberDir}
195        #IGCM_sys_Cp -r COMP/ PARAM/ ${StartDir}/${MemberDir}
196        IGCM_sys_Cd ${SUBMIT_DIR}/${StartDir}/${MemberDir}
197        ln -s ../../COMP
198        ln -s ../../PARAM
199        ln -s ../../POST
200        IGCM_sys_Cd ${RUN_DIR}
201        IGCM_sys_Cp config.card run.card.init ${SUBMIT_DIR}/${StartDir}/${MemberDir}
202        IGCM_sys_Cp Job_${ensemble_Ens_PERTURB_NAME} ${SUBMIT_DIR}/${StartDir}/${MemberDir}/${JobName}
203
204        # Dump command to be lauched in Qsub.txt
205        echo "cd ${StartDir}/${MemberDir}/ ;"  >> Qsub.txt
206        echo "sed -e 's/PeriodNb=.*/PeriodNb=48/' ${JobName} > ${JobName}.tmp ; " >> Qsub.txt
207        echo "mv ${JobName}.tmp ${JobName} ; " >> Qsub.txt
208        echo "ccc_msub ${JobName} ; cd -"      >> Qsub.txt
209
210      fi
211
212    # * Update files : config.card, Job_, COMP/comp.card
213      IGCM_ensemble_CastFilesUpdate ${DateBegin} ${DateEnd} ${RestartDate}
214
215    # * Apply noise on restart file
216      IGCM_ensemble_CastPerturbFile
217
218      (( i = i + 1 ))
219    done
220
221  # - Next DateBegin
222    echo "$DateBegin  $PeriodLengthInDays"
223    case ${ensemble_Ens_PERTURB_PERIODICITE} in
224    *[Yy]|*[Mm])
225      (( DateBegin = $( IGCM_date_AddDaysToGregorianDate ${DateBegin} ${PeriodLengthInDays} ) ))
226      ;;
227    esac
228    echo "New DateBegin = $DateBegin"
229    echo "========================================================================"
230  done
231  IGCM_debug_PopStack "IGCM_ensemble_CastPeriodicStarts"
232}
233
234function IGCM_ensemble_CastNonPeriodicStarts
235{
236  IGCM_debug_PushStack "IGCM_ensemble_CastNonPeriodicStarts"
237
238  #.. Manage non periodic starts => Loop over DateBegin ..
239  #   ==========================
240
241  [ ${CastNonPeriodicStart} = false ] && return
242
243  echo
244  IGCM_debug_Print 1 "Manage non periodic starts"
245
246  DateNum=0
247# ... Loop over ensemble_Ens_PERTURB_NONPERIODIC ...
248  echo ">${DateNum}<"
249  echo ">${#ensemble_Ens_PERTURB_NONPERIODIC[*]}<"
250  while [ $DateNum -lt ${#ensemble_Ens_PERTURB_NONPERIODIC[*]} ] ; do
251    DateBegin=${ensemble_Ens_PERTURB_NONPERIODIC[$DateNum]}
252    Duree=${ensemble_Ens_PERTURB_DUREE_NONPERIODIC[$DateNum]}
253    echo ">${DateBegin}<"
254    echo ">${Duree}<"
255
256  # - Determine number of day(s) in DUREE_NONPERIODIC
257    IGCM_date_GetYearMonth ${DateBegin} year month
258    DureeLengthInDays=$(( $( IGCM_date_DaysInCurrentPeriod ${DateBegin} ${Duree} ) - 1 ))
259
260  # - Build directory name
261    echo "========================================================================"
262    echo "ensemble_Ens_PERTURB_NAME = ${ensemble_Ens_PERTURB_NAME}"
263    IGCM_ensemble_CastDirectoryName ${ensemble_Ens_PERTURB_NAME} ${Duree} $year $month $StartDir
264
265  # - Determine DateEnd
266    (( DateEnd = $( IGCM_date_AddDaysToGregorianDate ${DateBegin} ${DureeLengthInDays} ) ))
267    echo ">${DateEnd}<"
268    echo "tout va bien 1"
269
270  # - Determine RestartDate
271    (( Offset = -1 ))
272    (( RestartDate = $( IGCM_date_AddDaysToGregorianDate ${DateBegin} ${Offset} ) ))
273
274    IGCM_debug_Print 2 "${DateBegin} => ${DateEnd} : ${StartDir}"
275
276  # -  Does $StartDir already exist ?
277    #echo "tout va bien 2" ${StartDir}
278    if [ ! -d ${SUBMIT_DIR}/${StartDir} ] ; then
279      echo "create dir"
280      IGCM_sys_Mkdir ${SUBMIT_DIR}/${StartDir}
281      IGCM_sys_Cd ${SUBMIT_DIR}/${StartDir}
282      ln -s ../../.resol .
283      ln -s ../../.libmpi .
284      IGCM_sys_Cd ${RUN_DIR}
285      echo "${DateBegin} ${DateEnd} ${StartDir}" >> CreatedDir.txt
286    fi
287    PeriodDateEnd=$( grep -m1 ${StartDir} CreatedDir.txt | cut -f2 -d\  )
288
289  # - Create directory in which to store new restart files if it does'nt already exist
290    RestartDir=${STORAGE}/IGCM_IN/${config_UserChoices_TagName}/${StartDir}
291    IGCM_sys_MkdirArchive ${RestartDir}
292
293  # - Loop over members
294    i=0
295    while [ $i -lt ${ensemble_Ens_PERTURB_MEMBER} ] ; do
296      MemberDir="${StartDir}${Alpha[$i]}"
297      IGCM_debug_Print 3 "${MemberDir}"
298
299      JobName="Job_${MemberDir}"
300
301    # * Create directory if it doesn't exist and copy files
302      if [ ! -d  ${SUBMIT_DIR}/${StartDir}/${MemberDir} ] ; then
303        IGCM_sys_Mkdir ${SUBMIT_DIR}/${StartDir}/${MemberDir}
304        #IGCM_sys_Cp -r COMP/ PARAM/ ${StartDir}/${MemberDir}
305        IGCM_sys_Cd ${SUBMIT_DIR}/${StartDir}/${MemberDir}
306        ln -s ../../COMP
307        ln -s ../../PARAM
308        ln -s ../../POST
309        IGCM_sys_Cd ${RUN_DIR}
310        IGCM_sys_Cp config.card run.card.init ${SUBMIT_DIR}/${StartDir}/${MemberDir}
311        IGCM_sys_Cp Job_${ensemble_Ens_PERTURB_NAME} ${SUBMIT_DIR}/${StartDir}/${MemberDir}/Job_${MemberDir}
312
313        # Dump command to be lauched in Qsub.txt
314        echo "cd ${StartDir}/${MemberDir}/ ;"  >> Qsub.txt
315        echo "sed -e 's/PeriodNb=.*/PeriodNb=48/' ${JobName} > ${JobName}.tmp ; " >> Qsub.txt
316        echo "mv ${JobName}.tmp ${JobName} ; " >> Qsub.txt
317        echo "ccc_msub ${JobName} ; cd -"      >> Qsub.txt
318      fi
319
320    # * Update files : config.card, Job_, COMP/comp.card
321      echo "${PeriodDateEnd} ? ${DateEnd}"
322      if [ ${PeriodDateEnd} -gt ${DateEnd} ] ; then
323        DateEnd=${PeriodDateEnd}
324      fi
325      IGCM_ensemble_CastFilesUpdate ${DateBegin} ${DateEnd} ${RestartDate}
326
327    # * Apply noise on restart file
328      IGCM_ensemble_CastPerturbFile
329      (( i = i + 1 ))
330    done
331
332    (( DateNum = DateNum + 1 ))
333  done
334  IGCM_debug_PopStack "IGCM_ensemble_CastNonPeriodicStarts"
335}
336
337function IGCM_ensemble_CastFilesUpdate
338{
339  IGCM_debug_PushStack "IGCM_ensemble_CastFilesUpdate"
340
341#  set -vx
342
343  # Debug Print :
344  echo
345  IGCM_debug_Print 1 "IGCM_ensemble_CastFilesUpdate :"
346
347  HumanDateBegin=$(   IGCM_date_ConvertFormatToHuman ${1} )
348  HumanDateEnd=$(     IGCM_date_ConvertFormatToHuman ${2} )
349  HumanRestartDate=$( IGCM_date_ConvertFormatToHuman ${3} )
350  # ==> config.card
351  IGCM_card_WriteOption ${SUBMIT_DIR}/${StartDir}/${MemberDir}/config.card Ensemble EnsembleRun 'y'
352  IGCM_card_WriteOption ${SUBMIT_DIR}/${StartDir}/${MemberDir}/config.card Ensemble EnsembleName ${ensemble_Ens_PERTURB_NAME}
353  IGCM_card_WriteOption ${SUBMIT_DIR}/${StartDir}/${MemberDir}/config.card Ensemble EnsembleDate ${StartDir}
354
355  IGCM_card_WriteOption ${SUBMIT_DIR}/${StartDir}/${MemberDir}/config.card UserChoices JobName   ${MemberDir}
356  IGCM_card_WriteOption ${SUBMIT_DIR}/${StartDir}/${MemberDir}/config.card UserChoices DateBegin ${HumanDateBegin}
357  IGCM_card_WriteOption ${SUBMIT_DIR}/${StartDir}/${MemberDir}/config.card UserChoices DateEnd   ${HumanDateEnd}
358  IGCM_card_WriteOption ${SUBMIT_DIR}/${StartDir}/${MemberDir}/config.card Restarts OverRule "n"
359
360  for Comp in ${config_ListOfComponents[*]} ; do
361    IGCM_card_WriteOption ${SUBMIT_DIR}/${StartDir}/${MemberDir}/config.card ${Comp} Restart "y"
362    IGCM_card_WriteOption ${SUBMIT_DIR}/${StartDir}/${MemberDir}/config.card ${Comp} RestartDate ${HumanRestartDate}
363    IGCM_card_WriteOption ${SUBMIT_DIR}/${StartDir}/${MemberDir}/config.card ${Comp} RestartJobName ${ensemble_Ens_PERTURB_INITFROM}
364    IGCM_card_WriteOption ${SUBMIT_DIR}/${StartDir}/${MemberDir}/config.card ${Comp} RestartPath ${ensemble_Ens_PERTURB_INITPATH}
365  done
366  # ==> Job
367  # SD : to be encapsulated in a sys function
368  sed -e "s/\(#MSUB -r \)[A-Za-z0-9]*\( *#.*\)/\1 ${MemberDir} \2/"           \
369      -e "s/\(#MSUB -o Script_Output_\)[A-Za-z0-9]*\(\.*\)/\1${MemberDir}\2/" \
370      -e "s/\(#MSUB -e Script_Output_\)[A-Za-z0-9]*\(\.*\)/\1${MemberDir}\2/" \
371      ${SUBMIT_DIR}/${StartDir}/${MemberDir}/Job_${MemberDir} > Job_${MemberDir}.tmp
372  IGCM_sys_Mv Job_${MemberDir}.tmp ${SUBMIT_DIR}/${StartDir}/${MemberDir}/Job_${MemberDir}
373
374  IGCM_debug_PopStack "IGCM_ensemble_CastFilesUpdate"
375}
376
377function IGCM_ensemble_CastDirectoryName
378{
379  IGCM_debug_PushStack "IGCM_ensemble_CastDirectoryName"
380
381  #.. Debug Print ..
382  echo
383  IGCM_debug_Print 1 "IGCM_ensemble_CastDirectoryName :"
384  echo
385
386#  set -vx
387
388  Name=$1
389  Duree=$2
390  year=$3
391  month=$4
392
393  # - Build directory name
394  case ${Duree} in
395    *Y|*y)
396      siecle="$( echo $year | cut -c1-2 )"
397      siecle=$( (( $siecle - 18 )) )
398      StartYear="${siecle}$( echo $year | cut -c3-4 )"
399      StartDir="${Name}${StartYear}"
400      ;;
401    *M|*m)
402      echo $month
403      siecle="$( echo $year | cut -c1-2 )"
404      siecle=$( (( $siecle - 18 )) )
405      StartYear="${siecle}$( echo $year | cut -c3-4 )"
406      StartMonth="${AlphaMonth[ (( $month - 1 )) ]}"
407      StartDir="${Name}${StartYear}${StartMonth}"
408      ;;
409  esac
410
411  IGCM_debug_PopStack "IGCM_ensemble_CastDirectoryName"
412}
413
414function IGCM_ensemble_CastPerturbFile
415{
416  IGCM_debug_PushStack "IGCM_ensemble_CastPerturbFile"
417
418  #.. Debug Print ..
419  echo
420  IGCM_debug_Print 1 "IGCM_ensemble_CastPerturbFile :"
421
422  #.. FileIn ? => RestartDate ..
423  DirIn="${ensemble_Ens_PERTURB_INITPATH}/${ensemble_Ens_PERTURB_INITFROM}/${PerturbComp}/Restart"
424  FileIn="${PerturbComp}_${ensemble_Ens_PERTURB_INITFROM}_${RestartDate}_${PerturbFile}.nc"
425  DirOut="${RestartDir}/${MemberDir}/${PerturbComp}/Restart"
426#  FileOut="${ensemble_Ens_PERTURB_INITFROM}_${RestartDate}_${PerturbFile}_${MemberDir}.nc"
427#  FileOut=${FileIn}
428  FileOut="${MemberDir}_${RestartDate}_${PerturbFile}.nc"
429  IGCM_debug_Print 1 "FileIn  = ${DirIn}/${FileIn}"
430  IGCM_debug_Print 1 "FileOut = ${DirOut}/${FileOut}"
431
432  # * Create member restart directory
433  IGCM_sys_TestDirArchive ${DirOut}
434  RET=$?
435  if [ $RET -gt 0 ] ; then
436    IGCM_sys_MkdirArchive ${DirOut}
437  fi
438
439  IGCM_sys_TestFileArchive ${DirOut}/${FileOut}
440  RET=$?
441  if [ $RET -gt 0 ] ; then
442
443#    ${PerturbExe} ${DirIn}/${FileIn} ${DirOut}/${FileOut} \
444#                  ${PerturbVar} ${PerturbAmp}
445#    IGCM_sys_Cp ${DirIn}/${FileIn} ${DirOut}/${FileOut}
446#    IGCM_sys_Chmod u+w ${DirOut}/${FileOut}
447#    ${PerturbExe} ${DirOut}/${FileOut} ${PerturbVar} ${PerturbAmp}
448
449    IGCM_sys_Get ${DirIn}/${FileIn} ${RUN_DIR}/${FileOut}
450    IGCM_sys_Chmod 744 ${RUN_DIR}/${FileOut}
451
452    IGCM_debug_Print 1 "${PerturbExe} ${RUN_DIR}/${FileOut} ${PerturbVar} ${PerturbAmp}"
453    echo
454    ${PerturbExe} ${RUN_DIR}/${FileOut} ${PerturbVar} ${PerturbAmp} > /dev/null 2>&1
455
456    if [ $? -ne 0 ] ; then
457      echo "Abend $( basename ${PerturbExe} )"
458      return
459    fi
460
461    IGCM_sys_Put_Out ${RUN_DIR}/${FileOut} ${DirOut}/ 744
462
463#    IGCM_sys_Cd ${DirOut}
464#    ln -s ${DirIn}/${ensemble_Ens_PERTURB_INITFROM}_${RestartDate}_flxat.nc ${MemberDir}_${RestartDate}_flxat.nc
465#    IGCM_sys_Cd -
466    IGCM_sys_RshArchive "cd ${DirOut} ; ln -s ${DirIn}/${ensemble_Ens_PERTURB_INITFROM}_${RestartDate}_flxat.nc ${MemberDir}_${RestartDate}_flxat.nc"
467  fi
468
469  #.. Update config.card..
470  IGCM_card_WriteOption ${SUBMIT_DIR}/${StartDir}/${MemberDir}/config.card ${PerturbComp} Restart "y"
471  IGCM_card_WriteOption ${SUBMIT_DIR}/${StartDir}/${MemberDir}/config.card ${PerturbComp} RestartDate    ${HumanRestartDate}
472  IGCM_card_WriteOption ${SUBMIT_DIR}/${StartDir}/${MemberDir}/config.card ${PerturbComp} RestartJobName ${MemberDir}
473  IGCM_card_WriteOption ${SUBMIT_DIR}/${StartDir}/${MemberDir}/config.card ${PerturbComp} RestartPath    ${RestartDir}/
474
475  IGCM_debug_PopStack "IGCM_ensemble_CastPerturbFile"
476}
Note: See TracBrowser for help on using the repository browser.