Ignore:
Timestamp:
10/07/11 13:59:22 (13 years ago)
Author:
acosce
Message:

Add Stomate to LMDZORINCA config

Location:
CONFIG/LMDZORINCA/tags/LMDZORINCA_v2/EXP_NMHC_AER
Files:
2 added
5 edited

Legend:

Unmodified
Added
Removed
  • CONFIG/LMDZORINCA/tags/LMDZORINCA_v2/EXP_NMHC_AER/COMP/lmdz.card

    r1557 r1565  
    1515[InitialStateFiles] 
    1616# start and startphy are the clim files computed with CLIM_CREATE (create_etat0_limit) job 
    17 List=   () 
     17List=   (/dmnfs/cont003/p86cozic/IGCM_OUT/LMDZORINCA/NMHC_AER/LOI_AR5WERF_E0_2006/ATM/Restart/LOI_AR5WERF_E0_2006_20051231_restart.nc, start.nc),\ 
     18        (/dmnfs/cont003/p86cozic/IGCM_OUT/LMDZORINCA/NMHC_AER/LOI_AR5WERF_E0_2006/ATM/Restart/LOI_AR5WERF_E0_2006_20051231_restartphy.nc, startphy.nc) 
    1819 
    1920[BoundaryFiles] 
  • CONFIG/LMDZORINCA/tags/LMDZORINCA_v2/EXP_NMHC_AER/COMP/orchidee.card

    r1064 r1565  
    44[UserChoices] 
    55NEWHYDROL=n 
    6 sechiba_LEVEL=10 
     6sechiba_LEVEL=11 
     7VEGET_UPDATE=1Y 
     8LAND_COVER_CHANGE=y 
    79 
    810[InitialStateFiles] 
    9 List=   (${R_INIT}/SRF/LMDZORINCA/soils_param.nc, .), \ 
    10         (${R_INIT}/SRF/LMDZORINCA/routing.nc, .) 
     11List=   (${R_INIT}/SRF/LMDZORINCA/soils_param.nc,      . ), \ 
     12        (${R_INIT}/SRF/LMDZORINCA/routing.nc,          . ) 
    1113 
    1214[BoundaryFiles] 
    13 List=   () 
    14 ListNonDel= (${R_BC}/SRF/LMDZORINCA/PFTmap.20C3M.nc, PFTmap.nc) 
     15List= (${BC_SRF_PFTmap}, PFTmap.nc) 
     16ListNonDel= () 
     17#Tous les 5 ans : 
     18#ListNonDel= (${BC_SRF_PFTmap}, PFTmap.nc) 
     19 
     20[SmoothFiles] 
     21List=() 
    1522 
    1623[ParametersFiles] 
     
    2229 
    2330[OutputText] 
    24 List=   (out_orchidee) 
     31List=   (out_orchidee, orchidee.def) 
    2532 
    2633[OutputFiles] 
    27 List=   (sechiba_history.nc, ${R_OUT_SRF_O_M}/${PREFIX}_1M_sechiba_history.nc, Post_1M_sechiba_history) \ 
    28         (sechiba_out_2.nc,   ${R_OUT_SRF_O_M}/${PREFIX}_1M_sechiba_out2.nc,    NONE) \ 
    29         (watchout.nc,        ${R_OUT_SRF_O_M}/${PREFIX}_1M_watchout.nc,        NONE) 
     34List=   (sechiba_history.nc, ${R_OUT_SRF_O_M}/${PREFIX}_1M_sechiba_history.nc, Post_1M_sechiba_history), \ 
     35        (sechiba_out_2.nc, ${R_OUT_SRF_O_H}/${PREFIX}_HF_sechiba_out_2.nc, Post_HF_sechiba_out_2), \ 
     36        (orchidee_watchout.nc, ${R_OUT_SRF_O_M}/${PREFIX}_1M_watchout.nc, NONE) 
    3037 
    3138[Post_1M_sechiba_history] 
    32 Patches = () 
    33 GatherWithInternal = (lon, lat, veget, time_counter, time_counter_bnds, Areas,Contfrac) 
    34 TimeSeriesVars2D = (nobiofrac, alb_nir, alb_vis, bqsb, evap, fluxlat, fluxsens, gqsb, netrad, qair, rain, runoff, snow, snownobio, snowf, subli, tair, temp_sol, tsol_max, tsol_min, drainage) 
    35 ChunckJob2D = NONE 
    36 TimeSeriesVars3D = (lai, maxvegetfrac, vegetfrac, CO2FLUX) 
    37 ChunckJob3D = NONE 
     39Patches= () 
     40GatherWithInternal= (lon, lat, veget, time_counter, time_counter_bnds, Areas, Contfrac) 
     41TimeSeriesVars2D= (nobiofrac, alb_nir, alb_vis, bqsb, evap, fluxlat, fluxsens, gqsb, netrad, qair, rain, runoff, snow, snownobio, snowf, subli, tair, temp_sol, tsol_max, tsol_min, drainage,mrsos, mrso, mrros, mrro, prveg, evspsblveg, evspsblsoi, tran, treeFrac, grassFrac, cropFrac, baresoilFrac, residualFrac) 
     42ChunckJob2D= NONE 
     43TimeSeriesVars3D= (lai, maxvegetfrac, vegetfrac, CO2FLUX, ptn, nee) 
     44ChunckJob3D= NONE 
     45Seasonal= ON 
     46 
     47[Post_HF_sechiba_out_2] 
     48Patches= () 
     49GatherWithInternal= (lon, lat, veget, time_counter, time_counter_bnds, Areas, Contfrac) 
     50TimeSeriesVars2D= (mrsos, mrro) 
     51ChunckJob2D= 50Y 
     52TimeSeriesVars3D= (ptn) 
     53ChunckJob3D= 10Y 
     54Seasonal= OFF 
  • CONFIG/LMDZORINCA/tags/LMDZORINCA_v2/EXP_NMHC_AER/COMP/orchidee.driver

    r858 r1565  
    2424    RESOL_SRF=ALL 
    2525 
     26    for frequency in ${config_SRF_WriteFrequency} ; do 
     27        case ${frequency} in 
     28            HF|hf) SRF_ok_hf=y ;; 
     29        esac 
     30    done 
     31 
    2632    IGCM_debug_PopStack "SRF_Initialize" 
     33} 
     34 
     35#----------------------------------------------------------------- 
     36function SRF_PeriodStart 
     37{ 
     38 
     39    IGCM_debug_PushStack "SRF_PeriodStart" 
     40    if [ ${year} -gt 2004 ] ; then 
     41 
     42        PFTmap="${R_BC}/SRF/LMDZORINCA/ORCHIDEE_PFTmap_2006to2100_AR5_RCP8.5_MESSAGE.nc" 
     43        if [ ${year} -eq 2005 ] ; then 
     44            IGCM_sys_Rm -f PFTmap.nc 
     45        fi 
     46 
     47    else 
     48        PFTmap="${R_BC}/SRF/LMDZORINCA/PFTmap_1850to2005_AR5_LUHa.rc2.nc" 
     49    fi 
     50 
     51    BC_SRF_PFTmap=${PFTmap} 
     52 
     53    IGCM_debug_PopStack "SRF_PeriodStart" 
     54 
    2755} 
    2856 
     
    3462    typeset SECHIBA_WRITE_STEP 
    3563 
    36     case ${config_SRF_WriteFrequency} in 
     64    SRF_WriteFrequency=$( echo ${config_SRF_WriteFrequency} | sed -e 's/\([0-9]*[yYmMdDs]\).*/\1/' )  
     65    case ${SRF_WriteFrequency} in 
    3766        *Y|*y)  
    38             WriteInYears=$( echo ${1} | awk -F '[yY]' "{print ${config_SRF_WriteFrequency}}" ) 
    39             PeriodLengthInYears=$( echo ${1} | awk -F '[yY]' "{print ${config_UserChoices_PeriodLength}}" ) 
     67            WriteInYears=$( echo ${SRF_WriteFrequency} | awk -F '[yY]' '{print $1}' ) 
     68            PeriodLengthInYears=$( echo ${config_UserChoices_PeriodLength} | awk -F '[yY]' '{print $1}' ) 
    4069            (( SECHIBA_WRITE_STEP = PeriodLengthInDays * WriteInYears / PeriodLengthInYears * 86400 )) ;; 
    41         *M|*m)  
    42             WriteInMonths=$( echo ${1} | awk -F '[mM]' "{print ${config_SRF_WriteFrequency}}" ) 
     70        1M) 
    4371            case ${config_UserChoices_PeriodLength} in 
    4472            *Y|*y) 
    45                 PeriodLengthInYears=$( echo ${1} | awk -F '[yY]' "{print ${config_UserChoices_PeriodLength}}" ) 
     73                SECHIBA_WRITE_STEP=-1. 
     74                ;; 
     75            *M|*m) 
     76                SECHIBA_WRITE_STEP=-1. 
     77                ;; 
     78            *) 
     79                (( SECHIBA_WRITE_STEP = $( IGCM_date_DaysInMonth $year $month ) * 86400 )) 
     80                ;; 
     81            esac 
     82            ;; 
     83        *M|*m)  
     84            WriteInMonths=$( echo ${SRF_WriteFrequency} | awk -F '[mM]' '{print $1}' ) 
     85            case ${config_UserChoices_PeriodLength} in 
     86            *Y|*y) 
     87                PeriodLengthInYears=$( echo ${config_UserChoices_PeriodLength} | awk -F '[yY]' '{print $1}' ) 
    4688                (( SECHIBA_WRITE_STEP = PeriodLengthInDays * 86400 / PeriodLengthInYears / 12  )) 
    4789                ;; 
    4890            *M|*m) 
    49                 PeriodLengthInMonths=$( echo ${1} | awk -F '[mM]' "{print ${config_UserChoices_PeriodLength}}" ) 
     91                PeriodLengthInMonths=$( echo ${config_UserChoices_PeriodLength} | awk -F '[mM]' '{print $1}' ) 
    5092                (( SECHIBA_WRITE_STEP = PeriodLengthInDays * WriteInMonths  * 86400 / PeriodLengthInMonths  )) 
    5193                ;; 
     
    60102            (( SECHIBA_WRITE_STEP = 86400 )) ;; 
    61103        *s) 
    62             WriteInSeconds=$( echo ${1} | awk -F '[s]' "{print ${config_SRF_WriteFrequency}}" ) 
     104            WriteInSeconds=$( echo ${SRF_WriteFrequency} | awk -F '[s]' '{print $1}' ) 
    63105            (( SECHIBA_WRITE_STEP = WriteInSeconds )) ;; 
    64106        *)  
    65             IGCM_debug_Exit "SRF_Update " ${config_SRF_WriteFrequency} " invalid WriteFrequency : choose in 1Y, 1M, 5D, 1D."  
     107            IGCM_debug_Exit "SRF_Update " ${SRF_WriteFrequency} " invalid WriteFrequency : choose in 1Y, 1M, 5D, 1D."  
    66108            IGCM_debug_Verif_Exit ;; 
    67109    esac 
    68110 
     111    # Dans le cas de test : 
     112#    if [ ${year} -eq 1979 ] ; then 
     113#       ORCHIDEE_sed VEGET_REINIT y 
     114#       ORCHIDEE_sed VEGET_YEAR 131 
     115#    fi 
     116#     
     117#    # Dans le cas historique : 
     118#    if [ ${year} -eq 2005 ] ; then 
     119#       ORCHIDEE_sed VEGET_REINIT y 
     120#       ORCHIDEE_sed VEGET_YEAR 1 
     121##      ORCHIDEE_sed VEGET_UPDATE '0Y' 
     122#    fi 
     123 
     124#     # Dans le cas sans changement de carte aprÚs 2005  
     125#     if [ ${year} -gt 2005 ] ; then 
     126#       ORCHIDEE_sed VEGET_UPDATE '0Y' 
     127#     fi 
     128  
    69129    ORCHIDEE_sed HYDROL_CWRR ${orchidee_UserChoices_NEWHYDROL} 
    70130 
     
    72132    ORCHIDEE_sed SECHIBA_HISTLEVEL ${orchidee_UserChoices_sechiba_LEVEL} 
    73133 
    74 ##    if [ ${year} -eq 1950 ] ; then 
    75 #       ORCHIDEE_sed SECHIBA_HISTFILE2 y 
    76 #       ORCHIDEE_sed SECHIBA_HISTLEVEL2 10 
    77 #       ORCHIDEE_sed WRITE_STEP2 86400.0 
    78 ##    fi 
     134    # Outputs HF in HISTFILE2 if required 
     135    [ X${SRF_ok_hf} = Xy ] && ORCHIDEE_sed SECHIBA_HISTFILE2 y 
     136    ORCHIDEE_sed SECHIBA_HISTLEVEL2 1 
     137    ORCHIDEE_sed WRITE_STEP2 10800.0 
    79138 
    80139    if ( [ ${CumulPeriod} -ne 1 ] || [ "${config_SRF_Restart}" != "n" ] ) ; then 
    81140        ORCHIDEE_sed SECHIBA_restart_in sechiba_rest_in.nc 
     141    fi 
     142 
     143# Desactaivate river routing for high resolution or very low resolution 
     144    if ( [ ${RESOL_ATM_3D} = "280x280x39-ZAf" ] || [ ${RESOL_ATM_3D} = "360x180x19-ZIn" ] || [ ${RESOL_ATM_3D} = "32x24x11" ] ) ; then 
     145        ORCHIDEE_sed  RIVER_ROUTING n 
     146    fi 
     147 
     148# Parameters set in orchidee.card 
     149    if [ ! X${orchidee_UserChoices_VEGET_UPDATE} = X ] ; then 
     150        ORCHIDEE_sed VEGET_UPDATE   ${orchidee_UserChoices_VEGET_UPDATE} 
     151    fi 
     152    if [ ! X${orchidee_UserChoices_LAND_COVER_CHANGE} = X ] ; then 
     153        ORCHIDEE_sed LAND_COVER_CHANGE ${orchidee_UserChoices_LAND_COVER_CHANGE} 
    82154    fi 
    83155 
  • CONFIG/LMDZORINCA/tags/LMDZORINCA_v2/EXP_NMHC_AER/PARAM/orchidee.def

    r858 r1565  
    248248STOMATE_HISTLEVEL = 10 
    249249# default = 10 
     250 
     251#-------------------------------------------------------------------- 
     252# STOMATE_IPCC_OUTPUT_FILE 
     253# This file is going to be created by the model 
     254#    and will contain the output from the model. 
     255#    This file is a truly COADS compliant netCDF file. 
     256#    It will be generated by the hist software from 
     257#    the IOIPSL package. 
     258# Name of file in which STOMATE's output is going 
     259# to be written 
     260STOMATE_IPCC_OUTPUT_FILE = stomate_ipcc_history.nc 
     261# default = stomate_ipcc_history.nc 
     262 
     263# STOMATE_IPCC_HIST_DT 
     264# Time step of the STOMATE IPCC history file 
     265# STOMATE IPCC history time step (d) 
     266STOMATE_IPCC_HIST_DT = -1 
     267# default = 0. 
    250268 
    251269# Write frequency for output files (SECHIBA in seconds et 
     
    572590 
    573591# Total depth of soil reservoir 
    574 HYDROL_SOIL_DEPTH = 2. 
     592HYDROL_SOIL_DEPTH = 4. 
    575593# default = 2. 
     594 
     595# Root profile 
     596# Default values were defined for 2 meters soil depth. 
     597# For 4 meters soil depth, you may use those ones : 
     598# 5., .4, .4, 1., .8, .8, 1., 1., .8, 4., 1., 4., 1. 
     599HYDROL_HUMCSTE = 5., .4, .4, 1., .8, .8, 1., 1., .8, 4., 1., 4., 1. 
     600# default =  5., .8, .8, 1., .8, .8, 1., 1., .8, 4., 4., 4., 4. 
    576601 
    577602# Initial restart deep soil moisture if not found in restart 
     
    682707# Add from Nathalie - the 28 of March 2006 - advice from Fred Hourdin 
    683708# By PFT. 
    684 RVEG_PFT = 1., 0.5, 0.5, 1., 1., 1., 1., 1., 1., 1., 0.5, 1., 0.5 
     709RVEG_PFT = 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1., 1. 
    685710# default = 1. 
    686711 
     
    703728# Year of the land_use vegetation map readed 
    704729# year off the pft map 
    705 # If LAND_USE (11 = 1860 - 1850 +1 for PFTmap.20C3M.nc)  
    706 VEGET_YEAR = 151 
     730# If LAND_USE (11 = 1860 - 1850 +1 for PFTmap.20C3M.nc, 1 for PFTmap_IPCC_2000.nc)  
     731VEGET_YEAR = 1 
    707732# default = 282 
     733 
     734# booleen to indicate that a new LAND USE file will be used (since 1.9.5 version). 
     735# The parameter is used to bypass veget_year count  
     736# and reinitialize it with VEGET_YEAR parameter. 
     737# Then it is possible to change LAND USE file. 
     738# If LAND_USE 
     739VEGET_REINIT = n 
     740# default = n 
    708741 
    709742# Update vegetation frequency (since 2.0 version) 
    710743# The veget datas will be update each this time step. 
    711744# If LAND_USE 
    712 VEGET_UPDATE = 0Y 
     745VEGET_UPDATE = XY 
    713746# default = 1Y 
    714747 
     
    718751# deforestation.                                 
    719752# If LAND_USE 
    720 LAND_COVER_CHANGE = n 
     753LAND_COVER_CHANGE = x 
    721754# default = y 
    722755 
  • CONFIG/LMDZORINCA/tags/LMDZORINCA_v2/EXP_NMHC_AER/config.card

    r1557 r1565  
    4545ATM= (lmdz, LMDZ4) 
    4646SRF= (orchidee, orchidee_1_9_5) 
     47SBG= (stomate, ORCHIDEE_1_9_5) 
    4748CHM= (inca, INCA3)  
    4849 
     
    5455ATM= (gcm.e, gcm.e) 
    5556SRF= ("", "") 
     57SBG= ("", "") 
    5658CHM= (inca.dat, inca.dat) 
    5759 
     
    6062[Restarts] 
    6163#D- If you want a GENERAL RULE FOR RESTARTS, put this flag to 'y' 
    62 OverRule=y 
     64OverRule=n 
    6365#D- Last day of the experience used as restart 
    6466RestartDate=2005-12-31 
     
    7981#D- If you want to produce time series, this flag determines 
    8082#D- frequency of post-processing submission 
    81 TimeSeriesFrequency=1Y 
     83TimeSeriesFrequency=NONE 
    8284#D- If you want to produce seasonal average, this flag determines 
    8385#D- the period of this average 
    84 SeasonalFrequency=5Y 
     86SeasonalFrequency=NONE 
    8587 
    8688#======================================================================== 
     
    113115OldName="" 
    114116 
     117#======================================================================== 
     118#D-- SBG - STOMATE 
     119[SBG] 
     120WriteFrequency="1M" 
     121# If config_Restarts_OverRule == 'n' next 4 params are read 
     122Restart=n 
     123#-- Last day of the experience used as restart for this component if Restart=y 
     124RestartDate=2000-01-31 
     125# Define restart simulation name for this component 
     126RestartJobName=EXP00 
     127RestartPath=${ARCHIVE}/IGCM_OUT/LMDZOR/DEVT/amip 
     128#-- Old component name for restart (if empty, use new name) 
     129OldName= 
    115130 
    116131#======================================================================== 
Note: See TracChangeset for help on using the changeset viewer.