source: CMIP6/OMIP/CM61-LR-OMIP1-CO2PIND-R1/DRIVER/opa9.driver @ 4818

Last change on this file since 4818 was 4818, checked in by cetlod, 4 years ago

Save new OMIP experiments
File size: 8.5 KB
Line 
1#!/bin/ksh
2#-----------------------------------------------------------------
3function OCE_Initialize
4{
5    IGCM_debug_PushStack "OCE_Initialize"
6
7    RESOL_OCE=${config_UserChoices_ExpType%%_*}
8    exptype=${config_UserChoices_ExpType##*/}
9
10    IGCM_debug_Print 1 " "
11    IGCM_debug_Print 1 "RESOL_OCE        : ${RESOL_OCE}"
12    IGCM_debug_Print 1 "FORCING          : ${exptype}"
13    IGCM_debug_Print 1 "EXPERIMENT TYPE  : ${opa9_UserChoices_forcing}"
14
15# Local function to find namelists parameters
16supergrep () {
17    grep "^ *$1 *=" $2 | sed -e "s% *\!.*%%"
18}
19
20    OPA_RDT=$( supergrep rn_rdt ${SUBMIT_DIR}/PARAM/NAMELIST/${RESOL_OCE}/namelist_${opa9_UserChoices_forcing}_${exptype}_cfg | sed 's/ *rn_rdt *=//' | sed 's/\. *,//' )
21
22    # Period Length In Days between DateBegin and first day of calendar 0001 01 01
23    # /!\ Needed by OPA namelist to compute file names /!\
24    (( DaysSinceJC = $( IGCM_date_DaysSinceJC ${DateBegin} ) + 1 ))
25
26    (( OPA_NPDT_DAY  =  86400 / OPA_RDT ))
27    (( OPA_NPDT_YEAR  =  OPA_NPDT_DAY * 365     ))
28
29    IGCM_debug_Print 1 " "
30    IGCM_debug_Print 1 "OPA_RDT        ${OPA_RDT}"
31    IGCM_debug_Print 1 "OPA_NPDT_DAY   ${OPA_NPDT_DAY}"
32    IGCM_debug_Print 1 "OPA_NPDT_YEAR  ${OPA_NPDT_YEAR}"
33
34
35    IGCM_debug_PopStack "OCE_Initialize"
36}
37
38#-----------------------------------------------------------------
39function OCE_Update
40{
41    IGCM_debug_PushStack "OCE_Update"
42
43
44    V1D_ENABLE=".FALSE."
45    V5D_ENABLE=".FALSE."
46    V1M_ENABLE=".FALSE."
47    V1Y_ENABLE=".FALSE."
48
49    for frequency in ${config_OCE_WriteFrequency} ; do
50        case ${frequency} in
51            ( 1D|1d ) V1D_ENABLE=".TRUE." ;;
52            ( 5D|5d ) V5D_ENABLE=".TRUE." ;;
53            ( 1M|1m ) V1M_ENABLE=".TRUE." ;;
54            ( *[yY] ) V1Y_ENABLE=".TRUE." ;;
55        esac
56    done
57
58    if [ "${exptype}" = "ia"   ] ; then
59       NbFileInter=${#opa9_UserChoices_ListForcing[*]}
60       (( LengthForcing=${opa9_UserChoices_YearForcingEnd} - ${opa9_UserChoices_YearForcingInit} + 1 ))
61       (( DiffForcing=abs(${opa9_UserChoices_YearForcingEnd} - ${year} )))
62       (( ModuloForcing=${DiffForcing} % ${LengthForcing} ))
63      if [ ${ModuloForcing} -eq 0 ] ; then
64          (( year_forcing=${opa9_UserChoices_YearForcingEnd} ))
65       else
66          if [ ${year} -ge ${opa9_UserChoices_YearForcingInit} ] && [ ${year} -le ${opa9_UserChoices_YearForcingEnd} ] ; then
67             year_forcing=${year}
68          elif [ ${year} -lt ${opa9_UserChoices_YearForcingInit} ] ; then
69             (( year_forcing=${opa9_UserChoices_YearForcingEnd} - ${ModuloForcing} ))
70          else
71             (( year_forcing=${opa9_UserChoices_YearForcingInit} + ${ModuloForcing} - 1 ))
72          fi
73       fi
74       (( year_forcing_m1 = ${year_forcing} - 1))
75       (( year_forcing_p1 = ${year_forcing} + 1))
76       (( i = 0 ))
77       while [ $i -lt ${NbFileInter} ] ; do
78          eval file=${opa9_UserChoices_ListForcing[$i]}
79          IGCM_sys_Get ${opa9_UserChoices_R_FORCING}/${file}_${year_forcing_m1}.nc  ${file}_y${year_m1}.nc
80          IGCM_sys_Get ${opa9_UserChoices_R_FORCING}/${file}_${year_forcing}.nc     ${file}_y${year}.nc
81          IGCM_sys_Get ${opa9_UserChoices_R_FORCING}/${file}_${year_forcing_p1}.nc  ${file}_y${year_p1}.nc
82          (( i = i + 1 ))
83       done
84
85    fi
86
87    # Period Length In Days between DateBegin and first day of calendar 0001 01 01
88    # Needed by OPA namelist to compute correct file names
89    (( PeriodDaysSinceJC =  $( IGCM_date_DaysSinceJC ${PeriodDateBegin} ) + 1  ))
90
91    ## Count of number of time steps from begin of OPA job
92    (( OPA_NIT000 = ( PeriodDaysSinceJC - DaysSinceJC ) * OPA_NPDT_DAY + 1 ))
93    (( OPA_NDT   = PeriodLengthInDays * OPA_NPDT_DAY ))
94    (( OPA_NITEND = OPA_NIT000 + OPA_NDT - 1))
95
96    ## Forcing to write one restart file at the end of period
97    OPA_NSTOCK="${OPA_NITEND}"
98
99    ## Verification of number of time steps per day
100    (( NB_SEC_DAY_MODEL = OPA_NPDT_DAY * OPA_RDT ))
101    (( NB_SEC_DAY = 60 * 60 * 24 ))
102
103    if [ ${NB_SEC_DAY_MODEL} -ne ${NB_SEC_DAY} ]
104        then
105        echo " NB_SEC_DAY_MODEL "${NB_SEC_DAY_MODEL} "NB_SEC_DAY" ${NB_SEC_DAY} "OPA_NPDT_DAY" ${OPA_NPDT_DAY} "OPA_RDT" ${OPA_RDT}
106        echo " VERIFY OPA_NPDT_DAY in opa9.card "
107        exit
108    fi
109
110   ##-- Restart configuration
111    if ( [ "${CumulPeriod}" -eq 1 ] && [ "${config_OCE_Restart}" = "n" ] ) ; then
112      OPA_LRSTAR=.FALSE.
113      OPA_NRSTDT=0
114      #echo "NO OPA RESTART"
115    elif ( [ "${CumulPeriod}" -eq 1 ] && [ "${config_OCE_Restart}" = "y" ] ) ; then
116      OPA_LRSTAR=.TRUE.
117      OPA_NRSTDT=0
118      #echo "OPA RESTART"
119    else
120      OPA_LRSTAR=.TRUE.
121      OPA_NRSTDT=2
122      #echo "OPA RESTART"
123    fi
124
125   
126    OPA_NMSH=0
127    if [ "${opa9_UserChoices_mesh_mask}" = "y" ]; then
128        OPA_NMSH=1
129        IGCM_card_WriteOption ${SUBMIT_DIR}/COMP/opa9.card UserChoices mesh_mask "n"
130    fi
131
132    ## nleapy configuration
133    case ${config_UserChoices_CalendarType} in
134        ( leap|gregorian) OPA_NLEAPY=1  ;;
135        ( noleap        ) OPA_NLEAPY=0  ;;
136        ( 360d          ) OPA_NLEAPY=30 ;;
137        ( *             ) OPA_NLEAPY=30 ;;
138    esac
139
140
141   typeset -r PRECIS=8
142   NEMO_END=$( echo $( awk "BEGIN { printf \"%0${PRECIS}d\",${OPA_NITEND} }" ) ) 
143
144
145   IGCM_comp_modifyNamelist blocker    namelist_cfg cn_exp    ${config_UserChoices_JobName}
146   IGCM_comp_modifyNamelist blocker    namelist_cfg nn_it000  ${OPA_NIT000}
147   IGCM_comp_modifyNamelist blocker    namelist_cfg nn_itend  ${OPA_NITEND}
148   IGCM_comp_modifyNamelist blocker    namelist_cfg ln_rstart ${OPA_LRSTAR}
149   IGCM_comp_modifyNamelist blocker    namelist_cfg nn_stock  ${OPA_NSTOCK}
150   IGCM_comp_modifyNamelist blocker    namelist_cfg nn_rstctl ${OPA_NRSTDT}
151   IGCM_comp_modifyNamelist blocker    namelist_cfg nn_date0  ${PeriodDateBegin}
152   IGCM_comp_modifyNamelist nonblocker namelist_cfg nn_msh    ${OPA_NMSH}
153   IGCM_comp_modifyNamelist blocker    namelist_cfg nn_leapy  ${OPA_NLEAPY} 
154
155   if [ X"$( echo ${config_UserChoices_ExpType} | grep CMIP6 )" != "X" ] ; then
156      OPA_LQSR_RGB=.FALSE.
157      OPA_LQSR_BIO=.TRUE.
158   else
159      OPA_LQSR_RGB=DEFAULT
160      OPA_LQSR_BIO=DEFAULT
161   fi
162   IGCM_comp_modifyNamelist blocker    namelist_cfg ln_qsr_rgb  ${OPA_LQSR_RGB}
163   IGCM_comp_modifyNamelist blocker    namelist_cfg ln_qsr_bio  ${OPA_LQSR_BIO}
164
165   IGCM_debug_Print 1 'Informations into iodef.xml : V1D_ENABLE V5D_ENABLE V1M_ENABLE V1Y_ENABLE '
166   IGCM_debug_Print 1 ${V1D_ENABLE}  ${V5D_ENABLE} ${V1M_ENABLE} ${V1Y_ENABLE}
167
168   IGCM_comp_modifyXmlFile nonblocker file_def_nemo-opa.xml 1d_opa enabled ${V1D_ENABLE}
169   IGCM_comp_modifyXmlFile nonblocker file_def_nemo-opa.xml 5d_opa enabled ${V5D_ENABLE}
170   IGCM_comp_modifyXmlFile nonblocker file_def_nemo-opa.xml 1m_opa enabled ${V1M_ENABLE}
171   IGCM_comp_modifyXmlFile nonblocker file_def_nemo-opa.xml 1y_opa enabled ${V1Y_ENABLE}
172
173   IGCM_comp_modifyXmlFile nonblocker file_def_nemo-opa.xml 1d_opa output_level ${opa9_UserChoices_OutputLevel}
174   IGCM_comp_modifyXmlFile nonblocker file_def_nemo-opa.xml 5d_opa output_level ${opa9_UserChoices_OutputLevel}
175   IGCM_comp_modifyXmlFile nonblocker file_def_nemo-opa.xml 1m_opa output_level ${opa9_UserChoices_OutputLevel}
176   IGCM_comp_modifyXmlFile nonblocker file_def_nemo-opa.xml 1y_opa output_level ${opa9_UserChoices_OutputLevel}
177
178   IGCM_comp_modifyXmlFile force context_nemo.xml ref_year NONE ${InitYear} 
179
180
181   #Long Name as global attribute (if LongName is not empty)
182   if [ ! "X${config_UserChoices_LongName}" = "X" ] ; then
183      listfile=$(ls file_def_nemo*.xml)
184      echo "<variable id=\"LongName\" type=\"string\">${config_UserChoices_LongName}</variable>" > add.tmp
185      for file in ${listfile}
186        do
187            cp ${file} ${file}.tmp
188            sed -e "/<file id/r add.tmp" \
189                ${file}.tmp > ${file}
190            rm ${file}.tmp
191        done
192        rm add.tmp
193   fi
194
195# update domain_def.xml file
196   case ${RESOL_OCE} in
197       ( *ORCA025* ) 
198        sed -e 's/<zoom_domain id="1point" ibegin="139" jbegin="119"/<zoom_domain id="1point" ibegin="600" jbegin="600"/' domain_def_nemo.xml > domain_def.xml.tmp 
199         IGCM_sys_Mv domain_def.xml.tmp domain_def_nemo.xml  ;
200   esac
201
202   IGCM_debug_PopStack "OCE_Update"
203}
204
205#-----------------------------------
206function OCE_Finalize
207{
208    IGCM_debug_PushStack "OCE_Finalize"
209
210    if [ -f date.file ] ; then
211        # Prefix use in opa.card AND in lim2.card :
212        DATE_OPA=$( cat date.file | \
213            sed "s/\ ${config_UserChoices_JobName}_[0-9]*[a-z]_\([0-9]*_[0-9]*\)_\ */\1/g" )
214        ###DATE_OPA=$( $DATE_OPABF | cut -c1-30 )
215        MainPrefix=${config_UserChoices_JobName}_1d_${DATE_OPA}
216        SecondPrefix=${config_UserChoices_JobName}_5d_${DATE_OPA}
217    fi
218
219    echo FINALIZE OCE !!!
220
221    IGCM_debug_PopStack "OCE_Finalize"
222}
223
Note: See TracBrowser for help on using the repository browser.