Changeset 1233 for CONFIG/IPSLCM

Timestamp:

12/30/10 10:07:18 (14 years ago)

Author:

mmaipsl

Message:

First CO2 interactive Control scripts.
Please, see ​https://forge.ipsl.jussieu.fr/igcmg/wiki/IPSLCM5ACo2 for informations.
And follow links in ​https://forge.ipsl.jussieu.fr/igcmg/wiki/IPSLCM5ACo2#Testsetruns
for example of use.

Location:

CONFIG/IPSLCM/IPSLCM5A/CpiControl

Files:

• COMP/lmdz.card (modified) (1 diff)
• COMP/lmdz.driver (modified) (3 diffs)
• config.card (modified) (1 diff)

CONFIG/IPSLCM/IPSLCM5A/CpiControl/COMP/lmdz.card

@@ -27 +27 @@
 LMDZ_NMC_daily=y
 LMDZ_NMC_hf=y
+
+# CO2_ppm init
+co2_init=284.725
 
 [InitialStateFiles]

CONFIG/IPSLCM/IPSLCM5A/CpiControl/COMP/lmdz.driver

@@ -88 +88 @@
 {
     IGCM_debug_PushStack "ATM_Update"
+
+    typeset CO2SBG CO2MBG CO2LU CO2_ppm CO2_ppm_prec fCO2_ff
+    typeset PreviousLinelog LastPeriodDateBegin LastPeriodDateEnd LastDatesPeriod LastPREFIX
 
     case ${config_UserChoices_PeriodLength} in
@@ -135 +138 @@
     [ ${CumulPeriod} -gt 1 ] && ( IGCM_sys_Get ${R_OUT_ATM_D}/${config_UserChoices_JobName}_Bands_${RESOL_ATM_3D}_${NUM_PROC_ATM}prc.dat_${RefPeriod} Bands_${RESOL_ATM_3D}_${NUM_PROC_ATM}prc.dat ; IGCM_sys_Chmod u+w Bands_${RESOL_ATM_3D}_${NUM_PROC_ATM}prc.dat )
 
+     ##-- GHG forcing for the current year
+
+    if ( ${FirstInitialize} ) ; then
+        
+        fCO2_ff=0
+        CO2SRF=0
+        CO2LU=0
+        CO2MBG=0
+        CO2_ppm=${lmdz_UserChoices_co2_init}
+        CO2_ppm_prec=${lmdz_UserChoices_co2_init}
+        
+        IGCM_sys_Cp ${SUBMIT_DIR}/co2.log.init ${SUBMIT_DIR}/co2.log
+        IGCM_sys_Chmod u+w ${SUBMIT_DIR}/co2.log
+
+        echo "${CumulPeriod} ${PeriodDateBegin} ${PeriodDateEnd} ${CO2_ppm} ${CO2_ppm_prec} ${fCO2_ff} ${CO2MBG} ${CO2SRF} ${CO2LU}" |   \
+             gawk '{printf("  %11d | %15s | %15s | %22.15g | %20.15g | %11.8g | %11.8g | %12.8g | %12.8g \n", \
+                              $1,$2,$3,$4,$5,$6,$7,$8,$9)}' >> ${SUBMIT_DIR}/co2.log
+
+    else
+
+#       get the fossile fuel file :
+        fCO2_ff=$( gawk "{if (match(\$0,\"${year}\ *${month}\")) {print \$5}}" ${SUBMIT_DIR}/PARAM/CMIP5_gridcar_CO2_emissions_fossil_fuel_Andres_1751-2007_monthly_SC_grid1x1.txt )
+        fCO2_ff=$( echo "${fCO2_ff} / 1000." | bc -l  )
+#       For Control, no fossil fuel emission :
+        fCO2_ff=0
+        echo "Fossil Fuel for control run  :" ${fCO2_ff}
+
+        # Get Previous line in co2.log
+        PreviousLinelog=$( tail -1 ${SUBMIT_DIR}/co2.log )
+
+        # Get Previous period output PREFIX
+        LastPeriodDateBegin=$( echo "${PreviousLinelog[*]}" | gawk -F '|' '{gsub(" ",""); print $2}' )
+        LastPeriodDateEnd=$( echo "${PreviousLinelog[*]}" | gawk -F '|' '{gsub(" ",""); print $3}' )
+        eval LastDatesPeriod=${LastPeriodDateBegin}_${LastPeriodDateEnd} > /dev/null 2>&1
+        # Prefix for save files of this period
+        eval LastPREFIX=${config_UserChoices_JobName}_${LastDatesPeriod}  > /dev/null 2>&1
+
+        # Ocean
+        
+        IGCM_sys_TestFileArchive ${R_OUT_OCE}/Debug/${LastPREFIX}_ocean.output
+        if [ $? = 0 ] ; then
+
+            if [ $DRYRUN -gt 1 ]; then
+                IGCM_sys_Cp ${R_OUT_OCE}/Debug/${LastPREFIX}_ocean.output ./
+            else
+                IGCM_sys_Get ${R_OUT_OCE}/Debug/${LastPREFIX}_ocean.output ./
+            fi
+            IGCM_comp_PrepareDeletedFiles ./${LastPREFIX}_ocean.output
+            CO2MBG=$( ${SUBMIT_DIR}/COMP/lmdz_analyse_pisces_out.awk ${LastPREFIX}_ocean.output )
+            CO2MBG=$( echo ${CO2MBG} | sed -e "s/ *//" )
+        else
+            echo "Error : file ${R_OUT_OCE}/Debug/${LastPREFIX}_ocean.output not found."
+            exit 1
+        fi
+
+        # Surface
+        IGCM_sys_TestFileArchive ${R_OUT_SRF}/Debug/${LastPREFIX}_out_orchidee_0000
+        if [ $? = 0 ] ; then
+            if [ $DRYRUN -gt 1 ]; then
+                IGCM_sys_Cp ${R_OUT_SRF}/Debug/${LastPREFIX}_out_orchidee_0000 ./
+            else
+                IGCM_sys_Get ${R_OUT_SRF}/Debug/${LastPREFIX}_out_orchidee_0000 ./
+            fi
+            IGCM_comp_PrepareDeletedFiles ./${LastPREFIX}_out_orchidee_0000
+            set +A CO2SBG -- $( ${SUBMIT_DIR}/COMP/lmdz_analyse_stomate_out.awk ${LastPREFIX}_out_orchidee_0000 )
+        else
+            echo "Error : file ${R_OUT_OCE}/Debug/${LastPREFIX}_ocean.output not found."
+            exit 1
+        fi
+        CO2SRF=$( echo ${CO2SBG[0]} | sed -e "s/ *//" )
+        CO2LU=$( echo ${CO2SBG[1]} | sed -e "s/ *//" )
+
+        # Previous CO2
+        CO2_ppm_prec=$( echo "${PreviousLinelog[*]}" | gawk -F '|' '{gsub(" ",""); print $4}' )
+
+        # Formula
+        CO2_ppm=$( echo "${CO2_ppm_prec} + (${fCO2_ff} + ${CO2MBG} + ${CO2SRF} + ${CO2LU}) / 2.12" | bc -l )
+
+        # Save CO2 values in ExeCpuLog variable contents 5 fields
+        echo "${CumulPeriod} ${PeriodDateBegin} ${PeriodDateEnd} ${CO2_ppm} ${CO2_ppm_prec} ${fCO2_ff} ${CO2MBG} ${CO2SRF} ${CO2LU}" |   \
+             gawk '{printf("  %11d | %15s | %15s | %22.15g | %20.15g | %11.8g | %11.8g | %12.8g | %12.8g \n", \
+                              $1,$2,$3,$4,$5,$6,$7,$8,$9)}' >> ${SUBMIT_DIR}/co2.log
+
+    fi
+
     ## physiq.def parameters modified from initial physiq_L${RESOL_ATM_Z}
     LMDZ_sed physiq.def OK_instan   ${OK_instan}
@@ -147 +235 @@
     LMDZ_sed physiq.def ok_hfCOSP   ${LMDZ_COSP_hf}
     LMDZ_sed physiq.def ok_histNMC  "${LMDZ_NMC_monthly}, ${LMDZ_NMC_daily}, ${LMDZ_NMC_hf}"
+#   For Control, we won't modify ATM CO2 :
+#    LMDZ_sed physiq.def co2_ppm     ${CO2_ppm}
 
     ## gcm.def parameters : no change since gcm.def_${RESOL_ATM_3D} is used and already modified

CONFIG/IPSLCM/IPSLCM5A/CpiControl/config.card

@@ -9 +9 @@
 [UserChoices]
 #============================
-JobName=EXP00
+JobName=CpiControl00
 #----- Short Name of Experiment
-ExperimentName=piControl
+ExperimentName=CpiControl
 #----- DEVT TEST PROD
 SpaceName=DEVT
-LongName="IPSLCM5A CMIP5 DEVT phase piControl example with large outputs"
-TagName=IPSLCM5A
+LongName="IPSLCM5A CMIP5 DEVT phase CpiControl example with large outputs"
+TagName=IPSLCM5A_C
 #============================
 #-- leap, noleap, 360d