source: trunk/tools/netCDF_scripts/make_1mm @ 43

#!/bin/sh
# -*- Mode: sh -*-
#
#  Script to build 1mm and annual mean from time serie
#
#  Uses nco tools
#
usage() {
echo
echo "usage : make_1mm [-d] [-g] [-v <field1>,<field2>...] [-t <record name>] [-g <grid>] [-m] <file> <nyears>"
echo
}
#
#
# --------------------------------------------------------------
#  E. Guilyardi - LODYC -  July 1999 - Nov 2007
#
# 
# Inits
#

while getopts g:dhv:t:m name
do
    case $name in
        d) set -x ;;
        g) type="$OPTARG";;
        v) fields="$OPTARG" ;;
        t) timef="$OPTARG" ;;
        m) nomean=1;;
        h) usage 
           echo
           exit 1
           ;;
    esac
done

#
# so that $1 is root
#
shift `expr $OPTIND - 1`
file=$1

enam=`echo $file | awk -F\_ '{print $1}'`
ave=`echo $file | awk -F\_ '{print $2}'`
date1=`echo $file | awk -F\_ '{print $3}'`
year1=`echo $date1 | sed 's/\(.*\)../\1/'`
date2=`echo $file | awk -F\_ '{print $4}'`
year2=`echo $date2 | sed 's/\(.*\)../\1/'`
ndates=$2

root=${enam}_${ave}_${date1}_${date2}

type=`echo $file | awk -F\_ '{print $5}'`
type2=`echo $file | awk -F\_ '{print $6}'`
type3=`echo $file | awk -F\_ '{print $7}'`
type4=`echo $file | awk -F\_ '{print $8}'`
type=${type}_${type2}_${type3}_${type4}
type=`echo $type | sed 's/\(.*\)\.nc.*/\1/'`


type=${type:-"auto"}
nomean=${nomean:-0}

echo ""
echo "Processing: $enam $ave $year1 $year2 $ndates $type"

case $type in
    t30) long="nav_lon" 
            lati="nav_lat"  
           time="time_counter" ;;
    lmdz*) long="lon" 
            lati="lat"  
            time="time_counter" ;;
    reg1) long="lon" 
            lati="lat"  
            time="time" ;;
    t42) long="lon" 
            lati="lat"  
            time="time" ;;
    t62) long="lon" 
           lati="lat"  
           time="time" ;;
    umat) long="longitude" 
            lati="latitude"  
            time="t" ;;
    umau) long="longitude" 
            lati="latitude"  
            time="t" ;;
    reg1mtrop) long="longitude" 
            lati="latitude"  
            time="t" ;;
    *) long="nav_lon" 
       lati="nav_lat"  
       time=`ncdump -h ${file} | grep UNLIMITED | awk '{print $1}'` ;;
esac 
time=${timef:-$time}
fields=${fields:-"-x -v $long,$lati"}
          
#
case $ave in
    1m) # output file name
        month1=`echo $date1 | sed 's/.*\(..\)/\1/'`
        month2=`echo $date2 | sed 's/.*\(..\)/\1/'`

        aveoum=1mm
        ;;

    *) echo " ERROR : case ave=$ave not relevant "
       exit 1
       ;;
esac

#
PATH=$PATH:/usr/local/netcdf/bin:/u/rech/eee/reee217/bin
#
out1m=${enam}_${aveoum}_01_12_${year1}-${year2}_$type.nc
outmean=${enam}_${ndates}y_${year1}_${year2}_$type.nc
#
optg="-O"
grid=0
if [ $grid -eq 1 ] ; then
# get grid
    ncks -O -v $lati,$long ${file} $out1m || echo 'no grid'
    optg="-A"
fi

m_list="00 01 02 03 04 05 06 07 08 09 10 11"

# average every month

for m in ${m_list}
    do
        m1=`expr $m + 1`
#       echo "   averaging month $m1 "
        ncks -O -d $time,$m,,12 ${file} ${root}_${type}_mean_month.nc || exit
        ncra -O -F -d $time,1,$ndates ${root}_${type}_mean_month.nc ${root}_${type}_mean_month_${m}.nc || exit
        rm ${root}_${type}_mean_month.nc
done
    
# build monthly means
    
echo " Building ${enam}_${aveoum}_01_12_${year1}-${year2}_$type.nc"
    
ncrcat ${optg} ${root}_${type}_mean_month_??.nc $out1m

if [ $nomean -eq 0 ] ; then 

    echo " Building ${outmean}"
    
    ncra ${optg} $out1m $outmean

fi

rm -f ${root}_${type}_mean_*.nc