[ferret_users] RE: regrid problem

[Zibiao Zhang <zzhang@xxxxxxx>]

Thank you Jaison, that solved the problem. Key point is paying attention
when regrid to a different dimension system, and using g=grd attributes.

-----Original Message-----
From: Jaison Kurian [mailto:jaison@xxxxxxxxxxxxxxxxxx]
Sent: Tuesday, July 17, 2007 2:21 PM
To: Zibiao Zhang
Subject: RE: regrid problem



  The point to note is that your original grid contains X, Y and
Z axes, but the newly defined "grd" contains only X and Y axes.

  So either you have to do like

     define axis/x=.......  xax
     define axis/y=.......  yax
     define grid/like=temp[d=1]/x=xax/y=yax grd
     ! above line is same as
     ! define grid/x=xax/y=yax/z=temp[d=1] grd
     let temperature = temp[d=1,g=grd]

  Or
     define axis/x=.......  xax
     define axis/y=.......  yax
     let temperature = temp[d=1,gx=xax,gy=yax]


If this fix/es worked for you, please do send a summary (a fresh
mail, discussing relevant points) to the mailing list.

Regards,

Jaison




On Tue, 17 Jul 2007, Zibiao Zhang wrote:

> Hi Jaison,
>
> It worked! thanks for your great help, I did not know there is a
difference
> between the two:
> let temperature=temp[g=grd]
> and
> let temperature = temp[d=1,gx=xax,gy=yax] ! if 3d_ave.nc is d=1
>
> thanks
> Zibiao
> -----Original Message-----
> From: Jaison Kurian [mailto:jaison@xxxxxxxxxxxxxxxxxx]
> Sent: Tuesday, July 17, 2007 1:38 PM
> To: Zibiao Zhang
> Subject: RE: regrid problem
>
>
>
>  After the script fails out, just use a "show data" command to make
> sure that either 3d_ave.nc is the only file opened or it is the
> default dataset. Also make sure that the variable temp exists in
> 3d_ave.nc.
>
>  If you pass above check, and still it fails, then instead of
>
>
>
> try
>      l
>      save/file=test.nc/append temperature
>
> Regards,
>
> Jaison
>
>
>
> On Tue, 17 Jul 2007, Zibiao Zhang wrote:
>
> > Hi Jaison,
> > I tried this already, they were the same.
> > Zibiao
> >
> > -----Original Message-----
> > From: Jaison Kurian [mailto:jaison@xxxxxxxxxxxxxxxxxx]
> > Sent: Tuesday, July 17, 2007 1:23 PM
> > To: Zibiao Zhang
> > Subject: Re: regrid problem
> >
> >
> > Hi Zibiao,
> >             Instead of
> >
> >     list/clobber/file=test.nc/format="CDF" temperature
> >
> > first remove test.nc manually, then use
> >
> >     save/file=test.nc/append temperature
> >
> >
> > please Let me know whether this works or not.....
> >
> > Jaison
> >
> >
> >
> >
> >
> > On Tue, 17 Jul 2007, Zibiao Zhang wrote:
> >
> > > Hi Jaison,
> > > I have a data file I called 3d_ave.nc has the grid as following:
> > >
> > > Grid GGR1
> > > Longitude200_450 	longitud 	251mi	82.08W(-82.08)	62.08W(-62.08)
> > > Latitude800_1050	Latitude		251   i	32.838N	47.96N
> > > Depth_1-		Depth(m)	10     i-	0		200
> > > normal 			T
> > >
> > > I'd like to resample the data point to 81x81, instead of the original
> > > 251x251.
> > > I used the following scripts:
> > >
> > > can data/all;
> > > use 3d_ave
> > > def axis/x=82.08w:62.08w:0.25 xax;
> > > def axis/y=32.83n:47.96n:0.19 yax;
> > > def grid/x=xax/y=yax grd
> > > let temperature=temp[g=grd]	! temp is the varialbe in file 3d_ave
> > > list temperature
> > > list/clobber/file=test.nc/format="CDF" temperature
> > >
> > > the result from "list temperature" looks fine, however when I save to
> the
> > > file test.nc, and list the temperature from the file, it has error:
> > >
> > > ** ERROR: variable unknown or not in data set: TEMP[G=grd]
> > >
> > > Can you tell me what is the problem?
> > >
> > > thank you very much,
> > > Zibiao
> > >
> > >
> > >
> >
> >
>
> --
> This message has been scanned for viruses and
> dangerous content by MailScanner, and is
> believed to be clean.
>
>
>

--
This message has been scanned for viruses and
dangerous content by MailScanner, and is
believed to be clean.