| About Triangulate [message #58579] |
Fri, 08 February 2008 00:44  |
qunjie0111
Messages: 9
Registered: September 2007
|
Junior Member |
|
|
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(fids,'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!
|
|
|
|
| Re: About Triangulate [message #58659 is a reply to message #58579] |
Fri, 08 February 2008 09:18  |
wgallery
Messages: 32
Registered: December 1998
|
Member |
|
|
On Feb 8, 3:44 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(fids,'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!
When executing the procedure:
Triangulate, x, y
the x and y parameters are rearranged (reordered.) Therefore, x and y
must be named variables (not expressions as you have). Try:
x = lon[0,index_theta]
y = lat[0,index_theta]
TRIANGULATE, x, y,
Tri,sphere=s,fvalue=f,/degrees
Good luck,
Bill Gallery
|
|
|
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
