#!/bin/sh

set -ax 

#first create list of ensemble members *
#run grmean with the list
grmean -s ensemble.sa.list -o sa.tmphag.mean.200803.monthly.grib 
#fix the date
force_grib_date_mon sa.tmphag.mean.200803.monthly.grib 2008030100
#run "grmap -g0" to create ctl, fix "0 to 2" in tdef of ctl, run "gribmap -0 -i"


#you can use grvar_ver2 to create standard deviation file from the list of grib files

grvar_ver2 -s ensemble.sa.list -o sa.tmphag.stdev.200803.monthly.grib 
force_grib_date_mon sa.tmphag.stdev.200803.monthly.grib 2008030100
