Changeset 257 for trunk

Timestamp:

07/26/07 19:35:23 (17 years ago)

Author:

smasson

Message:

allow interpolation from ORCA2 whithout mask

Location:

trunk/SRC

Files:

• Documentation/xmldoc/savesaxo.sh (modified) (1 diff)
• Interpolation/compute_fromirr_bilinear_weigaddr.pro (modified) (1 diff)

trunk/SRC/Documentation/xmldoc/savesaxo.sh

@@ -99 +99 @@
 #
 docompileall=1
-dohtml=1
+dohtml=0
 dooldtest=0
 doidlwavedoc=0

trunk/SRC/Interpolation/compute_fromirr_bilinear_weigaddr.pro

@@ -69 +69 @@
   jpjo = (size(olonin, /dimensions))[1]
 ; mask check
-  IF n_elements(omsk) EQ 1 AND omsk[0] EQ -1 THEN omsk = replicate(1b, jpio, jpjo)
+  IF n_elements(omsk) EQ 1 AND omsk[0] EQ -1 THEN BEGIN
+    omsk = replicate(1b, jpio, jpjo)
+   IF jpio EQ 182 AND jpjo EQ 149 THEN BEGIN
+      lontmp = (olonin[1:180, *] + 3600) MOD 360
+      lattmp = olat[1:180, *]
+      bad = abs(abs(lontmp - shift(lontmp, 1, 0)) - 180) LT 176 $
+            OR  abs(abs(lontmp - shift(lontmp, 0, 1)) - 180) LT 176 $
+            OR  abs(abs(lontmp - shift(lontmp, 1, 1)) - 180) LT 176 $
+            OR  abs(abs(lattmp - shift(lattmp, 1, 0)) - 180) LT 176 $
+            OR  abs(abs(lattmp - shift(lattmp, 0, 1)) - 180) LT 176 $
+            OR  abs(abs(lattmp - shift(lattmp, 1, 1)) - 180) LT 176
+      bad = bad AND lattmp LT 50
+      omsk[1:180, 1:148] = 1b - bad[*, 1:148]
+      omsk[0, *] = omsk[180, *]
+      omsk[181, *] = omsk[1, *]
+    ENDIF ELSE omsk = replicate(1b, jpio, jpjo)
+  ENDIF
   IF n_elements(amsk) EQ 1 AND amsk[0] EQ -1 THEN amsk = replicate(1b, jpia, jpja)
   IF n_elements(omsk) NE jpio*jpjo THEN BEGIN