averagefield_ensemble.cgi

#!/bin/bash
. ./init.cgi
# average an ensemble of fields
echo "Content-type: text/html"
echo

export DIR=`pwd`
. ./getargs.cgi

# no ensemble member selection yet

. ./queryfield.cgi
. ./myvinkhead.cgi "Ensemble mean" "$kindname $climfield" "nofollow,index"

ensfile=$file
file=`echo $file | sed -e 's/%%%/ave/' -e 's/%%/ave/' -e 's/\+\+\+/ave/' -e 's/\+\+/ave/'`
if [ $ensfile = $file ]; then
    echo "Error: input file is not an ensemble."
    . ./myvinkfoot.cgi
    exit
fi
file=data/`basename $file`

if [ -s $file ]; then
  echo "Using cached data<p>"
else
  echo "Computing ensemble mean ...<p>"
  echo "bin/averagefield_ensemble $ensfile mean $file" >> log/log
  ### echo bin/averagefield_ensemble $ensfile mean $file
  tmpfile=data/aap$$.nc
  ( bin/averagefield_ensemble $ensfile mean $tmpfile ) 2>&1
  mv $tmpfile $file
  if [ ! -s $file ]; then
    echo "Something went wrong in the averaging routine."
    echo "Most likely, the fields are too large. Please first select the region you are interested in and then average."
    . ./myvinkfoot.cgi
    exit
  fi
fi

eval `bin/getunits.sh $file`
ENSEMBLE=""
kindname="mean $kindname"
FORM_field=data/mean_`basename $FORM_field .info`.$EMAIL.info
cat > $FORM_field <<EOF
$file
NPERYEAR=$NPERYEAR
UNITS=$UNITS
$kindname
$climfield
EOF

. ./select.cgi