Re: problem with my script

[Ansley Manke <ansley.b.manke@xxxxxxxx>]

Hi Du,
Even though the @VAR transformation removes the dependence on time, Ferret
retains the information about the underlying grid of each variable. Because
they have different time ranges, the variables VAR, VAR1, and VAR2 cannot be
added. To combine them, you need to explicitly put them on a grid with
an abstract
or compressed time axis. You can use the RESHAPE function to do this.
(The function
COMPRESSL would also work.)

DEFINE AXIS/t=1:1:1 time1
LET var1_0 = RESHAPE(var1, t[gt=time1])

And likewise for the other variables.

Ansley Manke

du hui wrote:

> Hi, the following problem confused me I can plot "var1/n1", "var2/n2"
> and "var/n" , but when I tried to plot "var1/n1+var2/n2" it did not
> work. Can anyone give me some hints?
> Thanks a lot. The followint is the script and the error message
> go com.ttest.u200.jnl
> USE mm5_lon180.ctrl.nc
> USE indwp_control.nc
> use daily_reanal
> let n1=792
> let n2=731
> let n=730
> let var1=UNEW[d=1,k=5,T=@var]
> let var2=U[d=2,k=5,T=@var]
> let var=u[d=3,k=10,T=@var]
> let t1=UNEW[d=1,k=5,l=1:792@ave]
> let t2=U[d=2,k=5,l=1:731@ave]
> let to=u[d=3,k=10,l=1:730@ave]
> let se1=(var1/n1+var/n)^0.5
> let tt1=(t1-to)/se1
> let se2=(var2/n2+var/n)^0.5
> let tt2=(t2-to)/se2
>
> FILL/X=47.5E:177.5E/Y=22.50S:22.50N (var1/n1+var/n)^0.5
> *** NOTE: Ambiguous coordinates on X axis: (VAR1/N1+VAR/N)^0.5
> **ERROR: illegal limits: VAR is not in the range T=31-DEC-4134
> 00:00:15-JAN-4187 00:00
> Axis extremes are T=31-DEC-1989 12:00:31-DEC-1991 12:00
> CONTOUR/FILL/X=47.5E:177.5E/Y=22.50S:22.50N (var1/n1+var/n)^0.5
> Command file, command group, or REPEAT execution aborted
>
> _________________________________________________________________
> 与联机的朋友进行交流，请使用 MSN Messenger: http://messenger.msn.com/cn