[gradsusr] Problems converting grib2ctl
tschonholz
tschonholz at smn.gov.ar
Mon Jul 2 13:44:41 EDT 2018
Hi everyone! I´m posting this message because I´m finding a hard time
using grib2ctl
I use this bash to generate my ctl and idx on a server, but the original
ctl and idx are generated on different server than mine.
#!/bin/sh
#*****Script para procesar el analisis
today=$(date +%Y%m%d)
echo ${today}
#chequeo que exista la salida del WRF
if [ -d /nas-r5c/WRF/${today}_00/ ];
then
#si ocurre, genero el ctl
echo ***Proceso las salidas grib del modelo**
cd /home/mgatto/
#rm /home/mgatto/modelos/WRF/data/gribwrf.ctl
perl grib2ctl.pl -verf /nas-r5c/WRF/${today}_00/UPPOUT/wrfprs_d01.%f3
/home/mgatto/modelos/WRF/data/indexwrf00.idx
>/home/mgatto/modelos/WRF/data/gribwrf00.ctl
gribmap -i /home/mgatto/modelos/WRF/data/gribwrf00.ctl
else
exit
fi
--
doing this, i get this ctl:
dset /nas-r5c/WRF/20180702_00/UPPOUT/wrfprs_d01.%f3
index /home/mgatto/modelos/WRF/data/indexwrf00.idx
undef 9.999E+20
title /nas-r5c/WRF/20180702_00/UPPOUT/wrfprs_d01.000
* produced by grib2ctl v0.9.13
* command line options: -verf
/nas-r5c/WRF/20180702_00/UPPOUT/wrfprs_d01.%f3
/home/mgatto/modelos/WRF/data/indexwrf00.idx
dtype grib 255
options template
pdef 999 1249 lccr -54.387000 -94.331000 1 1 -35.000000 -35.000000
-65.000000 4000 4000
xdef 999 linear -94.331000 0.0624474869364198
ydef 1249 linear -54.387000 0.0363636363636364
tdef 49 linear 00Z02jul2018 1hr
* z has 28 levels, for prs
zdef 28 levels
1000 975 950 925 900 875 850 825 800 775 750 700 650 600 550 500 450 400
350 300 275 250 225 200 175 150 100 50
vars 130
and all the variables.
The problem is that the original ctl says this:
dset ^wrfprs_d01.%f3
index ^wrfprs_d01.idx
options template
undef 9.999E+20
title wrfprs_d01.00
* produced by grib2ctl v0.9.12.5p45
dtype grib 255
pdef 999 1249 lcc -54.387000 -94.331000 1 1 -35.000000 -35.000000
-65.000000 4000 4000
xdef 999 linear -94.331000 0.0624474869364198
ydef 1249 linear -56.853000 0.0363636363636364
tdef 49 linear 00Z02Jul2018 1hr
* z has 28 levels, for prs
zdef 28 levels
1000 975 950 925 900 875 850 825 800 775 750 700 650 600 550 500 450 400
350 300 275 250 225 200 175 150 100 50
vars 130
which means that I lose the south of my area because I don´t know why
it's not doing the same than the original.
I believe it has to do with my grib2ctl version, does anyone know what
happens or could help me please?
Thanks a lot!
Tamara Schonholz- Servicio Meteorologico Nacional-Argentina
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://gradsusr.org/pipermail/gradsusr/attachments/20180702/10bdb2e8/attachment.html>
More information about the gradsusr
mailing list