Changeset 187


Ignore:
Timestamp:
10/22/07 16:17:50 (17 years ago)
Author:
acosce
Message:

ACo : adapt IPSL_ESM_V1 to changes in libIGCM (platine, compatibility)

adapt IPSL_ESM_V1 to run PDVTLV2

Location:
CONFIG/trunk/IPSL_ESM_V1/EXP00
Files:
1 added
17 edited

Legend:

Unmodified
Added
Removed
  • CONFIG/trunk/IPSL_ESM_V1/EXP00/COMP/inca.card

    r103 r187  
     1[Compatibility] 
     2libIGCM=1.0 
     3 
    14[UserChoices] 
    25 
    36[InitialStateFiles] 
    4 List=   (${R_INIT}/CHM/start_AER_9671.nc, start_trac.nc) 
     7List=   (${R_INIT}/CHM/${config_UserChoices_TagName}/AER/INCA${RESOL_CHM}/etat0_aerosol.nc, start_trac.nc) 
    58         
    69 
     
    1215 
    1316ListNonDel= (${MODIPSL}/modeles/INCA3/SIMULATIONS/AER/inca.dat , .       )\ 
    14         (${R_INIT}/CHM/AER/rhvlmdz_96x71.txt             ,   rhv.dat      ) \ 
    15         (${R_INIT}/CHM/AER/wthlmdz_96x71.txt             ,   wth.dat      ) \ 
    16         (${R_INIT}/CHM/AER/clylmdz_96x71.txt             ,   cly.dat      ) \ 
    17         (${R_INIT}/CHM/AER/rhvEC11.09092002.txt          ,   rhvEC.txt    ) \ 
    18         (${R_INIT}/CHM/AER/wthEC11.09092002.txt          ,   wthEC.txt    ) \ 
    19         (${R_INIT}/CHM/AER/clyEC11.09092002.txt          ,   clyEC.txt    ) \ 
    20         (${R_INIT}/CHM/AER/sflx_lmdz_phy.couple_july07.nc   ,   sflx.nc      ) \ 
    21         (${R_INIT}/CHM/CH4/landuse_96x71.nc              ,   landuse.nc   ) \ 
    22         (${R_INIT}/CHM/CH4/npp_96x71.nc                  ,   npp.nc           ) \ 
    23         (${R_INIT}/CHM/CH4/oxydants_96x71.nc             ,   oxydants.nc  ) \ 
    24         (${R_INIT}/CHM/AER/obs_aerinput.prn              ,   .   ) 
     17        (${R_INIT}/CHM/${config_UserChoices_TagName}/AER/INCA${RESOL_CHM}/rhvlmdz.txt             ,   rhv.dat      ) \ 
     18        (${R_INIT}/CHM/${config_UserChoices_TagName}/AER/INCA${RESOL_CHM}/wthlmdz.txt             ,   wth.dat      ) \ 
     19        (${R_INIT}/CHM/${config_UserChoices_TagName}/AER/INCA${RESOL_CHM}/clylmdz.txt             ,   cly.dat      ) \ 
     20        (${R_INIT}/CHM/${config_UserChoices_TagName}/AER/rhvEC11.09092002.txt          ,   rhvEC.txt    ) \ 
     21        (${R_INIT}/CHM/${config_UserChoices_TagName}/AER/wthEC11.09092002.txt          ,   wthEC.txt    ) \ 
     22        (${R_INIT}/CHM/${config_UserChoices_TagName}/AER/clyEC11.09092002.txt          ,   clyEC.txt    ) \ 
     23        (${R_INIT}/CHM/${config_UserChoices_TagName}/AER/sflx_lmdz_phy.couple_july07.nc   ,   sflx.nc      ) \ 
     24        (${R_INIT}/CHM/${config_UserChoices_TagName}/CH4/INCA${RESOL_CHM}/landuse.nc              ,   landuse.nc   ) \ 
     25        (${R_INIT}/CHM/${config_UserChoices_TagName}/CH4/INCA${RESOL_CHM}/npp.nc                  ,   npp.nc          ) \ 
     26        (${R_INIT}/CHM/${config_UserChoices_TagName}/CH4/INCA${RESOL_CHM}/oxydants.nc             ,   oxydants.nc  ) \ 
     27        (${R_INIT}/CHM/${config_UserChoices_TagName}/AER/obs_aerinput.prn              ,   .   ) 
    2528 
    2629[ParametersFiles] 
  • CONFIG/trunk/IPSL_ESM_V1/EXP00/COMP/inca.driver

    r103 r187  
    88    IGCM_debug_PushStack "CHM_Initialize" 
    99 
     10    RESOL_CHM=$( echo $RESOL_ATM | awk "-FD" '{print $2}' ) 
     11    echo "resolution : " $RESOL_CHM 
    1012    echo "INCA Tag : LMDZ-INCA" ${config_UserChoices_TagName} 
    1113 
  • CONFIG/trunk/IPSL_ESM_V1/EXP00/COMP/lim.card

    r103 r187  
     1[Compatibility] 
     2libIGCM=1.0 
     3 
    14[UserChoices] 
    25 
     
    2326 
    2427[OutputFiles] 
    25 List=   (${PREFIX_MO}_icemod.nc, ${R_OUT_ICE_O_M}/${PREFIX}_1M_icemod.nc, Post_1M_icemod), \ 
    26         (${PREFIX_5D}_icemod.nc, ${R_OUT_ICE_O_D}/${PREFIX}_5D_icemod.nc, NONE), \ 
    27         (${PREFIX_DA}_icemod.nc, ${R_OUT_ICE_O_D}/${PREFIX}_1D_icemod.nc, NONE) 
     28List=   (${PREFIX_NWRITE}_${DATE_OPA}_icemod.nc, ${R_OUT_ICE_NWRITE}/${PREFIX}_${WF1}_icemod.nc, Post_1M_icemod) 
    2829 
    2930[Post_1M_icemod] 
  • CONFIG/trunk/IPSL_ESM_V1/EXP00/COMP/lim.driver

    r103 r187  
    2121    ##-- MUST BE DONE AFTER OCE_UPDATE 
    2222 
    23     (( ICE_NSTOCK = ORCA_NSTOCK )) 
    24     ## I PREFER (( ICE_NSTOCK = ORCA_NWRITE )) 
    25     ## BUT      (( ICE_NSTOCK = ORCA_NSTOCK )) is in AA_job 
     23    (( ICE_NSTOCK = ORCA_NWRITE )) 
    2624 
    2725    ##-- Restart configuration 
  • CONFIG/trunk/IPSL_ESM_V1/EXP00/COMP/lmdz.card

    r103 r187  
     1[Compatibility] 
     2libIGCM=1.0 
     3 
    14[UserChoices] 
    25 
     
    2730            (${R_BC}/ATM/${config_UserChoices_TagName}/${RESOL_ATM}/HISTORIQUE/so4.run2000.cdf, .),\ 
    2831            (${R_BC}/ATM/${config_UserChoices_TagName}/${RESOL_ATM}/HISTORIQUE/so4.run.nat.cdf, .),\ 
    29             (${R_BC}/ATM/${config_UserChoices_TagName}/${RESOL_ATM}/flux_iceberg, .),\ 
    3032            (${R_BC}/ATM/${config_UserChoices_TagName}/ISCCP/tautab.formatted, .),\ 
    3133            (${R_BC}/ATM/${config_UserChoices_TagName}/ISCCP/invtau.formatted, .) 
     
    3537        (${SUBMIT_DIR}/PARAM/gcm.def_${RESOL_ATM}, gcm.def), \ 
    3638        (${SUBMIT_DIR}/PARAM/run.def, .), \ 
    37         (${SUBMIT_DIR}/PARAM/traceur.def, .)  
     39        (${SUBMIT_DIR}/PARAM/traceur.def, .), \ 
     40        (${SUBMIT_DIR}/PARAM/Bands_${BandsResol}_${NUM_PROC_ATM}prc.dat, .) 
    3841 
    3942[RestartFiles] 
     
    4649[OutputFiles] 
    4750List=   (histmth.nc,      ${R_OUT_ATM_O_M}/${PREFIX}_1M_histmth.nc,      Post_1M_histmth), \ 
    48         (histISCCP.nc,    ${R_OUT_ATM_O_M}/${PREFIX}_1M_histISCCP.nc,    NONE),\ 
     51        (histISCCP.nc,    ${R_OUT_ATM_O_M}/${PREFIX}_1M_histISCCP.nc,    NONE), \ 
    4952        (histNMC.nc,      ${R_OUT_ATM_O_M}/${PREFIX}_1M_histNMC.nc,      NONE), \ 
    50         (dynzon.nc,       ${R_OUT_ATM_O_M}/${PREFIX}_1M_dynzon.nc,       Post_1M_dynzon), \ 
    51         (histday.nc,      ${R_OUT_ATM_O_D}/${PREFIX}_1D_histday.nc,      Post_1D_histday), \ 
     53        (dynzon.nc,       ${R_OUT_ATM_O_M}/${PREFIX}_1M_dynzon.nc,       NONE), \ 
     54        (histday.nc,      ${R_OUT_ATM_O_D}/${PREFIX}_1D_histday.nc,      NONE), \ 
    5255        (histhf.nc,       ${R_OUT_ATM_O_H}/${PREFIX}_HF_histhf.nc,       NONE), \ 
    5356        (histhf3d.nc,     ${R_OUT_ATM_O_H}/${PREFIX}_HF_histhf3d.nc,     NONE), \ 
     
    5861GatherWithInternal = (lon, lat, presnivs, time_counter) 
    5962MonitoringVars= (bils, nettop, precip, tsol) 
    60 TimeSeriesVars= (bils, cldh, cldl, cldm, cldq, cldt, evap, iwp, lwp, LWdnSFC, LWdnSFCclr, LWupSFCclr, LWdn200, LWdn200clr, LWup200, LWup200clr, precip, prw, pluc, psol, q2m, sens, slp, snow, snow_mass, soll0, soll, sols0, sols, SWdnSFC, SWupSFC, SWup200, SWup200clr, SWdnTOA, SWupTOA, SWdn200, SWdn200clr, SWdnSFCclr, SWupSFCclr, SWupTOAclr, t2m, taux_ter, taux_sic, taux_lic, taux_oce, tauy_ter, tauy_sic, tauy_lic, tauy_oce, lat_ter, lat_sic, lat_lic, lat_oce, topl0, topl, tops0, tops, topsad, topsai, tsol, tsol_ter, tsol_lic, tsol_oce, tsol_sic , u10m, u500, v10m, v500, w500, pourc_ter, pourc_oce, pourc_sic, pourc_lic, geop, ovap, rhum, rneb, pres, temp, vitu, vitv, vitw) 
     63TimeSeriesVars= (bils, cldh, cldl, cldm, cldq, cldt, evap, iwp, lwp, LWdnSFC, LWdnSFCclr, LWupSFCclr, LWdn200, LWdn200clr, LWup200, LWup200clr, precip, prw, pluc, psol, q2m, sens, slp, snowf, soll0, soll, sols0, sols, SWdnSFC, SWupSFC, SWup200, SWup200clr, SWdnTOA, SWupTOA, SWdn200, SWdn200clr, SWdnSFCclr, SWupSFCclr, SWupTOAclr, t2m, taux_ter, taux_sic, taux_lic, taux_oce, tauy_ter, tauy_sic, tauy_lic, tauy_oce, lat_ter, lat_sic, lat_lic, lat_oce, topl0, topl, tops0, tops, tsol, tsol_ter, tsol_lic, tsol_oce, tsol_sic , u10m, u500, v10m, v500, w500, pourc_ter, pourc_oce, pourc_sic, pourc_lic) 
    6164 
    6265[Post_1M_histISCCP] 
  • CONFIG/trunk/IPSL_ESM_V1/EXP00/COMP/lmdz.driver

    r103 r187  
    66 
    77    RESOL_ATM=$( echo $RESOL | awk "-Fx" '{print $2}' ) 
     8 
     9    case ${RESOL_ATM} in 
     10        LMD7245)   BandsResol=72x45x19   ;; 
     11        LMD9671)   BandsResol=96x71x19   ;; 
     12        LMD14496)  BandsResol=144x96x19  ;; 
     13        LMD144142) BandsResol=144x142x19 ;; 
     14    esac 
    815 
    916    ##--Frequency purpose .... 
     
    5057    IGCM_debug_PushStack "ATM_Update" 
    5158 
    52 ecritphy=${PeriodLengthInDays} 
     59    ecritphy=${PeriodLengthInDays} 
    5360 
    54 ##-- Remise ou non a zero de la date initiale de LMDZ pour le fichier run.def 
    55 if [ ${CumulPeriod} -eq 1 ] ; then 
    56      RAZ_DATE=1 
    57 else 
    58      RAZ_DATE=0 
    59 fi 
     61    ##-- Remise ou non a zero de la date initiale de LMDZ pour le fichier run.def 
     62    if [ ${CumulPeriod} -eq 1 ] ; then 
     63        RAZ_DATE=1 
     64    else 
     65        RAZ_DATE=0 
     66    fi 
    6067 
    61 ## Mise en forme du fichier physiq.def 
    62 sed -e "s/OK_instan=.*/OK_instan=${OK_instan}/"   \ 
    63     -e "s/OK_journe=.*/OK_journe=${OK_journe}/"   \ 
    64     -e "s/OK_mensuel=.*/OK_mensuel=${OK_mensuel}/" \ 
     68    ## Mise en forme du fichier physiq.def 
     69    sed -e "s/OK_instan=.*/OK_instan=${OK_instan}/"   \ 
     70        -e "s/OK_journe=.*/OK_journe=${OK_journe}/"   \ 
     71        -e "s/OK_mensuel=.*/OK_mensuel=${OK_mensuel}/" \ 
    6572        physiq.def > physiq.def.tmp 
    6673    IGCM_sys_Mv physiq.def.tmp physiq.def 
    6774 
    68 ## Mise en forme du fichier gcm.def 
    69 sed  -e "s/${PAT_ecritphy}/ecritphy=${ecritphy}/" \ 
    70      -e "s/${PAT_day_step}/day_step=${day_step}/" \ 
    71      -e "s/${PAT_iperiod}/iperiod=${iperiod}/" \ 
    72      -e "s/${PAT_iphysiq}/iphysiq=${iphysiq}/" \ 
     75    ## Mise en forme du fichier gcm.def 
     76    sed  -e "s/${PAT_ecritphy}/ecritphy=${ecritphy}/" \ 
     77        -e "s/${PAT_day_step}/day_step=${day_step}/" \ 
     78        -e "s/${PAT_iperiod}/iperiod=${iperiod}/" \ 
     79        -e "s/${PAT_iphysiq}/iphysiq=${iphysiq}/" \ 
    7380        gcm.def > gcm.def.tmp 
    7481    IGCM_sys_Mv gcm.def.tmp gcm.def 
    7582 
    76 ## Mise en forme du fichier run.def 
    77 sed -e "s/_dayref_/${InitDay}/" \ 
    78     -e "s/_anneeref_/${InitYear}/" \ 
    79     -e "s/_nday_/${PeriodLengthInDays}/" \ 
    80     -e "s/_raz_date_/${RAZ_DATE}/" \ 
     83    ## Mise en forme du fichier run.def 
     84    sed -e "s/_dayref_/${InitDay}/" \ 
     85        -e "s/_anneeref_/${InitYear}/" \ 
     86        -e "s/_nday_/${PeriodLengthInDays}/" \ 
     87        -e "s/_raz_date_/${RAZ_DATE}/" \ 
    8188        run.def > run.def.tmp 
    8289    IGCM_sys_Mv run.def.tmp run.def 
    8390 
    84 if ( ${FirstInitialize} ) ; then 
     91    if ( ${FirstInitialize} ) ; then 
    8592 
    86    if ( [ "${config_Restarts_OverRule}" = "n" ] && [ "${config_ATM_Restart}" = "n" ] ); then 
    87            IGCM_sys_Cp ${R_EXE}/create_etat0_limit.e ${RUN_DIR}/. 
    88            ./create_etat0_limit.e 
    89    fi 
     93        if ( [ "${config_Restarts_OverRule}" = "n" ] && [ "${config_ATM_Restart}" = "n" ] ); then 
     94            if [ ${DRYRUN} -le 2 ] ; then 
     95                IGCM_sys_Cp ${R_EXE}/create_etat0_limit.e ${RUN_DIR}/. 
     96                if [ ${DRYRUN} -le 1 ] ; then 
     97                    ./create_etat0_limit.e 
     98                fi 
     99            else 
     100                echo "EXECUTION of ./create_etat0_limit.e simulated" 
     101                echo "EXECUTION of ./create_etat0_limit.e simulated for DRYRUN = " $DRYRUN >> stack 
     102            fi 
     103        fi 
    90104 
    91 fi  
     105    fi  
    92106 
    93107    IGCM_debug_PopStack "ATM_Update" 
  • CONFIG/trunk/IPSL_ESM_V1/EXP00/COMP/oasis.card

    r103 r187  
     1[Compatibility]  
     2libIGCM=1.0  
     3          
    14[UserChoices] 
    25 
     
    1013            (${R_BC}/CPL/${config_UserChoices_TagName}/${RESOL}/masks.nc,                 masks.nc),\ 
    1114            (${R_BC}/CPL/${config_UserChoices_TagName}/${RESOL}/areas.nc,                 areas.nc),\ 
    12             (${R_BC}/CPL/${config_UserChoices_TagName}/${RESOL}/mozaic.wa2o.i4.r8,        wa2o.flx),\ 
    13             (${R_BC}/CPL/${config_UserChoices_TagName}/${RESOL}/mozaic.wo2a.i4.r8,        wo2a.tsg),\ 
    14             (${R_BC}/CPL/${config_UserChoices_TagName}/${RESOL}/mozaic.wa2o.runoff.i4.r8, wa2o.run),\ 
    15             (${R_BC}/CPL/${config_UserChoices_TagName}/${RESOL}/mozaic.wa2o.calvin.i4.r8, wa2o.cal) 
     15            (${R_BC}/CPL/${config_UserChoices_TagName}/${RESOL}/mozaic.wa2o,        wa2o.flx),\ 
     16            (${R_BC}/CPL/${config_UserChoices_TagName}/${RESOL}/mozaic.wo2a,        wo2a.tsg),\ 
     17            (${R_BC}/CPL/${config_UserChoices_TagName}/${RESOL}/mozaic.wa2o.runoff, wa2o.run),\ 
     18            (${R_BC}/CPL/${config_UserChoices_TagName}/${RESOL}/mozaic.wa2o.calvin, wa2o.cal) 
    1619 
    1720[ParametersFiles] 
     
    2730 
    2831[OutputFiles] 
    29 List=   (cpl_oce_tau.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_cpl_oce_tau.nc,    Post_cpl_oce_tau),\ 
    30         (cpl_oce_flx.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_cpl_oce_flx.nc,    Post_cpl_oce_flx),\ 
    31         (cpl_oce_sst.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_cpl_oce_sst.nc,    Post_cpl_oce_sst) 
     32List=   (cpl_oce_tau.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_cpl_oce_tau.nc,    NONE),\ 
     33        (cpl_oce_flx.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_cpl_oce_flx.nc,    NONE),\ 
     34        (cpl_oce_sst.nc,    ${R_OUT_CPL_O_M}/${PREFIX}_cpl_oce_sst.nc,    NONE) 
    3235 
    3336[Post_cpl_oce_tau] 
     
    4548MonitoringVars = 
    4649TimeSeriesVars = 
    47  
    48 [Post_cpl_atm_sst] 
    49 Patches = () 
    50 MonitoringVars = 
    51 TimeSeriesVars = 
    52  
    53 [Post_cpl_atm_tauflx] 
    54 Patches = () 
    55 MonitoringVars = 
    56 TimeSeriesVars = 
  • CONFIG/trunk/IPSL_ESM_V1/EXP00/COMP/oasis.driver

    r136 r187  
    11#!/bin/ksh 
    2  
    32#----------------------------------- 
    43CPL_Initialize () 
     
    98    echo INITIALIZE CPL !!! 
    109 
    11 #D- Default number of processor for lmdz and oasis 
    12 NUM_PROC_ATM=1 ; (( NUM_PROC_OASIS = NUM_PROC_ATM + 2 )) 
    13 #D- Number of processors used for lmdz and oasis coupler computed from PBS variable 
    14 [ -z "${PBS_NUM_PROC_TOT}" ] || { (( NUM_PROC_ATM = PBS_NUM_PROC_TOT - 1 )) ; (( NUM_PROC_OASIS = PBS_NUM_PROC_TOT + 1 )) ; } 
     10#D- MPI context used to launch the coupled model 
     11LIBMPI=$(cat ${SUBMIT_DIR}/../.libmpi) 
     12echo LIBMPI=${LIBMPI} 
    1513 
    16 echo PBS_NUM_PROC_TOT=${PBS_NUM_PROC_TOT} 
    17 echo NUM_PROC_ATM=${NUM_PROC_ATM} 
    18 echo NUM_PROC_OASIS=${NUM_PROC_OASIS} 
     14    IGCM_sys_build_run_file ${LIBMPI} 
    1915 
    2016    IGCM_debug_PopStack "CPL_Initialize" 
    21  
    2217 
    2318} 
     
    2621CPL_Update () 
    2722{ 
     23 
     24    IGCM_debug_PushStack "CPL_Update" 
    2825 
    2926############ Update Time and Counters ############## 
     
    3633sed -e "s/<initial_date>/${PeriodDateBegin}/" \ 
    3734    -e "s/<total_time>/${D_JOB_SEC}/" \ 
     35    -e "s/<lib_mpi>/${LIBMPI}/" \ 
    3836    -e "s/<nproc_atm>/${NUM_PROC_ATM}/g" \ 
    3937        namcouple > namcouple.tmp 
     
    4644    eval FileToBeDeleted[${#FileToBeDeleted[@]}]=anaisout 
    4745 
     46    IGCM_debug_PopStack "CPL_Update" 
    4847} 
    4948 
     
    5251{ 
    5352 
    54     IGCM_debug_PushStack "CPL_Update" 
     53    IGCM_debug_PushStack "CPL_Finalize" 
    5554 
    5655    echo FINALIZE CPL !!! 
    5756 
    58     IGCM_debug_PopStack "CPL_Update" 
     57    IGCM_debug_PopStack "CPL_Finalize" 
    5958 
    6059} 
  • CONFIG/trunk/IPSL_ESM_V1/EXP00/COMP/opa.card

    r103 r187  
     1[Compatibility] 
     2libIGCM=1.0 
     3 
    14[UserChoices] 
    25OPA_NPDT_JOURS=15 
     
    2629 
    2730[OutputFiles] 
    28 List=   (${PREFIX_MO}_grid_T.nc, ${R_OUT_OCE_O_M}/${PREFIX}_1M_grid_T.nc, Post_1M_grid_T),\ 
    29         (${PREFIX_MO}_grid_U.nc, ${R_OUT_OCE_O_M}/${PREFIX}_1M_grid_U.nc, Post_1M_grid_U) ,\ 
    30         (${PREFIX_MO}_grid_V.nc, ${R_OUT_OCE_O_M}/${PREFIX}_1M_grid_V.nc, Post_1M_grid_V),\ 
    31         (${PREFIX_MO}_grid_W.nc, ${R_OUT_OCE_O_M}/${PREFIX}_1M_grid_W.nc, Post_1M_grid_W),\ 
    32         (${PREFIX_MO}_diaznl.nc, ${R_OUT_OCE_O_M}/${PREFIX}_1M_diaznl.nc, Post_1M_diaznl),\ 
    33         (${PREFIX_MO}_S_tren.nc, ${R_OUT_OCE_O_M}/${PREFIX}_1M_S_tren.nc, NONE),\ 
    34         (${PREFIX_MO}_T_tren.nc, ${R_OUT_OCE_O_M}/${PREFIX}_1M_T_tren.nc, NONE),\ 
    35         (${PREFIX_MO}_trends.nc, ${R_OUT_OCE_O_M}/${PREFIX}_1M_trends.nc, NONE),\ 
    36         (${PREFIX_MO}_diagap.nc, ${R_OUT_OCE_O_M}/${PREFIX}_1M_diagap.nc, NONE),\ 
    37         (${PREFIX_5D}_grid_T.nc, ${R_OUT_OCE_O_D}/${PREFIX}_5D_grid_T.nc, NONE),\ 
    38         (${PREFIX_5D}_grid_U.nc, ${R_OUT_OCE_O_D}/${PREFIX}_5D_grid_U.nc, NONE),\ 
    39         (${PREFIX_5D}_grid_V.nc, ${R_OUT_OCE_O_D}/${PREFIX}_5D_grid_V.nc, NONE),\ 
    40         (${PREFIX_5D}_grid_W.nc, ${R_OUT_OCE_O_D}/${PREFIX}_5D_grid_W.nc, NONE),\ 
    41         (${PREFIX_5D}_diaznl.nc, ${R_OUT_OCE_O_D}/${PREFIX}_5D_diaznl.nc, NONE),\ 
    42         (${PREFIX_DA}_grid_T.nc, ${R_OUT_OCE_O_D}/${PREFIX}_1D_grid_T.nc, NONE),\ 
    43         (${PREFIX_DA}_grid_U.nc, ${R_OUT_OCE_O_D}/${PREFIX}_1D_grid_U.nc, NONE),\ 
    44         (${PREFIX_DA}_grid_V.nc, ${R_OUT_OCE_O_D}/${PREFIX}_1D_grid_V.nc, NONE),\ 
    45         (${PREFIX_DA}_grid_W.nc, ${R_OUT_OCE_O_D}/${PREFIX}_1D_grid_W.nc, NONE),\ 
    46         (${PREFIX_DA}_diaznl.nc, ${R_OUT_OCE_O_D}/${PREFIX}_1D_diaznl.nc, NONE) 
     31List=   (${PREFIX_NWRITE}_${DATE_OPA}_grid_T.nc, ${R_OUT_OCE_NWRITE}/${PREFIX}_${WF1}_grid_T.nc, Post_1M_grid_T),\ 
     32        (${PREFIX_NWRITE}_${DATE_OPA}_grid_U.nc, ${R_OUT_OCE_NWRITE}/${PREFIX}_${WF1}_grid_U.nc, Post_1M_grid_U) ,\ 
     33        (${PREFIX_NWRITE}_${DATE_OPA}_grid_V.nc, ${R_OUT_OCE_NWRITE}/${PREFIX}_${WF1}_grid_V.nc, Post_1M_grid_V),\ 
     34        (${PREFIX_NWRITE}_${DATE_OPA}_grid_W.nc, ${R_OUT_OCE_NWRITE}/${PREFIX}_${WF1}_grid_W.nc, Post_1M_grid_W),\ 
     35        (${PREFIX_NWRITE}_${DATE_OPA}_diaznl.nc, ${R_OUT_OCE_NWRITE}/${PREFIX}_${WF1}_diaznl.nc, Post_1M_diaznl),\ 
     36        (${PREFIX_NWRITE}_${DATE_OPA}_S_tren.nc, ${R_OUT_OCE_NWRITE}/${PREFIX}_${WF1}_S_tren.nc, NONE),\ 
     37        (${PREFIX_NWRITE}_${DATE_OPA}_T_tren.nc, ${R_OUT_OCE_NWRITE}/${PREFIX}_${WF1}_T_tren.nc, NONE),\ 
     38        (${PREFIX_NWRITE}_${DATE_OPA}_trends.nc, ${R_OUT_OCE_NWRITE}/${PREFIX}_${WF1}_trends.nc, NONE),\ 
     39        (${PREFIX_NWRITE}_${DATE_OPA}_diagap.nc, ${R_OUT_OCE_NWRITE}/${PREFIX}_${WF1}_diagap.nc, NONE),\ 
     40        (${PREFIX_NWRIHF}_${DATE_OPA}_grid_T.nc, ${R_OUT_OCE_NWRIHF}/${PREFIX}_${WF2}_grid_T.nc, NONE),\ 
     41        (${PREFIX_NWRIHF}_${DATE_OPA}_grid_U.nc, ${R_OUT_OCE_NWRIHF}/${PREFIX}_${WF2}_grid_U.nc, NONE),\ 
     42        (${PREFIX_NWRIHF}_${DATE_OPA}_grid_V.nc, ${R_OUT_OCE_NWRIHF}/${PREFIX}_${WF2}_grid_V.nc, NONE) 
    4743 
    4844[Post_1M_grid_T] 
  • CONFIG/trunk/IPSL_ESM_V1/EXP00/COMP/opa.driver

    r103 r187  
    6060    IGCM_debug_PushStack "OCE_Update" 
    6161 
     62    NbFreq=$( echo ${config_OCE_WriteFrequency} | wc -w ) 
     63 
    6264    ##--Write Frequency Purpose .... 
    63     for frequency in ${config_OCE_WriteFrequency} ; do 
     65    frequency=$( echo ${config_OCE_WriteFrequency} | awk "-F " '{print $1}' ) 
     66    factor=$( echo ${frequency} | sed -e "s/[yYmMdD]//" ) 
     67    case ${frequency} in 
     68        1Y|1y) 
     69            (( ORCA_NWRITE = ORCA_NPDT_JOUR * factor * $( IGCM_date_DaysInYear  ${year} ) ))          ; 
     70            PREFIX_NWRITE=${config_UserChoices_JobName}_${factor}y ; 
     71            R_OUT_OCE_NWRITE=${R_OUT_OCE_O_Y} ; 
     72            R_OUT_ICE_NWRITE=${R_OUT_ICE_O_Y} ; 
     73            WF1=${factor}Y ;; 
     74        1M|1m) 
     75            (( ORCA_NWRITE = ORCA_NPDT_JOUR * factor * $( IGCM_date_DaysInMonth ${year} ${month} ) )) ; 
     76            PREFIX_NWRITE=${config_UserChoices_JobName}_${factor}m ; 
     77            R_OUT_OCE_NWRITE=${R_OUT_OCE_O_M} ; 
     78            R_OUT_ICE_NWRITE=${R_OUT_ICE_O_M} ; 
     79            WF1=${factor}M ;; 
     80        *D|*d) 
     81            (( ORCA_NWRITE = ORCA_NPDT_JOUR * factor  ))  ; 
     82            PREFIX_NWRITE=${config_UserChoices_JobName}_${factor}d ; 
     83            R_OUT_OCE_NWRITE=${R_OUT_OCE_O_D} ; 
     84            R_OUT_ICE_NWRITE=${R_OUT_ICE_O_D} ; 
     85            WF1=${factor}D ;; 
     86        *) 
     87            (( ORCA_NWRITE = 0 ))                                                            ;; 
     88    esac 
     89 
     90    ##--Write Frequency Purpose .... concerning the second opa file type. 
     91    if [ ${NbFreq} -eq 2 ] ; then 
     92        frequency=$( echo ${config_OCE_WriteFrequency} | awk "-F " '{print $2}' ) 
     93        factor=$( echo ${frequency} | sed -e "s/[yYmMdD]//" ) 
    6494        case ${frequency} in 
    65             5D|5d) (( ORCA_NWRIHF=ORCA_NPDT_JOUR*5 )) ;; 
    66             1D|1d) (( ORCA_NWRIHF=ORCA_NPDT_JOUR )) ;; 
    67             *)     (( ORCA_NWRIHF=0 )) ;; 
     95            *Y|*y) 
     96                (( ORCA_NWRIHF = ORCA_NPDT_JOUR * factor * $( IGCM_date_DaysInYear  ${year} ) ))          ; 
     97                PREFIX_NWRIHF=${config_UserChoices_JobName}_${factor}y ; 
     98                R_OUT_OCE_NWRIHF=${R_OUT_OCE_O_Y} ; 
     99                WF2=${factor}Y ;; 
     100            *M|*m) 
     101                (( ORCA_NWRIHF = ORCA_NPDT_JOUR * factor * $( IGCM_date_DaysInMonth ${year} ${month} ) )) ; 
     102                PREFIX_NWRIHF=${config_UserChoices_JobName}_${factor}m ; 
     103                R_OUT_OCE_NWRIHF=${R_OUT_OCE_O_M} ; 
     104                WF2=${factor}M ;; 
     105            *D|*d) 
     106                (( ORCA_NWRIHF = ORCA_NPDT_JOUR * factor ))  ; 
     107                PREFIX_NWRIHF=${config_UserChoices_JobName}_${factor}d ; 
     108                R_OUT_OCE_NWRIHF=${R_OUT_OCE_O_D} ; 
     109                WF2=${factor}D ;; 
     110            *) 
     111                (( ORCA_NWRIHF = 0 ))                                                            ;; 
    68112        esac 
    69     done 
    70  
    71     (( ORCA_NWRITE = ORCA_NPDT_JOUR * PeriodLengthInDays )) 
     113    else 
     114        ORCA_NWRIHF=0 
     115    fi 
    72116 
    73117    # Period Length In Days between DateBegin and first day of calendar 0001 01 01 
     
    150194        DATE_OPA=$( cat date.file | \ 
    151195            sed "s/\ ${config_UserChoices_JobName}_[0-9]*[a-z]_\([0-9]*_[0-9]*\)_output\ */\1/g" ) 
    152         PREFIX_DA=${config_UserChoices_JobName}_1d_${DATE_OPA} 
    153         PREFIX_5D=${config_UserChoices_JobName}_5d_${DATE_OPA} 
    154         PREFIX_MO=${config_UserChoices_JobName}_1m_${DATE_OPA} 
     196 
     197        MainPrefix=${config_UserChoices_JobName}_1d_${DATE_OPA} 
     198        SecondPrefix=${config_UserChoices_JobName}_5d_${DATE_OPA} 
    155199    fi 
    156200 
  • CONFIG/trunk/IPSL_ESM_V1/EXP00/COMP/orchidee.card

    r103 r187  
     1[Compatibility] 
     2libIGCM=1.0 
     3 
    14[UserChoices] 
    25 
     
    1821[OutputText] 
    1922List=   (out_orchidee) 
    20 # avec la // : out_orchidee_* 
    2123 
    2224[OutputFiles] 
    2325List=   (sechiba_out.nc, ${R_OUT_SRF_O_M}/${PREFIX}_1M_sechiba_history.nc, Post_1M_sechiba_history) 
    24 #       \ 
    25 #       (stomate_out.nc, ${R_OUT_SRF_O_M}/${PREFIX}_1M_stomate_history.nc, Post_1M_stomate_history) 
    2626 
    2727[Post_1M_sechiba_history] 
     
    3030MonitoringVars = 
    3131TimeSeriesVars = (alb_nir, alb_vis, bqsb, evap, fluxlat, fluxsens, gqsb, netrad, qair, rain, runoff, snow, subli, tair, tsol_max, tsol_min, drainage) 
    32  
    33 [Post_1M_stomate_history] 
    34 Patches = () 
    35 GatherWithInternal = 
    36 MonitoringVars = 
    37 TimeSeriesVars = 
  • CONFIG/trunk/IPSL_ESM_V1/EXP00/PARAM/gcm.def_LMD9671

    r103 r187  
    4646##  Flag  pour la convection (1 pour LMD, 2 pour Tiedtke, 3 KE, 4 KE vect) 
    4747iflag_con=3 
     48## Facteur multiplication des precip convectives dans KE 
     49cvl_corr=1.0 
    4850## longitude en degres du centre du zoom                                  
    4951clon=0. 
  • CONFIG/trunk/IPSL_ESM_V1/EXP00/PARAM/namcouple_ORCA2xLMD9671

    r103 r187  
    2929# 
    3030 $CHANNEL 
    31            MPI2 
     31           <lib_mpi> 
    3232<nproc_atm>  <nproc_atm> 
    33331  1                
     
    477477 INT=1 
    478478# Interpolation method ou parametres mozaic 
    479 wa2o.run 94    4   4 
     479wa2o.run 94    4   158 
    480480# change units from m/s to kg/s, change of sign 
    481481# weights convert from kg/s to kg/m^2/s 
     
    497497 INT=1 
    498498# Interpolation method ou parametres mozaic 
    499 wa2o.run 94    4   4 
     499wa2o.run 94    4   158 
    500500# change units from m/s to kg/s, change of sign 
    501501# weights convert from kg/s to kg/m^2/s 
  • CONFIG/trunk/IPSL_ESM_V1/EXP00/PARAM/orchidee.def

    r103 r187  
    1 # 
    2 #************************************************************************** 
    3 #                    Namelist for ORCHIDEE 
    4 #************************************************************************** 
    51# 
    62# 
    7 #************************************************************************** 
    8 #          OPTIONS NOT SET 
    9 #************************************************************************** 
     3# Parameter file for IPSLCM4_v2 configuration 
     4# See comments : http://forge.ipsl.jussieu.fr/orchidee/ 
    105# 
    11 # 
    12 #************************************************************************** 
    13 #          Management of display in the run of ORCHIDEE 
    14 #************************************************************************** 
    15  
    16 # Model chatting level 
    17 # level of online diagnostics in STOMATE (0-4) 
    18 # With this variable, you can determine how much online information STOMATE 
    19 #  gives during the run. 0 means virtually no info. 
    20 BAVARD = 1 
    21 # default = 1 
    22  
    23 # Flag for debug information 
    24 # This option allows to switch on the output of debug 
    25 #         information without recompiling the code. 
    26 DEBUG_INFO = n 
    27 #default = n 
    28  
    29 # ORCHIDEE will print more messages 
    30 # This flag permits to print more debug messages in the run. 
    31 LONGPRINT = n 
    32 #default = n 
    33  
    34 #--------------------------------------------------------------------- 
    35  
    36 # Should the output follow the ALMA convention 
    37 # If this logical flag is set to true the model 
    38 #  will output all its data according to the ALMA  
    39 #  convention. It is the recommended way to write 
    40 #  data out of ORCHIDEE. 
    41 ALMA_OUTPUT = n 
    42 # default = n 
    43  
    44 # To reset the time coming from SECHIBA restart file 
    45 # This option allows the model to override the time 
    46 #  found in the restart file of SECHIBA with the time 
    47 #  of the first call. That is the restart time of the GCM. 
    48 SECHIBA_reset_time = n 
    49 # default = n 
    50  
    51 #************************************************************************** 
    52 #          Files : incoming / forcing / restart /output 
    53 #************************************************************************** 
    54 # Ancillary files : 
    55 #--------------------------------------------------------------------- 
    56  
    57 # Name of file from which the vegetation map is to be read 
    58 # If !IMPOSE_VEG 
    59 # If LAND_USE  
    60 #   default = pft_new.nc 
    61 #   The name of the file to be opened to read a vegetation 
    62 #   map (in pft) is to be given here.  
    63 # If !LAND_USE 
    64 #   default = ../surfmap/carteveg5km.nc 
    65 #   The name of the file to be opened to read the vegetation 
    66 #   map is to be given here. Usualy SECHIBA runs with a 5kmx5km 
    67 #   map which is derived from the IGBP one. We assume that we have 
    68 #   a classification in 87 types. This is Olson modified by Viovy. 
    69 VEGETATION_FILE = carteveg5km.nc 
    70  
    71  
    72 # Name of file from which the bare soil albedo 
    73 # If !IMPOSE_AZE 
    74 # The name of the file to be opened to read the soil types from  
    75 #  which we derive then the bare soil albedos. This file is 1x1  
    76 #  deg and based on the soil colors defined by Wilson and Henderson-Seller. 
    77 SOILALB_FILE = soils_param.nc 
    78 # default = ../surfmap/soils_param.nc 
    79  
    80 # Name of file from which soil types are read 
    81 # If !IMPOSE_VEG 
    82 # The name of the file to be opened to read the soil types.  
    83 #  The data from this file is then interpolated to the grid of 
    84 #  of the model. The aim is to get fractions for sand loam and 
    85 #  clay in each grid box. This information is used for soil hydrology 
    86 #  and respiration. 
    87 SOILTYPE_FILE = soils_param.nc 
    88 # default = ../surfmap/soils_param.nc 
    89  
    90 # Name of file from which the reference 
    91 # The name of the file to be opened to read 
    92 #  temperature is read 
    93 #  the reference surface temperature. 
    94 #  The data from this file is then interpolated 
    95 #  to the grid of the model. 
    96 #  The aim is to get a reference temperature either 
    97 #  to initialize the corresponding prognostic model 
    98 #  variable correctly (ok_dgvm = TRUE) or to impose it 
    99 #  as boundary condition (ok_dgvm = FALSE) 
    100 REFTEMP_FILE = reftemp.nc 
    101 # default = reftemp.nc 
    102  
    103 # Forcing file name 
    104 # Name of file containing the forcing data 
    105 # This is the name of the file which should be opened 
    106 # for reading the forcing data of the dim0 model. 
    107 # The format of the file has to be netCDF and COADS 
    108 # compliant. Cabauw.nc, islscp_for.nc, WG_cru.nc 
    109 FORCING_FILE = islscp_for.nc 
    110 # default = islscp_for.nc 
    111  
    112 # Input and output restart file for the driver 
    113 #--------------------------------------------------------------------- 
    114  
    115 # Name of restart to READ for initial conditions 
    116 # This is the name of the file which will be opened 
    117 #  to extract the initial values of all prognostic 
    118 #  values of the model. This has to be a netCDF file. 
    119 #  Not truly COADS compliant. NONE will mean that 
    120 #  no restart file is to be expected. 
    121 RESTART_FILEIN = NONE 
    122 # default = NONE 
    123  
    124 # Name of restart files to be created by the driver 
    125 # This variable give the  name for 
    126 #  the restart file. The restart software within 
    127 #  IOIPSL will add .nc if needed 
    128 RESTART_FILEOUT = driver_rest_out.nc 
    129 # default = driver_rest_out.nc 
    130  
    131  
    132 # Input and output restart file for SECHIBA : 
    133 #--------------------------------------------------------------------- 
    134  
    135 # Name of restart to READ for initial conditions 
    136 # This is the name of the file which will be opened 
    137 #  to extract the initial values of all prognostic 
    138 #  values of the model. This has to be a netCDF file. 
    139 #  Not truly COADS compliant. NONE will mean that 
    140 #  no restart file is to be expected. 
    141 SECHIBA_restart_in = _start_sech_ 
    142 # default = NONE 
    143  
    144 # Name of restart files to be created by SECHIBA 
    145 # This variable give the name for the restart files.  
    146 #  The restart software within IOIPSL will add .nc if needed. 
    147 SECHIBA_rest_out = sechiba_rest.nc 
    148 # default = sechiba_rest_out.nc 
    149  
    150 # Input and output restart file for STOMATE : 
    151 #--------------------------------------------------------------------- 
    152  
    153 # Name of restart to READ for initial conditions of STOMATE 
    154 # If STOMATE_OK_STOMATE || STOMATE_WATCHOUT 
    155 # This is the name of the file which will be opened of STOMATE 
    156 #   to extract the initial values of all prognostic values of STOMATE. 
    157 STOMATE_RESTART_FILEIN = NONE 
    158 # default = NONE 
    159  
    160 # Name of restart files to be created by STOMATE 
    161 # If STOMATE_OK_STOMATE || STOMATE_WATCHOUT 
    162 # This is the name of the file which will be opened 
    163 #        to write the final values of all prognostic values 
    164 #        of STOMATE. 
    165 STOMATE_RESTART_FILEOUT = stomate_rest_out.nc 
    166 # default = stomate_restart.nc 
    167  
    168 # Forcing files for TESTSTOMATE and FORCESOIL 
    169 #--------------------------------------------------------------------- 
    170  
    171 # Name of STOMATE's forcing file 
    172 # Name that will be given to STOMATE's offline forcing file 
    173 STOMATE_FORCING_NAME = stomate_forcing.nc 
    174 #default = NONE 
    175  
    176 # Size of STOMATE forcing data in memory (MB) 
    177 # This variable determines how many 
    178 #  forcing states will be kept in memory. 
    179 #  Must be a compromise between memory 
    180 #  use and frequeny of disk access. 
    181 STOMATE_FORCING_MEMSIZE = 50 
    182 # default = 50 
    183  
    184 # Name of STOMATE's carbon forcing file 
    185 # Name that will be given to STOMATE's carbon offline forcing file 
    186 STOMATE_CFORCING_NAME = stomate_Cforcing.nc 
    187 # default = NONE 
    188  
    189  
    190 # Produced forcing file name (SECHIBA puis STOMATE) : 
    191 #--------------------------------------------------------------------- 
    192  
    193 # ORCHIDEE will write out its forcing to a file 
    194 # This flag allows to write to a file all the variables 
    195 #  which are used to force the land-surface. The file  
    196 #  has exactly the same format than a normal off-line forcing 
    197 #  and thus this forcing can be used for forcing ORCHIDEE. 
    198 ORCHIDEE_WATCHOUT = n 
    199 # default = n 
    200  
    201 # Filenane for the ORCHIDEE forcing file 
    202 # If ORCHIDEE_WATCHOUT 
    203 # This is the name of the file in which the 
    204 #  forcing used here will be written for later use.  
    205 WATCHOUT_FILE = orchidee_watchout.nc 
    206 # default = orchidee_watchout.nc 
    207  
    208 # ORCHIDEE will write out with this frequency 
    209 # If ORCHIDEE_WATCHOUT 
    210 # This flag indicates the frequency of the write of the variables.  
    211 DT_WATCHOUT = 1800 
    212 # default = dt 
    213  
    214 # STOMATE does minimum service 
    215 # set to TRUE if you want STOMATE to read 
    216 #  and write its start files and keep track 
    217 #  of longer-term biometeorological variables. 
    218 #  This is useful if OK_STOMATE is not set, 
    219 #  but if you intend to activate STOMATE later. 
    220 #  In that case, this run can serve as a  
    221 #  spinup for longer-term biometeorological 
    222 #  variables. 
    223 STOMATE_WATCHOUT = n 
    224 # default = n 
    225  
    226 # Output file name (SECHIBA and STOMATE) : 
    227 #--------------------------------------------------------------------- 
    228 # Name of file in which the output is going 
    229 # This file is going to be created by the model 
    230 #  to be written 
    231 #  and will contain the output from the model. 
    232 #  This file is a truly COADS compliant netCDF file. 
    233 #  It will be generated by the hist software from 
    234 #  the IOIPSL package. 
    235 OUTPUT_FILE = sechiba_history.nc 
    236 # default = cabauw_out.nc 
    237  
    238 # Flag to switch on histfile 2 for SECHIBA (hi-frequency ?) 
    239 # This Flag switch on the second SECHIBA writing for hi (or low)  
    240 #  frequency writing. This second output is optional and not written 
    241 #  by default. 
    242 SECHIBA_HISTFILE2 = FALSE 
    243 # default  = FALSE 
    244  
    245 # Name of file in which the output number 2 is going 
    246 #   to be written 
    247 # If SECHIBA_HISTFILE2 
    248 # This file is going to be created by the model 
    249 #   and will contain the output 2 from the model. 
    250 SECHIVA_OUTPUT_FILE2 = sechiba_out_2.nc 
    251 # default  = sechiba_out_2.nc 
    252  
    253 # Name of file in which STOMATE's output is going to be written 
    254 # This file is going to be created by the model 
    255 #  and will contain the output from the model. 
    256 #  This file is a truly COADS compliant netCDF file. 
    257 #  It will be generated by the hist software from 
    258 #  the IOIPSL package. 
    259 STOMATE_OUTPUT_FILE = stomate_history.nc 
    260 # default = stomate_history.nc 
    261  
    262 # Write levels for outputs files (number of variables) : 
    263 #--------------------------------------------------------------------- 
    264  
    265 # SECHIBA history output level (0..10) 
    266 # Chooses the list of variables in the history file.  
    267 #  Values between 0: nothing is written; 10: everything is  
    268 #  written are available More details can be found on the web under documentation. 
    269 #  web under documentation. 
    270 SECHIBA_HISTLEVEL = 5 
    271 # default = 5 
    272  
    273 # SECHIBA history 2 output level (0..10) 
    274 # If SECHIBA_HISTFILE2 
    275 # Chooses the list of variables in the history file.  
    276 #   Values between 0: nothing is written; 10: everything is  
    277 #   written are available More details can be found on the web under documentation. 
    278 #   web under documentation. 
    279 # First level contains all ORCHIDEE outputs. 
    280 SECHIBA_HISTLEVEL2 = 1 
    281 # default = 1 
    282  
    283 # STOMATE history output level (0..10) 
    284 #  0: nothing is written; 10: everything is written 
    285 STOMATE_HISTLEVEL = 10 
    286 # default = 10 
    287  
    288 # Write frequency for output files (SECHIBA in seconds et 
    289 # STOMATE in days) : 
    290 #--------------------------------------------------------------------- 
    291 # Frequency in seconds at which to WRITE output 
    292 # This variables gives the frequency the output of 
    293 #  the model should be written into the netCDF file. 
    294 #  It does not affect the frequency at which the 
    295 #  operations such as averaging are done. 
    296 WRITE_STEP = 86400.0 
    297 # default = 86400.0 
    298  
    299 # Frequency in seconds at which to WRITE output 
    300 # If SECHIBA_HISTFILE2 
    301 # This variables gives the frequency the output 2 of 
    302 #   the model should be written into the netCDF file. 
    303 #   It does not affect the frequency at which the 
    304 #   operations such as averaging are done. 
    305 #   That is IF the coding of the calls to histdef 
    306 #   are correct ! 
    307 WRITE_STEP2 = 1800.0 
    308 # default = 1800.0 
    309  
    310 # STOMATE history time step (d) 
    311 # Time step of the STOMATE history file 
    312 # Care : this variable must be higher than DT_SLOW 
    313 STOMATE_HIST_DT = 10. 
    314 # default = 10 
    315  
    316 #************************************************************************** 
    317 #                             Area location 
    318 #************************************************************************** 
    319 #  The model will use the smalest regions from 
    320 #  region specified here and the one of the forcing file. 
    321  
    322 # Western limit of region 
    323 # Western limit of the region we are  
    324 #  interested in. Between -180 and +180 degrees 
    325 LIMIT_WEST = -180. 
    326 # default = -180. 
    327  
    328 # Eastern limit of region 
    329 # Eastern limit of the region we are 
    330 #  interested in. Between -180 and +180 degrees 
    331 LIMIT_EAST = 180. 
    332 # default = 180. 
    333  
    334 # Northern limit of region 
    335 # Northern limit of the region we are 
    336 #  interested in. Between +90 and -90 degrees 
    337 LIMIT_NORTH = 90. 
    338 # default = 90. 
    339  
    340 # Southern limit of region 
    341 # Southern limit of the region we are 
    342 #  interested in. Between 90 and -90 degrees 
    343 LIMIT_SOUTH = -90. 
    344 # default = -90. 
    345  
    346 ##************************************************************************** 
    347 #                       Simulation parameters 
    348 #************************************************************************** 
    349  
    350 # method of forcing 
    351 # A method is proposed by which the first atmospheric 
    352 #   level is not directly forced by observations but 
    353 #   relaxed with a time constant towards observations. 
    354 #   For the moment the methods tends to smooth too much 
    355 #   the diurnal cycle and introduces a time shift. 
    356 #   A more sophisticated method is needed. 
    357 RELAXATION = n 
    358 # default = n 
    359  
    360 # Time constant of the relaxation layer RELAXATION 
    361 # The time constant associated to the atmospheric 
    362 #  conditions which are going to be computed 
    363 #  in the relaxed layer. To avoid too much 
    364 #  damping the value should be larger than 1000. 
    365 RELAX_A = 1000. 
    366 # default = 1000.0 
    367  
    368 # Height at which T and Q are given 
    369 # The atmospheric variables (temperature and specific 
    370 #  humidity) are measured at a specific level. 
    371 #  The height of this level is needed to compute 
    372 #  correctly the turbulent transfer coefficients. 
    373 #  Look at the description of the forcing 
    374 #  DATA for the correct value. 
    375 HEIGHT_LEV1 = 2.0 
    376 # default  = 2.0 
    377  
    378 # Height at which the wind is given 
    379 # The height at which wind is needed to compute 
    380 #  correctly the turbulent transfer coefficients. 
    381 HEIGHT_LEVW = 10.0 
    382 # default  = 10.0 
    383  
    384 #--------------------------------------------------------------------- 
    385 # Weather generator or not : 
    386 #--------------------------------------------------------------------- 
    387  
    388 # Allow weather generator to create data. 
    389 # This flag allows the forcing-reader to generate 
    390 #  synthetic data if the data in the file is too sparse 
    391 #  and the temporal resolution would not be enough to 
    392 #  run the model. 
    393 ALLOW_WEATHERGEN = n 
    394 # default = n 
    395  
    396 # North-South Resolution 
    397 # If ALLOW_WEATHERGEN 
    398 # North-South Resolution of the region we are 
    399 #  interested in. In degrees 
    400 MERID_RES = 2. 
    401 # default = 2. 
    402  
    403 # East-West Resolution 
    404 # If ALLOW_WEATHERGEN 
    405 # East-West Resolution of the region we are 
    406 #  interested in. In degrees 
    407 ZONAL_RES = 2. 
    408 # default = 2. 
    409  
    410 # Use prescribed values 
    411 # If ALLOW_WEATHERGEN 
    412 # If this is set to 1, the weather generator 
    413 #   uses the monthly mean values for daily means. 
    414 #   If it is set to 0, the weather generator 
    415 #   uses statistical relationships to derive daily 
    416 #   values from monthly means. 
    417 IPPREC = 0 
    418 # default = 0 
    419  
    420 # Interpolation  or not IF split is larger than 1 
    421 # Choose IF you wish to interpolate linearly or not. 
    422 NO_INTER = y 
    423 INTER_LIN = n 
    424 # default : 
    425 #  NO_INTER = y 
    426 #  INTER_LIN = n 
    427  
    428 # Exact monthly precipitation 
    429 # If ALLOW_WEATHERGEN 
    430 # If this is set to y, the weather generator 
    431 #   will generate pseudo-random precipitations 
    432 #   whose monthly mean is exactly the prescribed one. 
    433 #   In this case, the daily precipitation (for rainy 
    434 #   days) is constant (that is, some days have 0 precip, 
    435 #   the other days have precip = Precip_month/n_precip, 
    436 #   where n_precip is the prescribed number of rainy days 
    437 #   per month). 
    438 WEATHGEN_PRECIP_EXACT = n 
    439 # default = n 
    440  
    441 # Calling frequency of weather generator (s) 
    442 # Determines how often the weather generator 
    443 #  is called (time step in s). Should be equal 
    444 #  to or larger than Sechiba's time step (say, 
    445 #  up to 6 times Sechiba's time step or so). 
    446 DT_WEATHGEN = 1800. 
    447 # default = 1800. 
    448  
    449 # Conserve net radiation in the forcing 
    450 # When the interpolation is used the net radiation 
    451 #  provided by the forcing is not conserved anymore. 
    452 #  This should be avoided and thus this option should 
    453 #  be TRUE (y). 
    454 #  This option is not used for short-wave if the 
    455 #  time-step of the forcing is longer than an hour. 
    456 #  It does not make sense to try and reconstruct 
    457 #  a diurnal cycle and at the same time conserve the  
    458 #  incoming solar radiation. 
    459 NETRAD_CONS = y 
    460 # default = y 
    461  
    462 # Write weather from generator into a forcing file 
    463 # This flag makes the weather generator dump its 
    464 #  generated weather into a forcing file which can 
    465 #  then be used to get the same forcing on different 
    466 #  machines. This only works correctly if there is 
    467 #  a restart file (otherwise the forcing at the first 
    468 #  time step is slightly wrong). 
    469 DUMP_WEATHER = n 
    470 # défault = n 
    471  
    472 # Name of the file that contains 
    473 #  the weather from generator 
    474 # If DUMP_WEATHER 
    475 DUMP_WEATHER_FILE = weather_dump.nc 
    476 # default = 'weather_dump.nc' 
    477  
    478 # Dump weather data on gathered grid 
    479 # If 'y', the weather data are gathered 
    480 #  for all land points. 
    481 # If DUMP_WEATHER 
    482 DUMP_WEATHER_GATHERED = y 
    483 # default = y 
    484  
    485  
    486 # Read Orbital Parameters 
    487  
    488 # Eccentricity Effect 
    489 # Use prescribed values 
    490 # IF ALLOW_WEATHERGEN 
    491 ECCENTRICITY = 0.016724 
    492 # default = 0.016724 
    493  
    494 # Longitude of perihelie 
    495 # Use prescribed values 
    496 # If ALLOW_WEATHERGEN 
    497 PERIHELIE = 102.04 
    498 # default = 102.04 
    499  
    500 # Use prescribed values 
    501 # If ALLOW_WEATHERGEN 
    502 OBLIQUITY = 23.446 
    503 # default = 23.446 
    504  
    505 #************************************************************************** 
    506 # length of simulation : 
    507 #--------------------------------------------------------------------- 
    508 # Length of the integration in time. 
    509 # Length of integration. By default the entire length 
    510 #        of the forcing is used. The FORMAT of this date can 
    511 #        be either of the following : 
    512 # n   : time step n within the forcing file 
    513 # nS  : n seconds after the first time-step in the file 
    514 # nD  : n days after the first time-step 
    515 # nM  : n month after the first time-step (year of 365 days) 
    516 # nY  : n years after the first time-step (year of 365 days) 
    517 #        Or combinations : 
    518 # nYmM: n years and m month 
    519 TIME_LENGTH = default 
    520 # default = depend on the time length and the number of time step in forcing file 
    521 #         = itau_len = itau_fin-itau_dep 
    522  
    523  
    524 # split time step : 
    525 #--------------------------------------------------------------------- 
    526  
    527 # Splits the timestep imposed by the forcing 
    528 # With this value the time step of the forcing 
    529 #  will be devided. In principle this can be run 
    530 #  in explicit mode but it is strongly suggested 
    531 #  to use the implicit method so that the 
    532 #  atmospheric forcing has a smooth evolution. 
    533 SPLIT_DT = 12 
    534 # default = 12 
    535  
    536 #  Time in the forcing file at which the model is started. 
    537 #  This time give the point in time at which the model 
    538 #  should be started.  
    539 #  If exists, the date of the restart file is use. 
    540 #  The FORMAT of this date can be either of the following : 
    541 #  n   : time step n within the forcing file 
    542 #  nS  : n seconds after the first time-step in the file 
    543 #  nD  : n days after the first time-step 
    544 #  nM  : n month after the first time-step (year of 365 days) 
    545 #  nY  : n years after the first time-step (year of 365 days) 
    546 #      Or combinations : 
    547 #  nYmM: n years and m month 
    548 TIME_SKIP = 0 
    549 # default = 0 
    550  
    551 # Number of time steps per year for carbon spinup 
    552 FORCESOIL_STEP_PER_YEAR = 12 
    553 # default = 12 
    554  
    555 # ??? 
    556 FORCESOIL_NB_YEAR = 1 
    557 # default = 1 
    558  
    559 # Spread the precipitation. 
    560 # Spread the precipitaiton over n steps of the splited forcing time step.  
    561 #  This is ONLY applied if the forcing time step has been splited (SPLIT_DT). 
    562 #  If the value indicated is greater than SPLIT_DT, SPLIT_DT is used for it. 
    563 SPRED_PREC = 1 
    564 # default = 1 
    565  
    566  
    567  
    568 #--------------------------------------------------------------------- 
    569 # Parametrization : 
    570 #--------------------------------------------------------------------- 
    571  
    572 # Activate STOMATE? 
    573 # set to TRUE if STOMATE is to be activated 
    574 STOMATE_OK_STOMATE = n 
    575 # default = n 
    576  
    577 # Activate DGVM? 
    578 # set to TRUE if Dynamic Vegetation DGVM is to be activated 
    579 STOMATE_OK_DGVM = n 
    580 # default = n 
    581  
    582 # Activate CO2? 
    583 # set to TRUE if photosynthesis is to be activated 
    584 STOMATE_OK_CO2 = y 
    585 # default = n 
    586  
    587 # Flag to force the value of atmospheric CO2 for vegetation. 
    588 # If this flag is set to true, the ATM_CO2 parameter is used 
    589 #  to prescribe the atmospheric CO2. 
    590 # This Flag is only use in couple mode. 
    591 FORCE_CO2_VEG = FALSE 
    592 # default = FALSE 
    593  
    594 # Value for atm CO2. 
    595 # If FORCE_CO2_VEG (in not forced mode) 
    596 # Value to prescribe the atm CO2. 
    597 #  For pre-industrial simulations, the value is 286.2 . 
    598 #  348. for 1990 year. 
    599 ATM_CO2 = 350. 
    600 # default = 350. 
    601  
    602  
    603 # Index of grid point for online diagnostics 
    604 # This is the index of the grid point which will be used for online diagnostics. 
    605 STOMATE_DIAGPT = 1 
    606 # default = 1 
    607  
    608 # constant tree mortality 
    609 # If yes, then a constant mortality is applied to trees.  
    610 #  Otherwise, mortality is a function of the trees'  
    611 #  vigour (as in LPJ). 
    612 LPJ_GAP_CONST_MORT = y 
    613 # default = y 
    614  
    615 # no fire allowed 
    616 # With this variable, you can allow or not 
    617 #  the estimation of CO2 lost by fire 
    618 FIRE_DISABLE = n 
    619 # default = n 
    620  
    621  
    622 # parameters describing the surface (vegetation + soil) : 
    623 #--------------------------------------------------------------------- 
    624 # 
    625 # Should the vegetation be prescribed 
    626 # This flag allows the user to impose a vegetation distribution 
    627 #  and its characterisitcs. It is espacially interesting for 0D 
    628 #  simulations. On the globe it does not make too much sense as 
    629 #  it imposes the same vegetation everywhere 
    630 IMPOSE_VEG = n 
    631 # default = n 
    632  
    633 # Flag to use old "interpolation" of vegetation map. 
    634 # IF NOT IMPOSE_VEG 
    635 #  If you want to recover the old (ie orchidee_1_2 branch)  
    636 #   "interpolation" of vegetation map. 
    637 SLOWPROC_VEGET_OLD_INTERPOL = n 
    638 # default = n 
    639  
    640 # Vegetation distribution within the mesh (0-dim mode) 
    641 # If IMPOSE_VEG 
    642 # The fraction of vegetation is read from the restart file. If 
    643 #  it is not found there we will use the values provided here. 
    644 SECHIBA_VEG__01 = 0.2 
    645 SECHIBA_VEG__02 = 0.0 
    646 SECHIBA_VEG__03 = 0.0 
    647 SECHIBA_VEG__04 = 0.0 
    648 SECHIBA_VEG__05 = 0.0 
    649 SECHIBA_VEG__06 = 0.0 
    650 SECHIBA_VEG__07 = 0.0 
    651 SECHIBA_VEG__08 = 0.0 
    652 SECHIBA_VEG__09 = 0.0 
    653 SECHIBA_VEG__10 = 0.8 
    654 SECHIBA_VEG__11 = 0.0 
    655 SECHIBA_VEG__12 = 0.0 
    656 SECHIBA_VEG__13 = 0.0 
    657 # default = 0.2, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.8, 0.0, 0.0, 0.0 
    658  
    659 # Maximum vegetation distribution within the mesh (0-dim mode) 
    660 # If IMPOSE_VEG 
    661 # The fraction of vegetation is read from the restart file. If 
    662 #  it is not found there we will use the values provided here. 
    663 SECHIBA_VEGMAX__01 = 0.2 
    664 SECHIBA_VEGMAX__02 = 0.0 
    665 SECHIBA_VEGMAX__03 = 0.0 
    666 SECHIBA_VEGMAX__04 = 0.0 
    667 SECHIBA_VEGMAX__05 = 0.0 
    668 SECHIBA_VEGMAX__06 = 0.0 
    669 SECHIBA_VEGMAX__07 = 0.0 
    670 SECHIBA_VEGMAX__08 = 0.0 
    671 SECHIBA_VEGMAX__09 = 0.0 
    672 SECHIBA_VEGMAX__10 = 0.8 
    673 SECHIBA_VEGMAX__11 = 0.0 
    674 SECHIBA_VEGMAX__12 = 0.0 
    675 SECHIBA_VEGMAX__13 = 0.0 
    676 # default = 0.2, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.8, 0.0, 0.0, 0.0 
    677  
    678 # LAI for all vegetation types (0-dim mode) 
    679 # If IMPOSE_VEG 
    680 # The maximum LAI used in the 0dim mode. The values should be found 
    681 #  in the restart file. The new values of LAI will be computed anyway 
    682 #  at the end of the current day. The need for this variable is caused 
    683 #  by the fact that the model may stop during a day and thus we have not 
    684 #  yet been through the routines which compute the new surface conditions. 
    685 SECHIBA_LAI__01 = 0. 
    686 SECHIBA_LAI__02 = 8. 
    687 SECHIBA_LAI__03 = 8. 
    688 SECHIBA_LAI__04 = 4. 
    689 SECHIBA_LAI__05 = 4.5 
    690 SECHIBA_LAI__06 = 4.5 
    691 SECHIBA_LAI__07 = 4. 
    692 SECHIBA_LAI__08 = 4.5 
    693 SECHIBA_LAI__09 = 4. 
    694 SECHIBA_LAI__10 = 2. 
    695 SECHIBA_LAI__11 = 2. 
    696 SECHIBA_LAI__12 = 2. 
    697 SECHIBA_LAI__13 = 2. 
    698 # default = 0., 8., 8., 4., 4.5, 4.5, 4., 4.5, 4., 2., 2., 2., 2. 
    699  
    700 # Height for all vegetation types (m) 
    701 # If IMPOSE_VEG 
    702 # The height used in the 0dim mode. The values should be found 
    703 #  in the restart file. The new values of height will be computed anyway 
    704 #  at the end of the current day. The need for this variable is caused 
    705 #  by the fact that the model may stop during a day and thus we have not 
    706 #  yet been through the routines which compute the new surface conditions. 
    707 SLOWPROC_HEIGHT__01 = 0. 
    708 SLOWPROC_HEIGHT__02 = 50. 
    709 SLOWPROC_HEIGHT__03 = 50. 
    710 SLOWPROC_HEIGHT__04 = 30. 
    711 SLOWPROC_HEIGHT__05 = 30. 
    712 SLOWPROC_HEIGHT__06 = 30. 
    713 SLOWPROC_HEIGHT__07 = 20. 
    714 SLOWPROC_HEIGHT__08 = 20. 
    715 SLOWPROC_HEIGHT__09 = 20. 
    716 SLOWPROC_HEIGHT__10 = .2 
    717 SLOWPROC_HEIGHT__11 = .2 
    718 SLOWPROC_HEIGHT__12 = .4 
    719 SLOWPROC_HEIGHT__13 = .4 
    720 # default = 0., 30., 30., 20., 20., 20., 15., 15., 15., .5, .6, 1.0, 1.0 
    721  
    722  
    723 # Fraction of the 3 soil types (0-dim mode) 
    724 # If IMPOSE_VEG 
    725 # Determines the fraction for the 3 soil types 
    726 #  in the mesh in the following order : sand loam and clay. 
    727 SOIL_FRACTIONS__01 = 0.28 
    728 SOIL_FRACTIONS__02 = 0.52 
    729 SOIL_FRACTIONS__03 = 0.20 
    730 # default = 0.28, 0.52, 0.20 
    731  
    732 # Temperature used for the initial guess of LAI 
    733 # If there is no LAI in the restart file, we may need 
    734 #  a temperature that is used to guess the initial LAI. 
    735 SLOWPROC_LAI_TEMPDIAG = 280. 
    736 # default = 280. 
    737  
    738 # Soil level (m) used for canopy development  
    739 # If STOMATE is not activated. 
    740 # The temperature at this soil depth is used to determine the LAI when 
    741 #   STOMATE is not activated. 
    742 SECHIBA_ZCANOP = 0.5 
    743 # default = 0.5 
    744  
    745 # Fraction of other surface types within the mesh (0-dim mode) 
    746 # If IMPOSE_VEG 
    747 # The fraction of ice, lakes, etc. is read from the restart file. If 
    748 #  it is not found there we will use the values provided here. 
    749 #  For the moment, there is only ice. 
    750 SECHIBA_FRAC_NOBIO = 0.0 
    751 # default = 0.0 
    752  
    753 # Fraction of the clay fraction (0-dim mode) 
    754 # If IMPOSE_VEG 
    755 # Determines the fraction of clay in the grid box. 
    756 CLAY_FRACTION = 0.2 
    757 # default = 0.2 
    758  
    759 # Should the surface parameters be prescribed 
    760 # This flag allows the user to impose the surface parameters 
    761 #  (Albedo Roughness and Emissivity). It is espacially interesting for 0D 
    762 #  simulations. On the globe it does not make too much sense as 
    763 #  it imposes the same vegetation everywhere 
    764 IMPOSE_AZE = n 
    765 # default = n 
    766  
    767 # Emissivity of the surface for LW radiation 
    768 # If IMPOSE_AZE 
    769 # The surface emissivity used for compution the LE emission 
    770 #  of the surface in a 0-dim version. Values range between  
    771 #  0.97 and 1.. The GCM uses 0.98. 
    772 CONDVEG_EMIS = 1.0 
    773 # default = 1.0 
    774  
    775 # SW visible albedo for the surface 
    776 # If IMPOSE_AZE 
    777 # Surface albedo in visible wavelengths to be used  
    778 #  on the point if a 0-dim version of SECHIBA is used.  
    779 #  Look at the description of the forcing data for  
    780 #  the correct value. 
    781 CONDVEG_ALBVIS = 0.25 
    782 # default = 0.25 
    783  
    784 # SW near infrared albedo for the surface 
    785 # If IMPOSE_AZE 
    786 # Surface albedo in near infrared wavelengths to be used  
    787 #  on the point if a 0-dim version of SECHIBA is used.  
    788 #  Look at the description of the forcing data for  
    789 #  the correct value. 
    790 CONDVEG_ALBNIR = 0.25 
    791 # default = 0.25 
    792  
    793 # Average method for z0 
    794 # If this flag is set to true (y) then the neutral Cdrag 
    795 #  is averaged instead of the log(z0). This should be 
    796 #  the prefered option. We still wish to keep the other 
    797 #  option so we can come back if needed. If this is 
    798 #  desired then one should set Z0CDRAG_AVE = n 
    799 Z0CDRAG_AVE = y 
    800 # default = y 
    801  
    802 # Surface roughness (m) 
    803 # If IMPOSE_AZE 
    804 # Surface rougness to be used on the point if a 0-dim version 
    805 #  of SECHIBA is used. Look at the description of the forcing   
    806 #  data for the correct value. 
    807 CONDVEG_Z0 = 0.15 
    808 # default = 0.15_stnd 
    809  
    810 # Height to be added to the height of the first level (m) 
    811 # If IMPOSE_AZE 
    812 # ORCHIDEE assumes that the atmospheric level height is counted 
    813 #  from the zero wind level. Thus to take into account the roughness 
    814 #  of tall vegetation we need to correct this by a certain fraction 
    815 #  of the vegetation height. This is called the roughness height in 
    816 #  ORCHIDEE talk. 
    817 ROUGHHEIGHT = 0.0 
    818 # default = 0.0 
    819  
    820 # The snow albedo used by SECHIBA 
    821 # This option allows the user to impose a snow albedo. 
    822 #  Default behaviour is to use the model of snow albedo 
    823 #  developed by Chalita (1993). 
    824 CONDVEG_SNOWA = default 
    825 # default = use the model of snow albedo developed by Chalita 
    826  
    827 # Switch bare soil albedo dependent (if TRUE) on soil wetness 
    828 # If TRUE, the model for bare soil albedo is the old formulation. 
    829 #  Then it depend on the soil dry or wetness. If FALSE, it is the  
    830 #  new computation that is taken, it is only function of soil color. 
     6STOMATE_OK_CO2=TRUE 
     7# STOMATE_OK_STOMATE is not set 
     8# STOMATE_OK_DGVM is not set 
     9# STOMATE_WATCHOUT is not set 
     10SECHIBA_restart_in=_start_sech_ 
     11SECHIBA_rest_out=sechiba_rest.nc 
     12SECHIBA_reset_time=y 
     13# SECHIBA_reset_time is not set 
     14OUTPUT_FILE=sechiba_out.nc 
     15WRITE_STEP=2592000 
     16SECHIBA_HISTLEVEL=5 
     17STOMATE_OUTPUT_FILE=stomate_history.nc 
     18STOMATE_HIST_DT=10. 
     19STOMATE_HISTLEVEL=0 
     20SECHIBA_DAY=0.0 
     21SECHIBA_ZCANOP=0.5 
     22DT_SLOW=86400. 
     23# IMPOSE_VEG is not set 
     24VEGETATION_FILE=carteveg5km.nc 
     25DIFFUCO_LEAFCI=233. 
     26CONDVEG_SNOWA=default 
     27# IMPOSE_AZE is not set 
     28SOILALB_FILE=soils_param.nc 
     29SOILTYPE_FILE=soils_param.nc  
     30ENERBIL_TSURF=280. 
     31HYDROL_SNOW=0.0 
     32HYDROL_SNOWAGE=0.0 
     33HYDROL_SNOWICE=0.0 
     34HYDROL_SNOWICEAGE=0.0 
     35HYDROL_HDRY=1.0 
     36HYDROL_HUMR=1.0 
     37HYDROL_BQSB=default 
     38HYDROL_GQSB=0.0 
     39HYDROL_DSG=0.0 
     40HYDROL_DSP=default 
     41HYDROL_QSV=0.0 
     42HYDROL_OK_HDIFF=n 
     43HYDROL_TAU_HDIFF=1800. 
     44THERMOSOIL_TPRO=280. 
     45RIVER_ROUTING=y 
     46ROUTING_FILE=routing.nc 
     47LAI_MAP=y 
     48LAI_FILE=lai2D.nc 
     49SECHIBA_QSINT=0.02 
    83150ALB_BARE_MODEL = FALSE 
    832 # default = FALSE 
    833  
    834 # Initial snow mass if not found in restart 
    835 # The initial value of snow mass if its value is not found 
    836 #   in the restart file. This should only be used if the model is  
    837 #   started without a restart file. 
    838 HYDROL_SNOW = 0.0 
    839 # default = 0.0 
    840  
    841  
    842 # Initial snow age if not found in restart 
    843 # The initial value of snow age if its value is not found 
    844 #  in the restart file. This should only be used if the model is  
    845 #  started without a restart file. 
    846 HYDROL_SNOWAGE = 0.0 
    847 # default = 0.0 
    848  
    849 # Initial snow amount on ice, lakes, etc. if not found in restart 
    850 # The initial value of snow if its value is not found 
    851 #  in the restart file. This should only be used if the model is  
    852 #  started without a restart file. 
    853 HYDROL_SNOW_NOBIO = 0.0 
    854 # default = 0.0 
    855  
    856 # Initial snow age on ice, lakes, etc. if not found in restart 
    857 # The initial value of snow age if its value is not found 
    858 #  in the restart file. This should only be used if the model is  
    859 #  started without a restart file. 
    860 HYDROL_SNOW_NOBIO_AGE = 0.0 
    861 # default = 0.0 
    862  
    863 # Initial dry soil height if not found in restart for ORCHIDEE_1.3 to 1.5 Tags only. 
    864 # The initial value of dry soil height if its value is not found 
    865 # in the restart file. This should only be used if the model is  
    866 # started without a restart file. 
    867 HYDROL_HDRY = 0.0 
    868 # default = 0.0 
    869  
    870 # Initial soil moisture stress if not found in restart 
    871 # The initial value of soil moisture stress if its value is not found 
    872 #  in the restart file. This should only be used if the model is  
    873 #  started without a restart file. 
    874 HYDROL_HUMR = 1.0 
    875 # default = 1.0 
    876  
    877 # Total depth of soil reservoir 
    878 HYDROL_SOIL_DEPTH = 2. 
    879 # default = 2. 
    880  
    881 # Initial restart deep soil moisture if not found in restart 
    882 # The initial value of deep soil moisture if its value is not found 
    883 #  in the restart file. This should only be used if the model is  
    884 #  started without a restart file. Default behaviour is a saturated soil. 
    885 HYDROL_BQSB = default 
    886 # default = Maximum quantity of water (Kg/M3) * Total depth of soil reservoir = 150. * 2 
    887  
    888 # Initial upper soil moisture if not found in restart 
    889 # The initial value of upper soil moisture if its value is not found 
    890 #  in the restart file. This should only be used if the model is  
    891 #  started without a restart file. 
    892 HYDROL_GQSB = 0.0 
    893 # default = 0.0 
    894  
    895 # Initial upper reservoir depth if not found in restart 
    896 # The initial value of upper reservoir depth if its value is not found 
    897 #  in the restart file. This should only be used if the model is  
    898 #  started without a restart file. 
    899 HYDROL_DSG = 0.0 
    900 # default = 0.0 
    901  
    902 # Initial dry soil above upper reservoir if not found in restart 
    903 # The initial value of dry soil above upper reservoir if its value  
    904 #  in the restart file. This should only be used if the model is  
    905 #  started without a restart file. The default behaviour 
    906 #  is to compute it from the variables above. Should be OK most of  
    907 #  the time. 
    908 HYDROL_DSP = default 
    909 # default = Total depth of soil reservoir - HYDROL_BQSB / Maximum quantity of water (Kg/M3) = 0.0 
    910  
    911 # Initial water on canopy if not found in restart 
    912 # The initial value of moisture on canopy if its value  
    913 #  in the restart file. This should only be used if the model is  
    914 #  started without a restart file. 
    915 HYDROL_QSV = 0.0 
    916 # default = 0.0 
    917  
    918 # Soil moisture on each soil tile and levels 
    919 # The initial value of mc if its value is not found 
    920 #  in the restart file. This should only be used if the model is  
    921 #  started without a restart file. 
    922 HYDROL_MOISTURE_CONTENT = 0.3 
    923 # default = 0.3 
    924  
    925 # US_NVM_NSTM_NSLM 
    926 # The initial value of us (relative moisture) if its value is not found 
    927 #  in the restart file. This should only be used if the model is  
    928 #  started without a restart file. 
    929 US_INIT = 0.0 
    930 # default = 0.0 
    931  
    932 # Coefficient for free drainage at bottom 
    933 # The initial value of free drainage if its value is not found 
    934 #  in the restart file. This should only be used if the model is  
    935 #  started without a restart file. 
    936 FREE_DRAIN_COEF = 1.0, 1.0, 1.0 
    937 # default = 1.0, 1.0, 1.0 
    938  
    939 # Bare soil evap on each soil if not found in restart 
    940 # The initial value of bare soils evap if its value is not found 
    941 #  in the restart file. This should only be used if the model is  
    942 #  started without a restart file. 
    943 EVAPNU_SOIL = 0.0 
    944 # default = 0.0 
    945  
    946  
    947 # Initial temperature if not found in restart 
    948 # The initial value of surface temperature if its value is not found 
    949 #  in the restart file. This should only be used if the model is  
    950 #  started without a restart file. 
    951 ENERBIL_TSURF = 280. 
    952 # default = 280. 
    953  
    954 # Initial Soil Potential Evaporation 
    955 # The initial value of soil potential evaporation if its value  
    956 #  is not found in the restart file. This should only be used if 
    957 #  the model is started without a restart file.  
    958 ENERBIL_EVAPOT = 0.0 
    959 # default = 0.0 
    960  
    961 # Initial soil temperature profile if not found in restart 
    962 # The initial value of the temperature profile in the soil if  
    963 #   its value is not found in the restart file. This should only  
    964 #   be used if the model is started without a restart file. Here 
    965 #   we only require one value as we will assume a constant  
    966 #   throughout the column. 
    967 THERMOSOIL_TPRO = 280. 
    968 # default = 280. 
    969  
    970 # Initial leaf CO2 level if not found in restart 
    971 # The initial value of leaf_ci if its value is not found 
    972 #  in the restart file. This should only be used if the model is 
    973 #  started without a restart file. 
    974 DIFFUCO_LEAFCI = 233. 
    975 # default = 233. 
    976  
    977  
    978 # Keep cdrag coefficient from gcm. 
    979 # Set to .TRUE. if you want q_cdrag coming from GCM. 
    980 #  Keep cdrag coefficient from gcm for latent and sensible heat fluxes. 
    981 #  TRUE if q_cdrag on initialization is non zero (FALSE for off-line runs). 
     51PERCENT_THROUGHFALL_PFT = 30., 30., 30., 30., 30., 30., 30., 30., 30., 30., 30., 30., 30. 
     52RVEG_PFT = .5, .5, .5, .5, .5, .5, .5, .5, .5, .5, .5, .5, .5 
    98253CDRAG_FROM_GCM = .TRUE. 
    983 # default =  IF q_cdrag == 0 ldq_cdrag_from_gcm = .FALSE. ELSE .TRUE. 
    984  
    985  
    986 # Artificial parameter to increase or decrease canopy resistance 
    987 # Add from Nathalie - the 28 of March 2006 - advice from Fred Hourdin 
    988 # By PFT. 
    989 RVEG_PFT = .5, .5, .5, .5, .5, .5, .5, .5, .5, .5, .5, .5, .5 
    990 # default = 1. 
    991  
    992  
    993 # Interception reservoir coefficient for ORCHIDEE_1.3 to 1.5 Tags only. 
    994 # Transforms leaf area index into size of interception reservoir 
    995 #  for slowproc_derivvar or stomate. 
    996 SECHIBA_QSINT = 0.02 
    997 # default = 0.1 
    998  
    999 #************************************************************************** 
    1000 # LAI 
    1001 #************************************************************************** 
    1002  
    1003 # Read the LAI map 
    1004 # It is possible to read a 12 month LAI map which will 
    1005 #  then be interpolated to daily values as needed. 
    1006 #  If n => type_of_lai (constant_veg.f90) 
    1007 #     - mean    : lai(ji,jv) = undemi * (llaimax(jv) + llaimin(jv)) 
    1008 #     - inter   : llaimin(jv) + tempfunc(stempdiag(ji,lcanop)) * (llaimax(jv) - llaimin(jv)) 
    1009 LAI_MAP = y 
    1010 # default = n 
    1011  
    1012 # Name of file from which the vegetation map is to be read 
    1013 # If LAI_MAP 
    1014 # The name of the file to be opened to read the LAI 
    1015 #  map is to be given here. Usualy SECHIBA runs with a 5kmx5km 
    1016 #  map which is derived from a Nicolas VIOVY one.  
    1017 LAI_FILE = lai2D.nc 
    1018 # default = ../surfmap/lai2D.nc 
    1019  
    1020 # Flag to use old "interpolation" of LAI 
    1021 # If LAI_MAP 
    1022 # If you want to recover the old (ie orchidee_1_2 branch)  
    1023 # "interpolation" of LAI map. 
    1024 SLOWPROC_LAI_OLD_INTERPOL = n 
    1025 # default = n 
    1026  
    1027 #************************************************************************** 
    1028 # LAND_USE 
    1029 #************************************************************************** 
    1030  
    1031 # Read a land_use vegetation map 
    1032 # pft values are needed, max time axis is 293 
    1033 LAND_USE = n 
    1034 # default = n 
    1035  
    1036 # Year of the land_use vegetation map readed 
    1037 # year off the pft map 
    1038 # If LAND_USE 
    1039 VEGET_YEAR = 282 
    1040 # default = 282 
    1041  
    1042 # Update vegetation frequency 
    1043 # The veget datas will be update each this time step. 
    1044 # If LAND_USE 
    1045 VEGET_LENGTH = 1Y 
    1046 # default = 1Y 
    1047  
    1048 # treat land use modifications 
    1049 # With this variable, you can use a Land Use map 
    1050 # to simulate anthropic modifications such as    
    1051 # deforestation.                                 
    1052 # If LAND_USE 
    1053 LAND_COVER_CHANGE = n 
    1054 # default = n 
    1055  
    1056 #************************************************************************** 
    1057  
    1058 # agriculture allowed? 
    1059 # With this variable, you can determine 
    1060 #  whether agriculture is allowed 
    1061 AGRICULTURE = y 
    1062 # default = y 
    1063  
    1064 # herbivores allowed? 
    1065 # With this variable, you can activate herbivores  
    1066 HERBIVORES = n 
    1067 # default = n 
    1068  
    1069 # treat expansion of PFTs across a grid cell? 
    1070 # With this variable, you can determine 
    1071 #  whether we treat expansion of PFTs across a 
    1072 #  grid cell. 
    1073 TREAT_EXPANSION = n 
    1074 # default = n 
    1075  
    1076 #************************************************************************** 
    1077  
    1078 # Time within the day simulated 
    1079 # This is the time spent simulating the current day. This variable is 
    1080 #  prognostic as it will trigger all the computations which are 
    1081 #  only done once a day. 
    1082 SECHIBA_DAY = 0.0 
    1083 # default = 0.0 
    1084  
    1085 # Time step of STOMATE and other slow processes 
    1086 # Time step (s) of regular update of vegetation 
    1087 #  cover, LAI etc. This is also the time step 
    1088 #  of STOMATE. 
    1089 DT_SLOW = 86400. 
    1090 # default = un_jour = 86400. 
    1091  
    1092 #************************************************************************** 
    1093  
    1094 # Allows to switch on the multilayer hydrology of CWRR 
    1095 # This flag allows the user to decide if the vertical 
    1096 #  hydrology should be treated using the multi-layer  
    1097 #  diffusion scheme adapted from CWRR by Patricia de Rosnay. 
    1098 #  by default the Choisnel hydrology is used. 
    1099 HYDROL_CWRR = n 
    1100 # default = n 
    1101  
    1102 # do horizontal diffusion? 
    1103 # If TRUE, then water can diffuse horizontally between 
    1104 #  the PFTs' water reservoirs. 
    1105 HYDROL_OK_HDIFF = n 
    1106 # default = n 
    1107   
    1108  
    1109 # time scale (s) for horizontal diffusion of water 
    1110 # If HYDROL_OK_HDIFF 
    1111 # Defines how fast diffusion occurs horizontally between 
    1112 #  the individual PFTs' water reservoirs. If infinite, no 
    1113 #  diffusion. 
    1114 HYDROL_TAU_HDIFF = 86400. 
    1115 # default = 86400. 
    1116  
    1117 # Percent of precip that is not intercepted by the canopy (only for TAG 1.6). 
    1118 # During one rainfall event, PERCENT_THROUGHFALL% of the incident rainfall 
    1119 #  will get directly to the ground without being intercepted. 
    1120 PERCENT_THROUGHFALL = 30. 
    1121 # default = 30. 
    1122  
    1123 # Percent by PFT of precip that is not intercepted by the canopy (since TAG 1.8). 
    1124 # During one rainfall event, PERCENT_THROUGHFALL_PFT% of the incident rainfall 
    1125 #  will get directly to the ground without being intercepted, for each PFT.. 
    1126 PERCENT_THROUGHFALL_PFT = 30., 30., 30., 30., 30., 30., 30., 30., 30., 30., 30., 30., 30. 
    1127 # default = 30., 30., 30., 30., 30., 30., 30., 30., 30., 30., 30., 30., 30. 
    1128  
    1129  
    1130 # Decides if we route the water or not 
    1131 # This flag allows the user to decide if the runoff 
    1132 #  and drainage should be routed to the ocean 
    1133 #  and to downstream grid boxes. 
    1134 RIVER_ROUTING = y 
    1135 # default = n 
    1136  
    1137 # Name of file which contains the routing information 
    1138 # The file provided here should allow the routing module to 
    1139 #  read the high resolution grid of basins and the flow direction  
    1140 #  from one mesh to the other. 
    1141 ROUTING_FILE = routing.nc 
    1142 # default = routing.nc 
    1143  
    1144 # Time step of th routing scheme 
    1145 # If RIVER_ROUTING 
    1146 # This values gives the time step in seconds of the routing scheme.  
    1147 #   It should be multiple of the main time step of ORCHIDEE. One day 
    1148 #   is a good value. 
    1149 ROUTING_TIMESTEP = 86400 
    1150 # default = 86400 
    1151  
    1152 # Number of rivers  
    1153 # If RIVER_ROUTING 
    1154 # This parameter chooses the number of largest river basins 
    1155 #  which should be treated as independently as rivers and not 
    1156 #  flow into the oceans as diffusion coastal flow. 
    1157 ROUTING_RIVERS = 50 
    1158 # default = 50 
    1159  
    1160 # Should we compute an irrigation flux  
    1161 # This parameters allows the user to ask the model 
    1162 #  to compute an irigation flux. This performed for the 
    1163 #  on very simple hypothesis. The idea is to have a good 
    1164 #  map of irrigated areas and a simple function which estimates 
    1165 #  the need to irrigate. 
    1166 DO_IRRIGATION = n 
    1167 # default = n 
    1168  
    1169 # Name of file which contains the map of irrigated areas 
    1170 # If IRRIGATE 
    1171 # The name of the file to be opened to read the field 
    1172 #  with the area in m^2 of the area irrigated within each 
    1173 #  0.5 0.5 deg grid box. The map currently used is the one 
    1174 #  developed by the Center for Environmental Systems Research  
    1175 #  in Kassel (1995). 
    1176 IRRIGATION_FILE = irrigated.nc 
    1177 # default = irrigated.nc 
    1178  
    1179 # Should we include floodplains  
    1180 # This parameters allows the user to ask the model 
    1181 #  to take into account the flood plains and return  
    1182 #  the water into the soil moisture. It then can go  
    1183 #  back to the atmopshere. This tried to simulate  
    1184 #  internal deltas of rivers. 
    1185 DO_FLOODPLAINS = n 
    1186 # default = n 
    1187  
    1188 #************************************************************************** 
  • CONFIG/trunk/IPSL_ESM_V1/EXP00/PARAM/physiq.def

    r103 r187  
    4343# parametres outputs 
    4444# 
     45# Flag pour hf : ok_hf=y/n 
     46ok_hf=n 
     47# Flag pour les variables de sorties : type_run=CLIM/AMIP 
     48# Pour ENSEMBLES mettre type_run=CLIM 
     49type_run=CLIM 
    4550#niveau de sortie "hf" lev_histhf avec 
    4651# - lev_histhf=3 => defaut 
     
    6469R_incl = 23.441 
    6570solaire = 1365. 
    66 co2_ppm = 280. 
     71co2_ppm = 348. 
    6772#RCO2 = co2_ppm * 1.0e-06  * 44.011/28.97 
    6873#RCO2 = 348. * 1.0e-06  * 44.011/28.97 
    6974#RCO2 =   5.286789092164308E-04 
    7075#RCO2 = 425.43e-06 
    71 CH4_ppb = 790. 
     76CH4_ppb = 1650. 
    7277#RCH4 = 1.65E-06* 16.043/28.97 
    7378#RCH4 =   9.137366240938903E-07 
    74 N2O_ppb = 270. 
     79N2O_ppb = 306. 
    7580#RN2O = 306.E-09* 44.013/28.97 
    7681#RN2O =    4.648939592682085E-07 
    77 CFC11_ppt = 0. 
     82CFC11_ppt = 280. 
    7883#RCFC11 = 280.E-12* 137.3686/28.97 
    7984#RCFC11 =    1.327690990680013E-09 
    80 CFC12_ppt = 0. 
     85CFC12_ppt = 484. 
    8186#RCFC12 = 484.E-12* 120.9140/28.97 
    8287#RCFC12 =    2.020102726958923E-09 
  • CONFIG/trunk/IPSL_ESM_V1/EXP00/PARAM/run.def

    r103 r187  
    1919## periode de stockage fichier histmoy (en jour) 
    2020periodav=1. 
     21## activation du calcul d equilibrage de charge 
     22adjust=n 
  • CONFIG/trunk/IPSL_ESM_V1/EXP00/config.card

    r138 r187  
    33# 
    44#======================================================================== 
     5#D-- Compatibility - 
     6[Compatibility] 
     7libIGCM=1.0 
    58#D-- UserChoices - 
    69[UserChoices] 
    710#============================ 
    811#-- (<8 chars MAX for JobName) 
    9 JobName=ESM_00 
     12JobName=ESM_100 
    1013LongName="script for ipsl_esm_v1 config" 
    1114TagName=IPSL_ESM_V1 
     
    2730#-- Total Number of Processors (minimum is 2 for a coupled configuration) 
    2831#JobNumProcTot=2 
    29 JobNumProcTot=6 
    30 #============================  
    31 #-- Run Options : here for oasis coupling  
    32 JobRunOptions='"-np 1 -max_np "${NUM_PROC_OASIS}' 
     32JobNumProcTot=4 
     33#============================ 
     34#-- Run Options : here for oasis coupling 
     35JobRunOptions='-f' 
    3336#============================ 
    3437# config.card configuration options : 
    35 R_INIT=/dmnfs/p86cozic/IGCM/INIT 
    36 R_BC=/dmnfs/p86cozic/IGCM/BC 
    37 #PeriodNb=2 
     38R_INIT=/dmnfs/p86cozic/IGCM_NEW/INIT 
     39R_BC=/dmnfs/p86cozic/IGCM_NEW/BC 
    3840#============================ 
    3941 
     
    4345#D- For each component, Name of component, Tag of component 
    4446ATM= (lmdz, LMDZ.4) 
    45 SRF= (orchidee, ORCHIDEE.1.3) 
     47SRF= (orchidee, ORCHIDEE.1.9) 
    4648OCE= (opa, OPA8.2) 
    4749ICE= (lim, LIM.1) 
     
    5254#D-- Executable - 
    5355[Executable] 
    54 Name=oasis 
     56Name=run_file 
    5557#D- For each component, Real name of executable, Name of executable for oasis 
    5658ATM= (gcm.e, lmdz.x) 
    57 #-- IPSL-CM4_IPCC =>? mod.def 
    5859SRF= ("", "") 
    5960OCE= (opa, opa.xx) 
     
    6970#D- Last day of the experience used as restart 
    7071RestartDate=1869-12-30 
    71 #D- Define restart simulation name (=> JOB_OS) 
     72#D- Define restart simulation name 
    7273RestartJobName=CD1 
    73 #D- Path Server Group Login (=> PSGL) 
    74 RestartPath=/u/rech/ces/rces452/IGCM_OUT/IPSLCM4_v1_OASIS3 
     74#D- Path Server Group Login 
     75RestartPath=${ARCHIVE}/IGCM_OUT/IPSL_ESM_V1 
    7576# Attention login depend de la machine 
    7677 
     
    102103# Define restart simulation name 
    103104RestartJobName=2L18 
    104 RestartPath=${ARCHIVE}/p86denv/SORTIES_CPL_IPSL 
     105RestartPath=${ARCHIVE}/IGCM_OUT/IPSL_ESM_V1 
    105106# Old component name for restart (if empty, use new name) 
    106107OldName= 
     
    115116# Define restart simulation name 
    116117RestartJobName=2L18 
    117 RestartPath=${ARCHIVE}/p86denv/SORTIES_CPL_IPSL 
     118RestartPath=${ARCHIVE}/IGCM_OUT/IPSL_ESM_V1 
    118119# Old component name for restart (if empty, use new name) 
    119120OldName= 
     
    128129# Define restart simulation name 
    129130RestartJobName=2L18            
    130 RestartPath=${ARCHIVE}/p86denv/SORTIES_CPL_IPSL 
     131RestartPath=${ARCHIVE}/IGCM_OUT/IPSL_ESM_V1 
    131132# Old component name for restart (if empty, use new name) 
    132133OldName= 
     
    141142# Define restart simulation name 
    142143RestartJobName=2L18 
    143 RestartPath=${ARCHIVE}/p86denv/SORTIES_CPL_IPSL 
     144RestartPath=${ARCHIVE}/IGCM_OUT/IPSL_ESM_V1 
    144145#-- Old component name for restart (if empty, use new name) 
    145146OldName= 
     
    154155# Define restart simulation name 
    155156RestartJobName=2L18 
    156 RestartPath=${ARCHIVE}/p86denv/SORTIES_CPL_IPSL 
     157RestartPath=${ARCHIVE}/IGCM_OUT/IPSL_ESM_V1 
    157158# Old component name for restart (if empty, use new name) 
    158159OldName= 
     
    167168#D-- Define restart simulation name 
    168169RestartJobName=Test_nvSc_AER 
    169 RestartPath=/dmnfs11/p86cozic/IGCM_OUT/AER 
     170RestartPath=${ARCHIVE}/IGCM_OUT/IPSL_ESM_V1 
    170171#D-- Old component name for restart (if empty, use new name) 
    171172OldName="" 
Note: See TracChangeset for help on using the changeset viewer.