Changeset 240 for trunk/SRC/ToBeReviewed/GRILLE
- Timestamp:
- 03/28/07 14:17:24 (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SRC/ToBeReviewed/GRILLE/domdef.pro
r232 r240 248 248 IF (dom[0] EQ -1) THEN BEGIN 249 249 IF keyword_set(findalways) THEN BEGIN 250 print, 'WARNING, empty T points box... we get the neighnors to define a new box...'250 ras = report('WARNING, empty T points box... we get the neighnors to define a new box...') 251 251 neig1 = neighbor(lon1, lat1, glamt, gphit, sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) 252 252 neig2 = neighbor(lon2, lat2, glamt, gphit, sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) … … 284 284 CASE 1 OF 285 285 (where(gridtype eq 'T'))[0] NE -1 OR (where(gridtype EQ 'W'))[0] NE -1:BEGIN 286 print, 'WARNING, empty U points box... we use the same index as T points...'286 ras = report('WARNING, empty U points box... we use the same index as T points...') 287 287 firstxu = firstxt & lastxu = lastxt & nxu = nxt 288 288 firstyu = firstyt & lastyu = lastyt & nyu = nyt 289 289 END 290 290 ELSE:BEGIN 291 print, 'WARNING, empty U points box... we get the neighnors to define a new box...'291 ras = report('WARNING, empty U points box... we get the neighnors to define a new box...') 292 292 neig1 = neighbor(lon1, lat1, glamu, gphiu, sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) 293 293 neig2 = neighbor(lon2, lat2, glamu, gphiu, sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) … … 327 327 CASE 1 OF 328 328 (where(gridtype eq 'T'))[0] NE -1 OR (where(gridtype EQ 'W'))[0] NE -1:BEGIN 329 print, 'WARNING, empty V points box... we use the same index as T points...'329 ras = report('WARNING, empty V points box... we use the same index as T points...') 330 330 firstxv = firstxt & lastxv = lastxt & nxv = nxt 331 331 firstyv = firstyt & lastyv = lastyt & nyv = nyt 332 332 END 333 333 (where(gridtype eq 'U'))[0] NE -1:BEGIN 334 print, 'WARNING, empty V points box... we use the same index as U points...'334 ras = report('WARNING, empty V points box... we use the same index as U points...') 335 335 firstxv = firstxu & lastxv = lastxu & nxv = nxu 336 336 firstyv = firstyu & lastyv = lastyu & nyv = nyu 337 337 END 338 338 ELSE:BEGIN 339 print, 'WARNING, empty V points box... we get the neighnors to define a new box...'339 ras = report('WARNING, empty V points box... we get the neighnors to define a new box...') 340 340 neig1 = neighbor(lon1, lat1, glamv, gphiv, sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) 341 341 neig2 = neighbor(lon2, lat2, glamv, gphiv, sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) … … 375 375 CASE 1 OF 376 376 (where(gridtype eq 'T'))[0] NE -1 OR (where(gridtype EQ 'W'))[0] NE -1:BEGIN 377 print, 'WARNING, empty F points box... we use the same index as T points...'377 ras = report('WARNING, empty F points box... we use the same index as T points...') 378 378 firstxf = firstxt & lastxf = lastxt & nxf = nxt 379 379 firstyf = firstyt & lastyf = lastyt & nyf = nyt 380 380 END 381 381 (where(gridtype eq 'U'))[0] NE -1:BEGIN 382 print, 'WARNING, empty F points box... we use the same index as U points...'382 ras = report('WARNING, empty F points box... we use the same index as U points...') 383 383 firstxf = firstxu & lastxf = lastxu & nxf = nxu 384 384 firstyf = firstyu & lastyf = lastyu & nyf = nyu 385 385 END 386 386 (where(gridtype eq 'V'))[0] NE -1:BEGIN 387 print, 'WARNING, empty F points box... we use the same index as V points...'387 ras = report('WARNING, empty F points box... we use the same index as V points...') 388 388 firstxf = firstxv & lastxf = lastxv & nxf = nxv 389 389 firstyf = firstyv & lastyf = lastyv & nyf = nyv 390 390 END 391 391 ELSE:BEGIN 392 print, 'WARNING, empty F points box... we get the neighnors to define a new box...'392 ras = report('WARNING, empty F points box... we get the neighnors to define a new box...') 393 393 neig1 = neighbor(lon1, lat1, glamf, gphif, sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) 394 394 neig2 = neighbor(lon2, lat2, glamf, gphif, sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) … … 497 497 IF (dom[0] EQ -1) THEN BEGIN 498 498 IF keyword_set(findalways) THEN BEGIN 499 print, 'WARNING, empty T points box... we get the neighnors to define a new box...'499 ras = report('WARNING, empty T points box... we get the neighnors to define a new box...') 500 500 neig1 = neighbor(lon1, lat1, glamt[fstx:lstx, *], gphit[fstx:lstx, *], sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) 501 501 neig2 = neighbor(lon2, lat2, glamt[fstx:lstx, *], gphit[fstx:lstx, *], sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) … … 528 528 CASE 1 OF 529 529 (where(gridtype eq 'T'))[0] NE -1 OR (where(gridtype EQ 'W'))[0] NE -1:BEGIN 530 print, 'WARNING, empty U points box... we use the same index as T points...'530 ras = report( 'WARNING, empty U points box... we use the same index as T points...') 531 531 firstyu = firstyt & lastyu = lastyt & nyu = nyt 532 532 END 533 533 ELSE:BEGIN 534 print, 'WARNING, empty U points box... we get the neighnors to define a new box...'534 ras = report('WARNING, empty U points box... we get the neighnors to define a new box...') 535 535 neig1 = neighbor(lon1, lat1, glamu[fstx:lstx, *], gphiu[fstx:lstx, *], sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) 536 536 neig2 = neighbor(lon2, lat2, glamu[fstx:lstx, *], gphiu[fstx:lstx, *], sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) … … 566 566 CASE 1 OF 567 567 (where(gridtype eq 'T'))[0] NE -1 OR (where(gridtype EQ 'W'))[0] NE -1:BEGIN 568 print, 'WARNING, empty V points box... we use the same index as T points...'568 ras = report('WARNING, empty V points box... we use the same index as T points...') 569 569 firstyv = firstyt & lastyv = lastyt & nyv = nyt 570 570 END 571 571 (where(gridtype eq 'U'))[0] NE -1:BEGIN 572 print, 'WARNING, empty V points box... we use the same index as U points...'572 ras = report('WARNING, empty V points box... we use the same index as U points...') 573 573 firstyv = firstyu & lastyv = lastyu & nyv = nyu 574 574 END 575 575 ELSE:BEGIN 576 print, 'WARNING, empty V points box... we get the neighnors to define a new box...'576 ras = report('WARNING, empty V points box... we get the neighnors to define a new box...') 577 577 neig1 = neighbor(lon1, lat1, glamv[fstx:lstx, *], gphiv[fstx:lstx, *], sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) 578 578 neig2 = neighbor(lon2, lat2, glamv[fstx:lstx, *], gphiv[fstx:lstx, *], sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) … … 608 608 CASE 1 OF 609 609 (where(gridtype eq 'T'))[0] NE -1 OR (where(gridtype EQ 'W'))[0] NE -1:BEGIN 610 print, 'WARNING, empty F points box... we use the same index as T points...'610 ras = report('WARNING, empty F points box... we use the same index as T points...') 611 611 firstyf = firstyt & lastyf = lastyt & nyf = nyt 612 612 END 613 613 (where(gridtype eq 'U'))[0] NE -1:BEGIN 614 print, 'WARNING, empty F points box... we use the same index as U points...'614 ras = report('WARNING, empty F points box... we use the same index as U points...') 615 615 firstyf = firstyu & lastyf = lastyu & nyf = nyu 616 616 END 617 617 (where(gridtype eq 'V'))[0] NE -1:BEGIN 618 print, 'WARNING, empty F points box... we use the same index as V points...'618 ras = report('WARNING, empty F points box... we use the same index as V points...') 619 619 firstyf = firstyv & lastyf = lastyv & nyf = nyv 620 620 END 621 621 ELSE:BEGIN 622 print, 'WARNING, empty F points box... we get the neighnors to define a new box...'622 ras = report('WARNING, empty F points box... we get the neighnors to define a new box...') 623 623 neig1 = neighbor(lon1, lat1, glamf[fstx:lstx, *], gphif[fstx:lstx, *], sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) 624 624 neig2 = neighbor(lon2, lat2, glamf[fstx:lstx, *], gphif[fstx:lstx, *], sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) … … 665 665 IF (dom[0] EQ -1) THEN BEGIN 666 666 IF keyword_set(findalways) THEN BEGIN 667 print, 'WARNING, empty T points box... we get the neighnors to define a new box...'667 ras = report('WARNING, empty T points box... we get the neighnors to define a new box...') 668 668 neig1 = neighbor(lon1, lat1, glamt[*, fsty:lsty], gphit[*, fsty:lsty], sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) 669 669 neig2 = neighbor(lon2, lat2, glamt[*, fsty:lsty], gphit[*, fsty:lsty], sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) … … 696 696 CASE 1 OF 697 697 (where(gridtype eq 'T'))[0] NE -1 OR (where(gridtype EQ 'W'))[0] NE -1:BEGIN 698 print, 'WARNING, empty U points box... we use the same index as T points...'698 ras = report('WARNING, empty U points box... we use the same index as T points...') 699 699 firstxu = firstxt & lastxu = lastxt & nxu = nxt 700 700 END 701 701 ELSE:BEGIN 702 print, 'WARNING, empty U points box... we get the neighnors to define a new box...'702 ras = report('WARNING, empty U points box... we get the neighnors to define a new box...') 703 703 neig1 = neighbor(lon1, lat1, glamu[*, fsty:lsty], gphiu[*, fsty:lsty], sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) 704 704 neig2 = neighbor(lon2, lat2, glamu[*, fsty:lsty], gphiu[*, fsty:lsty], sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) … … 734 734 CASE 1 OF 735 735 (where(gridtype eq 'T'))[0] NE -1 OR (where(gridtype EQ 'W'))[0] NE -1:BEGIN 736 print, 'WARNING, empty V points box... we use the same index as T points...'736 ras = report('WARNING, empty V points box... we use the same index as T points...') 737 737 firstxv = firstxt & lastxv = lastxt & nxv = nxt 738 738 END 739 739 (where(gridtype eq 'U'))[0] NE -1:BEGIN 740 print, 'WARNING, empty V points box... we use the same index as U points...'740 ras = report('WARNING, empty V points box... we use the same index as U points...') 741 741 firstxv = firstxu & lastxv = lastxu & nxv = nxu 742 742 END 743 743 ELSE:BEGIN 744 print, 'WARNING, empty V points box... we get the neighnors to define a new box...'744 ras = report('WARNING, empty V points box... we get the neighnors to define a new box...') 745 745 neig1 = neighbor(lon1, lat1, glamv[*, fsty:lsty], gphiv[*, fsty:lsty], sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) 746 746 neig2 = neighbor(lon2, lat2, glamv[*, fsty:lsty], gphiv[*, fsty:lsty], sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) … … 776 776 CASE 1 OF 777 777 (where(gridtype eq 'T'))[0] NE -1 OR (where(gridtype EQ 'W'))[0] NE -1:BEGIN 778 print, 'WARNING, empty F points box... we use the same index as T points...'778 ras = report('WARNING, empty F points box... we use the same index as T points...') 779 779 firstxf = firstxt & lastxf = lastxt & nxf = nxt 780 780 END 781 781 (where(gridtype eq 'U'))[0] NE -1:BEGIN 782 print, 'WARNING, empty F points box... we use the same index as U points...'782 ras = report('WARNING, empty F points box... we use the same index as U points...') 783 783 firstxf = firstxu & lastxf = lastxu & nxf = nxu 784 784 END 785 785 (where(gridtype eq 'V'))[0] NE -1:BEGIN 786 print, 'WARNING, empty F points box... we use the same index as V points...'786 ras = report('WARNING, empty F points box... we use the same index as V points...') 787 787 firstxf = firstxv & lastxf = lastxv & nxf = nxv 788 788 END 789 789 ELSE:BEGIN 790 print, 'WARNING, empty F points box... we get the neighnors to define a new box...'790 ras = report('WARNING, empty F points box... we get the neighnors to define a new box...') 791 791 neig1 = neighbor(lon1, lat1, glamf[*, fsty:lsty], gphif[*, fsty:lsty], sphere = keyword_set(key_onearth) * keyword_set(key_irregular)) 792 792 neig2 = neighbor(lon2, lat2, glamf[*, fsty:lsty], gphif[*, fsty:lsty], sphere = keyword_set(key_onearth) * keyword_set(key_irregular))
Note: See TracChangeset
for help on using the changeset viewer.