[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