Re: Problem about TRIANGULATE [message #58575] |
Fri, 08 February 2008 02:05  |
Spon
Messages: 178 Registered: September 2007
|
Senior Member |
|
|
On Feb 8, 8:49 am, qunjie0...@163.com wrote:
> HALLO:
> I wanted to interpolate the irregularly data (satellite data) to
> regularly_gridded data,and used TRIANGULATE and TRIGRID function ,but
> after run the .pro ,there was a problem:
>
> Attempt to store into an expression: <FLOAT Array[912]>.
> Execution halted at: HALOE_TRIGRID 176 E:\haloe\haloe_trigrid.pro
> $MAIN$
> MY CODE:
> fidr=ncdf_open('h:\haloe\HALOE_Version19_sr_2004_03.nc')
> varid=ncdf_varid(fidr,'latitude')
> ncdf_varget,fidr,varid,lat
>
> varid=ncdf_varid(fidr,'longitude')
> ncdf_varget,fidr,varid,lon
>
> varid= ncdf_varid(fidr,'theta')
> ncdf_varget,fidr,varid,theta
> ncdf_close,fidr
>
> index_theta=where(theta[0,*] ne 1.00000e+024)
> TRIANGULATE, lon[0,index_theta], lat[0,index_theta],
> Tri,sphere=s,fvalue=f,/degrees
> grid_data = trigrid(f,[2.5,2.5],[0,-80,357.5,80],sphere=s,/degrees)
> .
> .
> .
> What's the problem? Please help me!
> thanks!
Have a look at:
http://groups.google.com/group/comp.lang.idl-pvwave/browse_t hread/thread/1ca56b9c2bcad51/58a5fe20a6d16341?lnk=gst&q= Attempt+to+store+into+an+expression+triangulate#58a5fe20a6d1 6341
It boils down to having to define the x, y and fvalue variables before
passing them to triangulate so they can be sorted properly.
Take care,
Chris
|
|
|