Changeset 1316 for CONFIG/IPSLCM/IPSLCM5A/CpiControl/COMP

Timestamp:

02/22/11 10:36:38 (13 years ago)

Author:

cetlod

Message:

Improvment of interactive scripts for IPSLCM5A_C configuration

Location:

CONFIG/IPSLCM/IPSLCM5A/CpiControl/COMP

Files:

• lmdz.driver (modified) (2 diffs)
• lmdz_analyse_stomate_out.awk (modified) (5 diffs)
• pisces.driver (modified) (3 diffs)

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

@@ -141 +141 @@
 
     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", \
+        ##--Initialization of fluxes to an undefined value at the first run
+        UndefinedValueCo2=-9999
+        InitPeriodCo2=0
+        InitDateBeginCo2=${UndefinedValueCo2}
+        InitYearCo2=$(( ${year} - 1 )) 
+        InitDateEndCo2=${InitYearCo2}1231
+        fCO2_ff=${UndefinedValueCo2} 
+        CO2SBG=${UndefinedValueCo2}
+        CO2LU=${UndefinedValueCo2}
+        CO2MBG=${UndefinedValueCo2}
+        CO2_ppm_prec=${UndefinedValueCo2}
+        # Get the initial value of atmosp. pco2 
+        CO2_ppm=${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
+    
+        # Save CO2 values in ExeCpuLog variable contents 5 fields
+        echo "${InitPeriodCo2} ${InitDateBeginCo2} ${InitDateEndCo2} ${CO2_ppm} ${CO2_ppm_prec} ${fCO2_ff} ${CO2MBG} ${CO2SBG} ${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
+    
+        # Get the path where the log file co2.log is
+        PathCO2log=${SUBMIT_DIR} 
+        # Get Previous line in co2.log
+        PreviousLinelog=$( tail -1 ${PathCO2log}/co2.log )
 
     else
-
-#       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
+        # Get restart line in co2.log
+        PathCO2log=${SUBMIT_DIR}
+        PreviousLinelog=$( tail -1 ${PathCO2log}/co2.log ) 
+        # Get the value of atmosp. pco2 in co2.log
+        CO2_ppm=$( echo "${PreviousLinelog[*]}" | gawk -F '|' '{gsub(" ",""); print $4}' )
 
     fi
@@ -306 +263 @@
     [ ${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}
 
+    typeset CO2SBG_3V CO2SBG CO2MBG CO2SRF CO2LU CO2_ppm CO2_ppm_prec fCO2_ff
+
+    # For Control, no fossil fuel emission :
+    fCO2_ff=0
+    echo "Fossil Fuel for control run  :" ${fCO2_ff}
+
+    # Get the value of ocean carbon flux
+    CO2MBG=$( ${SUBMIT_DIR}/COMP/lmdz_analyse_pisces_out.awk ocean.output )
+    CO2MBG=$( echo ${CO2MBG} | sed -e "s/ *//" )
+    echo "Ocean carbon flux  :" ${CO2MBG}
+
+    # Get the value of land fluxes
+    set +A CO2SBG_3V -- $( ${SUBMIT_DIR}/COMP/lmdz_analyse_stomate_out.awk out_orchidee_0000 )
+    CO2SRF=$( echo ${CO2SBG_3V[0]} | sed -e "s/ *//" )
+    CO2LU=$( echo ${CO2SBG_3V[1]} | sed -e "s/ *//" )
+    CO2SBG=$( echo ${CO2SBG_3V[2]} | sed -e "s/ *//" )
+    echo "Land carbon flux (NEP)         :" ${CO2SRF}
+    echo "LU flux          (FLUC)        :" ${CO2LU}
+    echo "Total Land carbon flux (NBP)   :" ${CO2SBG}
+
+    # Previous CO2
+    CO2_ppm_prec=$( echo "${PreviousLinelog[*]}" | gawk -F '|' '{gsub(" ",""); print $4}' )
+    echo "Previous CO2  :" ${CO2_ppm_prec}
+    
+    # Formula
+    CO2_ppm=$( echo "${CO2_ppm_prec} + (${fCO2_ff} + ${CO2MBG} + ${CO2SBG}) / 2.12" | bc -l )
+    echo "New CO2  :" ${CO2_ppm}
+
+    # Save CO2 values in ExeCpuLog variable contents 5 fields
+    echo "${CumulPeriod} ${PeriodDateBegin} ${PeriodDateEnd} ${CO2_ppm} ${CO2_ppm_prec} ${fCO2_ff} ${CO2MBG} ${CO2SBG} ${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
+
     echo FINALIZE ATM !
 

CONFIG/IPSLCM/IPSLCM5A/CpiControl/COMP/lmdz_analyse_stomate_out.awk

@@ -25 +25 @@
 #==========================
 BEGIN {
-#  print "traitement de " ARGV[1]
+ # print "traitement de " ARGV[1]
 
   nbarg=ARGC
@@ -54 +54 @@
   counterLU=0
 
+  fluxNBP_found=0
+  counterNBP=0
+
   CO2file[1]=""
   LUfile[1]=""
+  NBPfile[1]=""
 }
 
@@ -79 +83 @@
     CO2file[counterCO2]=felts[2]
   }
-  else if (match($0, ".*GLOBAL net cflux_prod_total_sum.*")) {
+  else if (match($0, ".*GLOBAL net_cflux_prod_monthly.*")) {
     myprint("fluxLU_found = 1")
     fluxLU_found=1
@@ -87 +91 @@
     myprint("nb  :" nb)
     for (elt in felts) {
+        myprint(elt "- elt  :" felts[elt])
+    }
+    myprint("elts_co2  :" felts[2])
+
+    LUfile[counterLU]=felts[2]
+  }
+
+  else if (match($0, ".*GLOBAL net_biosp_prod_monthly.*")) {
+    myprint("fluxNBP_found = 1")
+    fluxNBP_found=1
+    counterNBP=counterNBP+1
+
+    nb=split($0,felts, "=")
+    myprint("nb  :" nb)
+    for (elt in felts) {
         myprint(elt "- elt  :" felts[elt])
     }
-    myprint("elts_lu  :" felts[2])
+    myprint("elts_nbp  :" felts[2])
 
-    LUfile[counterLU]=felts[2]
+    NBPfile[counterNBP]=felts[2]
   }
 }
@@ -106 +125 @@
             myprint("value  :" CO2file[i])
         }
-        
-        if ( fluxLU_found == 1 ) {
-            myprint("LU value for " file)
-            myprint("counterLU  :" counterLU)
-            for (i=1; i<=counterLU; i++) {  
-                print LUfile[i]
-                myprint("value  :" LUfile[i])
+
+        if ( fluxLU_found == 1 ) {
+            myprint("LU value for " file)
+            myprint("counterLU  :" counterLU)
+            for (i=1; i<=counterLU; i++) {
+                print LUfile[i]
+                myprint("value  :" LUfile[i])
+            }
+
+            if ( fluxNBP_found == 1 ) {
+               myprint("NBP value for " file)
+               myprint("counterNBP  :" counterNBP)
+               for (i=1; i<=counterNBP; i++) {  
+                    print NBPfile[i]
+                    myprint("value  :" NBPfile[i])
+               }
+            
+               exit_value=0
+               exit 0
             }
-            
+            else {
+                 myprint("ERROR : NO NBP value for " file)
+                 print "0003"
+                 exit 1
+            }
+
             exit_value=0
             exit 0
-        }
-        else {
-            myprint("ERROR : NO LU value for " file)
-            print "0002"
-            exit 1
-        }
-
-        exit_value=0
-        exit 0
+        }
+        else {
+            myprint("ERROR : NO LU value for " file)
+            print "0002"
+            print "0003"
+            exit 1
+        }
     }
     else {
-        myprint("ERROR : NO CO2 value for " file)
-        print "0001"
-        print "0002"
-        exit 1
-    }
-    
+         myprint("ERROR : NO CO2 value for " file)
+         print "0001"
+         print "0002"
+         print "0003"
+         exit 1
+       }
   }
   else {

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

@@ -23 +23 @@
     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} )
 
@@ -60 +61 @@
     fi
 
+    ##-- GHG forcing for the current year. Normally constant
+    PathCO2log=${SUBMIT_DIR}
+    PreviousLinelog=$( tail -1 ${PathCO2log}/co2.log )
+    IPCC_PIS_CO2=$( echo "${PreviousLinelog[*]}" | gawk -F '|' '{gsub(" ",""); print $4}' )
+
     ##-- Update namelist_top and namelist_pisces 
 
@@ -68 +74 @@
     IGCM_sys_Mv namelist_top.tmp namelist_top
 
-    sed -e "s%${PAT_PIS_LRIVER}%          ln_river=${PISCES_LRIVER}%"                \
+    sed -e "s%${PAT_PIS_ATCCO2}%          atcco2=${IPCC_PIS_CO2}%"                \
+        -e "s%${PAT_PIS_LRIVER}%          ln_river=${PISCES_LRIVER}%"                \
         namelist_pisces > namelist_pisces.tmp