Index: trunk/SRC/ToBeReviewed/GRILLE/domdef.pro
===================================================================
--- trunk/SRC/ToBeReviewed/GRILLE/domdef.pro	(revision 232)
+++ trunk/SRC/ToBeReviewed/GRILLE/domdef.pro	(revision 240)
@@ -248,5 +248,5 @@
       IF (dom[0] EQ -1) THEN BEGIN
         IF keyword_set(findalways) THEN BEGIN
-          print, 'WARNING, empty T points box... we get the neighnors to define a new box...'
+          ras = report('WARNING, empty T points box... we get the neighnors to define a new box...')
           neig1 = neighbor(lon1, lat1, glamt, gphit, sphere = keyword_set(key_onearth) * keyword_set(key_irregular))
           neig2 = neighbor(lon2, lat2, glamt, gphit, sphere = keyword_set(key_onearth) * keyword_set(key_irregular))
@@ -284,10 +284,10 @@
             CASE 1 OF
               (where(gridtype eq 'T'))[0] NE -1 OR (where(gridtype EQ 'W'))[0] NE -1:BEGIN
-                print, 'WARNING, empty U points box... we use the same index as T points...'
+                ras = report('WARNING, empty U points box... we use the same index as T points...')
                 firstxu = firstxt & lastxu = lastxt & nxu = nxt
                 firstyu = firstyt & lastyu = lastyt & nyu = nyt
               END
               ELSE:BEGIN
-                print, 'WARNING, empty U points box... we get the neighnors to define a new box...'
+                ras = report('WARNING, empty U points box... we get the neighnors to define a new box...')
                 neig1 = neighbor(lon1, lat1, glamu, gphiu, sphere = keyword_set(key_onearth) * keyword_set(key_irregular))
                 neig2 = neighbor(lon2, lat2, glamu, gphiu, sphere = keyword_set(key_onearth) * keyword_set(key_irregular))
@@ -327,15 +327,15 @@
             CASE 1 OF
               (where(gridtype eq 'T'))[0] NE -1 OR (where(gridtype EQ 'W'))[0] NE -1:BEGIN
-                print, 'WARNING, empty V points box... we use the same index as T points...'
+                ras = report('WARNING, empty V points box... we use the same index as T points...')
                 firstxv = firstxt & lastxv = lastxt & nxv = nxt
                 firstyv = firstyt & lastyv = lastyt & nyv = nyt
               END
               (where(gridtype eq 'U'))[0] NE -1:BEGIN
-                print, 'WARNING, empty V points box... we use the same index as U points...'
+                ras = report('WARNING, empty V points box... we use the same index as U points...')
                 firstxv = firstxu & lastxv = lastxu & nxv = nxu
                 firstyv = firstyu & lastyv = lastyu & nyv = nyu
               END
               ELSE:BEGIN
-                print, 'WARNING, empty V points box... we get the neighnors to define a new box...'
+                ras = report('WARNING, empty V points box... we get the neighnors to define a new box...')
                 neig1 = neighbor(lon1, lat1, glamv, gphiv, sphere = keyword_set(key_onearth) * keyword_set(key_irregular))
                 neig2 = neighbor(lon2, lat2, glamv, gphiv, sphere = keyword_set(key_onearth) * keyword_set(key_irregular))
@@ -375,20 +375,20 @@
             CASE 1 OF
               (where(gridtype eq 'T'))[0] NE -1 OR (where(gridtype EQ 'W'))[0] NE -1:BEGIN
-                print, 'WARNING, empty F points box... we use the same index as T points...'
+                ras = report('WARNING, empty F points box... we use the same index as T points...')
                 firstxf = firstxt & lastxf = lastxt & nxf = nxt
                 firstyf = firstyt & lastyf = lastyt & nyf = nyt
               END
               (where(gridtype eq 'U'))[0] NE -1:BEGIN
-                print, 'WARNING, empty F points box... we use the same index as U points...'
+                ras = report('WARNING, empty F points box... we use the same index as U points...')
                 firstxf = firstxu & lastxf = lastxu & nxf = nxu
                 firstyf = firstyu & lastyf = lastyu & nyf = nyu
               END
               (where(gridtype eq 'V'))[0] NE -1:BEGIN
-                print, 'WARNING, empty F points box... we use the same index as V points...'
+                ras = report('WARNING, empty F points box... we use the same index as V points...')
                 firstxf = firstxv & lastxf = lastxv & nxf = nxv
                 firstyf = firstyv & lastyf = lastyv & nyf = nyv
               END
               ELSE:BEGIN
-                print, 'WARNING, empty F points box... we get the neighnors to define a new box...'
+                ras = report('WARNING, empty F points box... we get the neighnors to define a new box...')
                 neig1 = neighbor(lon1, lat1, glamf, gphif, sphere = keyword_set(key_onearth) * keyword_set(key_irregular))
                 neig2 = neighbor(lon2, lat2, glamf, gphif, sphere = keyword_set(key_onearth) * keyword_set(key_irregular))
@@ -497,5 +497,5 @@
           IF (dom[0] EQ -1) THEN BEGIN
             IF keyword_set(findalways) THEN BEGIN
-              print, 'WARNING, empty T points box... we get the neighnors to define a new box...'
+              ras = report('WARNING, empty T points box... we get the neighnors to define a new box...')
               neig1 = neighbor(lon1, lat1, glamt[fstx:lstx, *], gphit[fstx:lstx, *], sphere = keyword_set(key_onearth) * keyword_set(key_irregular))
               neig2 = neighbor(lon2, lat2, glamt[fstx:lstx, *], gphit[fstx:lstx, *], sphere = keyword_set(key_onearth) * keyword_set(key_irregular))
@@ -528,9 +528,9 @@
                 CASE 1 OF
                   (where(gridtype eq 'T'))[0] NE -1 OR (where(gridtype EQ 'W'))[0] NE -1:BEGIN
-                    print, 'WARNING, empty U points box... we use the same index as T points...'
+                    ras = report( 'WARNING, empty U points box... we use the same index as T points...')
                     firstyu = firstyt & lastyu = lastyt & nyu = nyt
                   END
                   ELSE:BEGIN
-                    print, 'WARNING, empty U points box... we get the neighnors to define a new box...'
+                    ras = report('WARNING, empty U points box... we get the neighnors to define a new box...')
                     neig1 = neighbor(lon1, lat1, glamu[fstx:lstx, *], gphiu[fstx:lstx, *], sphere = keyword_set(key_onearth) * keyword_set(key_irregular))
                     neig2 = neighbor(lon2, lat2, glamu[fstx:lstx, *], gphiu[fstx:lstx, *], sphere = keyword_set(key_onearth) * keyword_set(key_irregular))
@@ -566,13 +566,13 @@
                 CASE 1 OF
                   (where(gridtype eq 'T'))[0] NE -1 OR (where(gridtype EQ 'W'))[0] NE -1:BEGIN
-                    print, 'WARNING, empty V points box... we use the same index as T points...'
+                    ras = report('WARNING, empty V points box... we use the same index as T points...')
                     firstyv = firstyt & lastyv = lastyt & nyv = nyt
                   END
                   (where(gridtype eq 'U'))[0] NE -1:BEGIN
-                    print, 'WARNING, empty V points box... we use the same index as U points...'
+                    ras = report('WARNING, empty V points box... we use the same index as U points...')
                     firstyv = firstyu & lastyv = lastyu & nyv = nyu
                   END
                   ELSE:BEGIN
-                    print, 'WARNING, empty V points box... we get the neighnors to define a new box...'
+                    ras = report('WARNING, empty V points box... we get the neighnors to define a new box...')
                     neig1 = neighbor(lon1, lat1, glamv[fstx:lstx, *], gphiv[fstx:lstx, *], sphere = keyword_set(key_onearth) * keyword_set(key_irregular))
                     neig2 = neighbor(lon2, lat2, glamv[fstx:lstx, *], gphiv[fstx:lstx, *], sphere = keyword_set(key_onearth) * keyword_set(key_irregular))
@@ -608,17 +608,17 @@
                 CASE 1 OF
                   (where(gridtype eq 'T'))[0] NE -1 OR (where(gridtype EQ 'W'))[0] NE -1:BEGIN
-                    print, 'WARNING, empty F points box... we use the same index as T points...'
+                    ras = report('WARNING, empty F points box... we use the same index as T points...')
                     firstyf = firstyt & lastyf = lastyt & nyf = nyt
                   END
                   (where(gridtype eq 'U'))[0] NE -1:BEGIN
-                    print, 'WARNING, empty F points box... we use the same index as U points...'
+                    ras = report('WARNING, empty F points box... we use the same index as U points...')
                     firstyf = firstyu & lastyf = lastyu & nyf = nyu
                   END
                   (where(gridtype eq 'V'))[0] NE -1:BEGIN
-                    print, 'WARNING, empty F points box... we use the same index as V points...'
+                    ras = report('WARNING, empty F points box... we use the same index as V points...')
                     firstyf = firstyv & lastyf = lastyv & nyf = nyv
                   END
                   ELSE:BEGIN
-                    print, 'WARNING, empty F points box... we get the neighnors to define a new box...'
+                    ras = report('WARNING, empty F points box... we get the neighnors to define a new box...')
                     neig1 = neighbor(lon1, lat1, glamf[fstx:lstx, *], gphif[fstx:lstx, *], sphere = keyword_set(key_onearth) * keyword_set(key_irregular))
                     neig2 = neighbor(lon2, lat2, glamf[fstx:lstx, *], gphif[fstx:lstx, *], sphere = keyword_set(key_onearth) * keyword_set(key_irregular))
@@ -665,5 +665,5 @@
           IF (dom[0] EQ -1) THEN BEGIN
             IF keyword_set(findalways) THEN BEGIN
-              print, 'WARNING, empty T points box... we get the neighnors to define a new box...'
+              ras = report('WARNING, empty T points box... we get the neighnors to define a new box...')
               neig1 = neighbor(lon1, lat1, glamt[*, fsty:lsty], gphit[*, fsty:lsty], sphere = keyword_set(key_onearth) * keyword_set(key_irregular))
               neig2 = neighbor(lon2, lat2, glamt[*, fsty:lsty], gphit[*, fsty:lsty], sphere = keyword_set(key_onearth) * keyword_set(key_irregular))
@@ -696,9 +696,9 @@
                 CASE 1 OF
                   (where(gridtype eq 'T'))[0] NE -1 OR (where(gridtype EQ 'W'))[0] NE -1:BEGIN
-                    print, 'WARNING, empty U points box... we use the same index as T points...'
+                    ras = report('WARNING, empty U points box... we use the same index as T points...')
                     firstxu = firstxt & lastxu = lastxt & nxu = nxt
                   END
                   ELSE:BEGIN
-                    print, 'WARNING, empty U points box... we get the neighnors to define a new box...'
+                    ras = report('WARNING, empty U points box... we get the neighnors to define a new box...')
                     neig1 = neighbor(lon1, lat1, glamu[*, fsty:lsty], gphiu[*, fsty:lsty], sphere = keyword_set(key_onearth) * keyword_set(key_irregular))
                     neig2 = neighbor(lon2, lat2, glamu[*, fsty:lsty], gphiu[*, fsty:lsty], sphere = keyword_set(key_onearth) * keyword_set(key_irregular))
@@ -734,13 +734,13 @@
                 CASE 1 OF
                   (where(gridtype eq 'T'))[0] NE -1 OR (where(gridtype EQ 'W'))[0] NE -1:BEGIN
-                    print, 'WARNING, empty V points box... we use the same index as T points...'
+                    ras = report('WARNING, empty V points box... we use the same index as T points...')
                     firstxv = firstxt & lastxv = lastxt & nxv = nxt
                   END
                   (where(gridtype eq 'U'))[0] NE -1:BEGIN
-                    print, 'WARNING, empty V points box... we use the same index as U points...'
+                    ras = report('WARNING, empty V points box... we use the same index as U points...')
                     firstxv = firstxu & lastxv = lastxu & nxv = nxu
                   END
                   ELSE:BEGIN
-                    print, 'WARNING, empty V points box... we get the neighnors to define a new box...'
+                    ras = report('WARNING, empty V points box... we get the neighnors to define a new box...')
                     neig1 = neighbor(lon1, lat1, glamv[*, fsty:lsty], gphiv[*, fsty:lsty], sphere = keyword_set(key_onearth) * keyword_set(key_irregular))
                     neig2 = neighbor(lon2, lat2, glamv[*, fsty:lsty], gphiv[*, fsty:lsty], sphere = keyword_set(key_onearth) * keyword_set(key_irregular))
@@ -776,17 +776,17 @@
                 CASE 1 OF
                   (where(gridtype eq 'T'))[0] NE -1 OR (where(gridtype EQ 'W'))[0] NE -1:BEGIN
-                    print, 'WARNING, empty F points box... we use the same index as T points...'
+                    ras = report('WARNING, empty F points box... we use the same index as T points...')
                     firstxf = firstxt & lastxf = lastxt & nxf = nxt
                   END
                   (where(gridtype eq 'U'))[0] NE -1:BEGIN
-                    print, 'WARNING, empty F points box... we use the same index as U points...'
+                    ras = report('WARNING, empty F points box... we use the same index as U points...')
                     firstxf = firstxu & lastxf = lastxu & nxf = nxu
                   END
                   (where(gridtype eq 'V'))[0] NE -1:BEGIN
-                    print, 'WARNING, empty F points box... we use the same index as V points...'
+                    ras = report('WARNING, empty F points box... we use the same index as V points...')
                     firstxf = firstxv & lastxf = lastxv & nxf = nxv
                   END
                   ELSE:BEGIN
-                    print, 'WARNING, empty F points box... we get the neighnors to define a new box...'
+                    ras = report('WARNING, empty F points box... we get the neighnors to define a new box...')
                     neig1 = neighbor(lon1, lat1, glamf[*, fsty:lsty], gphif[*, fsty:lsty], sphere = keyword_set(key_onearth) * keyword_set(key_irregular))
                     neig2 = neighbor(lon2, lat2, glamf[*, fsty:lsty], gphif[*, fsty:lsty], sphere = keyword_set(key_onearth) * keyword_set(key_irregular))
