Changeset 10 for branches/AllPostFred/AA_create_se

Timestamp:

05/30/08 10:10:34 (16 years ago)

Author:

mmaipsl

Message:

MM: All frequencies in post-treatement.

File:

• branches/AllPostFred/AA_create_se (modified) (4 diffs)

branches/AllPostFred/AA_create_se

@@ -139 +139 @@
 DateBegin=${DateBegin:=$( IGCM_date_ConvertFormatToGregorian ${config_UserChoices_DateBegin} )}
 #
-PeriodDateEnd=${PeriodDateEnd:=$( IGCM_date_ConvertFormatToGregorian ${config_UserChoices_DateEnd} )}
+DateEnd=${DateEnd:=$( IGCM_date_ConvertFormatToGregorian ${config_UserChoices_DateEnd} )}
 #
 freq=$( echo ${config_Post_SeasonalFrequency} |  sed -e "s/[yY]//" )
 #
-an_fin=$( expr $( IGCM_date_ConvertGregorianDateToJulian $PeriodDateEnd ) / 1000 )
+an_fin=$( expr $( IGCM_date_ConvertGregorianDateToJulian $DateEnd ) / 1000 )
 an_deb=$( IGCM_date_YearDigit $(( an_fin - freq + 1 )) )
 an_fin=$( IGCM_date_YearDigit ${an_fin} )
@@ -227 +227 @@
         FlagDir=$( echo ${file} | awk -F "_" '{print $1}' )
         case ${FlagDir} in
-            1Y)    FreqDir=YE  ;;
-            1M)    FreqDir=MO  ;;
-            5D|1D) FreqDir=DA  ;;
+            *Y)    FreqDir=YE  ;;
+            *M)    FreqDir=MO  ;;
+            5D|1D|*D) FreqDir=DA  ;;
             HF)    FreqDir=HF  ;;
             INS)   FreqDir=INS ;;
@@ -286 +286 @@
           IGCM_debug_Verif_Exit_Post
           ;;
+        # For a Multi-Year simulation :
+        *Y|*y)
+          PeriodLengthInYears=$( echo ${config_UserChoices_PeriodLength} | awk -F '[yY]' '{print $1}' )
+          (( an = an_deb ))
+          # Initialize array
+          set -A liste_file_se
+          set -A liste_file_se_tmp
+          while [ ${an} -le ${an_fin} ] ; do
+              moisL=$( IGCM_date_DaysInMonth ${an} ${mois} )
+              #
+              an4=$(printf "%04i\n" ${an} )
+              #
+              (( diy = 0 ))
+              (( ian = 0 ))
+              while [ ${ian} -lt ${PeriodLengthInYears} ] ; do
+                  #
+                  ian4=$(printf "%04i\n" $(( an + ian )) )
+                  (( diy = diy + $( IGCM_date_DaysInYear $ian4 ) ))
+                  (( ian = ian + 1 ))
+              done
+              afin=$( IGCM_date_AddDaysToGregorianDate ${an4}0101 $(( diy - 1 )) )
+              date=${an4}0101_${afin}
+              set +A liste_file_se ${liste_file_se[*]} ${file_path}${config_UserChoices_JobName}_${date}_${file}.nc
+              set +A liste_file_se_tmp ${liste_file_se_tmp[*]} ${config_UserChoices_JobName}_${date}_${file}.nc
+              (( an = an + PeriodLengthInYears ))
+          done
+          #
+          IGCM_sys_Get /l liste_file_se[*] ${RUN_DIR_PATH}
+          eval ExitFlag=false
+          IGCM_debug_Verif_Exit_Post
+          
+          # Apply IGCM_Patch if needed
+          if [ X$( eval echo \${LISTE_PATCHES_${file}[0]} ) !=  X${NULL_STR} ]; then
+              for file_t in $( eval echo ${liste_file_se_tmp[*]} ); do
+                  for Patch in $( eval echo \${LISTE_PATCHES_${file}[*]} ); do
+                      IGCM_${Patch} ${file_t}
+                  done
+              done
+          fi
+
+          for mois in  01 02 03 04 05 06 07 08 09 10 11 12 ; do
+          #
+              (( nb = ${mois} - 1 ))
+          #
+              OUT_SE[${nb}]=${config_UserChoices_JobName}_SE${mois}_${an_deb}_${an_fin}_${file}.nc
+              IGCM_sys_ncra -O -d time_counter,${nb},,12 ${liste_file_se_tmp[*]} ${OUT_SE[${nb}]}
+          done
+          IGCM_sys_Rm ${liste_file_se_tmp[*]}
+          eval ExitFlag=false
+          IGCM_debug_Verif_Exit_Post
+          ;;
         # For Monthly simulation :
         1M|1m)
@@ -322 +373 @@
               fi
               
-              IGCM_sys_ncra -O ${liste_file_se_tmp[*]} ${OUT_SE[${nb}]}
-              IGCM_sys_Rm ${liste_file_se_tmp[*]}
-              eval ExitFlag=false
+              IGCM_sys_ncra -O ${liste_file_se_tmp} ${OUT_SE[${nb}]}
+              MENAGE_LOCAL ${liste_file_se_tmp}
+              #
               IGCM_debug_Verif_Exit_Post
           done
+          ;;
+
+        # For Multi-Monthly simulation :
+        *M|*m)
+          PeriodLengthInMonths=$( echo ${config_UserChoices_PeriodLength} | awk -F '[mM]' '{print $1}' )
+          # Must be a divisor of 12 !!! e.g. 2,3,4 and 6 only
+          if ( [ $PeriodLengthInMonths -ne 2 ] && \
+               [ $PeriodLengthInMonths -ne 3 ] && \
+               [ $PeriodLengthInMonths -ne 4 ] && \
+               [ $PeriodLengthInMonths -ne 6 ] ) ; then
+              echo "create_se.job can't work with this monthly Period : ",$PeriodLengthInMonths,"M !!"
+              echo "It must be a divisor of twelve."
+              echo "We will stop there."
+              IGCM_debug_Verif_Exit_Post
+          else
+              #
+              (( NbSet = 12 / PeriodLengthInMonths ))
+              (( moisd = 01 ))
+              (( Mset = 1 )) 
+              while [ ${Mset} -le ${NbSet} ] ; do
+                  moisf=$(printf "%02i\n" $(( moisd + PeriodLengthInMonths - 1 )) )
+                  moisd=$(printf "%02i\n" ${moisd} )
+                  #
+                  (( an = an_deb ))
+                  # Initialize array
+                  set -A liste_file_se
+                  set -A liste_file_se_tmp
+                  while [ ${an} -le ${an_fin} ] ; do
+                      moisL=$( IGCM_date_DaysInMonth ${an} ${moisf} )
+                      #
+                      an4=$(printf "%04i\n" ${an} )
+                      #
+                      date=${an4}${moisd}01_${an4}${moisf}${moisL}
+                      set +A liste_file_se ${liste_file_se[*]} ${file_path}${config_UserChoices_JobName}_${date}_${file}.nc
+                      set +A liste_file_se_tmp ${liste_file_se_tmp[*]} ${config_UserChoices_JobName}_${date}_${file}.nc
+                      (( an = an + 1 ))
+                  done
+                  #
+                  IGCM_sys_Get /l liste_file_se[*] ${RUN_DIR_PATH}
+                  eval ExitFlag=false
+                  IGCM_debug_Verif_Exit_Post
+              
+                  # Apply IGCM_Patch if needed
+                  if [ X$( eval echo \${LISTE_PATCHES_${file}[0]} ) !=  X${NULL_STR} ]; then
+                      for file_t in $( eval echo ${liste_file_se_tmp[*]} ); do
+                          for Patch in $( eval echo \${LISTE_PATCHES_${file}[*]} ); do
+                              IGCM_${Patch} ${file_t}
+                          done
+                      done
+                  fi
+        
+                  (( imois = moisd ))
+                  (( i = 0 ))
+                  while [ ${imois} -le ${moisf} ] ; do
+                      #
+                      (( nb = ${imois} - 1 ))
+                      #
+                      OUT_SE[${nb}]=${config_UserChoices_JobName}_SE${imois}_${an_deb}_${an_fin}_${file}.nc
+                      IGCM_sys_ncra -O -d time_counter,${i},,${PeriodLengthInMonths} ${liste_file_se_tmp[*]} ${OUT_SE[${nb}]}
+                      (( imois = imois + 1 ))
+                      (( i = i + 1 ))
+                  done
+
+                  IGCM_sys_Rm ${liste_file_se_tmp[*]}
+                  eval ExitFlag=false
+                  IGCM_debug_Verif_Exit_Post
+                  moisf=$(printf "%02i\n" $(( moisf + 1 )) )
+
+                  (( Mset = Mset + 1 ))
+              done
+          fi
           ;;
         *D|*d)