Changeset 236 for trunk/SRC/Interpolation/inquad.pro
- Timestamp:
- 03/21/07 18:18:44 (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SRC/Interpolation/inquad.pro
r232 r236 261 261 ; points out of the cells 262 262 IF (where(chtest EQ 0))[0] NE -1 THEN BEGIN 263 IF NOT keyword_set(noprint) THEN print, 'Points number '+strjoin(strtrim(where(chtest EQ 0), 1), ', ')+' are out of the grid' 263 IF NOT keyword_set(noprint) THEN BEGIN 264 ras = report('Points number '+strjoin(strtrim(where(chtest EQ 0), 1), ', ')+' are out of the grid') 265 ENDIF 264 266 stop 265 267 ENDIF 266 268 ; points in more than one cell 267 269 IF (where(chtest GT 1))[0] NE -1 THEN BEGIN 268 IF NOT keyword_set(noprint) THEN print, 'Points number '+strjoin(strtrim(where(chtest GT 1), 1), ', ')+' are in more than one cell' 270 IF NOT keyword_set(noprint) THEN BEGIN 271 ras = report('Points number '+strjoin(strtrim(where(chtest GT 1), 1), ', ')+' are in more than one cell') 272 ENDIF 269 273 stop 270 274 ENDIF … … 277 281 CASE 1 OF 278 282 found[0] EQ -1:BEGIN 279 IF NOT keyword_set(noprint) THEN print, 'The point is out of the cells' 283 IF NOT keyword_set(noprint) THEN BEGIN 284 ras = report('The point is out of the cells') 285 ENDIF 280 286 IF keyword_set(onsphere) THEN BEGIN 281 287 ; restore old !map... … … 289 295 END 290 296 n_elements(found) GT ntofind:BEGIN 291 IF NOT keyword_set(noprint) THEN print, 'The point is in more than one cell' 297 IF NOT keyword_set(noprint) THEN BEGIN 298 ras = report('The point is in more than one cell') 299 ENDIF 292 300 END 293 301 ELSE:
Note: See TracChangeset
for help on using the changeset viewer.