source: CONFIG/IPSLCM/IPSLCM5/branches/IPSLCM5_v2/EXP00/COMP/lmdz.driver @ 716

Last change on this file since 716 was 697, checked in by mafoipsl, 15 years ago

Parameters included to prepare 365 days calendar simulation.
File size: 6.8 KB
Line 
1#!/bin/ksh
2#-----------------------------------------------------------------
3function ATM_Initialize
4{
5    IGCM_debug_PushStack "ATM_Initialize"
6
7    RESOL_ATM=$( echo $RESOL | awk "-Fx" '{print $2}' | awk "-F-" '{print $1}')
8
9    [ -f ${SUBMIT_DIR}/../.resol ] && eval $(grep RESOL_ATM_3D ${SUBMIT_DIR}/../.resol) || RESOL_ATM_3D=96x95x19
10
11    RESOL_ATM_Z=$( echo ${RESOL_ATM_3D} | awk "-Fx" '{print $3}' )
12
13    ##-- Calendar type for LMDZ
14    case ${config_UserChoices_CalendarType} in
15        leap|gregorian)
16            CalendarTypeForLmdz=earth_366d;;
17        noleap)
18            CalendarTypeForLmdz=earth_365d;;
19        360d)
20            CalendarTypeForLmdz=earth_360d;;
21        *)
22            CalendarTypeForLmdz=earth_360d
23    esac
24
25    ##--Frequency purpose ....
26    ##--  Initialisation  ....
27    OK_instan=n
28    OK_journe=n
29    OK_mensuel=n
30    ok_hf=n
31
32    case ${config_UserChoices_PeriodLength} in
33        1Y|1y|1M|1m) OK_mensuel=y ;;
34        5D|5d|1D|1d) OK_journe=y ;;
35    esac
36
37    for frequency in ${config_ATM_WriteFrequency} ; do
38        case ${frequency} in
39            5D|5d|1D|1d) OK_journe=y ;;
40        esac
41        case ${frequency} in
42            HF|hf) ok_hf=y ;;
43        esac
44    done
45
46    ##--Variables used by LMDZ in physiq.def --
47    PAT_INST=$(        grep -v ^# ${SUBMIT_DIR}/PARAM/physiq.def_L${RESOL_ATM_Z}  | grep 'OK_instan='   )
48    PAT_JOUR=$(        grep -v ^# ${SUBMIT_DIR}/PARAM/physiq.def_L${RESOL_ATM_Z}  | grep 'OK_journe='   )
49    PAT_MOIS=$(        grep -v ^# ${SUBMIT_DIR}/PARAM/physiq.def_L${RESOL_ATM_Z}  | grep 'OK_mensuel='  )
50    PAT_HFRE=$(        grep -v ^# ${SUBMIT_DIR}/PARAM/physiq.def_L${RESOL_ATM_Z}  | grep 'ok_hf='       )
51    PAT_ecrit_mth=$(   grep -v ^# ${SUBMIT_DIR}/PARAM/physiq.def_L${RESOL_ATM_Z}  | grep 'ecrit_mth='   )
52    PAT_ecrit_ISCCP=$( grep -v ^# ${SUBMIT_DIR}/PARAM/physiq.def_L${RESOL_ATM_Z}  | grep 'ecrit_ISCCP=' )
53    ##--Variables used by LMDZ in gcm.def --
54    PAT_iphysiq=$(   grep -v ^# ${SUBMIT_DIR}/PARAM/gcm.def_${RESOL_ATM_3D} | grep iphysiq  )
55    PAT_iperiod=$(   grep -v ^# ${SUBMIT_DIR}/PARAM/gcm.def_${RESOL_ATM_3D} | grep iperiod  )
56    PAT_day_step=$(  grep -v ^# ${SUBMIT_DIR}/PARAM/gcm.def_${RESOL_ATM_3D} | grep day_step )
57
58    ## Read LMDZ_adjust option in ${compname}.card --> put adjust=y and save Bands_*.dat in ${SUBMIT_DIR}/PARAM
59    IGCM_card_DefineVariableFromOption ${SUBMIT_DIR}/COMP/${compname}.card UserChoices LMDZ_adjust
60    eval LMDZ_adjust=\${${compname}_UserChoices_LMDZ_adjust} > /dev/null 2>&1
61
62    ##-- This could be define in lmdz.card, inside section [UserChoices]
63    ##-- Otherwise we get the value in *.def
64    iperiod=$( echo ${PAT_iperiod} | awk -F= '{print $2}' )
65    iphysiq=$( echo ${PAT_iphysiq} | awk -F= '{print $2}' )
66    # day_step : number of steps per day (multiple of iperiod)
67    #(( day_step = 48 * iphysiq ))
68    day_step=$( echo ${PAT_day_step} | awk -F= '{print $2}' )
69
70    #[ -f ${SUBMIT_DIR}/PARAM/Bands_${RESOL_ATM_3D}_${NUM_PROC_ATM}prc.dat ] && IGCM_sys_Cp ${SUBMIT_DIR}/PARAM/Bands_${RESOL_ATM_3D}_${NUM_PROC_ATM}prc.dat .
71
72    IGCM_debug_PopStack "ATM_Initialize"
73}
74
75#-----------------------------------------------------------------
76function ATM_Update
77{
78    IGCM_debug_PushStack "ATM_Update"
79
80    LMDZ_ecrit_mth=${PeriodLengthInDays}.
81    LMDZ_ecrit_ISCCP=${PeriodLengthInDays}.
82    LMDZ_periodav=${PeriodLengthInDays}.
83
84    ##-- Remise ou non a zero de la date initiale de LMDZ pour le fichier run.def
85    if [ ${CumulPeriod} -eq 1 ] ; then
86        RAZ_DATE=1
87    else
88        RAZ_DATE=0
89    fi
90
91    ## Mise en forme du fichier physiq.def
92    sed -e "s/${PAT_INST}/OK_instan=${OK_instan}/"                 \
93        -e "s/${PAT_JOUR}/OK_journe=${OK_journe}/"                 \
94        -e "s/${PAT_MOIS}/OK_mensuel=${OK_mensuel}/"               \
95        -e "s/${PAT_HFRE}/ok_hf=${ok_hf}/"                         \
96        -e "s/${PAT_ecrit_mth}/ecrit_mth=${LMDZ_ecrit_mth}/"       \
97        -e "s/${PAT_ecrit_ISCCP}/ecrit_ISCCP=${LMDZ_ecrit_ISCCP}/" \
98        physiq.def > physiq.def.tmp
99    IGCM_sys_Mv physiq.def.tmp physiq.def
100
101    ## Mise en forme du fichier gcm.def
102    sed -e "s/${PAT_day_step}/day_step=${day_step}/"    \
103        -e "s/${PAT_iperiod}/iperiod=${iperiod}/"       \
104        -e "s/${PAT_iphysiq}/iphysiq=${iphysiq}/"       \
105        gcm.def > gcm.def.tmp
106    IGCM_sys_Mv gcm.def.tmp gcm.def
107
108    ## Mise en forme du fichier run.def
109    sed -e "s/_dayref_/${InitDay}/"                    \
110        -e "s/_anneeref_/${InitYear}/"                 \
111        -e "s/_calend_/${CalendarTypeForLmdz}/"        \
112        -e "s/_nday_/${PeriodLengthInDays}/"           \
113        -e "s/_raz_date_/${RAZ_DATE}/"                 \
114        -e "s/_periodav_/${LMDZ_periodav}/"            \
115        -e "s/_adjust_/${LMDZ_adjust}/"                \
116        run.def > run.def.tmp
117    IGCM_sys_Mv run.def.tmp run.def
118
119    ## Read ByPass_hgardfou_teta option in ${compname}.card --> divide teta* by 2 if [ $ByPass_hgardfou_teta = 1 ]
120    IGCM_card_DefineVariableFromOption ${SUBMIT_DIR}/COMP/${compname}.card UserChoices ByPass_hgardfou_teta
121    eval ByPass_hgardfou_teta=\${${compname}_UserChoices_ByPass_hgardfou_teta} > /dev/null 2>&1
122    if [ ${ByPass_hgardfou_teta} = y ] ; then
123        awk '{ if ($0 ~ /^teta.*=/) {split($0,a,"=") ; print a[1]"="a[2]/2"."} else print $0}' gcm.def > gcm.def.tmp
124        IGCM_sys_Mv gcm.def.tmp gcm.def
125        echo
126        IGCM_debug_Print 1 "ByPass_hgardfou_teta : ^teta*/2 in gcm.def"
127        echo
128        cat gcm.def
129        ByPass_hgardfou_teta=n
130        IGCM_card_WriteOption ${SUBMIT_DIR}/COMP/${compname}.card UserChoices ByPass_hgardfou_teta "${ByPass_hgardfou_teta}"
131    fi
132
133    ## Read ByPass_hgardfou_mats option in ${compname}.card --> purmats=y 2 if [ $ByPass_hgardfou_mats = 1 ]
134    IGCM_card_DefineVariableFromOption ${SUBMIT_DIR}/COMP/${compname}.card UserChoices ByPass_hgardfou_mats
135    eval ByPass_hgardfou_mats=\${${compname}_UserChoices_ByPass_hgardfou_mats} > /dev/null 2>&1
136    if [ ${ByPass_hgardfou_mats} = y ] ; then
137        sed -e "s/^purmats=.*/purmats=y/" gcm.def > gcm.def.tmp
138        IGCM_sys_Mv gcm.def.tmp gcm.def
139        echo
140        IGCM_debug_Print 1 "ByPass_hgardfou_mats : purmats=y in gcm.def"
141        echo
142        cat gcm.def
143        ByPass_hgardfou_mats=n
144        IGCM_card_WriteOption ${SUBMIT_DIR}/COMP/${compname}.card UserChoices ByPass_hgardfou_mats "${ByPass_hgardfou_mats}"
145    fi
146
147    if ( ${FirstInitialize} ) ; then
148
149        if ( [ "${config_Restarts_OverRule}" = "n" ] && [ "${config_ATM_Restart}" = "n" ] ); then
150            if [ ${DRYRUN} -le 2 ] ; then
151                IGCM_sys_Cp ${R_EXE}/create_etat0_limit.e ${RUN_DIR}/.
152                if [ ${DRYRUN} -le 1 ] ; then
153                    ./create_etat0_limit.e
154                fi
155            else
156                echo "EXECUTION of ./create_etat0_limit.e simulated"
157                echo "EXECUTION of ./create_etat0_limit.e simulated for DRYRUN = " $DRYRUN >> stack
158            fi
159        fi
160
161    fi
162
163    IGCM_debug_PopStack "ATM_Update"
164}
165
166#-----------------------------------
167function ATM_Finalize
168{
169    IGCM_debug_PushStack "ATM_Finalize"
170
171    [ ${LMDZ_adjust} = y ] && IGCM_sys_Cp Bands_${RESOL_ATM_3D}_${NUM_PROC_ATM}prc.dat ${SUBMIT_DIR}/PARAM
172
173    echo FINALIZE ATM !
174
175    IGCM_debug_PopStack "ATM_Finalize"
176}
Note: See TracBrowser for help on using the repository browser.