[gradsusr] Cross section of temperature

[Geraldo Deniro]

Hi,
I wanted to create a script to plot a cross section of temperature of 100hPa until 1000hPa of the CFSR (Climate Forecast System Reanalysis) 4x daily but I'm not getting, when i created the script to plot the cross section of the zonal wind, the script worked  perfectly but in this case it did not work. I really don't know how to do it, the script that was used to try to create is below but GrADS says:
''Data request error: Invalid grid coordinatesvarying dimension 2 decreases: 27 to 1error occured getting variables 'tmpprs'error occured at column 1''
The script:
'reinit''open pgrbanl.cdas.ctl'
'set lat -26''set lon 270 340''set lev 100 1000'
temp=1while(temp<=12)'set t ' temp'q time'
data=subwrd(result,3)day=substr(data,4,2)month=substr(data,6,3)year=substr(data,9,4)time=substr(data,1,3)
'set display color white''set mpdset brmap_hires''set grads off''set grid off''set zlog on''set gxout shaded''run rgbset.gs''set clevs 5 10 15 20 25 30 35 40''set ccols 45 43 42 41 21 23 25 27 28''d tmpprs-273''set gxout contour''d tmpprs-273''run cbar''draw title CFSR: CROSS SECTION TEMPERATURE\ '%day' '%month' '%year' '%time'''draw ylab PRESSURE LEVELS (hPa)''draw xlab LATITUDE''set display color white''printim TEMP_'%month'_'%day'_'%year'_'%time'.jpg x1600 y1400''c'
temp=temp+1endwhile'quit'
If somebody could help me, i'd be so thankful

Geraldo Deniro contact: geraldodeniro at ufrj.br Graduate Student of Meteorology Department of Meteorology Institute of Geosciences Federal University of Rio de Janeiro 		 	   		  
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20140802/df5d9eb8/attachment.html