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

Last change on this file since 1026 was 1026, checked in by sdipsl, 10 years ago

Addon for 3D perturbation mechanism. see #197
Cosmetics
SN
  • Property svn:keywords set to Revision Author Date
File size: 51.7 KB
Line 
1#!/bin/ksh
2
3#**************************************************************
4# Author: Sebastien Denvil, Sonia Labetoulle, Nicolas Lebas, Sebastien Nguyen
5# Contact: Nicolas.Lebas__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# >>> Date ensemble <<<
13# Author: Nicolas Lebas (adapted from Sonia Labetoulle)
14# Contact: Nicolas.Lebas__at__locean-ipsl.upmc.fr
15# IPSL (2014)
16#
17# >>> Add 3D perturbation maps to oceanic restart <<<
18# Author: Sebastien Nguyen
19# Contact: Sebastien.Nguyen__at__locean-ipsl.upmc.fr
20# IPSL (2014)
21#
22#**************************************************************
23
24# Read which ensemble type are active
25function IGCM_ensemble_Init
26{
27  IGCM_debug_PushStack "IGCM_ensemble_Init"
28
29  IGCM_card_DefineVariableFromOption ${F_CFG_ENS} Ens_PERTURB active
30  IGCM_card_DefineVariableFromOption ${F_CFG_ENS} Ens_DATE active
31  IGCM_card_DefineVariableFromOption ${F_CFG_ENS} Ens_PARAMETRIC active
32
33  IGCM_debug_Print 1 "Ens_PERTURB ACTIVE     = ${ensemble_Ens_PERTURB_active}"
34  IGCM_debug_Print 1 "Ens_DATE ACTIVE        = ${ensemble_Ens_DATE_active}"
35  IGCM_debug_Print 1 "Ens_PARAMETRIC ACTIVE  = ${ensemble_Ens_PARAMETRIC_active}"
36  echo ""
37
38  IGCM_debug_PopStack "IGCM_ensemble_Init"
39}
40
41# Set Alphanumerical variables ajust to member nb
42function IGCM_ensemble_SetAlpha
43{
44  IGCM_debug_PushStack "IGCM_ensemble_SetAlpha"
45
46  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
47  set -A AlphaMonth a b c d e f g h i j k l
48
49  IGCM_debug_PopStack "IGCM_ensemble_SetAlpha"
50}
51
52############### Perturb ENSEMBLE #################
53function IGCM_ensemble_CastInit
54{
55  IGCM_debug_PushStack "IGCM_ensemble_CastInit"
56
57  IGCM_sys_Mkdir ${RUN_DIR}
58
59  IGCM_sys_Cp ${SUBMIT_DIR}/config.card   ${RUN_DIR}
60  IGCM_sys_Cp ${SUBMIT_DIR}/ensemble.card ${RUN_DIR}
61  IGCM_sys_Cp ${SUBMIT_DIR}/Job_*         ${RUN_DIR}
62  IGCM_sys_Cp ${SUBMIT_DIR}/run.card.init ${RUN_DIR}
63  if [ -f ${SUBMIT_DIR}/Qsub.* ]; then
64    IGCM_sys_Cp ${SUBMIT_DIR}/Qsub.*        ${RUN_DIR}
65  fi
66  if [ -f ${SUBMIT_DIR}/Qclean.* ]; then
67    IGCM_sys_Cp ${SUBMIT_DIR}/Qclean.*      ${RUN_DIR}
68  fi
69
70  # Useful?
71  #if [ -f  ${SUBMIT_DIR}/CreatedDir.txt ] ; then
72  #  IGCM_sys_Cp ${SUBMIT_DIR}/CreatedDir.txt ${RUN_DIR}
73  #fi
74  # Useful?
75  #if [ -f  ${SUBMIT_DIR}/Qsub.sh ] ; then
76  #  IGCM_sys_Cp ${SUBMIT_DIR}/Qsub.sh ${RUN_DIR}
77  #fi
78  echo ${PWD}
79
80  IGCM_card_DefineVariableFromOption ${F_CFG_ENS} Ens_PERTURB active
81  IGCM_card_DefineVariableFromOption ${F_CFG_ENS} Ens_PERTURB NAME
82  IGCM_card_DefineVariableFromOption ${F_CFG_ENS} Ens_PERTURB BEGIN_INIT
83  IGCM_card_DefineVariableFromOption ${F_CFG_ENS} Ens_PERTURB END_INIT
84  IGCM_card_DefineVariableFromOption ${F_CFG_ENS} Ens_PERTURB PERIODICITY
85  IGCM_card_DefineArrayFromOption    ${F_CFG_ENS} Ens_PERTURB NONPERIODIC
86  IGCM_card_DefineVariableFromOption ${F_CFG_ENS} Ens_PERTURB LENGTH
87  IGCM_card_DefineArrayFromOption    ${F_CFG_ENS} Ens_PERTURB LENGTH_NONPERIODIC
88  IGCM_card_DefineVariableFromOption ${F_CFG_ENS} Ens_PERTURB MEMBER
89  IGCM_card_DefineArrayFromOption    ${F_CFG_ENS} Ens_PERTURB MEMBER_LIST
90  IGCM_card_DefineArrayFromOption    ${F_CFG_ENS} Ens_PERTURB MEMBER_NAMESLIST
91  IGCM_card_DefineVariableFromOption ${F_CFG_ENS} Ens_PERTURB MEMBER_INITFROM
92  IGCM_card_DefineVariableFromOption ${F_CFG_ENS} Ens_PERTURB MEMBER_INITPATH
93  IGCM_card_DefineArrayFromOption    ${F_CFG_ENS} Ens_PERTURB PERTURB_BIN
94  IGCM_card_DefineVariableFromOption ${F_CFG_ENS} Ens_PERTURB INITFROM
95  IGCM_card_DefineVariableFromOption ${F_CFG_ENS} Ens_PERTURB INITPATH
96  IGCM_card_DefineVariableFromOption config.card UserChoices JobName
97  IGCM_card_DefineVariableFromOption config.card UserChoices TagName
98  IGCM_card_DefineVariableFromOption config.card UserChoices CalendarType
99  IGCM_card_DefineArrayFromSection   config.card ListOfComponents
100
101  echo
102  IGCM_debug_Print 1 "[Ens_PERTURB]"
103  IGCM_debug_Print 1 "ACTIVE            = ${ensemble_Ens_PERTURB_active}"
104  IGCM_debug_Print 1 "NAME              = ${ensemble_Ens_PERTURB_NAME}"
105  IGCM_debug_Print 1 "BEGIN_INIT        = ${ensemble_Ens_PERTURB_BEGIN_INIT}"
106  IGCM_debug_Print 1 "END_INIT          = ${ensemble_Ens_PERTURB_END_INIT}"
107  IGCM_debug_Print 1 "PERIODICITY       = ${ensemble_Ens_PERTURB_PERIODICITY}"
108  IGCM_debug_Print 1 "NONPERIODIC       = ${ensemble_Ens_PERTURB_NONPERIODIC[*]}"
109  IGCM_debug_Print 1 "LENGTH             = ${ensemble_Ens_PERTURB_LENGTH}"
110  IGCM_debug_Print 1 "LENGTH_NONPERIODIC = ${ensemble_Ens_PERTURB_LENGTH_NONPERIODIC[*]}"
111  IGCM_debug_Print 1 "MEMBER            = ${ensemble_Ens_PERTURB_MEMBER}"
112  IGCM_debug_Print 1 "MEMBER_LIST       = ${ensemble_Ens_PERTURB_MEMBER_LIST[*]}"
113  IGCM_debug_Print 1 "MEMBER_NAMESLIST  = ${ensemble_Ens_PERTURB_MEMBER_NAMESLIST[*]}"
114  IGCM_debug_Print 1 "MEMBER_INITFROM   = ${ensemble_Ens_PERTURB_MEMBER_INITFROM}"
115  IGCM_debug_Print 1 "MEMBER_INITPATH   = ${ensemble_Ens_PERTURB_MEMBER_INITPATH}"
116  IGCM_debug_Print 1 "PERTURB_BIN       = ${ensemble_Ens_PERTURB_PERTURB_BIN[*]}"
117  IGCM_debug_Print 1 "INITFROM          = ${ensemble_Ens_PERTURB_INITFROM}"
118  IGCM_debug_Print 1 "INITPATH          = ${ensemble_Ens_PERTURB_INITPATH}"
119  IGCM_debug_Print 1 "JobName           = ${config_UserChoices_JobName}"
120  IGCM_debug_Print 1 "TagName           = ${config_UserChoices_TagName}"
121  IGCM_debug_Print 1 "CalendarType      = ${config_UserChoices_CalendarType}"
122  IGCM_debug_Print 1 "ListOfComponents  = ${config_ListOfComponents[*]}"
123
124  PerturbExe=${ensemble_Ens_PERTURB_PERTURB_BIN[0]}
125
126  case ${PerturbExe} in
127  AddNoise)
128    PerturbComp=${ensemble_Ens_PERTURB_PERTURB_BIN[1]}
129    PerturbFile=${ensemble_Ens_PERTURB_PERTURB_BIN[2]}
130    PerturbVar=${ensemble_Ens_PERTURB_PERTURB_BIN[3]}
131    PerturbAmp=${ensemble_Ens_PERTURB_PERTURB_BIN[4]}
132
133    IGCM_debug_Print 1 "PerturbExe  = ${PerturbExe}"
134    IGCM_debug_Print 1 "PerturbFile = ${PerturbFile}"
135    IGCM_debug_Print 1 "PerturbComp = ${PerturbComp}"
136    IGCM_debug_Print 1 "PerturbVar  = ${PerturbVar}"
137    IGCM_debug_Print 1 "PerturbAmp  = ${PerturbAmp}"
138    ;;
139  AddPertu3DOCE)
140    PerturbComp=${ensemble_Ens_PERTURB_PERTURB_BIN[1]}
141    PerturbFile=${ensemble_Ens_PERTURB_PERTURB_BIN[2]}
142    PerturbVar=${ensemble_Ens_PERTURB_PERTURB_BIN[3]}
143    PerturbMask=${ensemble_Ens_PERTURB_PERTURB_BIN[4]}
144
145    IGCM_card_DefineVariableFromOption ${F_CFG_ENS} Ens_PERTURB MASKPATH
146
147    IGCM_debug_Print 1 "PerturbExe  = ${PerturbExe}"
148    IGCM_debug_Print 1 "PerturbFile = ${PerturbFile}"
149    IGCM_debug_Print 1 "PerturbComp = ${PerturbComp}"
150    IGCM_debug_Print 1 "PerturbVar  = ${PerturbVar}"
151    IGCM_debug_Print 1 "PerturbMask = ${PerturbMask}"
152    IGCM_debug_Print 1 "MASK PATH   = ${ensemble_Ens_PERTURB_MASKPATH}"
153    ;;
154  esac
155
156  IGCM_ensemble_SetAlpha ${ensemble_Ens_PERTURB_MEMBER}
157
158  # A few checks Period case:
159
160  # ... Check PERIODICITY ...
161  case ${ensemble_Ens_PERTURB_PERIODICITY} in
162  NONE)
163    IGCM_debug_Print 1 "periodic start not active"
164    CastPeriodicStart=false
165    ;;
166  *[Yy]|*[Mm])
167    CastPeriodicStart=true
168    IGCM_debug_Print 1 "Periodic length : ${ensemble_Ens_PERTURB_PERIODICITY}" ;;
169  *)
170    IGCM_debug_Exit "IGCM_ensemble_CastInit ${ensemble_Ens_PERTURB_PERIODICITY} : invalid PERIODICITY"
171    IGCM_debug_Exit "Choose a value in *Y or *M"
172    IGCM_debug_Verif_Exit ;;
173  esac
174  # ... Check LENGTH ...
175  case ${ensemble_Ens_PERTURB_LENGTH} in
176  *[Yy]|*[Mm])
177    IGCM_debug_Print 1 "Periodic duration : ${ensemble_Ens_PERTURB_LENGTH}" ;;
178  *)
179    IGCM_debug_Exit "IGCM_ensemble_CastInit ${ensemble_Ens_PERTURB_LENGTH} invalid LENGTH"
180    IGCM_debug_Exit "Choose a value in choose in *Y or *M"
181    IGCM_debug_Verif_Exit ;;
182  esac
183
184  # A few checks for the Non-Periodic case:
185  DateNum=0
186  while [ ${DateNum} -lt ${#ensemble_Ens_PERTURB_NONPERIODIC[*]} ] ; do
187
188    # - Check LENGTH_NONPERIODIC
189    case ${ensemble_Ens_PERTURB_LENGTH_NONPERIODIC[${DateNum}]} in
190    _0_)
191      IGCM_debug_Print 1 "non-periodic start not active"
192      CastNonPeriodicStart=false
193      ;;
194    *[Yy]|*[Mm])
195      IGCM_debug_Print 1 "Non-periodic duration : ${ensemble_Ens_PERTURB_LENGTH_NONPERIODIC[${DateNum}]}"
196      CastNonPeriodicStart=true
197      ;;
198    *)
199      IGCM_debug_Exit "IGCM_ensemble_CastInit ${ensemble_Ens_PERTURB_LENGTH_NONPERIODIC[${DateNum}]} : invalid LENGTH"
200      IGCM_debug_Exit "choose in *Y or *M"
201      IGCM_debug_Verif_Exit ;;
202    esac
203    (( DateNum = DateNum + 1 ))
204  done
205
206  # A few checks for the MEMBER_LIST case:
207  case ${ensemble_Ens_PERTURB_MEMBER_LIST[0]} in
208  _0_)
209    IGCM_debug_Print 1 "list of perturbation maps not active"
210    CastMemberList=false
211    ;;
212  *)
213    if [ ${CastPeriodicStart} = "true" ] ; then
214      IGCM_debug_Exit "list of perturbation maps for periodic start not implemented, will stop execution"
215      IGCM_debug_Verif_Exit
216    elif [ ${CastNonPeriodicStart} = "true" ] ; then
217      IGCM_debug_Exit "list of perturbation maps for non periodic start not implemented, will stop execution"
218      IGCM_debug_Verif_Exit
219    fi
220
221    # test that MEMBER_NAMESLIST and MEMBER_LIST have the same size
222
223    if [ ${#ensemble_Ens_PERTURB_MEMBER_LIST[*]} -ne ${#ensemble_Ens_PERTURB_MEMBER_NAMESLIST[*]} ] ; then
224      IGCM_debug_Exit "number of elements in MEMBER_LIST and MEMBER_NAMESLIST differ"
225      IGCM_debug_Verif_Exit
226    fi
227
228    IGCM_debug_Print 1 "list of perturbation maps : ${ensemble_Ens_PERTURB_MEMBER_LIST[*]}"
229    IGCM_debug_Print 1 "list of members names : ${ensemble_Ens_PERTURB_MEMBER_NAMESLIST[*]}"
230    CastMemberList=true
231    ;;
232  esac
233
234#  IGCM_debug_Exit "fin du test MEMBER_LIST"
235#  IGCM_debug_Verif_Exit
236
237  # Need to know all the restart filename of the component we will apply the noise to
238  IGCM_card_DefineArrayFromOption config.card ListOfComponents ${PerturbComp}
239  eval compname=\${config_ListOfComponents_${PerturbComp}[0]} > /dev/null 2>&1
240
241  # Target the component's card we apply the noise to
242  card=${SUBMIT_DIR}/COMP/${compname}.card
243
244  # Read the restart file list. To be used later
245  IGCM_card_DefineArrayFromOption ${card} RestartFiles List
246  ListFilesName=${compname}_RestartFiles_List
247  eval FileName0=\${${ListFilesName}[0]} > /dev/null 2>&1
248  eval NbFiles=\${#${ListFilesName}[@]} > /dev/null 2>&1
249
250  # Check
251  IGCM_debug_Print 1 "Nb Restart Files  = ${NbFiles}"
252
253  IGCM_debug_PopStack "IGCM_ensemble_CastInit"
254}
255
256function IGCM_ensemble_CastPeriodicStarts
257{
258  IGCM_debug_PushStack "IGCM_ensemble_CastPeriodicStarts"
259
260  [ ${CastPeriodicStart} = false ] && return
261
262  echo
263  IGCM_debug_Print 1 "Manage periodic starts"
264
265#.. Manage periodic starts ..
266#   ======================
267
268# ... Loop over DateBegin ...
269  eval DateBegin=\${ensemble_Ens_PERTURB_BEGIN_INIT}
270
271  while [ ${DateBegin} -le ${ensemble_Ens_PERTURB_END_INIT} ] ; do
272    IGCM_date_GetYearMonth ${DateBegin} year month
273
274  # - Determine number of day(s) in PERIODICITY
275    PeriodLengthInDays=$( IGCM_date_DaysInCurrentPeriod ${DateBegin} ${ensemble_Ens_PERTURB_PERIODICITY} )
276
277  # - Determine number of day(s) in LENGTH
278    DureeLengthInDays=$(( $( IGCM_date_DaysInCurrentPeriod ${DateBegin} ${ensemble_Ens_PERTURB_LENGTH} ) - 1 ))
279
280  # - Determine DateEnd
281    (( DateEnd = $( IGCM_date_AddDaysToGregorianDate ${DateBegin} ${DureeLengthInDays} ) ))
282
283  # - Build directory name
284    IGCM_ensemble_CastDirectoryName ${ensemble_Ens_PERTURB_NAME} ${ensemble_Ens_PERTURB_PERIODICITY} $year $month $StartDir
285
286  # - Determine RestartDate
287    (( Offset = -1 ))
288    (( RestartDate = $( IGCM_date_AddDaysToGregorianDate ${DateBegin} ${Offset} ) ))
289
290    IGCM_debug_Print 2 "${DateBegin} => ${DateEnd} : ${StartDir}"
291    echo "${DateBegin} ${DateEnd} ${StartDir}" >> ${RUN_DIR}/CreatedDir.txt
292
293  # - Create directory for current DateBegin
294    if [ ! -d  ${StartDir} ] ; then
295      IGCM_sys_Mkdir ${SUBMIT_DIR}/${StartDir}
296      IGCM_sys_Cd ${SUBMIT_DIR}/${StartDir}
297      ln -s ../../.resol .
298      ln -s ../../.libmpi .
299      IGCM_sys_Cd ${RUN_DIR}
300    fi
301
302  # - Create directory to store modified restart files
303    RestartDir=${STORAGE}/IGCM_IN/${config_UserChoices_TagName}/${StartDir}
304    IGCM_sys_MkdirArchive ${RestartDir}
305
306  # - Loop over members
307    i=0
308    while [ $i -lt ${ensemble_Ens_PERTURB_MEMBER} ] ; do
309      MemberDir="${StartDir}${Alpha[$i]}"
310      echo
311      IGCM_debug_Print 3 "${MemberDir}"
312
313      JobName="Job_${MemberDir}"
314
315    # * Create directory if it doesn't exist and copy/link files
316      if [ ! -d  ${SUBMIT_DIR}/${StartDir}/${MemberDir} ] ; then
317        IGCM_sys_Mkdir ${SUBMIT_DIR}/${StartDir}/${MemberDir}
318        IGCM_sys_Cd ${SUBMIT_DIR}/${StartDir}/${MemberDir}
319        ln -s ../../COMP
320        ln -s ../../PARAM
321        ln -s ../../POST
322        ln -s ../../DRIVER
323        IGCM_sys_Cd ${RUN_DIR}
324        IGCM_sys_Cp config.card run.card.init ${SUBMIT_DIR}/${StartDir}/${MemberDir}
325        IGCM_sys_Cp Job_${config_UserChoices_JobName} ${SUBMIT_DIR}/${StartDir}/${MemberDir}/${JobName}
326
327        # Dump command to be lauched
328        echo "cd ${StartDir}/${MemberDir}/ ;"  >> ${RUN_DIR}/Qsub.${StartDir}.sh
329        echo "${SUBMIT} ${JobName} ; cd -"     >> ${RUN_DIR}/Qsub.${StartDir}.sh
330
331        echo "cd ${StartDir}/${MemberDir}/ ;"  >> ${RUN_DIR}/Qclean.month.${StartDir}.sh
332        echo "${libIGCM}/clean_month.job ; cd -"     >> ${RUN_DIR}/Qclean.month.${StartDir}.sh
333
334        echo "cd ${StartDir}/${MemberDir}/ ;"  >> ${RUN_DIR}/Qclean.year.${StartDir}.sh
335        echo "${libIGCM}/clean_year.job ; cd -"     >> ${RUN_DIR}/Qclean.year.${StartDir}.sh
336
337        # * Update files : config.card, Job_, COMP/comp.card
338        IGCM_ensemble_CastFilesUpdate ${DateBegin} ${DateEnd} ${RestartDate}
339
340        # * Apply noise on restart file
341        IGCM_ensemble_CastPerturbFile
342      fi
343
344      (( i = i + 1 ))
345    done
346
347    # Done. Save ${StartDir} submission text file
348    IGCM_sys_Cp ${RUN_DIR}/Qsub.${StartDir}.sh ${SUBMIT_DIR}
349    IGCM_sys_Cp ${RUN_DIR}/Qclean.month.${StartDir}.sh ${SUBMIT_DIR}
350    IGCM_sys_Cp ${RUN_DIR}/Qclean.year.${StartDir}.sh ${SUBMIT_DIR}
351
352  # - Next DateBegin
353    echo "$DateBegin  $PeriodLengthInDays"
354    case ${ensemble_Ens_PERTURB_PERIODICITY} in
355    *[Yy]|*[Mm])
356      (( DateBegin = $( IGCM_date_AddDaysToGregorianDate ${DateBegin} ${PeriodLengthInDays} ) ))
357      ;;
358    esac
359    echo "New DateBegin = $DateBegin"
360    echo "========================================================================"
361  done
362  IGCM_debug_PopStack "IGCM_ensemble_CastPeriodicStarts"
363}
364
365function IGCM_ensemble_CastNonPeriodicStarts
366{
367  IGCM_debug_PushStack "IGCM_ensemble_CastNonPeriodicStarts"
368
369  #.. Manage non periodic starts => Loop over DateBegin ..
370  #   ==========================
371
372  [ ${CastNonPeriodicStart} = false ] && return
373
374  echo
375  IGCM_debug_Print 1 "Manage non periodic starts"
376
377  DateNum=0
378# ... Loop over ensemble_Ens_PERTURB_NONPERIODIC ...
379  echo ">${DateNum}<"
380  echo ">${#ensemble_Ens_PERTURB_NONPERIODIC[*]}<"
381  while [ ${DateNum} -lt ${#ensemble_Ens_PERTURB_NONPERIODIC[*]} ] ; do
382    DateBegin=${ensemble_Ens_PERTURB_NONPERIODIC[${DateNum}]}
383    Duree=${ensemble_Ens_PERTURB_LENGTH_NONPERIODIC[${DateNum}]}
384    echo ">${DateBegin}<"
385    echo ">${Duree}<"
386
387  # - Determine number of day(s) in LENGTH_NONPERIODIC
388    IGCM_date_GetYearMonth ${DateBegin} year month
389    DureeLengthInDays=$(( $( IGCM_date_DaysInCurrentPeriod ${DateBegin} ${Duree} ) - 1 ))
390
391  # - Build directory name
392    echo "========================================================================"
393    echo "ensemble_Ens_PERTURB_NAME = ${ensemble_Ens_PERTURB_NAME}"
394    IGCM_ensemble_CastDirectoryName ${ensemble_Ens_PERTURB_NAME} ${Duree} $year $month $StartDir
395
396  # - Determine DateEnd
397    (( DateEnd = $( IGCM_date_AddDaysToGregorianDate ${DateBegin} ${DureeLengthInDays} ) ))
398    echo ">${DateEnd}<"
399    echo "tout va bien 1"
400
401  # - Determine RestartDate
402    (( Offset = -1 ))
403    (( RestartDate = $( IGCM_date_AddDaysToGregorianDate ${DateBegin} ${Offset} ) ))
404
405    IGCM_debug_Print 2 "${DateBegin} => ${DateEnd} : ${StartDir}"
406
407  # -  Does $StartDir already exist ?
408    #echo "tout va bien 2" ${StartDir}
409    if [ ! -d ${SUBMIT_DIR}/${StartDir} ] ; then
410      echo "create dir"
411      IGCM_sys_Mkdir ${SUBMIT_DIR}/${StartDir}
412      IGCM_sys_Cd ${SUBMIT_DIR}/${StartDir}
413      ln -s ../../.resol .
414      ln -s ../../.libmpi .
415      IGCM_sys_Cd ${RUN_DIR}
416      echo "${DateBegin} ${DateEnd} ${StartDir}" >> ${RUN_DIR}/CreatedDir.txt
417    fi
418    PeriodDateEnd=$( grep -m1 ${StartDir} ${RUN_DIR}/CreatedDir.txt | cut -f2 -d\  )
419
420  # - Create directory in which to store new restart files if it does'nt already exist
421    RestartDir=${STORAGE}/IGCM_IN/${config_UserChoices_TagName}/${StartDir}
422    IGCM_sys_MkdirArchive ${RestartDir}
423
424  # - Loop over members
425    i=0
426    while [ $i -lt ${ensemble_Ens_PERTURB_MEMBER} ] ; do
427      MemberDir="${StartDir}${Alpha[$i]}"
428      IGCM_debug_Print 3 "${MemberDir}"
429
430      JobName="Job_${MemberDir}"
431
432    # * Create directory if it doesn't exist and copy files
433      if [ ! -d  ${SUBMIT_DIR}/${StartDir}/${MemberDir} ] ; then
434        IGCM_sys_Mkdir ${SUBMIT_DIR}/${StartDir}/${MemberDir}
435        #IGCM_sys_Cp -r COMP/ PARAM/ ${StartDir}/${MemberDir}
436        IGCM_sys_Cd ${SUBMIT_DIR}/${StartDir}/${MemberDir}
437        ln -s ../../COMP
438        ln -s ../../PARAM
439        ln -s ../../POST
440        ln -s ../../DRIVER
441        IGCM_sys_Cd ${RUN_DIR}
442        IGCM_sys_Cp config.card run.card.init ${SUBMIT_DIR}/${StartDir}/${MemberDir}
443        IGCM_sys_Cp Job_${config_UserChoices_JobName} ${SUBMIT_DIR}/${StartDir}/${MemberDir}/Job_${MemberDir}
444
445        # Dump command to be lauched
446        echo "cd ${StartDir}/${MemberDir}/ ;"  >> Qsub.${StartDir}.sh
447        echo "${SUBMIT} ${JobName} ; cd -"     >> Qsub.${StartDir}.sh
448
449        echo "cd ${StartDir}/${MemberDir}/ ;"  >> ${RUN_DIR}/Qclean.month.${StartDir}.sh
450        echo "${libIGCM}/clean_month.job ; cd -"     >> ${RUN_DIR}/Qclean.month.${StartDir}.sh
451
452        echo "cd ${StartDir}/${MemberDir}/ ;"  >> ${RUN_DIR}/Qclean.year.${StartDir}.sh
453        echo "${libIGCM}/clean_year.job ; cd -"     >> ${RUN_DIR}/Qclean.year.${StartDir}.sh
454
455        # * Update files : config.card, Job_, COMP/comp.card
456        echo "${PeriodDateEnd} ? ${DateEnd}"
457        if [ ${PeriodDateEnd} -gt ${DateEnd} ] ; then
458          DateEnd=${PeriodDateEnd}
459        fi
460        IGCM_ensemble_CastFilesUpdate ${DateBegin} ${DateEnd} ${RestartDate}
461
462        # * Apply noise on restart file
463        IGCM_ensemble_CastPerturbFile
464      fi
465
466      (( i = i + 1 ))
467    done
468
469    # Done. Save ${StartDir} submission text file
470    IGCM_sys_Cp ${RUN_DIR}/Qsub.${StartDir}.sh ${SUBMIT_DIR}
471    IGCM_sys_Cp ${RUN_DIR}/Qclean.month.${StartDir}.sh ${SUBMIT_DIR}
472    IGCM_sys_Cp ${RUN_DIR}/Qclean.year.${StartDir}.sh ${SUBMIT_DIR}
473
474    (( DateNum = DateNum + 1 ))
475  done
476  IGCM_debug_PopStack "IGCM_ensemble_CastNonPeriodicStarts"
477}
478
479function IGCM_ensemble_CastMemberList
480{
481  IGCM_debug_PushStack "IGCM_ensemble_CastMemberList"
482
483  if [ ${CastMemberList} = false ] ; then
484    IGCM_debug_PopStack "IGCM_ensemble_CastMemberList"
485    return
486  fi
487
488  echo
489  IGCM_debug_Print 1 "Manage members list"
490
491#.. Manage members list ..
492#   ======================
493
494  # DateBegin
495  eval DateBegin=\${ensemble_Ens_PERTURB_BEGIN_INIT}
496
497  IGCM_date_GetYearMonth ${DateBegin} year month
498
499  # - Determine number of day(s) in LENGTH
500  DureeLengthInDays=$(( $( IGCM_date_DaysInCurrentPeriod ${DateBegin} ${ensemble_Ens_PERTURB_LENGTH} ) - 1 ))
501
502  # - Determine DateEnd
503  DateEnd = $( IGCM_date_AddDaysToGregorianDate ${DateBegin} ${DureeLengthInDays} )
504
505  # bad hack enforce yearly for parent directory name
506  # - Build directory name
507  IGCM_ensemble_CastDirectoryName ${ensemble_Ens_PERTURB_NAME} 1Y $year $month $StartDir
508
509  # - Determine RestartDate
510  (( Offset = -1 ))
511  RestartDate = $( IGCM_date_AddDaysToGregorianDate ${DateBegin} ${Offset} )
512
513  IGCM_debug_Print 2 "${DateBegin} => ${DateEnd} : ${StartDir}"
514  echo "${DateBegin} ${DateEnd} ${StartDir}" >> ${RUN_DIR}/CreatedDir.txt
515
516  # - Create directory for current DateBegin
517  if [ ! -d  ${StartDir} ] ; then
518    IGCM_sys_Mkdir ${SUBMIT_DIR}/${StartDir}
519    IGCM_sys_Cd ${SUBMIT_DIR}/${StartDir}
520    ln -s ../../.resol .
521    ln -s ../../.libmpi .
522    IGCM_sys_Cd ${RUN_DIR}
523  fi
524
525  # - Create directory to store modified restart files
526  RestartDir=${STORAGE}/IGCM_IN/${config_UserChoices_TagName}/${StartDir}
527  IGCM_sys_MkdirArchive ${RestartDir}
528
529  # - Loop over members
530  i=0
531  nbmember=${#ensemble_Ens_PERTURB_MEMBER_LIST[*]}
532  while [ $i -lt $nbmember ] ; do
533    MemberDir=${ensemble_Ens_PERTURB_MEMBER_NAMESLIST[${i}]}
534    MemberVec=${ensemble_Ens_PERTURB_MEMBER_LIST[${i}]}
535
536    JobName="Job_${MemberDir}"
537    echo
538    IGCM_debug_Print 3 "${MemberDir}"
539
540    # * Create directory if it doesn't exist and copy/link files
541    if [ ! -d  ${SUBMIT_DIR}/${StartDir}/${MemberDir} ] ; then
542      IGCM_sys_Mkdir ${SUBMIT_DIR}/${StartDir}/${MemberDir}
543      IGCM_sys_Cd ${SUBMIT_DIR}/${StartDir}/${MemberDir}
544      ln -s ../../COMP
545      ln -s ../../PARAM
546      ln -s ../../POST
547      ln -s ../../DRIVER
548      IGCM_sys_Cd ${RUN_DIR}
549      IGCM_sys_Cp config.card run.card.init ${SUBMIT_DIR}/${StartDir}/${MemberDir}
550      IGCM_sys_Cp Job_${config_UserChoices_JobName} ${SUBMIT_DIR}/${StartDir}/${MemberDir}/${JobName}
551
552      # Dump command to be lauched
553      echo "cd ${StartDir}/${MemberDir}/ ;"     >> ${RUN_DIR}/Qsub.${StartDir}.sh
554      echo "${SUBMIT} ${JobName} ; cd -"        >> ${RUN_DIR}/Qsub.${StartDir}.sh
555
556      echo "cd ${StartDir}/${MemberDir}/ ;"     >> ${RUN_DIR}/Qclean.month.${StartDir}.sh
557      echo "${libIGCM}/clean_month.job ; cd -"  >> ${RUN_DIR}/Qclean.month.${StartDir}.sh
558
559      echo "cd ${StartDir}/${MemberDir}/ ;"     >> ${RUN_DIR}/Qclean.year.${StartDir}.sh
560      echo "${libIGCM}/clean_year.job ; cd -"   >> ${RUN_DIR}/Qclean.year.${StartDir}.sh
561
562      # * Update files : config.card, Job_, COMP/comp.card
563      IGCM_ensemble_CastFilesUpdate ${DateBegin} ${DateEnd} ${RestartDate}
564
565      # * Apply noise on restart file
566      IGCM_ensemble_CastPerturbFile
567    fi
568
569    (( i = i + 1 ))
570  done
571
572  # Done. Save ${StartDir} submission text file
573  IGCM_sys_Cp ${RUN_DIR}/Qsub.${StartDir}.sh ${SUBMIT_DIR}
574  IGCM_sys_Cp ${RUN_DIR}/Qclean.month.${StartDir}.sh ${SUBMIT_DIR}
575  IGCM_sys_Cp ${RUN_DIR}/Qclean.year.${StartDir}.sh ${SUBMIT_DIR}
576
577  IGCM_debug_PopStack "IGCM_ensemble_CastMemberList"
578}
579
580function IGCM_ensemble_CastFilesUpdate
581{
582  IGCM_debug_PushStack "IGCM_ensemble_CastFilesUpdate"
583
584  # Debug Print :
585  echo
586  IGCM_debug_Print 1 "IGCM_ensemble_CastFilesUpdate :"
587
588  HumanDateBegin=$(   IGCM_date_ConvertFormatToHuman ${1} )
589  HumanDateEnd=$(     IGCM_date_ConvertFormatToHuman ${2} )
590  HumanRestartDate=$( IGCM_date_ConvertFormatToHuman ${3} )
591  # ==> config.card
592  IGCM_card_WriteOption ${SUBMIT_DIR}/${StartDir}/${MemberDir}/config.card Ensemble EnsembleRun 'y'
593  IGCM_card_WriteOption ${SUBMIT_DIR}/${StartDir}/${MemberDir}/config.card Ensemble EnsembleName ${ensemble_Ens_PERTURB_NAME}
594  IGCM_card_WriteOption ${SUBMIT_DIR}/${StartDir}/${MemberDir}/config.card Ensemble EnsembleDate ${StartDir}
595
596  IGCM_card_WriteOption ${SUBMIT_DIR}/${StartDir}/${MemberDir}/config.card UserChoices JobName   ${MemberDir}
597  IGCM_card_WriteOption ${SUBMIT_DIR}/${StartDir}/${MemberDir}/config.card UserChoices DateBegin ${HumanDateBegin}
598  IGCM_card_WriteOption ${SUBMIT_DIR}/${StartDir}/${MemberDir}/config.card UserChoices DateEnd   ${HumanDateEnd}
599  IGCM_card_WriteOption ${SUBMIT_DIR}/${StartDir}/${MemberDir}/config.card Restarts OverRule "n"
600
601  for comp in ${config_ListOfComponents[*]} ; do
602    IGCM_card_WriteOption ${SUBMIT_DIR}/${StartDir}/${MemberDir}/config.card ${comp} Restart "y"
603    IGCM_card_WriteOption ${SUBMIT_DIR}/${StartDir}/${MemberDir}/config.card ${comp} RestartDate ${HumanRestartDate}
604    IGCM_card_WriteOption ${SUBMIT_DIR}/${StartDir}/${MemberDir}/config.card ${comp} RestartJobName ${ensemble_Ens_PERTURB_INITFROM}
605    IGCM_card_WriteOption ${SUBMIT_DIR}/${StartDir}/${MemberDir}/config.card ${comp} RestartPath ${ensemble_Ens_PERTURB_INITPATH}
606  done
607
608  # ==> Job
609  sed -e "s/\(#.*\)${config_UserChoices_JobName}\( *#.*\)/\1${MemberDir} \2/"            \
610      -e "s/\(#.*Script_Output_\)${config_UserChoices_JobName}\(\.*\)/\1${MemberDir}\2/" \
611      -e "s/^PeriodNb=.*/PeriodNb=60/"                                                   \
612      ${SUBMIT_DIR}/${StartDir}/${MemberDir}/Job_${MemberDir} > Job_${MemberDir}.tmp
613  IGCM_sys_Mv Job_${MemberDir}.tmp ${SUBMIT_DIR}/${StartDir}/${MemberDir}/Job_${MemberDir}
614
615  IGCM_debug_PopStack "IGCM_ensemble_CastFilesUpdate"
616}
617
618function IGCM_ensemble_CastDirectoryName
619{
620  IGCM_debug_PushStack "IGCM_ensemble_CastDirectoryName"
621
622  #.. Debug Print ..
623  echo
624  IGCM_debug_Print 1 "IGCM_ensemble_CastDirectoryName :"
625  echo
626
627  Name=$1
628  Duree=$2
629  year=$3
630  month=$4
631
632  # - Build directory name
633  case ${Duree} in
634  *Y|*y)
635    siecle="$( echo $year | cut -c1-2 )"
636    siecle=$( (( $siecle - 18 )) )
637    StartYear="${siecle}$( echo $year | cut -c3-4 )"
638    StartDir="${Name}${StartYear}"
639    ;;
640  *M|*m)
641    echo $month
642    siecle="$( echo $year | cut -c1-2 )"
643    siecle=$( (( $siecle - 18 )) )
644    StartYear="${siecle}$( echo $year | cut -c3-4 )"
645    StartMonth="${AlphaMonth[ (( 10#${month} - 1 )) ]}"
646    StartDir="${Name}${StartYear}${StartMonth}"
647    ;;
648  esac
649
650  IGCM_debug_PopStack "IGCM_ensemble_CastDirectoryName"
651}
652
653function IGCM_ensemble_CastPerturbFile
654{
655  IGCM_debug_PushStack "IGCM_ensemble_CastPerturbFile"
656
657  typeset i i_ j
658  typeset -Z4 j4
659  typeset file_out file_out_
660
661  #.. Debug Print ..
662  echo
663  IGCM_debug_Print 1 "IGCM_ensemble_CastPerturbFile :"
664
665  #.. FileIn ? => RestartDate ..
666  DirIn="${ensemble_Ens_PERTURB_INITPATH}/${ensemble_Ens_PERTURB_INITFROM}/${PerturbComp}/Restart"
667  DirInTar="${ensemble_Ens_PERTURB_INITPATH}/${ensemble_Ens_PERTURB_INITFROM}/RESTART"
668  FileIn="${ensemble_Ens_PERTURB_INITFROM}_${RestartDate}_${PerturbFile}"
669  DirOut="${RestartDir}/${MemberDir}/${PerturbComp}/Restart"
670
671  # * Create member restart directory
672  IGCM_sys_TestDirArchive ${DirOut}
673  RET=$?
674  if [ $RET -gt 0 ] ; then
675    IGCM_sys_MkdirArchive ${DirOut}
676  fi
677
678  FileOut="${MemberDir}_${RestartDate}_${PerturbFile}"
679  IGCM_debug_Print 1 "FileIn  = ${DirIn}/${FileIn}"
680  IGCM_debug_Print 1 "FileOut = ${DirOut}/${FileOut}.nc"
681
682  IGCM_sys_TestFileArchive ${DirOut}/${FileOut}.nc
683  RET=$?
684  if [ $RET -gt 0 ] ; then
685
686    # * Look for the restart file we apply the noise to
687
688    # restart file list pertaining to the component we apply the noise to
689    # but not being the precise restart file we will apply the noise to
690    unset OtherFileInList
691    # generic restart filename list (like flxat, sstoc, restart, ...)
692    unset OtherGenericList
693
694    if ( [ X${FileName0} != X${NULL_STR} ] && [ X${FileName0} != XNONE ] ) ; then
695      (( i=0 ))
696      until [ $i -ge ${NbFiles} ]; do
697
698        (( i_ = i+1 ))
699        eval file_out_=\${${ListFilesName}[$i_]} > /dev/null 2>&1
700        eval file_out=${file_out_}
701
702        generic_restart_file_name_out=$( basename ${file_out} .nc )
703
704        if [ ! ${generic_restart_file_name_out} = ${PerturbFile} ] ; then
705          set +A OtherFileInList ${OtherFileInList[*]} ${PerturbComp}_${ensemble_Ens_PERTURB_INITFROM}_${RestartDate}_${generic_restart_file_name_out}\*.nc
706          set +A OtherGenericList ${OtherGenericList[*]} ${generic_restart_file_name_out}
707        fi
708
709        (( i=i+3 ))
710      done
711    fi
712    # How many restart files other than the one we will apply the noise to
713    NbOtherFiles=${#OtherGenericList[*]}
714
715    ##########################
716    # TO BE A FUNCTION BEGIN #
717    ##########################
718
719    if [ $( IGCM_sys_TestFileBuffer ${DirIn}/${FileIn}*.nc ; echo $? ) = 0 ] ; then
720      IGCM_debug_Print 3 "Buffered restart"
721      Buffered=true
722      Archived=false
723      Tared=false
724      nb_restart_file=$(IGCM_sys_CountFileBuffer ${DirIn}/${FileIn}_????.nc)
725    elif [ $( IGCM_sys_TestFileArchive ${DirIn}/${FileIn}*.nc ; echo $? ) = 0 ] ; then
726      IGCM_debug_Print 3 "Archived restart"
727      Buffered=false
728      Archived=true
729      Tared=false
730      nb_restart_file=$(IGCM_sys_CountFileArchive ${DirIn}/${FileIn}_????.nc)
731    else
732      IGCM_debug_Print 3 "Tared restart"
733      Buffered=false
734      Archived=false
735      Tared=true
736
737      # Look for the tar file we want if we did not found it already
738      for PotentialTarFile in $( find ${ensemble_Ens_PERTURB_INITPATH}/${ensemble_Ens_PERTURB_INITFROM}/RESTART -name "${ensemble_Ens_PERTURB_INITFROM}_*restart*.tar" -print ) ; do
739        IsMatching=$( echo ${PotentialTarFile##*/} | sed "s:_restart::" | sed "s:^${ensemble_Ens_PERTURB_INITFROM}_::" | sed "s:\.tar$::" | gawk -F_ -v restartdate=${RestartDate} '{if (($1 <= restartdate) && ($2 >= restartdate)) {print $1"_"$2}}' )
740        if [ ! X${IsMatching} = X ] ; then
741          TarFileFound=${PotentialTarFile}
742          break
743        fi
744      done
745
746      # Extract relevant restart files
747      IGCM_debug_Print 1 "tar xvf ${TarFileFound} ${PerturbComp}_${FileIn}*.nc ${OtherFileInList[*]}"
748      tar xvf ${TarFileFound} ${PerturbComp}_${FileIn}*.nc ${OtherFileInList[*]}
749      nb_restart_file=$( IGCM_sys_CountFileBuffer ${PerturbComp}_${FileIn}_????.nc )
750    fi
751
752    # Move around and perturb restart files so as to be able to start hindcast/forecast simulation members
753    if [ ${nb_restart_file} -gt 1 ] ; then
754      j=0
755      until [ $j -ge ${nb_restart_file} ]; do
756        j4=${j}
757        if [ X${Buffered} = Xtrue ] ; then
758          IGCM_sys_GetBuffer ${DirIn}/${FileIn}_${j4}.nc ${RUN_DIR}/${FileOut}_${j4}.nc
759
760          cd ${DirOut}
761          for generic in ${OtherGenericList[*]} ; do
762            ln -s ${DirIn}/${ensemble_Ens_PERTURB_INITFROM}_${RestartDate}_${generic}.nc ${MemberDir}_${RestartDate}_${generic}_${j4}.nc
763          done
764          cd -
765
766        elif [ X${Archived} = Xtrue ] ; then
767          IGCM_sys_Get ${DirIn}/${FileIn}_${j4}.nc ${RUN_DIR}/${FileOut}_${j4}.nc
768
769          for generic in ${OtherGenericList[*]} ; do
770            IGCM_sys_RshArchive "cd ${DirOut} ; ln -s ${DirIn}/${ensemble_Ens_PERTURB_INITFROM}_${RestartDate}_${generic}.nc ${MemberDir}_${RestartDate}_${generic}_${j4}.nc"
771          done
772
773        elif [ X${Tared} = Xtrue ] ; then
774          IGCM_debug_Print 2 "IGCM_sys_Mv ${PerturbComp}_${FileIn}_${j4}.nc ${RUN_DIR}/${FileOut}_${j4}.nc"
775          IGCM_sys_Mv ${PerturbComp}_${FileIn}_${j4}.nc ${RUN_DIR}/${FileOut}_${j4}.nc
776
777          for generic in ${OtherGenericList[*]} ; do
778            IGCM_sys_Mv ${PerturbComp}_${ensemble_Ens_PERTURB_INITFROM}_${RestartDate}_${generic}.nc ${DirOut}/${MemberDir}_${RestartDate}_${generic}_${j4}.nc
779          done
780
781        fi
782        (( j=j+1 ))
783      done
784    else
785      if [ X${Buffered} = Xtrue ] ; then
786        IGCM_sys_GetBuffer ${DirIn}/${FileIn}.nc ${RUN_DIR}/${FileOut}.nc
787
788        cd ${DirOut}
789        for generic in ${OtherGenericList[*]} ; do
790          ln -s ${DirIn}/${ensemble_Ens_PERTURB_INITFROM}_${RestartDate}_${generic}.nc ${MemberDir}_${RestartDate}_${generic}.nc
791        done
792        cd -
793
794      elif [ X${Archived} = Xtrue ] ; then
795        IGCM_sys_Get ${DirIn}/${FileIn}.nc ${RUN_DIR}/${FileOut}.nc
796
797        for generic in ${OtherGenericList[*]} ; do
798          IGCM_sys_RshArchive "cd ${DirOut} ; ln -s ${DirIn}/${ensemble_Ens_PERTURB_INITFROM}_${RestartDate}_${generic}.nc ${MemberDir}_${RestartDate}_${generic}.nc"
799        done
800
801      elif [ X${Tared} = Xtrue ] ; then
802        IGCM_debug_Print 2 "IGCM_sys_Mv ${PerturbComp}_${FileIn}.nc ${RUN_DIR}/${FileOut}.nc"
803        IGCM_sys_Mv ${PerturbComp}_${FileIn}.nc ${RUN_DIR}/${FileOut}.nc
804
805        for generic in ${OtherGenericList[*]} ; do
806          IGCM_debug_Print 2 "IGCM_sys_Mv ${PerturbComp}_${ensemble_Ens_PERTURB_INITFROM}_${RestartDate}_${generic}.nc ${DirOut}/${MemberDir}_${RestartDate}_${generic}.nc"
807          IGCM_sys_Mv ${PerturbComp}_${ensemble_Ens_PERTURB_INITFROM}_${RestartDate}_${generic}.nc ${DirOut}/${MemberDir}_${RestartDate}_${generic}.nc
808        done
809
810      fi
811    fi
812
813    ########################
814    # TO BE A FUNCTION END #
815    ########################
816
817# treat the perturbation on different components by looking at the executable name
818
819    case ${PerturbExe} in
820    (AddNoise)
821      IGCM_sys_Chmod 644 ${RUN_DIR}/${FileOut}.nc
822
823      IGCM_debug_Print 1 "${PerturbExe} ${RUN_DIR}/${FileOut}.nc ${PerturbVar} ${PerturbAmp}"
824      echo
825
826      ${PerturbExe} ${RUN_DIR}/${FileOut}.nc ${PerturbVar} ${PerturbAmp}  > /dev/null 2>&1
827      if [ $? -ne 0 ] ; then
828        IGCM_debug_Exit "Abend $( basename ${PerturbExe} )"
829        IGCM_debug_Verif_Exit
830      fi
831      IGCM_sys_Put_Out ${RUN_DIR}/${FileOut}.nc ${DirOut}/ 644
832      ;;
833    (AddPertu3DOCE)
834      # where to find the pattern we apply to the restart
835      PatternFile=${ensemble_Ens_PERTURB_MEMBER_INITPATH}/${ensemble_Ens_PERTURB_MEMBER_INITFROM}/${MemberVec}.nc
836
837      # where to find the land mask for the grid
838      MaskFile=${ensemble_Ens_PERTURB_MASKPATH}/${PerturbMask}
839
840      # if there is multiple restart files rebuild restart file
841      if [ ${nb_restart_file} -gt 1 ] ; then
842        IGCM_debug_Print 1 "rebuild files ${FileOut}_????.nc"
843        IGCM_sys_rebuild ${RUN_DIR}/${FileOut}.nc ${RUN_DIR}/${FileOut}_????.nc
844      fi
845
846      # there is now a single restart file
847      IGCM_sys_Chmod 644 ${RUN_DIR}/${FileOut}.nc
848
849      IGCM_debug_Print 1 "${PerturbExe} ${RUN_DIR}/${FileOut}.nc ${PerturbVar} ${PatternFile} ${MaskFile}"
850      echo
851
852      # add pattern to restart file on variable PerturbVar
853      ${PerturbExe} ${RUN_DIR}/${FileOut}.nc ${PerturbVar} ${PatternFile} ${MaskFile}  > /dev/null 2>&1
854      if [ $? -ne 0 ] ; then
855        IGCM_debug_Exit "Abend $( basename ${PerturbExe} )"
856        IGCM_debug_Verif_Exit
857      fi
858      IGCM_sys_Put_Out ${RUN_DIR}/${FileOut}.nc ${DirOut}/ 644
859     ;;
860    esac
861
862  fi
863
864  #.. Update config.card..
865  IGCM_card_WriteOption ${SUBMIT_DIR}/${StartDir}/${MemberDir}/config.card ${PerturbComp} Restart "y"
866  IGCM_card_WriteOption ${SUBMIT_DIR}/${StartDir}/${MemberDir}/config.card ${PerturbComp} RestartDate    ${HumanRestartDate}
867  IGCM_card_WriteOption ${SUBMIT_DIR}/${StartDir}/${MemberDir}/config.card ${PerturbComp} RestartJobName ${MemberDir}
868  IGCM_card_WriteOption ${SUBMIT_DIR}/${StartDir}/${MemberDir}/config.card ${PerturbComp} RestartPath    ${RestartDir}/
869  IGCM_card_WriteOption ${SUBMIT_DIR}/${StartDir}/${MemberDir}/config.card Ensemble EnsembleType "Ens_PERTURB"
870
871  IGCM_debug_PopStack "IGCM_ensemble_CastPerturbFile"
872}
873
874############### Date ENSEMBLE #################
875function IGCM_ensemble_DateInit
876{
877  IGCM_debug_PushStack "IGCM_ensemble_DateInit"
878
879  IGCM_sys_Mkdir ${RUN_DIR}
880
881  IGCM_sys_Cp ${SUBMIT_DIR}/config.card   ${RUN_DIR}
882  IGCM_sys_Cp ${SUBMIT_DIR}/ensemble.card ${RUN_DIR}
883  IGCM_sys_Cp ${SUBMIT_DIR}/Job_*         ${RUN_DIR}
884  IGCM_sys_Cp ${SUBMIT_DIR}/run.card.init ${RUN_DIR}
885  if [ -f ${SUBMIT_DIR}/Qsub.* ]; then
886    IGCM_sys_Cp ${SUBMIT_DIR}/Qsub.*        ${RUN_DIR}
887  fi
888  if [ -f ${SUBMIT_DIR}/Qclean.* ]; then
889    IGCM_sys_Cp ${SUBMIT_DIR}/Qclean.*      ${RUN_DIR}
890  fi
891
892  IGCM_card_DefineVariableFromOption ${F_CFG_ENS} Ens_DATE active
893  IGCM_card_DefineVariableFromOption ${F_CFG_ENS} Ens_DATE NAME
894  IGCM_card_DefineVariableFromOption ${F_CFG_ENS} Ens_DATE BEGIN_INIT
895  IGCM_card_DefineVariableFromOption ${F_CFG_ENS} Ens_DATE BEGIN_RESTART
896  IGCM_card_DefineVariableFromOption ${F_CFG_ENS} Ens_DATE END_INIT
897  IGCM_card_DefineVariableFromOption ${F_CFG_ENS} Ens_DATE PERIODICITY
898  IGCM_card_DefineArrayFromOption    ${F_CFG_ENS} Ens_DATE NONPERIODIC
899  IGCM_card_DefineArrayFromOption    ${F_CFG_ENS} Ens_DATE RESTART_NONPERIODIC
900  IGCM_card_DefineVariableFromOption ${F_CFG_ENS} Ens_DATE LENGTH
901  IGCM_card_DefineArrayFromOption    ${F_CFG_ENS} Ens_DATE LENGTH_NONPERIODIC
902  IGCM_card_DefineVariableFromOption ${F_CFG_ENS} Ens_DATE INITFROM
903  IGCM_card_DefineVariableFromOption ${F_CFG_ENS} Ens_DATE INITFROM_NONPERIODIC
904  IGCM_card_DefineVariableFromOption ${F_CFG_ENS} Ens_DATE INITPATH
905  IGCM_card_DefineVariableFromOption config.card UserChoices JobName
906  IGCM_card_DefineVariableFromOption config.card UserChoices TagName
907  IGCM_card_DefineVariableFromOption config.card UserChoices CalendarType
908  IGCM_card_DefineArrayFromSection   config.card ListOfComponents
909
910  echo
911  IGCM_debug_Print 1 "[Ens_DATE]"
912  IGCM_debug_Print 1 "ACTIVE               = ${ensemble_Ens_DATE_active}"
913  IGCM_debug_Print 1 "NAME                 = ${ensemble_Ens_DATE_NAME}"
914  IGCM_debug_Print 1 "BEGIN_INIT           = ${ensemble_Ens_DATE_BEGIN_INIT}"
915  IGCM_debug_Print 1 "END_INIT             = ${ensemble_Ens_DATE_END_INIT}"
916  IGCM_debug_Print 1 "PERIODICITY          = ${ensemble_Ens_DATE_PERIODICITY}"
917  IGCM_debug_Print 1 "BEGIN_RESTART        = ${ensemble_Ens_DATE_BEGIN_RESTART}"
918  IGCM_debug_Print 1 "NONPERIODIC          = ${ensemble_Ens_DATE_NONPERIODIC[*]}"
919  IGCM_debug_Print 1 "RESTART_NONPERIODIC  = ${ensemble_Ens_DATE_RESTART_NONPERIODIC[*]}"
920  IGCM_debug_Print 1 "LENGTH               = ${ensemble_Ens_DATE_LENGTH}"
921  IGCM_debug_Print 1 "LENGTH_NONPERIODIC   = ${ensemble_Ens_DATE_LENGTH_NONPERIODIC[*]}"
922  IGCM_debug_Print 1 "INITFROM             = ${ensemble_Ens_DATE_INITFROM}"
923  IGCM_debug_Print 1 "INITFROM_NONPERIODIC = ${ensemble_Ens_DATE_INITFROM_NONPERIODIC[*]}"
924  IGCM_debug_Print 1 "INITPATH             = ${ensemble_Ens_DATE_INITPATH}"
925  IGCM_debug_Print 1 "JobName              = ${config_UserChoices_JobName}"
926  IGCM_debug_Print 1 "TagName              = ${config_UserChoices_TagName}"
927  IGCM_debug_Print 1 "CalendarType         = ${config_UserChoices_CalendarType}"
928  IGCM_debug_Print 1 "ListOfComponents     = ${config_ListOfComponents[*]}"
929  echo ""
930
931  ensemble_Ens_DATE_MEMBER=1 # actually use only 1 member
932  IGCM_ensemble_SetAlpha ${ensemble_Ens_DATE_MEMBER}
933
934  IGCM_debug_Print 1 "Check args..."
935  DatePeriodicStart=false
936  DateNonPeriodicStart=false
937
938  # ... Check LENGTH ...
939  case ${ensemble_Ens_DATE_LENGTH} in
940  *[Yy]|*[Mm])
941    IGCM_debug_Print 1 "Default simulation duration : ${ensemble_Ens_DATE_LENGTH}" ;;
942  *)
943    IGCM_debug_Exit "IGCM_ensemble_DateInit ${ensemble_Ens_DATE_LENGTH} invalid LENGTH"
944    IGCM_debug_Exit "Choose a value in choose in *Y or *M"
945    IGCM_debug_Verif_Exit ;;
946  esac
947
948  # ***************************************
949  # A few checks Period case:
950  # ***************************************
951  # if all Periodic params are not filled: desactivate Periodic mode
952  totalPeriodArgs=4
953  periodFillArgs=0
954
955  if [[ X${ensemble_Ens_DATE_BEGIN_RESTART} != "X" ]]; then
956    (( periodFillArgs = periodFillArgs + 1 ))
957  fi
958
959  if [[ X${ensemble_Ens_DATE_BEGIN_INIT} != "X" ]]; then
960    (( periodFillArgs = periodFillArgs + 1 ))
961  fi
962
963  if [[ X${ensemble_Ens_DATE_END_INIT} != "X" ]]; then
964    (( periodFillArgs = periodFillArgs + 1 ))
965  fi
966
967  if [[ X${ensemble_Ens_DATE_PERIODICITY} != "X" ]]; then
968    (( periodFillArgs = periodFillArgs + 1 ))
969
970    # ... Check PERIODICITY ...
971    case ${ensemble_Ens_DATE_PERIODICITY} in
972    *[Yy]|*[Mm])
973      IGCM_debug_Print 1 "Periodic length : ${ensemble_Ens_DATE_PERIODICITY}" ;;
974    *)
975      IGCM_debug_Exit "IGCM_ensemble_DateInit ${ensemble_Ens_DATE_PERIODICITY} : invalid PERIODICITY"
976      IGCM_debug_Exit "Choose a value in *Y or *M"
977      IGCM_debug_Verif_Exit ;;
978    esac
979  fi # if periodicity
980
981  if [[ ${periodFillArgs} = ${totalPeriodArgs} ]]; then
982    DatePeriodicStart=true
983  else
984    if [[ ${periodFillArgs} = 0 ]]; then
985      IGCM_debug_Print 1 "Periodic start NOT ACTIVE"
986      DatePeriodicStart=false
987    else
988      IGCM_debug_Exit "IGCM_ensemble_DateInit missing arguments for Periodic mode!"
989      IGCM_debug_Exit "Get only ${periodFillArgs} on ${totalPeriodArgs} args. Check ${F_CFG_ENS} file."
990      IGCM_debug_Verif_Exit
991    fi
992  fi
993
994  # ***************************************
995  # A few checks for the Non-Periodic case:
996  # ***************************************
997  if [[ ${#ensemble_Ens_DATE_NONPERIODIC[*]} != ${#ensemble_Ens_DATE_RESTART_NONPERIODIC[*]} ]] ; then
998    IGCM_debug_Exit "IGCM_ensemble_DateInit: NONPERIODIC and RESTART_NONPERIODIC lists have different sizes"
999    IGCM_debug_Verif_Exit
1000  fi
1001
1002  if [[ ${#ensemble_Ens_DATE_NONPERIODIC[*]} > 0 ]] && [[ ${ensemble_Ens_DATE_NONPERIODIC[*]} != _0_ ]]; then
1003    DateNonPeriodicStart=true
1004
1005    # Use LENGTH if no NONPERIODIC_LENGTH given
1006    if [[ ${#ensemble_Ens_DATE_LENGTH_NONPERIODIC[*]} < ${#ensemble_Ens_DATE_NONPERIODIC[*]} ]] ; then
1007      IGCM_debug_Print 1 "WARNING: LENGTH_NONPERIODIC is not fill (or not correctly). Use LENGTH value '${ensemble_Ens_DATE_LENGTH}' for all NONPERIODIC runs"
1008      DateNum=0
1009      while [ ${DateNum} -lt ${#ensemble_Ens_DATE_NONPERIODIC[*]} ] ; do
1010        ensemble_Ens_DATE_LENGTH_NONPERIODIC[${DateNum}]=${ensemble_Ens_DATE_LENGTH}
1011        (( DateNum = DateNum + 1 ))
1012      done
1013    fi
1014
1015    # Use INITFROM if no NONPERIODIC_INITFROM given
1016    if [ ${#ensemble_Ens_DATE_INITFROM_NONPERIODIC[*]} -lt ${#ensemble_Ens_DATE_NONPERIODIC[*]} ] ; then
1017      IGCM_debug_Print 1 "WARNING: INITFROM_NONPERIODIC is not fill (or not correctly). Use INITFROM value '${ensemble_Ens_DATE_INITFROM}' for all NONPERIODIC runs"
1018      DateNum=0
1019      while [ ${DateNum} -lt ${#ensemble_Ens_DATE_NONPERIODIC[*]} ] ; do
1020        ensemble_Ens_DATE_INITFROM_NONPERIODIC[${DateNum}]=${ensemble_Ens_DATE_INITFROM}
1021        (( DateNum = DateNum + 1 ))
1022      done
1023    fi
1024  else
1025    IGCM_debug_Print 1 "Non-Periodic start NOT ACTIVE"
1026    DateNonPeriodicStart=false
1027  fi
1028
1029  if [[ ${DateNonPeriodicStart} = true ]]; then
1030    DateNum=0
1031    while [ ${DateNum} -lt ${#ensemble_Ens_DATE_NONPERIODIC[*]} ] ; do
1032      # - Check LENGTH_NONPERIODIC
1033      case ${ensemble_Ens_DATE_LENGTH_NONPERIODIC[${DateNum}]} in
1034      *[Yy]|*[Mm])
1035        IGCM_debug_Print 1 "Non-periodic duration ${DateNum}: ${ensemble_Ens_DATE_LENGTH_NONPERIODIC[${DateNum}]}"
1036        ;;
1037      *)
1038        IGCM_debug_Exit "IGCM_ensemble_DateInit ${ensemble_Ens_DATE_LENGTH_NONPERIODIC[${DateNum}]} : invalid NON PERIODIC LENGTH"
1039        IGCM_debug_Exit "choose in *Y or *M"
1040        IGCM_debug_Verif_Exit ;;
1041      esac
1042
1043      # - Check RESTART_NONPERIODIC
1044      case ${ensemble_Ens_DATE_RESTART_NONPERIODIC[${DateNum}]} in
1045      _0_)
1046        IGCM_debug_Exit "IGCM_ensemble_DateInit ${ensemble_Ens_DATE_RESTART_NONPERIODIC[${DateNum}]} : invalid NON PERIODIC RESTART"
1047        IGCM_debug_Verif_Exit ;;
1048      esac
1049
1050      (( DateNum = DateNum + 1 ))
1051    done
1052  fi # DateNonPeriodicStart = true
1053
1054  IGCM_debug_PopStack "IGCM_ensemble_DateInit"
1055}
1056
1057function IGCM_ensemble_DatePeriodicStarts
1058{
1059  IGCM_debug_PushStack "IGCM_ensemble_DatePeriodicStarts"
1060
1061  [ ${DatePeriodicStart} = false ] && return
1062
1063  echo
1064  IGCM_debug_Print 1 ">>>  MANAGE PERIODIC STARTS  <<<"
1065
1066#.. Manage periodic starts ..
1067#   ======================
1068
1069  # - Build directory name
1070  StartDir="${ensemble_Ens_DATE_NAME}"
1071
1072  # - Create directory for current DateBegin
1073  if [ ! -d  ${StartDir} ] ; then
1074    IGCM_sys_Mkdir ${SUBMIT_DIR}/${StartDir}
1075    IGCM_sys_Cd ${SUBMIT_DIR}/${StartDir}
1076    ln -s ../../.resol .
1077    ln -s ../../.libmpi .
1078    IGCM_sys_Cd ${RUN_DIR}
1079  fi
1080
1081  # ... Loop over DateBegin ...
1082  eval DateBegin=\${ensemble_Ens_DATE_BEGIN_INIT}
1083  eval RestartDate=\${ensemble_Ens_DATE_BEGIN_RESTART}
1084
1085  DateNum=0
1086  while [ ${DateBegin} -le ${ensemble_Ens_DATE_END_INIT} ] ; do
1087    IGCM_date_GetYearMonth ${DateBegin} year month
1088
1089    echo "========================================================================"
1090    echo "New DateBegin = $DateBegin"
1091
1092  # - Determine number of day(s) in PERIODICITY
1093    PeriodLengthInDays=$( IGCM_date_DaysInCurrentPeriod ${DateBegin} ${ensemble_Ens_DATE_PERIODICITY} )
1094
1095  # - Determine number of day(s) in LENGTH
1096    DureeLengthInDays=$(( $( IGCM_date_DaysInCurrentPeriod ${DateBegin} ${ensemble_Ens_DATE_LENGTH} ) - 1 ))
1097
1098  # - Determine DateEnd
1099    (( DateEnd = $( IGCM_date_AddDaysToGregorianDate ${DateBegin} ${DureeLengthInDays} ) ))
1100
1101    IGCM_debug_Print 2 "${DateBegin} => ${DateEnd} : ${StartDir}"
1102    echo "${DateBegin} ${DateEnd} ${StartDir}" >> ${RUN_DIR}/CreatedDir.txt
1103
1104  # - Loop over members (default =1 no member)
1105    i=0
1106    while [ $i -lt ${ensemble_Ens_DATE_MEMBER} ] ; do
1107      MemberDir="${ensemble_Ens_DATE_NAME}${DateNum}${Alpha[$i]}_per"
1108      echo
1109      IGCM_debug_Print 3 "${MemberDir}"
1110
1111      JobName="Job_${MemberDir}"
1112
1113    # * Create directory if it doesn't exist and copy/link files
1114      if [ ! -d  ${SUBMIT_DIR}/${StartDir}/${MemberDir} ] ; then
1115        IGCM_sys_Mkdir ${SUBMIT_DIR}/${StartDir}/${MemberDir}
1116        IGCM_sys_Cd ${SUBMIT_DIR}/${StartDir}/${MemberDir}
1117        ln -s ../../COMP
1118        ln -s ../../PARAM
1119        ln -s ../../POST
1120        ln -s ../../DRIVER
1121        IGCM_sys_Cd ${RUN_DIR}
1122        IGCM_sys_Cp config.card run.card.init ${SUBMIT_DIR}/${StartDir}/${MemberDir}
1123        IGCM_sys_Cp Job_${config_UserChoices_JobName} ${SUBMIT_DIR}/${StartDir}/${MemberDir}/${JobName}
1124
1125        # Dump command to be lauched
1126        echo "cd ${StartDir}/${MemberDir}/ ;"  >> ${RUN_DIR}/Qsub.${StartDir}.sh
1127        echo "${SUBMIT} ${JobName} ; cd -"     >> ${RUN_DIR}/Qsub.${StartDir}.sh
1128
1129        echo "cd ${StartDir}/${MemberDir}/ ;"  >> ${RUN_DIR}/Qclean.month.${StartDir}.sh
1130        echo "${libIGCM}/clean_month.job ; cd -"     >> ${RUN_DIR}/Qclean.month.${StartDir}.sh
1131
1132        echo "cd ${StartDir}/${MemberDir}/ ;"  >> ${RUN_DIR}/Qclean.year.${StartDir}.sh
1133        echo "${libIGCM}/clean_year.job ; cd -"     >> ${RUN_DIR}/Qclean.year.${StartDir}.sh
1134
1135        # * Update files : config.card, Job_, COMP/comp.card
1136        IGCM_ensemble_DateFilesUpdate ${DateBegin} ${DateEnd} ${RestartDate}
1137      fi
1138
1139      (( i = i + 1 ))
1140    done
1141
1142  # - Next DateBegin & RestartDate
1143    echo "$DateBegin  $PeriodLengthInDays"
1144    case ${ensemble_Ens_DATE_PERIODICITY} in
1145    *[Yy]|*[Mm])
1146      (( DateBegin =   $( IGCM_date_AddDaysToGregorianDate ${DateBegin}   ${PeriodLengthInDays} ) ))
1147      (( RestartDate = $( IGCM_date_AddDaysToGregorianDate ${RestartDate} ${PeriodLengthInDays} ) ))
1148      ;;
1149    esac
1150
1151    (( DateNum = DateNum + 1 )) # increment number of restart date
1152  done
1153
1154  # Done. Save ${StartDir} submission text file
1155  IGCM_sys_Cp ${RUN_DIR}/Qsub.${StartDir}.sh ${SUBMIT_DIR}
1156  IGCM_sys_Cp ${RUN_DIR}/Qclean.month.${StartDir}.sh ${SUBMIT_DIR}
1157  IGCM_sys_Cp ${RUN_DIR}/Qclean.year.${StartDir}.sh ${SUBMIT_DIR}
1158
1159  IGCM_debug_PopStack "IGCM_ensemble_DatePeriodicStarts"
1160}
1161
1162function IGCM_ensemble_DateNonPeriodicStarts
1163{
1164  IGCM_debug_PushStack "IGCM_ensemble_DateNonPeriodicStarts"
1165
1166  #.. Manage non periodic starts => Loop over DateBegin ..
1167  #   ==========================
1168
1169  [ ${DateNonPeriodicStart} = false ] && return
1170
1171  echo
1172  IGCM_debug_Print 1 ">>>  MANAGE NON PERIODIC STARTS  <<<"
1173
1174  # - Build directory name
1175  echo ""
1176  echo "========================================================================"
1177  echo "ensemble_Ens_DATE_NAME = ${ensemble_Ens_DATE_NAME}"
1178  StartDir="${ensemble_Ens_DATE_NAME}"
1179
1180  # -  Does $StartDir already exist ?
1181  if [ ! -d ${SUBMIT_DIR}/${StartDir} ] ; then
1182    IGCM_sys_Mkdir ${SUBMIT_DIR}/${StartDir}
1183    IGCM_sys_Cd ${SUBMIT_DIR}/${StartDir}
1184    ln -s ../../.resol .
1185    ln -s ../../.libmpi .
1186    IGCM_sys_Cd ${RUN_DIR}
1187  fi
1188
1189  DateNum=0
1190  # ... Loop over ensemble_Ens_DATE_NONPERIODIC ...
1191  while [ ${DateNum} -lt ${#ensemble_Ens_DATE_NONPERIODIC[*]} ] ; do
1192    DateBegin=${ensemble_Ens_DATE_NONPERIODIC[${DateNum}]}
1193    Duree=${ensemble_Ens_DATE_LENGTH_NONPERIODIC[${DateNum}]}
1194    RestartDate=${ensemble_Ens_DATE_RESTART_NONPERIODIC[${DateNum}]}
1195    InitFrom=${ensemble_Ens_DATE_INITFROM_NONPERIODIC[${DateNum}]}
1196
1197  # - Determine number of day(s) in LENGTH_NONPERIODIC
1198    IGCM_date_GetYearMonth ${DateBegin} year month
1199    DureeLengthInDays=$(( $( IGCM_date_DaysInCurrentPeriod ${DateBegin} ${Duree} ) - 1 ))
1200
1201  # - Determine DateEnd
1202    (( DateEnd = $( IGCM_date_AddDaysToGregorianDate ${DateBegin} ${DureeLengthInDays} ) ))
1203
1204    IGCM_debug_Print 2 "${DateBegin} => ${DateEnd} : ${StartDir}"
1205    echo "${DateBegin} ${DateEnd} ${StartDir}" >> ${RUN_DIR}/CreatedDir.txt
1206
1207    PeriodDateEnd=$( grep -m1 ${StartDir} ${RUN_DIR}/CreatedDir.txt | cut -f2 -d\  )
1208
1209  # - Loop over members
1210    i=0
1211    while [ $i -lt ${ensemble_Ens_DATE_MEMBER} ] ; do
1212      MemberDir="${ensemble_Ens_DATE_NAME}${DateNum}${Alpha[$i]}"
1213      IGCM_debug_Print 3 "${MemberDir}"
1214
1215      JobName="Job_${MemberDir}"
1216
1217    # * Create directory if it doesn't exist and copy files
1218      if [ ! -d  ${SUBMIT_DIR}/${StartDir}/${MemberDir} ] ; then
1219        IGCM_sys_Mkdir ${SUBMIT_DIR}/${StartDir}/${MemberDir}
1220        IGCM_sys_Cd ${SUBMIT_DIR}/${StartDir}/${MemberDir}
1221        ln -s ../../COMP
1222        ln -s ../../PARAM
1223        ln -s ../../POST
1224        ln -s ../../DRIVER
1225        IGCM_sys_Cd ${RUN_DIR}
1226        IGCM_sys_Cp config.card run.card.init ${SUBMIT_DIR}/${StartDir}/${MemberDir}
1227        IGCM_sys_Cp Job_${config_UserChoices_JobName} ${SUBMIT_DIR}/${StartDir}/${MemberDir}/Job_${MemberDir}
1228
1229        # Dump command to be lauched
1230        echo "cd ${StartDir}/${MemberDir}/ ;"  >> ${RUN_DIR}/Qsub.${StartDir}.sh
1231        echo "${SUBMIT} ${JobName} ; cd -"     >> ${RUN_DIR}/Qsub.${StartDir}.sh
1232
1233        echo "cd ${StartDir}/${MemberDir}/ ;"  >> ${RUN_DIR}/Qclean.month.${StartDir}.sh
1234        echo "${libIGCM}/clean_month.job ; cd -"     >> ${RUN_DIR}/Qclean.month.${StartDir}.sh
1235
1236        echo "cd ${StartDir}/${MemberDir}/ ;"  >> ${RUN_DIR}/Qclean.year.${StartDir}.sh
1237        echo "${libIGCM}/clean_year.job ; cd -"     >> ${RUN_DIR}/Qclean.year.${StartDir}.sh
1238
1239        # * Update files : config.card, Job_, COMP/comp.card
1240        IGCM_ensemble_DateFilesUpdate ${DateBegin} ${DateEnd} ${RestartDate} ${InitFrom}
1241      fi
1242
1243      (( i = i + 1 ))
1244    done
1245
1246    # Done. Save ${StartDir} submission text file
1247    IGCM_sys_Cp ${RUN_DIR}/Qsub.${StartDir}.sh ${SUBMIT_DIR}
1248    IGCM_sys_Cp ${RUN_DIR}/Qclean.month.${StartDir}.sh ${SUBMIT_DIR}
1249    IGCM_sys_Cp ${RUN_DIR}/Qclean.year.${StartDir}.sh ${SUBMIT_DIR}
1250
1251    (( DateNum = DateNum + 1 ))
1252  done
1253
1254  # Done. Save ${StartDir} submission text file
1255  IGCM_sys_Cp ${RUN_DIR}/Qsub.${StartDir}.sh ${SUBMIT_DIR}
1256  IGCM_sys_Cp ${RUN_DIR}/Qclean.month.${StartDir}.sh ${SUBMIT_DIR}
1257  IGCM_sys_Cp ${RUN_DIR}/Qclean.year.${StartDir}.sh ${SUBMIT_DIR}
1258
1259  IGCM_debug_PopStack "IGCM_ensemble_DateNonPeriodicStarts"
1260}
1261
1262function IGCM_ensemble_DateFilesUpdate
1263{
1264  IGCM_debug_PushStack "IGCM_ensemble_DateFilesUpdate"
1265
1266  # Debug Print :
1267  echo
1268  IGCM_debug_Print 1 "IGCM_ensemble_DateFilesUpdate :"
1269
1270  HumanDateBegin=$(   IGCM_date_ConvertFormatToHuman ${1} )
1271  HumanDateEnd=$(     IGCM_date_ConvertFormatToHuman ${2} )
1272  HumanRestartDate=$( IGCM_date_ConvertFormatToHuman ${3} )
1273  if [[ X${4} != "X" ]]; then
1274    initFrom=${4} # non periodic config (INITFROM could be different between members)
1275  else
1276    initFrom=${ensemble_Ens_DATE_INITFROM} # periodic (same INITFROM value)
1277  fi
1278
1279  # ==> config.card
1280  # [ENSEMBLE]
1281  IGCM_card_WriteOption ${SUBMIT_DIR}/${StartDir}/${MemberDir}/config.card Ensemble EnsembleRun 'y'
1282  IGCM_card_WriteOption ${SUBMIT_DIR}/${StartDir}/${MemberDir}/config.card Ensemble EnsembleName ${ensemble_Ens_DATE_NAME}
1283  IGCM_card_WriteOption ${SUBMIT_DIR}/${StartDir}/${MemberDir}/config.card Ensemble EnsembleDate ${HumanDateBegin}
1284  IGCM_card_WriteOption ${SUBMIT_DIR}/${StartDir}/${MemberDir}/config.card Ensemble EnsembleType "Ens_DATE"
1285
1286  # [UserChoices]
1287  IGCM_card_WriteOption ${SUBMIT_DIR}/${StartDir}/${MemberDir}/config.card UserChoices JobName   ${MemberDir}
1288  IGCM_card_WriteOption ${SUBMIT_DIR}/${StartDir}/${MemberDir}/config.card UserChoices DateBegin ${HumanDateBegin}
1289  IGCM_card_WriteOption ${SUBMIT_DIR}/${StartDir}/${MemberDir}/config.card UserChoices DateEnd   ${HumanDateEnd}
1290
1291  # [Restarts]
1292  IGCM_card_WriteOption ${SUBMIT_DIR}/${StartDir}/${MemberDir}/config.card Restarts OverRule "y"
1293  IGCM_card_WriteOption ${SUBMIT_DIR}/${StartDir}/${MemberDir}/config.card Restarts RestartDate ${HumanRestartDate}
1294  IGCM_card_WriteOption ${SUBMIT_DIR}/${StartDir}/${MemberDir}/config.card Restarts RestartJobName ${initFrom}
1295  IGCM_card_WriteOption ${SUBMIT_DIR}/${StartDir}/${MemberDir}/config.card Restarts RestartPath ${ensemble_Ens_DATE_INITPATH}
1296
1297    # [ATM/OCE/...]
1298  for comp in ${config_ListOfComponents[*]} ; do
1299    IGCM_card_WriteOption ${SUBMIT_DIR}/${StartDir}/${MemberDir}/config.card ${comp} Restart "n"
1300  done
1301  unset initFrom
1302
1303  # ==> Job
1304  sed -e "s/\(#.*\)${config_UserChoices_JobName}\( *#.*\)/\1${MemberDir} \2/"            \
1305      -e "s/\(#.*Script_Output_\)${config_UserChoices_JobName}\(\.*\)/\1${MemberDir}\2/" \
1306      -e "s/^PeriodNb=.*/PeriodNb=60/"                                                   \
1307      ${SUBMIT_DIR}/${StartDir}/${MemberDir}/Job_${MemberDir} > Job_${MemberDir}.tmp
1308  IGCM_sys_Mv Job_${MemberDir}.tmp ${SUBMIT_DIR}/${StartDir}/${MemberDir}/Job_${MemberDir}
1309
1310  IGCM_debug_PopStack "IGCM_ensemble_DateFilesUpdate"
1311}
1312
1313############### Parametric ENSEMBLE #################
Note: See TracBrowser for help on using the repository browser.