Changeset 1550 for CONFIG/IPSLCM/IPSLCM5-2/EXPALL/DRIVER/pisces.driver

Timestamp:

08/30/11 14:53:52 (13 years ago)

Author:

jgipsl

Message:

• Added choice of config.def file : set ConfType?=preind/actuel/annuel in lmdz.card. The file config.def_$ConfType? will be used. The differences are :

preind : aer_type = preind, solar and GHG pre-industrial values, pmagic=0 (as piControl, CpiControl?)
actuel : aer_type = actuel, solar and GHG present values, pmagic=0.008 (as EXP00)
annuel : aer_type = annuel, solar and GHG must be updated with annual values from file, pmagic=0.008 (as historical,Chistorical)

• For pisces, take same value for CO2 as for LMDZ except for CarbonCycle?=y (same treatment as before).

NB! I suppose that LMDZ ATM_Uptdate always is done before pisces MBG_Update.

• Added choice of output.def file : set OutLev?=low(as EXP00)/medium(as piControl)/high(as histrical) in lmdz.card. The output are now the same as original IPSLCM5A for each experience.

=> In EXPALL, all 5 experiences correspond to original IPSLCM5A experiences, including possibility of new physics in LMDZ as IPSLCM5B and carbon cycle as IPSLCM5C. More testing is needed.

File:

• CONFIG/IPSLCM/IPSLCM5-2/EXPALL/DRIVER/pisces.driver (modified) (3 diffs)

CONFIG/IPSLCM/IPSLCM5-2/EXPALL/DRIVER/pisces.driver

@@ -30 +30 @@
 {
     IGCM_debug_PushStack "MBG_Update"
+
+    # Local function to find parameters in LMDZ .def files
+    lmdzgrep () {
+        grep $1 $2 |grep -v "#" | awk "-F=" '{print $2}'
+    }
 
     ##-- Restart configuration
@@ -75 +80 @@
     IGCM_sys_Mv iodef.xml.tmp iodef.xml
 
-    ##-- CO2 forcing : update atcco2 value in namelist_pisces for some cases
+    ##-- CO2 forcing : update atcco2 value in namelist_pisces
     if [ "${piseces_UserChoices_CARBON_CYCLE}" = "y" ]; then
         # Special case for carbon cycle
@@ -86 +91 @@
             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}%"             \
+    else
+        # Take same CO2 as for LMDZ, in parameter file config.def.
+        # ATM_Update is done before MBG_Update and therefor config.def file contains the correct value for CO2.
+        echo 'Update atcco2 in namelist_pisces with same value as for LMDZ'
+        LMDZ_CO2=$( lmdzgrep co2_ppm config.def )
+        sed -e "s%${PAT_PIS_ATCCO2}%          atcco2=${LMDZ_CO2}%"             \
             namelist_pisces > namelist_pisces.tmp
         IGCM_sys_Mv namelist_pisces.tmp namelist_pisces