<HTML>
<HEAD>
<META content="text/html; charset=iso-8859-1" http-equiv=Content-Type>
<META content="OPENWEBMAIL" name=GENERATOR>
</HEAD>
<BODY bgColor=#ffffff>
Dear Theresa,
<br />
<br />try this piece of script below:
<br />
<br /><font size="2"> 'reinit'
<br /> 'open E:\test.ctl'
<br /> 'set t 5'
<br /> 'set lat 7'
<br /> 'set lon
107'
<br />* 'define r1= sst' !NOT NECESSARY, since you have already set a single point.
<br /> 'd sst'
<br />r1=subwrd(result,4)
<br /><font size="2">say r1
<br />'close 1'</font></font>
<br /><font size="2"> 'sdfopen sst.nc'
<br /> 'set t 5'
<br /> 'set lat 7'
<br /> 'set lon 107'
<br /> *'define r2=
sst' !NOT NECESSARY
<br />'d sst'
<br /> <font size="2">r2=subwrd(result,4)
<br /><font size="2">say r2
<br />
<br /></font></font></font><font size="2">'define aver= ('r1'+'r2')/2'
<br />
<br />'d aver'
<br /> Ar=subwrd(result,4)
<br /> say Ar</font>
<br />
<br />NOTE that, now, r1 and r2 are "script variables" (not "grads prompt variables"), since they were "captured" with "subwrd", and not with "define". So, r1 and r2 have to appear between quotes in the Ar definition. Note also that you can not use this kind of construction in the grads prompt terminal:
<br />
<br /><font size="2">'d sst'
<br /> <font size="2">r2=subwrd(result,4)
<br />
<br /></font></font>
<br />So, script variables CAN NOT be directly displayed in the grads graphic window: they exist only in the script memory.
<br />
<br />On the other hand, DEFINED variables exist in the "GrADS prompt memory" and in the "script memory". These defined variables can be directly displayed in the grads graphic window. Of course, to make computations with defined variables, you don't need to put quotes in the variables. Is it clear enough?
<br />
<br />Please, let me know if this doesn't work as you expect.
<br />
<br />Regards,
<br />Ricardo
<br />
<br /><font size="2"><b>On Tue, 5 Feb 2013 18:10:01 +0000, T. haghroosta wrote</b>
<br />>
Dear Ricardo,
<br />>
<br />> The script that you helped me to write, was
run correctly. Now I would like to calculate the average of one parameter in two
different files. I wrote this script, but it doesn't show me any outputs. Could
you please tell me what is the problem?
<br />>
<br />> 'reinit'
<br />>
'open E:\test.ctl'
<br />> 'set t 5'
<br />> 'set lat 7'
<br />> 'set lon
107'
<br />> 'define r1= sst'
<br />> 'close 1'
<br />> 'sdfopen sst.nc'
<br />>
'set t 5'
<br />> 'set lat 7'
<br />> 'set lon 107'
<br />> 'define r2=
sst'
<br />> 'define aver= (r1+r2)/2'
<br />> 'd aver'
<br />>
Ar=subwrd(result,4)
<br />> say Ar
<br />>
<br />> Regards,
<br />>
<br />> T.
<br />>
<br />> >
<br />> >
<br />> > Today's
Topics:
<br />> >
<br />> > 1. Re: Average in two different time and
locations (T. haghroosta)
<br />> >
<br />> >
<br />> >
----------------------------------------------------------------------
<br />>
>
<br />
<br!>
<br />
<br />
</br!></font>
</BODY>
</HTML>