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

Last change on this file since 711 was 711, checked in by sdipsl, 16 years ago
  • RESOL_ATM must not contain any z level information (as long as atlas repository has not been adapted to 39 levels)
File size: 4.5 KB
Line 
1# $Id$
2#
3#!/bin/ksh
4#-----------------------------------------------------------------
5function ATM_Initialize
6{
7    IGCM_debug_PushStack "ATM_Initialize"
8
9    RESOL_ATM=$( echo $RESOL | awk "-F-" '{print $1}')
10
11    [ -f ${SUBMIT_DIR}/../.resol ] && eval $(grep RESOL_ATM_3D ${SUBMIT_DIR}/../.resol) || RESOL_ATM_3D=96x95x19
12
13    RESOL_ATM_Z=$( echo $RESOL_ATM_3D | awk "-Fx" '{print $3}' )
14
15    ##- Create_etat0_limit version
16    IGCM_card_DefineVariableFromOption ${SUBMIT_DIR}/COMP/lmdz.card UserChoices CREATE
17    CREATE=${lmdz_UserChoices_CREATE}
18    echo create_etat0_limit version : ${CREATE}
19
20    ##- Default number of processor for lmdz
21    NUM_PROC_ATM=1
22    if [ X"${BATCH_NUM_PROC_TOT}" != X ] ; then
23        NUM_PROC_ATM=${BATCH_NUM_PROC_TOT}
24    fi
25 
26    ##--Frequency purpose ....
27    ##--  Initialisation  ....
28    OK_instan=n
29    OK_journe=n
30    OK_mensuel=n
31    ok_hf=n
32
33    case ${config_UserChoices_PeriodLength} in
34        1Y|1y|1M|1m) OK_mensuel=y ;;
35        5D|5d|1D|1d) OK_journe=y ;;
36    esac
37
38    for frequency in ${config_ATM_WriteFrequency} ; do
39        case ${frequency} in
40            5D|5d|1D|1d) OK_journe=y ;;
41        esac
42        case ${frequency} in
43            HF|hf) ok_hf=y ;;
44        esac
45    done
46
47    ##--Variables used by LMDZ in physiq.def --
48    PAT_INST=$(      grep -v ^# ${SUBMIT_DIR}/PARAM/physiq.def_L${RESOL_ATM_Z}  | grep 'OK_instan='  )
49    PAT_JOUR=$(      grep -v ^# ${SUBMIT_DIR}/PARAM/physiq.def_L${RESOL_ATM_Z}  | grep 'OK_journe='  )
50    PAT_MOIS=$(      grep -v ^# ${SUBMIT_DIR}/PARAM/physiq.def_L${RESOL_ATM_Z}  | grep 'OK_mensuel=' )
51    PAT_HFRE=$(      grep -v ^# ${SUBMIT_DIR}/PARAM/physiq.def_L${RESOL_ATM_Z}  | grep 'ok_hf='      )
52    PAT_ecrit_mth=$( grep -v ^# ${SUBMIT_DIR}/PARAM/physiq.def_L${RESOL_ATM_Z}  | grep 'ecrit_mth='  )
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    ecrit_mth=${PeriodLengthInDays}.
81
82    ##-- Remise ou non a zero de la date initiale de LMDZ pour le fichier run.def
83    if [ ${CumulPeriod} -eq 1 ] ; then
84        RAZ_DATE=1
85    else
86        RAZ_DATE=0
87    fi
88
89    ## Mise en forme du fichier physiq.def
90    sed -e "s/${PAT_INST}/OK_instan=${OK_instan}/"      \
91        -e "s/${PAT_JOUR}/OK_journe=${OK_journe}/"      \
92        -e "s/${PAT_MOIS}/OK_mensuel=${OK_mensuel}/"    \
93        -e "s/${PAT_HFRE}/ok_hf=${ok_hf}/"              \
94        -e "s/${PAT_ecrit_mth}/ecrit_mth=${ecrit_mth}/" \
95        physiq.def > physiq.def.tmp
96    IGCM_sys_Mv physiq.def.tmp physiq.def
97
98    ## Mise en forme du fichier gcm.def
99    sed -e "s/${PAT_day_step}/day_step=${day_step}/"    \
100        -e "s/${PAT_iperiod}/iperiod=${iperiod}/"       \
101        -e "s/${PAT_iphysiq}/iphysiq=${iphysiq}/"       \
102        gcm.def > gcm.def.tmp
103    IGCM_sys_Mv gcm.def.tmp gcm.def
104
105    ## Mise en forme du fichier run.def
106    sed -e "s/_dayref_/${InitDay}/"                    \
107        -e "s/_anneeref_/${InitYear}/"                 \
108        -e "s/_nday_/${PeriodLengthInDays}/"           \
109        -e "s/_raz_date_/${RAZ_DATE}/"                 \
110        -e "s/_adjust_/${LMDZ_adjust}/"                \
111        run.def > run.def.tmp
112    IGCM_sys_Mv run.def.tmp run.def
113
114    IGCM_debug_PopStack "ATM_Update"
115}
116
117#-----------------------------------
118function ATM_Finalize
119{
120    IGCM_debug_PushStack "ATM_Finalize"
121
122    [ ${LMDZ_adjust} = y ] && IGCM_sys_Cp Bands_${RESOL_ATM_3D}_${NUM_PROC_ATM}prc.dat ${SUBMIT_DIR}/PARAM
123
124    echo FINALIZE ATM !
125
126    IGCM_debug_PopStack "ATM_Finalize"
127}
Note: See TracBrowser for help on using the repository browser.