Changeset 307 for trunk

Timestamp:

11/05/07 11:18:32 (17 years ago)

Author:

smasson

Message:

minor bugfix in error test

File:

• trunk/SRC/Interpolation/compute_fromirr_bilinear_weigaddr.pro (modified) (2 diffs)

trunk/SRC/Interpolation/compute_fromirr_bilinear_weigaddr.pro

@@ -317 +317 @@
   a = -1 &  b = -1              ; free memory
 ; mask the weight to suppress the corner located on land
-  newaweig = newaweig*((omsk)[newaaddr])
+  newaweig = newaweig*(omsk[newaaddr])
 ; for cell with some land corner,
 ; we have to redistribute the weight on the remaining water corners
 ; weights normalization
   totalweig = total(newaweig, 1, /double)
-  IF abs(min(totalweig, max = ma)) LT 1.e-6 then stop
-  IF abs(1.d - ma) GT 1.e-6 then stop
+  IF min(totalweig, max = ma) LE 0.d then stop
+  IF ma GT 1.d then stop
   newaweig = newaweig/(replicate(1.d, 4)#totalweig)
 ; weights
@@ -329 +329 @@
   weig[*, awater] = temporary(newaweig)
 ; address
-  addr = dblarr(4, jpia*jpja)
+  addr = lonarr(4, jpia*jpja)
   addr[*, awater] = temporary(newaaddr)
 ;