source: CONFIG/LMDZOR/branches/LMDZ4OR_v3/CREATE/COMP/lmdz.driver @ 1066

Last change on this file since 1066 was 1066, checked in by lfairhead, 14 years ago

Variable LMDZ_ajust is not set for CREATE run and cannot be used

File size: 4.9 KB
Line 
1#!/bin/ksh
2#-----------------------------------------------------------------
3function ATM_Initialize
4{
5    IGCM_debug_PushStack "ATM_Initialize"
6
7    RESOL_ATM=${RESOL}
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    ##-- This could be define in lmdz.card, inside section [UserChoices]
59    ##-- Otherwise we get the value in *.def
60    iperiod=$( echo ${PAT_iperiod} | awk -F= '{print $2}' )
61    iphysiq=$( echo ${PAT_iphysiq} | awk -F= '{print $2}' )
62    # day_step : number of steps per day (multiple of iperiod)
63    #(( day_step = 48 * iphysiq ))
64    day_step=$( echo ${PAT_day_step} | awk -F= '{print $2}' )
65
66    IGCM_debug_PopStack "ATM_Initialize"
67}
68
69#-----------------------------------------------------------------
70function ATM_Update
71{
72    IGCM_debug_PushStack "ATM_Update"
73
74    case ${config_UserChoices_PeriodLength} in
75        *Y|*y) 
76               LMDZ_ecrit_mth=30.
77               LMDZ_ecrit_ISCCP=30.
78               LMDZ_periodav=30.
79               if [ "${config_UserChoices_CalendarType}" != "360d" ] ; then
80                  echo Change calendartype in config.card. PeriodLength=1Y allowed only for CalendarType=360d
81                  exit
82               fi
83               ;;
84        *)
85               LMDZ_ecrit_mth=${PeriodLengthInDays}.
86               LMDZ_ecrit_ISCCP=${PeriodLengthInDays}.
87               LMDZ_periodav=${PeriodLengthInDays}.
88               ;;
89    esac
90
91    ##-- Remise ou non a zero de la date initiale de LMDZ pour le fichier run.def
92    if [ ${CumulPeriod} -eq 1 ] ; then
93        RAZ_DATE=1
94    else
95        RAZ_DATE=0
96    fi
97
98    ## Mise en forme du fichier physiq.def
99    sed -e "s/${PAT_INST}/OK_instan=${OK_instan}/"                 \
100        -e "s/${PAT_JOUR}/OK_journe=${OK_journe}/"                 \
101        -e "s/${PAT_MOIS}/OK_mensuel=${OK_mensuel}/"               \
102        -e "s/${PAT_HFRE}/ok_hf=${ok_hf}/"                         \
103        -e "s/${PAT_ecrit_mth}/ecrit_mth=${LMDZ_ecrit_mth}/"       \
104        -e "s/${PAT_ecrit_ISCCP}/ecrit_ISCCP=${LMDZ_ecrit_ISCCP}/" \
105        physiq.def > physiq.def.tmp
106    IGCM_sys_Mv physiq.def.tmp physiq.def
107
108    echo cat
109    cat physiq.def
110
111    ## Mise en forme du fichier gcm.def
112    sed -e "s/${PAT_day_step}/day_step=${day_step}/"    \
113        -e "s/${PAT_iperiod}/iperiod=${iperiod}/"       \
114        -e "s/${PAT_iphysiq}/iphysiq=${iphysiq}/"       \
115        gcm.def > gcm.def.tmp
116    IGCM_sys_Mv gcm.def.tmp gcm.def
117
118    echo cat
119    cat gcm.def
120
121    ## Mise en forme du fichier run.def
122    sed -e "s/_dayref_/${InitDay}/"                    \
123        -e "s/_anneeref_/${InitYear}/"                 \
124        -e "s/_calend_/${CalendarTypeForLmdz}/"        \
125        -e "s/_nday_/${PeriodLengthInDays}/"           \
126        -e "s/_raz_date_/${RAZ_DATE}/"                 \
127        -e "s/_periodav_/${LMDZ_periodav}/"            \
128        -e "s/_adjust_/n/"                \
129        run.def > run.def.tmp
130    IGCM_sys_Mv run.def.tmp run.def
131
132    echo cat run
133    cat run.def
134
135    IGCM_debug_PopStack "ATM_Update"
136}
137
138#-----------------------------------
139function ATM_Finalize
140{
141    IGCM_debug_PushStack "ATM_Finalize"
142
143    echo FINALIZE ATM !
144
145    IGCM_debug_PopStack "ATM_Finalize"
146}
Note: See TracBrowser for help on using the repository browser.