[gradsusr] Request is completely outside file limits

[Dave Allured]

Tereza,

That is a good point, thank you.  For files that contain only a 
single level, such as in Simba's example, that is the best way, 
because you leave it up to Grads to automatically select the only 
available level.  Level numbers do not need to be entered into the 
script at all.

For multi-level files, set lev or equivalent is mandatory, of course.

--Dave

On 8/5/2011 10:39 AM, Tereza Cavazos wrote:
> Possibly it is better not to define any level, just the time, do not
> define the winds, make the difference directly:
>
> 'set t 26353 26712'
> 'd uwnd.2-uwnd.1'
> 'draw title..."
>
> On 8/4/2011 9:53 PM, Dave Allured wrote:
>> These two lines are missing file number references:
>>
>> 'define u=ave(uwnd,t=26353,t=26712)'
>> 'define u1=ave(uwnd,t=26353,t=26712)'
>>
>> Both refer to the same /default/ file; see the top of the page on
>> Grads Variables in the Users Guide.  So one of the variable
>> references is getting the error because you are indexing a level
>> number that goes with the opposite file.
>>
>> Insert file numbers .1 and .2 :
>>
>> 'define u=ave(uwnd.1,t=26353,t=26712)'
>> 'define u1=ave(uwnd.2,t=26353,t=26712)'
>>
>> But keep the "set lev" commands on the same lines as before.
>>
>> Also, when requesting help, please identify the line where the error
>> happens.  If it is not already in the error message, insert "say"
>> commands to locate the line with the problem.
>>
>> --Dave
>> 1 9:
>> On 8/4/20137 PM, Chidzambwa Simba wrote:
>>> Dear Users,
>>> I am computing average zonal wind shear for Dec73-Feb74 from NetCDF
>>> files with DJF months only for all years. I am looking at shear
>>> between 850mb (first file) and 200mb (second file). My grads script
>>> is as follows;
>>> /'reinit'
>>> 'sdfopen uwndDJF_1900_2008.cdf'
>>> 'sdfopen uwndDJF_200_1900_2008.cdf'
>>> 'set background 1'
>>> 'set annot 0 6'
>>> 'set map 0 1 8'
>>> 'set xlopts 0 4'
>>> 'set ylopts 0 4'
>>> 'set lat -60 0'
>>> 'set lon 100 180'
>>> 'set gxout contour'
>>> 'set lev 850'
>>> 'define u=ave(uwnd,t=26353,t=26712)'
>>> 'set lev 200'
>>> 'define u1=ave(uwnd,t=26353,t=26712)'
>>> 'd u1-u'
>>> 'draw title 1973-74 Average zonal DJF wind shear (m/s)'
>>> 'printim ZW_Map_DJF7374.jpeg x400 y300'
>>> /
>>> and when I query my files this is what I get;
>>> /ga->   q file 1
>>> File 1 :
>>> Descriptor: uwndDJF_1900_2008.cdf
>>> Binary: uwndDJF_1900_2008.cdf
>>> Type = Gridded
>>> Xsize = 41 Ysize = 31 Zsize = 1 Tsize = 39096 Esize = 1
>>> Number of Variables = 1
>>> uwnd 1 t,z,y,x 6-hourly U-wind on Pressure Levels
>>> ga->   q file 2
>>> File 2 :
>>> Descriptor: uwndDJF_200_1900_2008.cdf
>>> Binary: uwndDJF_200_1900_2008.cdf
>>> Type = Gridded
>>> Xsize = 41 Ysize = 31 Zsize = 1 Tsize = 39096 Esize = 1
>>> Number of Variables = 1
>>> uwnd 1 t,z,y,x 6-hourly U-wind on Pressure Levels/
>>> When I run my script I get an error message /" Data Request Warning:
>>> Request is completely outside file limits"/ and I don't understand
>>> why. The areal limits are correct and I have cross-checked them and
>>> the time-steps are 6-hourly daily wind components. May somebody
>>> assist me please?
>>> Regards
>>> Simba
>> _______________________________________________
>> gradsusr mailing list
>> gradsusr at gradsusr.org
>> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr