Changeset 1596 for branches/libIGCM/ORCA2_LIM/IGCM00/COMP/opa9.driver

Timestamp:

2009-08-10T12:11:36+02:00 (15 years ago)

Author:

flavoni

Message:

Changes in opa9.card and opa9.driver to have possibility to do interannual runs, see ticket : #524

File:

• branches/libIGCM/ORCA2_LIM/IGCM00/COMP/opa9.driver (modified) (2 diffs)

branches/libIGCM/ORCA2_LIM/IGCM00/COMP/opa9.driver

@@ -52 +52 @@
     (( DaysSinceJC = $( IGCM_date_DaysSinceJC ${DateBegin} ) + 1 ))
 
+    # Definition from opa9.card of List_jobsbc and List_runsbc used to do to interannual and climatological runs 
+   IGCM_card_DefineArrayFromOption ${SUBMIT_DIR}/COMP/opa9.card Interannual List_jobsbc
+   IGCM_card_DefineArrayFromOption ${SUBMIT_DIR}/COMP/opa9.card Interannual List_runsbc
+        set -A ListFormulationJobsbc -- \${opa9_Interannual_Listjobsc${opa9_Interannual_List_jobsbc}[*]}
+        set -A ListFormulationRunsbc -- \${opa9_Interannual_Listrunsbc${opa9_Interannual_List_runsbc}[*]}
+
     IGCM_debug_PopStack "OCE_Initialize"
 }
@@ -58 +64 @@
 function OCE_Update
 {
-###    set -vx
     IGCM_debug_PushStack "OCE_Update"
+
+    # Interannual run
+    # In this case job need to calculate previous and following year and give, and if every run needs to add "_y" before every year 
+    if [ X${opa9_UserChoices_Interannual_Run} = Xy ] ; then
+
+   typeset file fileo
+
+        eval NbFileInter=${#opa9_Interannual_List_jobsbc[*]}
+
+   (( i = 0 ))
+   while [ $i -lt ${NbFileInter} ] ; do
+       eval file=${opa9_Interannual_List_jobsbc[$i]}
+       eval jobsbc_file_ym1=${R_BC}/OCE/${config_UserChoices_TagName}/${opa9_UserChoices_OPA_version}/${file}_$(( year - 1 )).nc
+       eval jobsbc_file_y=${R_BC}/OCE/${config_UserChoices_TagName}/${opa9_UserChoices_OPA_version}/${file}_${year}.nc
+       eval jobsbc_file_yp1=${R_BC}/OCE/${config_UserChoices_TagName}/${opa9_UserChoices_OPA_version}/${file}_$(( year + 1 )).nc
+     # Copy of List_jobsbc in List_runsbc if this last one is empty 
+            if [ X${opa9_Interannual_List_runsbc[$i]} = X${NULL_STR} ] ; then
+         eval fileo=${opa9_Interannual_List_jobsbc[$i]}
+               eval runsbc_file_ym1=${file}_y$(( year - 1 )).nc
+               eval runsbc_file_y=${file}_y${year}.nc
+               eval runsbc_file_yp1=${file}_y$(( year + 1 )).nc
+       else
+         eval fileo=${opa9_Interannual_List_runsbc[$i]}
+               eval runsbc_file_ym1=${fileo}_y$(( year - 1 )).nc
+               eval runsbc_file_y=${fileo}_y${year}.nc
+               eval runsbc_file_yp1=${fileo}_y$(( year + 1 )).nc
+      fi
+
+       if [ X${Period} = X1 ] ; then
+      if [ ${month} -eq 01 ] ; then
+          IGCM_sys_Get ${jobsbc_file_ym1} ${runsbc_file_ym1}
+      fi
+      if [ ! -f ${runsbc_file_y} ] ; then
+          IGCM_sys_Get ${jobsbc_file_y} ${runsbc_file_y}
+      fi
+      if [ ! -f ${runsbc_file_yp1} ] ; then
+          IGCM_sys_Get ${jobsbc_file_yp1} ${runsbc_file_yp1}
+      fi
+       else 
+      if [ ! -f ${runsbc_file_yp1} ] ; then
+          IGCM_sys_Get ${jobsbc_file_yp1} ${runsbc_file_yp1}
+      fi
+       fi 
+
+       (( i = i + 1 ))
+   done
+    # End interannual
+    elif [ X${opa9_UserChoices_Interannual_Run} = Xn ] ; then
+    # Climatological run
+       typeset file fileo
+
+       eval NbFileInter=${#opa9_Interannual_List_jobsbc[*]}
+
+        (( i = 0 ))
+        while [ $i -lt ${NbFileInter} ] ; do
+            eval file=${opa9_Interannual_List_jobsbc[$i]}
+            eval jobsbc_file=${R_BC}/OCE/${config_UserChoices_TagName}/${opa9_UserChoices_OPA_version}/${file}.nc
+            eval NbFilerunsbc=${#opa9_Interannual_List_runsbc[*]}
+            if [ ${NbFilerunsbc} -eq 0  ] ; then
+                eval fileo=${opa9_Interannual_List_jobsbc[$i]}
+                eval runsbc_file=${file}.nc
+            else
+                eval fileo=${opa9_Interannual_List_runsbc[$i]}
+                eval runsbc_file=${fileo}.nc
+            fi
+            IGCM_sys_Get ${jobsbc_file} ${runsbc_file}
+            (( i = i + 1 ))
+   done
+    fi
+    # End climatological
 
     NbFreq=$( echo ${config_OCE_WriteFrequency} | wc -w )