Changeset 1353


Ignore:
Timestamp:
03/28/11 15:04:16 (13 years ago)
Author:
cetlod
Message:

Registration of used configuration files for CpiControl1

Location:
CMIP5/esmControl/CpiControl1
Files:
1 added
7 edited

Legend:

Unmodified
Added
Removed

  • CMIP5/esmControl/CpiControl1/COMP/lmdz.card

    r1320 r1353  
    8484Patches= () 
    8585GatherWithInternal = (lon, lat, presnivs, time_counter, time_counter_bnds, aire) 
    86 TimeSeriesVars2D = (bils, cldh, cldl, cldm, cldq, cldt, evap, flat, ftime_con, iwp, lwp, LWdnSFC, LWdnSFCclr, LWupSFCclr, LWdn200, LWdn200clr, LWup200, LWup200clr, precip, prw, pluc, psol, q2m, rh2m, sens, sicf, slp, snow, soll0, soll, sols0, sols, SWdnSFC, SWupSFC, SWup200, SWdnTOA, SWupTOA, SWdn200, SWdnSFCclr, SWupSFCclr, SWupTOAclr, t2m, t2m_max, t2m_min, taux, taux_ter, taux_sic, taux_lic, taux_oce, tauy, tauy_ter, tauy_sic, tauy_lic, tauy_oce, lat_ter, lat_sic, lat_lic, lat_oce, nettop, pbase, ptop, topl0, topl, tops0, tops, topswad, topswai, tsol, tsol_ter, tsol_lic, tsol_oce, tsol_sic , u10m, u500, v10m, v500, w500, pourc_ter, pourc_oce, pourc_sic, pourc_lic, wind10m, od550aer, od550lt1aer, absvisaer, loadoa, loadbc, loadso4, loaddust, loadss, sconcoa, sconcbc, sconcso4, sconcdust, sconcss, reffclwtop, cldncl, swtoaas_ant, swtoacs_ant, swtoacf_ant, swsrfas_ant, swsrfcs_ant, swsrfcf_ant, swtoaas_nat, swtoacs_nat, swtoacf_nat, swsrfas_nat, swsrfcs_nat, swsrfcf_nat, solswai, swtoacf_zero, swsrfcf_zero) 
     86TimeSeriesVars2D = (bils, cldh, cldl, cldm, cldq, cldt, evap, flat, ftime_con, iwp, lwp, LWdnSFC, LWdnSFCclr, LWupSFCclr, LWdn200, LWdn200clr, LWup200, LWup200clr, precip, prw, pluc, psol, q2m, rh2m, sens, sicf, slp, snow, soll0, soll, sols0, sols, SWdnSFC, SWupSFC, SWup200, SWdnTOA, SWupTOA, SWdn200, SWdnSFCclr, SWupSFCclr, SWupTOAclr, t2m, t2m_max, t2m_min, taux, taux_ter, taux_sic, taux_lic, taux_oce, tauy, tauy_ter, tauy_sic, tauy_lic, tauy_oce, lat_ter, lat_sic, lat_lic, lat_oce, nettop, pbase, ptop, topl0, topl, tops0, tops, topswad, topswai, tsol, tsol_ter, tsol_lic, tsol_oce, tsol_sic , u10m, u500, v10m, v500, w500, pourc_ter, pourc_oce, pourc_sic, pourc_lic, wind10m, od550aer, od550lt1aer, absvisaer, loadoa, loadbc, loadso4, loaddust, loadss, sconcoa, sconcbc, sconcso4, sconcdust, sconcss, reffclwtop, cldncl) 
    8787ChunckJob2D = NONE 
    88 TimeSeriesVars3D = (upwd, lwcon, iwcon, lcc3dcon, lcc3dstra, temp, theta, ovap, ovapinit, geop, vitu, vitv, vitw, pres, rneb, rnebcon, rhum, mc, paprs, scdnc, ec550aer, concoa, concbc, concso4, concss, concdust, reffclws, reffclwc, Aga) 
     88TimeSeriesVars3D = (upwd, lwcon, iwcon, lcc3dcon, lcc3dstra, temp, theta, ovap, ovapinit, geop, vitu, vitv, vitw, pres, rneb, rnebcon, rhum, mc, paprs, ec550aer, concoa, concbc, concso4, concss, concdust, reffclws, reffclwc, Aga) 
    8989ChunckJob3D = 100Y 
    9090Seasonal=ON 
    … …  
    120120Patches= () 
    121121GatherWithInternal = (lon, lat, presnivs, time_counter, time_counter_bnds) 
    122 TimeSeriesVars2D = (bils, cldh, cldl, cldm, cldq, cldt, evap, flat, precip, pluc, plul, prw, psol, q2m, sens, sicf, slp, soll, sols, snow, SWdnSFC, t2m, t2m_max, t2m_min, topl, tops, tsol, t700, u10m, v10m, wind10m, wind10max, w500) 
     122TimeSeriesVars2D = (bils, cldh, cldl, cldm, cldq, cldt, evap, flat, precip, pluc, plul, prw, psol, q2m, sens, sicf, slp, soll, sols, SWdnSFC, t2m, t2m_max, t2m_min, topl, tops, tsol, u10m, v10m, wind10m, wind10max) 
    123123ChunckJob2D = 200Y 
    124 TimeSeriesVars3D = (lwcon, iwcon, mc, rhum, rneb, temp, theta, ovap, geop, vitu, vitv, vitw, pres, paprs) 
     124TimeSeriesVars3D = (lwcon, iwcon, mc, rhum, temp, theta, ovap, geop, vitu, vitv, vitw, pres, paprs) 
    125125ChunckJob3D = OFF 
    126126Seasonal=OFF 
    … …  
    183183Patches= () 
    184184GatherWithInternal = (lon, lat, presnivs, time_counter, time_counter_bnds) 
    185 TimeSeriesVars2D = (precip, flat, sens, LWdnSFC, LWupSFC, SWdnSFC, SWupSFC, pluc, snow, LWdnSFCclr, SWdnSFCclr, LWupSFCclr, SWupSFCclr, cldt) 
     185TimeSeriesVars2D = (precip, flat, sens, LWdnSFC, LWupSFC, SWdnSFC, SWupSFC, pluc, snow, LWdnSFCclr, SWdnSFCclr, LWupSFCclr, cldt) 
    186186ChunckJob2D = OFF 
    187187TimeSeriesVars3D = () 

  • CMIP5/esmControl/CpiControl1/COMP/lmdz.driver

    r1316 r1353  
    139139 
    140140     ##-- GHG forcing for the current year 
    141  
    142     if ( ${FirstInitialize} ) ; then 
    143         ##--Initialization of fluxes to an undefined value at the first run 
    144         UndefinedValueCo2=-9999 
    145         InitPeriodCo2=0 
     141    if ( ${FirstInitialize} && [ X"${config_ATM_Restart}" = X"n" ]) ; then 
     142        UndefinedValueCo2=-9999 
     143        InitPeriodCo2=$(( ${CumulPeriod} - 1 )) 
    146144        InitDateBeginCo2=${UndefinedValueCo2} 
    147         InitYearCo2=$(( ${year} - 1 ))  
     145        InitYearCo2=$(( ${year} - 1 )) 
    148146        InitDateEndCo2=${InitYearCo2}1231 
    149         fCO2_ff=${UndefinedValueCo2}  
     147        fCO2_ff=${UndefinedValueCo2} 
    150148        CO2SBG=${UndefinedValueCo2} 
    151149        CO2LU=${UndefinedValueCo2} 
    152150        CO2MBG=${UndefinedValueCo2} 
    153151        CO2_ppm_prec=${UndefinedValueCo2} 
    154         # Get the initial value of atmosp. pco2  
    155152        CO2_ppm=${lmdz_UserChoices_co2_init} 
    156153 
    157154        IGCM_sys_Cp ${SUBMIT_DIR}/co2.log.init ${SUBMIT_DIR}/co2.log 
    158155        IGCM_sys_Chmod u+w ${SUBMIT_DIR}/co2.log 
    159      
    160         # Save CO2 values in ExeCpuLog variable contents 5 fields 
     156 
    161157        echo "${InitPeriodCo2} ${InitDateBeginCo2} ${InitDateEndCo2} ${CO2_ppm} ${CO2_ppm_prec} ${fCO2_ff} ${CO2MBG} ${CO2SBG} ${CO2LU}" |   \ 
    162158             gawk '{printf("  %11d | %15s | %15s | %22.15g | %20.15g | %11.8g | %11.8g | %12.8g | %12.8g \n", \ 
    163159                              $1,$2,$3,$4,$5,$6,$7,$8,$9)}' >> ${SUBMIT_DIR}/co2.log 
    164      
    165         # Get the path where the log file co2.log is 
    166         PathCO2log=${SUBMIT_DIR}  
    167         # Get Previous line in co2.log 
     160 
     161        PathCO2log=${SUBMIT_DIR} 
    168162        PreviousLinelog=$( tail -1 ${PathCO2log}/co2.log ) 
    169163 
    170164    else 
     165 
    171166        # Get restart line in co2.log 
    172         PathCO2log=${SUBMIT_DIR} 
    173         PreviousLinelog=$( tail -1 ${PathCO2log}/co2.log )  
    174         # Get the value of atmosp. pco2 in co2.log 
    175         CO2_ppm=$( echo "${PreviousLinelog[*]}" | gawk -F '|' '{gsub(" ",""); print $4}' ) 
     167 
     168        # NEW parameter in config.card (even for Overule all restarts, you must modify the line after this comment ): 
     169        # !! NEW for Carbone cycle !! Path where to find old co2.log file for restart CO2. 
     170        #SubmitRestartPath= 
     171        if ( ${FirstInitialize} && [ X"${config_ATM_Restart}" = X"y" ]) ; then 
     172            IGCM_card_DefineVariableFromOption ${SUBMIT_DIR}/config.card ATM SubmitRestartPath 
     173            PathCO2log=${config_ATM_SubmitRestartPath} 
     174            Date_r=$( IGCM_date_ConvertFormatToGregorian ${config_ATM_RestartDate} ) 
     175            PreviousLinelog=$( grep ${Date_r} ${PathCO2log}/co2.log ) 
     176 
     177            IGCM_sys_Cp ${SUBMIT_DIR}/co2.log.init ${SUBMIT_DIR}/co2.log 
     178            IGCM_sys_Chmod u+w ${SUBMIT_DIR}/co2.log 
     179 
     180            # Save Last Line of control in new co2.log file 
     181            InitPeriodCo2=0 
     182            InitDateBeginCo2=$( echo "${PreviousLinelog[*]}" | gawk -F '|' '{gsub(" ",""); print $2}' ) 
     183            InitDateEndCo2=$( echo "${PreviousLinelog[*]}" | gawk -F '|' '{gsub(" ",""); print $3}' ) 
     184            fCO2_ff=$( echo "${PreviousLinelog[*]}" | gawk -F '|' '{gsub(" ",""); print $6}' ) 
     185            CO2SBG=$( echo "${PreviousLinelog[*]}" | gawk -F '|' '{gsub(" ",""); print $8}' ) 
     186            CO2LU=$( echo "${PreviousLinelog[*]}" | gawk -F '|' '{gsub(" ",""); print $9}' ) 
     187            CO2MBG=$( echo "${PreviousLinelog[*]}" | gawk -F '|' '{gsub(" ",""); print $7}' ) 
     188            CO2_ppm_prec=$( echo "${PreviousLinelog[*]}" | gawk -F '|' '{gsub(" ",""); print $5}' ) 
     189            CO2_ppm=$( echo "${PreviousLinelog[*]}" | gawk -F '|' '{gsub(" ",""); print $4}' ) 
     190 
     191            echo "${InitPeriodCo2} ${InitDateBeginCo2} ${InitDateEndCo2} ${CO2_ppm} ${CO2_ppm_prec} ${fCO2_ff} ${CO2MBG} ${CO2SBG} ${CO2LU}" |   \ 
     192                  gawk '{printf("  %11d | %15s | %15s | %22.15g | %20.15g | %11.8g | %11.8g | %12.8g | %12.8g \n", \ 
     193                              $1,$2,$3,$4,$5,$6,$7,$8,$9)}' >> ${SUBMIT_DIR}/co2.log 
     194 
     195        else 
     196            PathCO2log=${SUBMIT_DIR} 
     197            PreviousLinelog=$( tail -1 ${PathCO2log}/co2.log ) 
     198            CO2_ppm=$( echo "${PreviousLinelog[*]}" | gawk -F '|' '{gsub(" ",""); print $4}' ) 
     199        fi 
    176200 
    177201    fi 
    … …  
    190214    LMDZ_sed physiq.def ok_histNMC  "${LMDZ_NMC_monthly}, ${LMDZ_NMC_daily}, ${LMDZ_NMC_hf}" 
    191215#   For Control, we won't modify ATM CO2 : 
    192 #    LMDZ_sed physiq.def co2_ppm     ${CO2_ppm} 
     216    LMDZ_sed physiq.def co2_ppm     ${CO2_ppm} 
    193217 
    194218    ## gcm.def parameters : no change since gcm.def_${RESOL_ATM_3D} is used and already modified 
    … …  
    263287    [ ${CumulPeriod} -le ${LMDZ_NbPeriod_adjust} ] && IGCM_sys_Put_Out Bands_${RESOL_ATM_3D}_${NUM_PROC_ATM}prc.dat ${R_OUT_ATM_D}/${config_UserChoices_JobName}_Bands_${RESOL_ATM_3D}_${NUM_PROC_ATM}prc.dat_${CumulPeriod} 
    264288 
     289#   get the fossile fuel file : 
     290 
    265291    typeset CO2SBG_3V CO2SBG CO2MBG CO2SRF CO2LU CO2_ppm CO2_ppm_prec fCO2_ff 
    266292 
    267     # For Control, no fossil fuel emission : 
    268293    fCO2_ff=0 
    269     echo "Fossil Fuel for control run  :" ${fCO2_ff} 
    270  
    271     # Get the value of ocean carbon flux 
     294 
    272295    CO2MBG=$( ${SUBMIT_DIR}/COMP/lmdz_analyse_pisces_out.awk ocean.output ) 
    273296    CO2MBG=$( echo ${CO2MBG} | sed -e "s/ *//" ) 
    274297    echo "Ocean carbon flux  :" ${CO2MBG} 
    275298 
    276     # Get the value of land fluxes 
    277299    set +A CO2SBG_3V -- $( ${SUBMIT_DIR}/COMP/lmdz_analyse_stomate_out.awk out_orchidee_0000 ) 
    278300    CO2SRF=$( echo ${CO2SBG_3V[0]} | sed -e "s/ *//" ) 
    279301    CO2LU=$( echo ${CO2SBG_3V[1]} | sed -e "s/ *//" ) 
    280302    CO2SBG=$( echo ${CO2SBG_3V[2]} | sed -e "s/ *//" ) 
    281     echo "Land carbon flux (NEP)        :" ${CO2SRF} 
    282     echo "LU flux          (FLUC)        :" ${CO2LU} 
    283     echo "Total Land carbon flux (NBP)  :" ${CO2SBG} 
     303    echo "Land carbon flux :" ${CO2SRF} 
     304    echo "LU flux          :" ${CO2LU} 
     305    echo "Total Land carbon flux :" ${CO2SBG} 
    284306 
    285307    # Previous CO2 
    286308    CO2_ppm_prec=$( echo "${PreviousLinelog[*]}" | gawk -F '|' '{gsub(" ",""); print $4}' ) 
    287309    echo "Previous CO2  :" ${CO2_ppm_prec} 
    288      
     310 
    289311    # Formula 
    290312    CO2_ppm=$( echo "${CO2_ppm_prec} + (${fCO2_ff} + ${CO2MBG} + ${CO2SBG}) / 2.12" | bc -l ) 
    … …  
    295317         gawk '{printf("  %11d | %15s | %15s | %22.15g | %20.15g | %11.8g | %11.8g | %12.8g | %12.8g \n", \ 
    296318                           $1,$2,$3,$4,$5,$6,$7,$8,$9)}' >> ${SUBMIT_DIR}/co2.log 
    297  
    298319    echo FINALIZE ATM ! 
    299320 

  • CMIP5/esmControl/CpiControl1/COMP/pisces.driver

    r1316 r1353  
    7878        namelist_pisces > namelist_pisces.tmp 
    7979 
    80     IGCM_sys_Mv namelist_pisces.tmp namelist_pisces 
     80    IGCM_sys_Mv namelist_pisces.tmp  namelist_pisces 
    8181 
    8282 

  • CMIP5/esmControl/CpiControl1/PARAM/iodef.xml

    r1248 r1353  
    349349 
    350350         <!-- table 2.2 --> 
    351           <field ref="uice_ipa"     name="iicevelu"   description="Ice velocity along i-axis at I-point (ice presence average)"  level="3"    /> 
    352           <field ref="vice_ipa"     name="iicevelv"   description="Ice velocity along j-axis at I-point (ice presence average)"  level="3"    /> 
    353           <field ref="ice_cover"    name="ileadfra"   description="sea_ice_area_fraction"                                        level="3"    /> 
    354           <field ref="icethic_cea"  name="iicethic"   description="sea_ice_thickness"                                            level="3"    /> 
     351          <field ref="uice_ipa"     name="iicevelu"   description="Ice velocity along i-axis at I-point (ice presence average)"  level="2"    /> 
     352          <field ref="vice_ipa"     name="iicevelv"   description="Ice velocity along j-axis at I-point (ice presence average)"  level="2"    /> 
     353          <field ref="ice_cover"    name="ileadfra"   description="sea_ice_area_fraction"                                        level="2"    /> 
     354          <field ref="icethic_cea"  name="iicethic"   description="sea_ice_thickness"                                            level="2"    /> 
    355355        </file>    
    356356 
    … …  
    435435        <field ref="ice_cover"         name="soicecov"    description="sea_ice_area_fraction"                                    /> 
    436436        <field ref="BLT"               name="blt"         description="Barrier Layer Thickness"                    level="3"     /> 
    437         <field ref="mld|dt|"           name="mld_dt02"    description="Mixed Layer Depth |dt|=0.2 (ref.10m)"       level="3"     /> 
     437        <field ref="mld|dt|"           name="mld_dt02"    description="Mixed Layer Depth |dt|=0.2 (ref.10m)"         level="3"     /> 
    438438 
    439439          </file>    

  • CMIP5/esmControl/CpiControl1/PARAM/physiq.def_L39

    r1283 r1353  
    4949flag_pres         =  2, 3, 7, 6, 6 
    5050name_pres         =  pres 
    51 # w500 for histday 
    52 flag_w500         =  1, 5, 3, 10, 10 
    53 name_w500         =  w500 
    54 # t700 for histday 
    55 flag_t700         =  1, 5, 3, 10, 10 
    56 name_t700         =  t700 
    5751# Parametres niveau de sorties differents fichiers  
    5852# 
    … …  
    6963# - lev_histday=4 => + champs sous-surfaces 
    7064# - lev_histday=5 => + champs F. Aires 
    71 lev_histday=2 
     65lev_histday=5 
    7266### lev_histmth=0-4, niveau de sortie fichier "histmth.nc" 
    7367# - lev_histmth=0 => pas de sorties lev_histmth.nc 

  • CMIP5/esmControl/CpiControl1/README

    r1055 r1353  
    44 
    55Experiment Name comes from CMIP5 DRS. See: http://cmip-pcmdi.llnl.gov/cmip5/docs/cmip5_data_reference_syntax.pdf 
    6 Output files are stored into : login/IGCM_OUT/IPSLCM5A/DEVT/pdControl/EXP00 directory 
     6Output files are stored into : login/IGCM_OUT/IPSLCM5A_C/PROD/CpiControl directory 
     7 
     8plot_co2.jnl is a ferret script to plot co2.log file - remove the | before running the script - 
    79 
    810 

  • CMIP5/esmControl/CpiControl1/config.card

    r1316 r1353  
    99[UserChoices] 
    1010#============================ 
    11 JobName=CpiControl00 
     11JobName=CpiControl01 
    1212#----- Short Name of Experiment 
    1313ExperimentName=CpiControl 
    1414#----- DEVT TEST PROD 
    15 SpaceName=DEVT 
    16 LongName="IPSLCM5A CMIP5 DEVT phase CpiControl example with large outputs" 
     15SpaceName=PROD 
     16LongName="IPSLCM5A_C CMIP5 PROD phase CpiControl : Fossil fuel equal 0  and increasing interactive atmospheric pCO2" 
    1717TagName=IPSLCM5A_C 
    1818#============================ 
    … …  
    2121#-- Experiment dates : Beginning and ending 
    2222#-- "YYYY-MM-DD" 
    23 DateBegin=1850-01-01 
    24 DateEnd=1899-12-31 
     23DateBegin=1800-01-01 
     24DateEnd=2099-12-31 
    2525#============================ 
    2626#-- 1Y, 1M, 5D, 1D Period Length of one trunk of simulation 
    … …  
    5454#D- If you want to restart all components from the same simulation, 
    5555#D-      put OveRule flag to 'y' and set the next 3 parameters 
    56 OverRule=n 
     56OverRule=y 
    5757#D- Last day of the experience used as restart for all components 
    58 RestartDate=1899-12-31 
     58RestartDate=1850-12-31 
    5959#D- Define restart simulation name for all components 
    60 RestartJobName=EXP00 
     60RestartJobName=CpiControl00 
    6161#D- Path Server Group Login 
    62 RestartPath=${ARCHIVE}/IGCM_OUT/IPSLCM5A/DEVT/piControl 
     62RestartPath=${ARCHIVE}/IGCM_OUT/IPSLCM5A_C/PROD/CpiControl 
    6363 
    6464#======================================================================== 
    … …  
    7474RestartJobName=EXP00 
    7575RestartPath=${ARCHIVE}/IGCM_OUT/IPSLCM5A/DEVT/piControl 
     76# !! NEW for Carbone cycle !! Path where to find old co2.log file for restart CO2. 
     77SubmitRestartPath=${SUBMIT_DIR}/../CpiControl0 
    7678# Old component name for restart (if empty, use new name) 
    7779OldName= 
    … …  
    193195RebuildFrequency=5Y 
    194196#D- Do we rebuild parallel output from archive (use NONE to use SCRATCHDIR as buffer) 
    195 RebuildFromArchive=true 
     197RebuildFromArchive=none 
    196198#D- If you want to produce time series, this flag determines 
    197199#D- frequency of post-processing submission (NONE if you don't want) 
    198 TimeSeriesFrequency=10Y 
     200TimeSeriesFrequency=5Y 
    199201#D- If you want to produce seasonal average, this flag determines 
    200202#D- the period of this average (NONE if you don't want) 
    201 SeasonalFrequency=10Y 
     203SeasonalFrequency=5Y 
    202204#D- Offset for seasonal average first start dates ; same unit as SeasonalFrequency 
    203205#D- Usefull if you do not want to consider the first X simulation's years 
Note: See TracChangeset for help on using the changeset viewer.