gfs plot

teresa tiberi crancyfera at GMAIL.COM
Thu May 22 12:48:34 EDT 2008


Hi,

I' trying to plot the gfs data from
URL="http://www.ftp.ncep.noaa.gov/data/nccf/com/gfs/prod/"

by this procedure:

 ${DIR}scriptgfs/wesley/get_inv.pl $URLFILE.idx | egrep "${WHAT}"\
   |${DIR}scriptgfs/wesley/get_grib.pl $URLFILE ${RUN_DIR}${HHMM}${FILE}

I put them together by

 /usr/local/bin/wgrib -d all $fileout -grib -append -o ${YYYYMMDDHH}.grb
and since +192,I made a second grb because resolution change from 1
degree to 2.5
 /usr/local/bin/wgrib -d all $fileout -grib -append -o ${YYYYMMDDHH}.192.grb

and after this I run
/home/grib/bin/grib2ctl.pl -verf ${YYYYMMDDHH}.grb > ${YYYYMMDDHH}.ctl
/usr/local/bin/gribmap -i ${YYYYMMDDHH}.ctl
/home/grib/bin/grib2ctl.pl -verf ${YYYYMMDDHH}.192.grb > ${YYYYMMDDHH}.ctl
/usr/local/bin/gribmap -i ${YYYYMMDDHH}.192.ctl

so I'm ready to open with grads, but they are rotated data and when I
ask a single value it is not rigth value,but it's "from other lat
lon".

could someone help me?
thanks

Teresa



More information about the gradsusr mailing list