Changeset 1353 for CMIP5/esmControl/CpiControl1/COMP
- Timestamp:
- 03/28/11 15:04:16 (13 years ago)
- Location:
- CMIP5/esmControl/CpiControl1/COMP
- Files:
-
- 3 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
Note: See TracChangeset
for help on using the changeset viewer.