gxout print

[Carlos Minjarez]

Thank you I figured out the problem, it was actually taking the values of
the grids around my interest region.

Thank you any way
carlos
> Why your X and Y ranges are not integers? It is my understanding that X
> and Y refer to grid points, so if you have a script that says something
> like 1300.5 you are referring to an interpolated value.  At least that
> is how I understand it.  I managed  to extract wind speeds simply  by
> referring to  grid points in  nested while-loops, more or less like you.
>
> Tuomo Lauri
>
> Carlos Minjarez wrote:
>> Has any body try to print the actual values in an ascii file.
>> I am trying but it is printing the interpolated values and I am just
>> interested in the natural grid value.
>> I will appreciate your help
>>
>> carlos
>>
>> 'reinit'
>> datadir = 'textdata'
>>
>> 'open ST4.ctl'
>> count = 1633
>> while(count <=1657)
>> 'set t ' count
>> *'set lon -111.531 -109.817'
>> 'set X -1300.5 -1298.5'
>> *'set lat 30.8284 32.4566'
>> 'set Y 90.5 92.5'
>> 'q dims'
>> val1 = sublin(result,5)
>>
>> val2 = subwrd(val1,6)
>> 'set gxout print'
>>
>> 'set prnopts %8.0f 5 1'
>>
>> *'set prnopts %10.3e 1 1'
>> 'd precip'
>> file = "ST4data"count"-"%val2%".txt"
>>
>> result1 = write(file, result)
>>
>> say count
>> count = count+1
>>
>>
>> endwhile
>>
>