Changeset 44 for trunk/tools

Timestamp:

01/31/08 07:59:54 (16 years ago)

Author:

ericg

Message:

Clean up of make_1mm

File:

• trunk/tools/netCDF_scripts/make_1mm (modified) (6 diffs)

trunk/tools/netCDF_scripts/make_1mm

@@ -2 +2 @@
 # -*- Mode: sh -*-
 #
-#  Script to build 1mm and annual mean from time serie
+#  Script to build 1mm and full record annual mean from monthly time serie
 #
 #  Uses nco tools
@@ -8 +8 @@
 usage() {
 echo
-echo "usage : make_1mm [-d] [-g] [-v <field1>,<field2>...] [-t <record name>] [-g <grid>] [-m] <file> <nyears>"
+echo "usage : make_1mm [-d] [-v <field1>,<field2>...] [-t <record name>] [-m] <file> <nyears>"
+echo "        -d : toggle debug mode"
+echo "        -v : specify fields to work on"
+echo "        -t : specify time record name (should be detected automatically)"
+echo "        -m : do not make full record time mean"
+echo "        <file>   : netcdf post-it format file"
+echo "        <nyears> : number of years in file"
 echo
 }
@@ -20 +26 @@
 #
 
-while getopts g:dhv:t:m name
+while getopts dhv:t:m name
 do
     case $name in
         d) set -x ;;
-        g) type="$OPTARG";;
         v) fields="$OPTARG" ;;
         t) timef="$OPTARG" ;;
@@ -36 +41 @@
 
 #
-# so that $1 is root
+# decode file name
 #
 shift `expr $OPTIND - 1`
@@ -65 +70 @@
 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 
+# detect time dimension
+
+time=`ncdump -h ${file} | grep UNLIMITED | awk '{print $1}'`
 time=${timef:-$time}
 fields=${fields:-"-x -v $long,$lati"}
@@ -112 +91 @@
 
 #
-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"