[NCEP.List.Nomads-Users] problem subset grib2 files

[Eric Altshuler]

Enzo,

One possible problem is that the longitude range of the subsetted grib2 data is offset by +/- 360 degrees from the value you input to your script as 'lon1'. For example, the grib2 data may have a longitude range of 270 to 320 instead of -90 to -40. The physical location on the earth is the same in both cases but grads will not recognize the 360 degree offset. If you input 'lon1' as -70, grads will return undefined value. Try using 290 (which equals -70 + 360) instead of -70 and see if that works.

Best regards,

Eric L. Altshuler
Assistant Research Scientist
Center for Ocean-Land-Atmosphere Studies
4041 Powder Mill Road, Suite 302
Calverton, MD 20705-3106
USA

E-mail: ela at cola.iges.org
Phone: (301) 902-1257
Fax: (301) 595-9793

----- Original Message -----
From: "Enzo Campetella" <efcam at speedy.com.ar>
To: "ncep list nomads-users" <ncep.list.nomads-users at noaa.gov>
Cc: "Jun Wang" <Jun.Wang at noaa.gov>, "Jordan Alpert" <Jordan.Alpert at noaa.gov>, "Wesley Ebisuzaki" <Wesley.Ebisuzaki at noaa.gov>
Sent: Friday, July 3, 2009 2:23:44 AM GMT -05:00 US/Canada Eastern
Subject: [NCEP.List.Nomads-Users] problem subset grib2 files

Hello!!!

My query is about problem using subset files in grib2 format / opendap.

I use Grads 2.0a5 client. Newer one doesnt work in Windows System.

I get problems when I want to make a meteogram or a graphics like script patchs at the end of  my mail.

When I work with global grib2 files, script run perfect and I can obtain the graphic Temperature vs. time.

But when I work using a subseting grib2 file, script doesnt run. In this case I get this legend: Cannot Plot Data - all undefined values / Data Request Warning: Request beyond file limits.

I tryed to change sentence set lat / set lon but I cant solve my problem. Why it run well using global files, and doesnt work using subset files?

Normally I work using this area: -90  - 40 lon / -60 -20 lat

Thanks, Enzo Campetella

Script eg:

'reinit'
'open c:\gblav\grillachica\gblav.ctl'

'q file 1'
res=sublin(result,5)
t2=subwrd(res,12)
'set t 1 't2

prompt 'Ingrese la latitud: '
pull lat1
prompt 'Ingrese la longitud: '
pull lon1

lati=lat1-0.1
lats=lat1+0.1
'set lat ' lat1
'set lon ' lon1
'set grads off'

'set parea 0.5 10.5 3 8'
'set gxout line'
'set grid on'
*'set xlab off'
'set ylab on'

'set cthick 10'
'd tmp2m-273.16'

'set cthick 2'
'd tmax2m-273.16'
'd tmin2m-273.16'

'set parea 0.5 10.5 1 3'

'set lat ' lati ' ' lats
'set gxout barb'
'set ccolor 2'
'set xyrev on'
'set ylab off'
'set grid vertical'
'd ugrd10m*1.92;vgrd10m*1.92'

say 'Pulse una tecla para ver humedad'
pull tecla

'set parea off'
'c'
'set grid on'
'set xlab on'
'set ylab on'
'set grads off'
'set lat ' lat1
'set lon ' lon1

'set cthick 10'
'd rh2m'




_______________________________________________
NCEP.List.Nomads-Users mailing list
NCEP.List.Nomads-Users at noaa.gov
https://lstsrv.ncep.noaa.gov/mailman/listinfo/ncep.list.nomads-users