Changeset 1535 for CONFIG/IPSLCM/IPSLCM5-2/piControl

Timestamp:

08/18/11 17:41:19 (13 years ago)

Author:

jgipsl

Message:

Modifications in working configuration :

Modifications in lmdz.card and lmdz.driver

• GHG forcing files should now be copied in files with standard names such as CO2.txt, CH4.txt, N2O.txt, CFC11.txt and CFC12.txt. Now the

lmdz.driver will update the GHG parameters in physiq.def as soon as one
of theses files are added in lmdz.card.

• removed CO2_1765_2005.txt from Chistorical/COMP/lmdz.card, the file is not used in this experience

Modifications in pisces.driver and pisces.card

• add option CARBON_CYCLE=y to be set in pisces.card section UserChoices? for experience Chistorical and CpiControl?. With this option, co2 will

be read from co2.log file as before.

• modified pisces.driver to read co2 from file CO2.txt if existing. As done in lmdz.driver
• all pisces.driver are now the same

Location:

CONFIG/IPSLCM/IPSLCM5-2/piControl/COMP

Files:

• lmdz.driver (modified) (1 diff)
• pisces.driver (modified) (3 diffs)

CONFIG/IPSLCM/IPSLCM5-2/piControl/COMP/lmdz.driver

@@ -134 +134 @@
 
     [ ${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 :
+    ##   If forcing file exist in run directory, read values for the current year and modify physiq.def
+    if [ -f SOLARANDVOLCANOES.txt ] ; then
+        IPCC_SOLAR=`grep Annee_${year} SOLARANDVOLCANOES.txt | awk -F= '{print $2}'`
+        LMDZ_sed physiq.def solaire     ${IPCC_SOLAR}
+    fi
+    if [ -f CO2.txt ] ; then
+        IPCC_CO2=`grep Annee_${year} CO2.txt | awk -F= '{print $2}'`
+        LMDZ_sed physiq.def co2_ppm     ${IPCC_CO2}
+    fi
+    if [ -f CH4.txt ] ; then
+        IPCC_CH4=`grep Annee_${year} CH4.txt | awk -F= '{print $2}'`
+        LMDZ_sed physiq.def CH4_ppb     ${IPCC_CH4}
+    fi
+    if [ -f N2O.txt ] ; then
+        IPCC_N2O=`grep Annee_${year} N2O.txt | awk -F= '{print $2}'`
+        LMDZ_sed physiq.def N2O_ppb     ${IPCC_N2O}
+    fi
+    if [ -f CFC11.txt ] ; then
+        IPCC_CFC11=`grep Annee_${year} CFC11.txt | awk -F= '{print $2}'`
+        LMDZ_sed physiq.def CFC11_ppt   ${IPCC_CFC11}
+    fi
+    if [ -f CFC12.txt ] ; then
+        IPCC_CFC12=`grep Annee_${year} CFC12.txt | awk -F= '{print $2}'`
+        LMDZ_sed physiq.def CFC12_ppt   ${IPCC_CFC12}
+    fi
 
     ## physiq.def parameters modified from initial physiq_L${RESOL_ATM_Z}

CONFIG/IPSLCM/IPSLCM5-2/piControl/COMP/pisces.driver

@@ -11 +11 @@
 }
 
-    ##--Variables used by OPA --
-
-
     ##--Variables used by MBG --
 
@@ -23 +20 @@
     PAT_TOP_LRSTTR=$( supergrep ln_rsttr          ${SUBMIT_DIR}/PARAM/namelist_top_${RESOL_OCE} )
     PAT_TOP_NRSTTR=$( supergrep nrsttr            ${SUBMIT_DIR}/PARAM/namelist_top_${RESOL_OCE} )
+    PAT_PIS_ATCCO2=$( supergrep atcco2            ${SUBMIT_DIR}/PARAM/namelist_pisces_${RESOL_OCE} )
     PAT_PIS_LRIVER=$( supergrep ln_river          ${SUBMIT_DIR}/PARAM/namelist_pisces_${RESOL_OCE} )
 
@@ -77 +75 @@
     IGCM_sys_Mv iodef.xml.tmp iodef.xml
 
+    ##-- CO2 forcing : update atcco2 value in namelist_pisces for some cases
+    if [ "${piseces_UserChoices_CARBON_CYCLE}" = "y" ]; then
+        # Special case for carbon cycle
+        # Activate this case by putting CARBON_CYCLE=y in pisces.card section UserChoices
+        echo 'Update atcco2 in namelist_pisces with value from file co2.log'
+        PathCO2log=${SUBMIT_DIR}
+        PreviousLinelog=$( tail -1 ${PathCO2log}/co2.log )
+        IPCC_PIS_CO2=$( echo "${PreviousLinelog[*]}" | gawk -F '|' '{gsub(" ",""); print $4}' )
+        sed -e "s%${PAT_PIS_ATCCO2}%          atcco2=${IPCC_PIS_CO2}%"             \
+            namelist_pisces > namelist_pisces.tmp
+        IGCM_sys_Mv namelist_pisces.tmp namelist_pisces
+    elif [ -f CO2.txt ] ; then
+        # This case is activated if the file CO2.txt is found in run directory, add file in BoundaryFiles in pisces.card or lmdz.card
+        echo 'Update atcco2 in namelist_pisces with annual value from file CO2.txt'
+        IPCC_PIS_CO2=`grep Annee_${year} CO2.txt | awk -F= '{print $2}'`
+        sed -e "s%${PAT_PIS_ATCCO2}%          atcco2=${IPCC_PIS_CO2}%"             \
+            namelist_pisces > namelist_pisces.tmp
+        IGCM_sys_Mv namelist_pisces.tmp namelist_pisces
+    fi
+
     ##-- Update namelist_top and namelist_pisces