[gradsusr] Request is completely outside file limits

[Tereza Cavazos]

Hi Dave,

Yes, that is a good point. Thanks !

Tereza

On 8/5/2011 11:32 AM, Dave Allured wrote:
> 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
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr