Changeset 1353
- Timestamp:
- 03/28/11 15:04:16 (14 years ago)
- Location:
- CMIP5/esmControl/CpiControl1
- Files:
-
- 1 added
- 7 edited
Legend:
- Unmodified
- Added
- Removed
-
CMIP5/esmControl/CpiControl1/COMP/lmdz.card
r1320 r1353 84 84 Patches= () 85 85 GatherWithInternal = (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)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) 87 87 ChunckJob2D = 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)88 TimeSeriesVars3D = (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) 89 89 ChunckJob3D = 100Y 90 90 Seasonal=ON … … 120 120 Patches= () 121 121 GatherWithInternal = (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)122 TimeSeriesVars2D = (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) 123 123 ChunckJob2D = 200Y 124 TimeSeriesVars3D = (lwcon, iwcon, mc, rhum, rneb,temp, theta, ovap, geop, vitu, vitv, vitw, pres, paprs)124 TimeSeriesVars3D = (lwcon, iwcon, mc, rhum, temp, theta, ovap, geop, vitu, vitv, vitw, pres, paprs) 125 125 ChunckJob3D = OFF 126 126 Seasonal=OFF … … 183 183 Patches= () 184 184 GatherWithInternal = (lon, lat, presnivs, time_counter, time_counter_bnds) 185 TimeSeriesVars2D = (precip, flat, sens, LWdnSFC, LWupSFC, SWdnSFC, SWupSFC, pluc, snow, LWdnSFCclr, SWdnSFCclr, LWupSFCclr, SWupSFCclr,cldt)185 TimeSeriesVars2D = (precip, flat, sens, LWdnSFC, LWupSFC, SWdnSFC, SWupSFC, pluc, snow, LWdnSFCclr, SWdnSFCclr, LWupSFCclr, cldt) 186 186 ChunckJob2D = OFF 187 187 TimeSeriesVars3D = () -
CMIP5/esmControl/CpiControl1/COMP/lmdz.driver
r1316 r1353 139 139 140 140 ##-- 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 )) 146 144 InitDateBeginCo2=${UndefinedValueCo2} 147 InitYearCo2=$(( ${year} - 1 )) 145 InitYearCo2=$(( ${year} - 1 )) 148 146 InitDateEndCo2=${InitYearCo2}1231 149 fCO2_ff=${UndefinedValueCo2} 147 fCO2_ff=${UndefinedValueCo2} 150 148 CO2SBG=${UndefinedValueCo2} 151 149 CO2LU=${UndefinedValueCo2} 152 150 CO2MBG=${UndefinedValueCo2} 153 151 CO2_ppm_prec=${UndefinedValueCo2} 154 # Get the initial value of atmosp. pco2155 152 CO2_ppm=${lmdz_UserChoices_co2_init} 156 153 157 154 IGCM_sys_Cp ${SUBMIT_DIR}/co2.log.init ${SUBMIT_DIR}/co2.log 158 155 IGCM_sys_Chmod u+w ${SUBMIT_DIR}/co2.log 159 160 # Save CO2 values in ExeCpuLog variable contents 5 fields 156 161 157 echo "${InitPeriodCo2} ${InitDateBeginCo2} ${InitDateEndCo2} ${CO2_ppm} ${CO2_ppm_prec} ${fCO2_ff} ${CO2MBG} ${CO2SBG} ${CO2LU}" | \ 162 158 gawk '{printf(" %11d | %15s | %15s | %22.15g | %20.15g | %11.8g | %11.8g | %12.8g | %12.8g \n", \ 163 159 $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} 168 162 PreviousLinelog=$( tail -1 ${PathCO2log}/co2.log ) 169 163 170 164 else 165 171 166 # 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 176 200 177 201 fi … … 190 214 LMDZ_sed physiq.def ok_histNMC "${LMDZ_NMC_monthly}, ${LMDZ_NMC_daily}, ${LMDZ_NMC_hf}" 191 215 # 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} 193 217 194 218 ## gcm.def parameters : no change since gcm.def_${RESOL_ATM_3D} is used and already modified … … 263 287 [ ${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} 264 288 289 # get the fossile fuel file : 290 265 291 typeset CO2SBG_3V CO2SBG CO2MBG CO2SRF CO2LU CO2_ppm CO2_ppm_prec fCO2_ff 266 292 267 # For Control, no fossil fuel emission :268 293 fCO2_ff=0 269 echo "Fossil Fuel for control run :" ${fCO2_ff} 270 271 # Get the value of ocean carbon flux 294 272 295 CO2MBG=$( ${SUBMIT_DIR}/COMP/lmdz_analyse_pisces_out.awk ocean.output ) 273 296 CO2MBG=$( echo ${CO2MBG} | sed -e "s/ *//" ) 274 297 echo "Ocean carbon flux :" ${CO2MBG} 275 298 276 # Get the value of land fluxes277 299 set +A CO2SBG_3V -- $( ${SUBMIT_DIR}/COMP/lmdz_analyse_stomate_out.awk out_orchidee_0000 ) 278 300 CO2SRF=$( echo ${CO2SBG_3V[0]} | sed -e "s/ *//" ) 279 301 CO2LU=$( echo ${CO2SBG_3V[1]} | sed -e "s/ *//" ) 280 302 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} 284 306 285 307 # Previous CO2 286 308 CO2_ppm_prec=$( echo "${PreviousLinelog[*]}" | gawk -F '|' '{gsub(" ",""); print $4}' ) 287 309 echo "Previous CO2 :" ${CO2_ppm_prec} 288 310 289 311 # Formula 290 312 CO2_ppm=$( echo "${CO2_ppm_prec} + (${fCO2_ff} + ${CO2MBG} + ${CO2SBG}) / 2.12" | bc -l ) … … 295 317 gawk '{printf(" %11d | %15s | %15s | %22.15g | %20.15g | %11.8g | %11.8g | %12.8g | %12.8g \n", \ 296 318 $1,$2,$3,$4,$5,$6,$7,$8,$9)}' >> ${SUBMIT_DIR}/co2.log 297 298 319 echo FINALIZE ATM ! 299 320 -
CMIP5/esmControl/CpiControl1/COMP/pisces.driver
r1316 r1353 78 78 namelist_pisces > namelist_pisces.tmp 79 79 80 IGCM_sys_Mv namelist_pisces.tmp namelist_pisces80 IGCM_sys_Mv namelist_pisces.tmp namelist_pisces 81 81 82 82 -
CMIP5/esmControl/CpiControl1/PARAM/iodef.xml
r1248 r1353 349 349 350 350 <!-- 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" /> 355 355 </file> 356 356 … … 435 435 <field ref="ice_cover" name="soicecov" description="sea_ice_area_fraction" /> 436 436 <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" /> 438 438 439 439 </file> -
CMIP5/esmControl/CpiControl1/PARAM/physiq.def_L39
r1283 r1353 49 49 flag_pres = 2, 3, 7, 6, 6 50 50 name_pres = pres 51 # w500 for histday52 flag_w500 = 1, 5, 3, 10, 1053 name_w500 = w50054 # t700 for histday55 flag_t700 = 1, 5, 3, 10, 1056 name_t700 = t70057 51 # Parametres niveau de sorties differents fichiers 58 52 # … … 69 63 # - lev_histday=4 => + champs sous-surfaces 70 64 # - lev_histday=5 => + champs F. Aires 71 lev_histday= 265 lev_histday=5 72 66 ### lev_histmth=0-4, niveau de sortie fichier "histmth.nc" 73 67 # - lev_histmth=0 => pas de sorties lev_histmth.nc -
CMIP5/esmControl/CpiControl1/README
r1055 r1353 4 4 5 5 Experiment 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 6 Output files are stored into : login/IGCM_OUT/IPSLCM5A_C/PROD/CpiControl directory 7 8 plot_co2.jnl is a ferret script to plot co2.log file - remove the | before running the script - 7 9 8 10 -
CMIP5/esmControl/CpiControl1/config.card
r1316 r1353 9 9 [UserChoices] 10 10 #============================ 11 JobName=CpiControl0 011 JobName=CpiControl01 12 12 #----- Short Name of Experiment 13 13 ExperimentName=CpiControl 14 14 #----- DEVT TEST PROD 15 SpaceName= DEVT16 LongName="IPSLCM5A CMIP5 DEVT phase CpiControl example with large outputs"15 SpaceName=PROD 16 LongName="IPSLCM5A_C CMIP5 PROD phase CpiControl : Fossil fuel equal 0 and increasing interactive atmospheric pCO2" 17 17 TagName=IPSLCM5A_C 18 18 #============================ … … 21 21 #-- Experiment dates : Beginning and ending 22 22 #-- "YYYY-MM-DD" 23 DateBegin=18 50-01-0124 DateEnd= 1899-12-3123 DateBegin=1800-01-01 24 DateEnd=2099-12-31 25 25 #============================ 26 26 #-- 1Y, 1M, 5D, 1D Period Length of one trunk of simulation … … 54 54 #D- If you want to restart all components from the same simulation, 55 55 #D- put OveRule flag to 'y' and set the next 3 parameters 56 OverRule= n56 OverRule=y 57 57 #D- Last day of the experience used as restart for all components 58 RestartDate=18 99-12-3158 RestartDate=1850-12-31 59 59 #D- Define restart simulation name for all components 60 RestartJobName= EXP0060 RestartJobName=CpiControl00 61 61 #D- Path Server Group Login 62 RestartPath=${ARCHIVE}/IGCM_OUT/IPSLCM5A /DEVT/piControl62 RestartPath=${ARCHIVE}/IGCM_OUT/IPSLCM5A_C/PROD/CpiControl 63 63 64 64 #======================================================================== … … 74 74 RestartJobName=EXP00 75 75 RestartPath=${ARCHIVE}/IGCM_OUT/IPSLCM5A/DEVT/piControl 76 # !! NEW for Carbone cycle !! Path where to find old co2.log file for restart CO2. 77 SubmitRestartPath=${SUBMIT_DIR}/../CpiControl0 76 78 # Old component name for restart (if empty, use new name) 77 79 OldName= … … 193 195 RebuildFrequency=5Y 194 196 #D- Do we rebuild parallel output from archive (use NONE to use SCRATCHDIR as buffer) 195 RebuildFromArchive= true197 RebuildFromArchive=none 196 198 #D- If you want to produce time series, this flag determines 197 199 #D- frequency of post-processing submission (NONE if you don't want) 198 TimeSeriesFrequency= 10Y200 TimeSeriesFrequency=5Y 199 201 #D- If you want to produce seasonal average, this flag determines 200 202 #D- the period of this average (NONE if you don't want) 201 SeasonalFrequency= 10Y203 SeasonalFrequency=5Y 202 204 #D- Offset for seasonal average first start dates ; same unit as SeasonalFrequency 203 205 #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.