Generating variable with several layers

[Erko Jakobson]

Hello

Thank you very much Heinz Juergen, I didn't know the function const, now the
problem is solved and I can go on with my calculations.


Erko



On Tue, 4 Sep 2007 18:45:14 +0200, Heinz Juergen Punge
<heinzjuergen.punge at ZMAW.DE> wrote:

>Hi Erko,
>are the variables Q_lev'z' defined on the other levels? Try this
>(assuming you have a ctl file with the 27 levels open):
>zi=1
>'set z 1 27'
>while zi<=27
>'QQ_lev'zi'=const(Q_lev'zi',0,-u)'
>  ...
>  zi=zi+1
>endwhile
>
>Maybe defval only works for lat/lon dimensions?
>You may alternatively want to try something like this :
>zi=1
>'set z 1 27'
>'QQ=const(Q_lev'zi',0)
>*i.e. a new variable that is zero everywhere
>while zi<=27
>  'set z 'zi
>  'q dims'
>  levi=subwrd(sublin(result,4),6)
>  'set z 1 27'
>  'QQ=QQ+const(maskout(Q_lev'zi',-abs(lev-'levi')),0,-u)'
>  zi=zi+1
>endwhile
>
>Best,
>Heinz Juergen
>
>
>Erko Jakobson wrote:
>> Hello again
>>
>> One solution I tried was function set defval. As it works only 2D, then I
>> fixed lon value and made a new 2D variable QQ, the variables are lat and z.
>> 'set defval QQ 3 10 2.46'
>> (here z=3, lat=10, new value for QQ is 2.46)
>> But it gives error "MODIFY DEFVAL Error:  Out of Range", and it gives this
>> error always when z<>1, with z=1 the function is working correctly.
>>
>> Hoping to get some new ideas to go on
>>
>> Erko Jakobson,
>> PhD student,
>> University of Tartu,
>> Estonia
>>
>>
>>
>>
>> On Mon, 3 Sep 2007 13:27:37 +0200, Erko Jakobson <erkokkk at HOT.EE> wrote:
>>
>>
>>>I have Q seasonal averages for 27 level separately (Q_lev1 to Q_lev27). For
>>>making the vertical profiles, I have to put them in only one variable Q,
>>>that consists 27 levels, so I can display the profiles:
>>>
>>>'set lat 55 90'
>>>'set lon 90'
>>>'set z 1 27'
>>>'d Q'
>>>
>>>It should be function like Q(z)=Q_lev'z', but for me it seems, that in Grads
>>>it is not so easy:(
>>>
>>>
>>>Thank you very much for any advice,
>>>
>>>Erko Jakobson,
>>>PhD student,
>>>University of Tartu,
>>>Estonia
>>
>>