source: CMIP6/CFMIP/CM61-LR-2xCO2-01/DRIVER/oasis.driver @ 3987

Last change on this file since 3987 was 3987, checked in by lfairhead, 6 years ago

Import of CFMIP 2xCO2 IPSL-CM6A-LR configuration

File size: 4.1 KB
Line 
1#!/bin/ksh
2#-----------------------------------
3CPL_Initialize ()
4{
5    IGCM_debug_PushStack "CPL_Initialize"
6 
7    echo INITIALIZE CPL !!!
8
9# Local function to find parameters in LMDZ .def files
10    lmdzgrep () {
11        grep $1 $2 |grep -v "#" | awk "-F=" '{print $2}' | awk "-F " '{print $1}'
12    }
13
14    #RESOL_CPL=$( echo $RESOL | awk "-F-" '{print $1}' )
15    RESOL_CPL=${RESOL_OCE}x${RESOL_ATM}
16    IGCM_debug_Print 1 "RESOL_CPL      : ${RESOL_CPL}"
17
18    CPL_OkDay="n"
19    CPL_OkMonth="n"
20
21    for CPL_Frequency in ${config_CPL_WriteFrequency} ; do
22        case ${CPL_Frequency} in
23            *D|*d) CPL_OkDay="y"   ;;
24            *M)    CPL_OkMonth="y" ;;
25        esac
26    done
27
28    IGCM_debug_PopStack "CPL_Initialize"
29
30}
31
32#-----------------------------------
33CPL_Update ()
34{
35    IGCM_debug_PushStack "CPL_Update"
36
37    ############ Update Time and Counters ##############
38    ##-- Calcul de la duree du job (en secondes)
39    (( D_JOB_SEC = 3600 * 24 * PeriodLengthInDays ))
40
41    ############ Update Other Parameters  ##############
42    ## Read Output_mode option in oasis.card --> Put "y" to have coupling variables in ouptut files
43    CPL_OutputMode="n"
44    IGCM_card_DefineVariableFromOption ${SUBMIT_DIR}/COMP/oasis.card UserChoices OutputMode
45    eval CPL_OutputMode=\${oasis_UserChoices_OutputMode} > /dev/null 2>&1
46    echo CPL_OutputMode = ${CPL_OutputMode}
47    CPL_OutputVar='EXPORTED'
48    if [ ${CPL_OutputMode} = y ] ; then
49        CPL_OutputVar='EXPOUT'
50    fi
51
52    FreqCoupling=${oasis_UserChoices_FreqCoupling:-86400}
53    FreqCouplingRoffCalv=${oasis_UserChoices_FreqCouplingRoffCalv:-86400}
54
55    ## Calculate lag for transfer of fields from atmosphere -> ocean.
56    ## LagAtm="1day in sec"/day_step * iphysiq/nsplit_phys
57    ## First get day_step, iphyisq and nsplit_phys from LMDZ parameter files
58    day_step=$(    lmdzgrep day_step gcm.def )
59    iphysiq=$(     lmdzgrep iphysiq gcm.def )
60    nsplit_phys=$( lmdzgrep nsplit_phys physiq.def )
61
62    IGCM_debug_Print 3 "Parametre de calcul Lag LMDZ recuperes " ${day_step} ${iphysiq} ${nsplit_phys}
63
64    [[ ${nsplit_phys} -eq 0 || "X${nsplit_phys}" = "X" ]] && nsplit_phys=1
65
66    IGCM_debug_Print 3 "Parametre de calcul Lag LMDZ mis a jour " ${day_step} ${iphysiq} ${nsplit_phys}
67                                                             
68    (( LagAtm = 86400 * iphysiq / day_step / nsplit_phys ))
69    (( LagAtmRoffCalv = $FreqCoupling + $LagAtm ))
70    (( LagOce = 86400 / OPA_NPDT_DAY ))
71
72    IGCM_debug_Print 1 "OPA_NPDT_DAY : ${OPA_NPDT_DAY}"
73    IGCM_debug_Print 1 "LagAtm        : ${LagAtm}"
74    IGCM_debug_Print 1 "LagAtmRoffCalv : ${LagAtmRoffCalv}"
75    IGCM_debug_Print 1 "LagOce        : ${LagOce}"
76
77    ## Use of Lucia ##
78    lucia_ok=""
79    if [ ${oasis_UserChoices_Lucia} = y ] ; then
80        lucia_ok="-1"
81    fi
82   
83    ############ Update Parameter Files   ##############
84    sed -e "s/<total_time>/${D_JOB_SEC}/" \
85        -e "s/<lucia_ok>/${lucia_ok}/" \
86        -e "s/<output_mode>/${CPL_OutputVar}/g" \
87        -e "s/<freq_coupling>/${FreqCoupling}/g" \
88        -e "s/<freq_coupling_roff_calv>/${FreqCouplingRoffCalv}/g" \
89        -e "s/<lag_atm>/${LagAtm}/g" \
90        -e "s/<lag_atm_roff_calv>/${LagAtmRoffCalv}/g" \
91        -e "s/<lag_oce>/${LagOce}/g" \
92        namcouple > namcouple.tmp
93    IGCM_sys_Mv namcouple.tmp namcouple
94    # For IPSLCM6 and use of OASIS-MCT
95    if [ ${CumulPeriod} -eq 1 ]; then
96        ncrename -v OIceFrac,OIceFrc sstoc.nc > /dev/null 2>&1
97    fi
98    # For use of XIOS
99    IGCM_comp_modifyXmlFile force iodef.xml using_oasis NONE true
100
101    IGCM_debug_PopStack "CPL_Update"
102}
103
104#-----------------------------------
105CPL_Finalize ()
106{
107    IGCM_debug_PushStack "CPL_Finalize"
108
109    if [ "${CPL_OutputMode}" = "y" ] ; then
110
111        [[ ${CPL_OkDay}   = "y" ]] && date >TEMPO_1D_cpl_atm.nc
112        [[ ${CPL_OkMonth} = "y" ]] && date >TEMPO_1M_cpl_atm.nc
113        [[ ${CPL_OkDay}   = "y" ]] && date >TEMPO_1D_cpl_oce.nc
114        [[ ${CPL_OkMonth} = "y" ]] && date >TEMPO_1M_cpl_oce.nc
115
116    fi
117
118    ## Use of Lucia ##
119    if [ ${oasis_UserChoices_Lucia} = y ] ; then
120        echo "Use of LUCIA"
121    #   To be changed
122    #   On Curie
123    #   /ccc/scratch/cont003/dsm/p86caub/LUCIA/lucia
124    #   On Ada
125    #   /linkhome/rech/psl/rpsl035/LUCIA/lucia
126    fi
127
128    echo FINALIZE CPL !!!
129
130    IGCM_debug_PopStack "CPL_Finalize"
131}
Note: See TracBrowser for help on using the repository browser.