If you can replace sph_scat with griddata, you should give that a try.
Haje
"Jeff Patrick" <jeff.patrick@computer.org> wrote in message
news:b368fde2.0406180703.2a86ddb2@posting.google.com...
> Hello all,
>
> The section of code that follows appears to work fine right up until
> it gets to the contour command where the program hangs. No errors are
> given and IDL must be shutdown from task manager and restarted. The
> only thing I've changed from some working code I have is the SPH_SCAT
> commands used to be calls to a custom krigging function. Any ideas
> what may be causing this to happen?
>
> Any suggestions are welcome,
>
> Jeff
>
> arrX = SPH_SCAT(arrLon, arrLat, arrXValue, BOUNDS=[1.,0.,360., 89.],
> GS=[2.,1.], BOUT=bout)
>
> arrY = SPH_SCAT(arrLon, arrLat, arrYValue, BOUNDS=[1.,0.,360., 89.],
> GS=[2.,1.], BOUT=bout)
>
> FOR j = 0.0, 89.0, 1.0 DO BEGIN
>
> FOR i = 0.0, 360.0, 2.0 DO BEGIN
>
> pro_ConvMagGeo, 89.0 - j, i - dblShift + 1, geo_lat, geo_lon
>
> IF j EQ 0.0 AND i EQ 0.0 THEN BEGIN
>
> arrGLat = [geo_lat]
>
> arrGLon = [geo_lon]
>
> ENDIF ELSE BEGIN
>
> arrGLat = [arrGLat, geo_lat]
>
> arrGLon = [arrGLon, geo_lon]
>
> ENDELSE
>
> ENDFOR
>
> ENDFOR
>
> sngH = SQRT(((arrX)^2) + ((arrY)^2))
>
> MAP_SET, intLat, intLon, 0, /SATELLITE, SCALE=(intZoom * 1E6),
> /ISOTROPIC, /HORIZON, E_HORIZON={FILL:0, COLOR:'000000'xL},
> XMARGIN=[0,17], YMARGIN=[0,0], /NOERASE
>
> hexTextCol = 'CCCCCC'xL ; gray text
>
> CONTOUR, sngH, arrGLon, arrGLat, /IRREGULAR, /ISOTROPIC, /CLOSED,
> /CELL_FILL, /OVERPLOT,LEVELS=[ 20, 40, 60, 80, 100, 120, 140, 160,
> 180, 200, 220, 240, 260, 280, 300, 320, 340, 360, 380, 400, 420, 440,
> 460, 480, 500, 520, 540, 560, 580, 600, 620, 640, 660, 680, 700, 720,
> 740, 760, 780, 800, 820, 840, 860, 880, 900, 920, 940, 960, 980, 1000]
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
