Changeset 3097 for TOOLS

Timestamp:

01/25/17 16:51:45 (8 years ago)

Author:

jgipsl

Message:

Added case to create yearly files

File:

• TOOLS/INTERP_NUDGE/interp_from_era.ksh (modified) (5 diffs)

TOOLS/INTERP_NUDGE/interp_from_era.ksh

@@ -13 +13 @@
 
 # gridfile : file grilles_gcm.nc containing longitudes and latitudes of the destination grid
-gridfile=/ccc/store/cont003/dsm/p86ghatt/IGCM_OUT/LMDZ/ELC-128x118x39/ATM/Output/Grid/ELC-128x118x39_grilles_gcm.nc
+gridfile=/ccc/store/cont003/dsm/p86ghatt/IGCM_OUT/LMDZ/ELI-144x142x79/ATM/Output/Grid/ELI-144x142x79_1979_grilles_gcm.nc
 
 # indir : path for directory containing the files to interpolate
 #indir=/u/rech/psl/rpsl376/ERAI/NETCDF/GLOBAL_075/4xdaily/AN_PL
-indir=/ccc/store/cont003/subipsl/subipsl/dmf_import/ERAI/NETCDF/GLOBAL_075/4xdaily/AN_PL
+indir=/ccc/work/cont003/subipsl/subipsl/ERAI/NETCDF/GLOBAL_075/4xdaily/AN_PL/
 # insuffix : suffix in input files
 insuffix=.aphei.GLOBAL_075.nc
@@ -26 +26 @@
 # Variables for output :
 # outdir : directory to store output files
-outdir=$WORKDIR/NUDGEFILES/ERAI-128x118x39
+outdir=$WORKDIR/NUDGEFILES/ERAI-144x142x79
 # outsuffix : suffix to output files
 outsuffix=_erai.nc
@@ -32 +32 @@
 # Time variables :
 # first_year : first year to interpolate
-first_year=2000
+first_year=1979
 # last_year : last year to interpolate
-last_year=2000
+last_year=2016
 
+# Frequency of output files : yearly or monthly
+# OutFreq='monthly' or OutFrequ='yearly'
+OutFreq='yearly'
+ 
 # rundir : temporary run directory
 rundir=`pwd`/rundir
@@ -99 +103 @@
 
 
+    # Do different treatement if monthly or yearly output files are wanted
+    if [[ $OutFreq = 'monthly' ]] ; then
     # 2 - Loop over all files and add a time step : 
     #     take the first time step from next coming month and add last in current month
@@ -139 +145 @@
     done
 
+    elif [[ $OutFreq = 'yearly' ]] ; then
+        # 3a. Concate the monthly files into yearly files
+        mkdir -p  ${outdir}/allyears
+        
+        year=$first_year
+        while [[ ${year} -le ${last_year} ]] ; do
+        echo Now merge year=$year for variable $var
+            cdo mergetime ${outdir}/${year}/${var}_${year}* ${outdir}/allyears/${var}_${year}${outsuffix}
+            year=`expr $year + 1`
+        done
 
+        # 3b.Add a last timestep in the end of the yearly files from next coming year
+            year=$first_year
+            while [[ ${year} -le ${last_year} ]] ; do
+              # Find next month and year
+              nextyear=`expr $year + 1`
+                
+              file=${outdir}/allyears/${var}_${year}${outsuffix}
+              nextfile=${outdir}/allyears/${var}_${nextyear}${outsuffix}
+              if [ -f ${nextfile} ] ; then
+                 echo First time step from file ${nextfile} will be added last in file ${file}
+                 cp -f $file file.nc
+                 cp -f $nextfile next.nc
+                 ncks -O -d TIME,0,0 next.nc tmp.nc
+                 ncrcat -O file.nc tmp.nc tmptot.nc
+                 mv tmptot.nc $file
+              else
+                 echo Next file ${nextfile} does not exist. Nothing will be done for ${file}
+              fi
+
+              # Update month and year
+              year=$nextyear
+          done
+    fi
     
 # End of loop variable list