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
