Changeset 259 for branches/ORCHIDEE_EXT/ORCHIDEE_OL/FLUXNET/Job_FLUXNET_validation

Timestamp:

2011-06-17T16:55:45+02:00 (13 years ago)

Author:

didier.solyga

Message:

Forget to add ENSEMBLE from the trunk version to the merge

File:

• branches/ORCHIDEE_EXT/ORCHIDEE_OL/FLUXNET/Job_FLUXNET_validation (modified) (5 diffs)

branches/ORCHIDEE_EXT/ORCHIDEE_OL/FLUXNET/Job_FLUXNET_validation

@@ -212 +212 @@
 
     IGCM_card_WriteOption ${New_SUBMIT_DIR}/COMP/spinup.card UserChoices DRIVER_NORESTART y
+    IGCM_card_WriteOption ${New_SUBMIT_DIR}/COMP/spinup.card UserChoices DRIVER_TIMELENGTH n
 
     eval IGCM_card_WriteOption ${New_SUBMIT_DIR}/COMP/spinup.card UserChoices duree_nostomate $( correct_duree ${fluxnet_SPINUP_duree_nostomate} ${TIME_YEAR} )
@@ -225 +226 @@
 
     # TimeSeries : ALMAOUTPUT for SECHIBA
-
-# WARNING : You can use this first option if you set CONSERVE=y in fluxnet.card
+    typeset option
+    for option in ${config_SubJobPost[*]} ; do
+        eval value=\${config_SubJobPost_${option}}
+        eval echo ${option} ${value}
+        if [ X${value} != X ] ; then
+            eval IGCM_card_WriteOption ${New_SUBMIT_DIR}/config.card SubJobPost ${option} \${config_SubJobPost_${option}}
+        fi
+    done
 #    eval IGCM_card_WriteOption ${New_SUBMIT_DIR}/config.card SRF TimeSeriesFrequency ${TIME_YEAR}Y
-    eval IGCM_card_WriteOption ${New_SUBMIT_DIR}/config.card SRF TimeSeriesFrequency NONE
+#    eval IGCM_card_WriteOption ${New_SUBMIT_DIR}/config.card SRF TimeSeriesFrequency NONE
     set -A sechiba_TimeSeriesVars2D -- nobiofrac SWnet LWnet Qh Qle Qg Qf DelSurfHeat DelColdCont Snowf Rainf Evap \
                                      Qs Qsb Qsm DelSoilMoist DelSWE DelIntercept AvgSurfT RadT Albedo SWE SoilMoist SoilWet SoilTemp PotEvap \
@@ -247 +254 @@
 
 
-    # modification on orchidee.def
-    IGCM_sys_Cp ${SUBMIT_DIR}/orchidee.def ${New_SUBMIT_DIR}/PARAM/orchidee.def
+    # modification on parameter files
+    IGCM_sys_Cp ${SUBMIT_DIR}/PARAM/sechiba.def ${New_SUBMIT_DIR}/PARAM/sechiba.def
+    IGCM_sys_Cp ${SUBMIT_DIR}/PARAM/stomate.def ${New_SUBMIT_DIR}/PARAM/stomate.def
+    IGCM_sys_Cp ${SUBMIT_DIR}/PARAM/driver.def ${New_SUBMIT_DIR}/PARAM/driver.def
 
     eval first_option=${fluxnet_SubJobParams[0]} > /dev/null 2>&1
@@ -255 +264 @@
         eval first_option=${spinup_SubJobParams[0]} > /dev/null 2>&1
             
-        typeset option option_
+        typeset option option_ PARAM_FILE
         for option in ${fluxnet_SubJobParams[*]} ; do
-            eval echo "${option}="\${fluxnet_SubJobParams_${option}}
-            eval ORCHIDEE_sed ${option} \${fluxnet_SubJobParams_${option}}
+            PARAM_FILE=$( echo ${option} | awk -F_ '{print $1}' )
+            echo "For parameter file ${PARAM_FILE}.def"
+            option_=$( echo ${option} | sed -e "s/${PARAM_FILE}_//" )
+            eval echo "${option_}="\${fluxnet_SubJobParams_${option}}
+            case ${PARAM_FILE} in
+                sechiba)
+                    eval SECHIBA_sed ${option_} \${fluxnet_SubJobParams_${option}}
+                    ;;
+                stomate)
+                    eval STOMATE_sed ${option_} \${fluxnet_SubJobParams_${option}}
+                    ;;
+                driver)
+                    eval DRIVER_sed ${option_} \${fluxnet_SubJobParams_${option}}
+                    ;;
+            esac
 
             if [ X${first_option} != X"Error:" ] ; then
@@ -271 +293 @@
 
     # Replace vegetation
-    eval ORCHIDEE_sed IMPOSE_VEG y
-
-    ipft=1
-    while [ $ipft -lt 14 ] ; do
-        if [ $ipft -lt 10 ] ; then
-            numpft=0${ipft}
-        else
-            numpft=${ipft}
-        fi
-        ORCHIDEE_sed SECHIBA_VEG__${numpft} ${fluxnet_FLUXNET_Sites[$(( index + ColumnPFT + ipft - 1 )) ]}
-        ORCHIDEE_sed SECHIBA_VEGMAX__${numpft} ${fluxnet_FLUXNET_Sites[$(( index + ColumnPFT + ipft - 1 )) ]}
-        ORCHIDEE_sed SECHIBA_LAI__${numpft} ${fluxnet_FLUXNET_Sites[$(( index + ColumnLAI + ipft - 1 )) ]}
-        let ipft=ipft+1
+    IGCM_card_WriteOption ${New_SUBMIT_DIR}/COMP/spinup.card UserChoices impose_veg y
+
+    (( iphys = 0 ))
+    while [ $iphys -lt ${fluxnet_FLUXNET_NbSitesParam} ] ; do
+        
+        # Detect if physical parameter has 
+        #  array description or line per PFT description in the parameter file
+        SearchParam=""
+        ${fluxnet_FLUXNET_CompSitesParam[${iphys}]}_grep ${fluxnet_FLUXNET_NameSitesParam[${iphys}]}
+
+        case ${SearchParam} in
+            undefined)
+                echo "ERROR : parameter ${fluxnet_FLUXNET_NameSitesParam[${iphys}]} doen't exist in ${fluxnet_FLUXNET_CompSitesParam[${iphys}]} parameter file !"
+                echo "You must correct NameSiteParam option in fluxnet.card file or add new parameter in component.def file"
+                echo "We must STOP here."
+                exit 1
+                ;;
+            line)
+                ipft=1
+                while [ $ipft -le ${NbPFTs} ] ; do
+                    if [ $ipft -lt 10 ] ; then
+                        numpft=0${ipft}
+                    else
+                        numpft=${ipft}
+                    fi
+            
+                    ${fluxnet_FLUXNET_CompSitesParam[${iphys}]}_sed ${fluxnet_FLUXNET_NameSitesParam[${iphys}]}__${numpft} ${fluxnet_FLUXNET_Sites[$(( index + 4 + ( iphys * NbPFTs ) + ipft - 1 )) ]}
+                    (( ipft = ipft + 1 ))
+                done
+                ;;
+            vector)
+                unset ValueSiteParams
+                ipft=0
+                while [ $ipft -lt ${NbPFTs} ] ; do
+                    ValueSiteParams[${ipft}]=${fluxnet_FLUXNET_Sites[$(( index + 4 + ( iphys * NbPFTs ) + ipft - 1 )) ]}
+                    (( ipft = ipft + 1 ))
+                done
+                ${fluxnet_FLUXNET_CompSitesParam[${iphys}]}_sed ${fluxnet_FLUXNET_NameSitesParam[${iphys}]} "${ValueSiteParams[*]}"
+                ;;
+        esac
+
+        (( iphys = iphys + 1 ))
     done