Changeset 816 for misc/GISTEMP/GISTEMP.pro

Timestamp:

08/27/07 20:08:59 (13 months ago)

Author:

mankoff

Message:

Misc logic checks for years with little or no data. Changed all location markers to pins from airplanes due to Google Earth bug (blue airplanes)

Files:

• misc/GISTEMP/GISTEMP.pro (modified) (5 diffs)

misc/GISTEMP/GISTEMP.pro

@@ -92 +92 @@
     s[i].seasons.data[0:n_elements(year)-1,0:3] = [[djf],[mam],[jja],[son]]
     s[i].annual = ann
-    ygd = where( s[i].years ne 0 )
-    y0 = min( s[i].years[ygd] )
-    s[i].years = SHIFT( s[i].years, y0-1880 )
-    s[i].months.data = SHIFT( s[i].months.data, y0-1880 )
-    s[i].seasons.data = SHIFT( s[i].seasons.data, y0-1880 )
-    s[i].annual = SHIFT( s[i].annual, y0-1880 )
+    ygd = where( s[i].years ne 0, count )
+    if count ne 0 then begin
+       y0 = min( s[i].years[ygd] )
+       s[i].years = SHIFT( s[i].years, y0-1880 )
+       s[i].months.data = SHIFT( s[i].months.data, y0-1880 )
+       s[i].seasons.data = SHIFT( s[i].seasons.data, y0-1880 )
+       s[i].annual = SHIFT( s[i].annual, y0-1880 )
+    endif
 
     print, i
@@ -266 +268 @@
 load_panoply
 
-if n_elements( where( s.years ne 0 ) ) LE 20 then MESSAGE, "Less than 20 data points at this station", /CONTINUE
+nyrs = n_elements(where( s.years ne 0 ) )
+if nyrs le 20 then MESSAGE, "Less than 20 data points at this station", /CONTINUE
 
 dd = s.months.data
 bad = where( dd eq 0 or dd eq 999.90, nbad, complement=gd, ncomplement=ngd )
+if n_elements(gd) le 2 then begin
+   MESSAGE, "Error... Not enough data to produce image", /CONTINUE
+   return
+end
 img = BYTSCL( dd, min=min(dd[gd]), max=max(dd[gd]), top=252 )+1 ; 1 to 253
 img[ bad ] = 254
@@ -561 +568 @@
 ; Airport or not (use airplane icon)
 printf, lun, '<Icon>'
-printf, lun, ghcn.a eq 'A'?  $
-        '<href>http://maps.google.com/mapfiles/kml/pal2/icon56.png</href>' : $
-        '<href>http://maps.google.com/mapfiles/kml/paddle/wht-blank.png</href>'
+;; Airplane icon shows up blue (Google Earth bug?), so only do pins...
+;; printf, lun, ghcn.a eq 'A'?  $
+;;         '<href>http://maps.google.com/mapfiles/kml/pal2/icon56.png</href>' : $
+;;         '<href>http://maps.google.com/mapfiles/kml/paddle/wht-blank.png</href>'
+printf, lun, '<href>http://maps.google.com/mapfiles/kml/paddle/wht-blank.png</href>'
 printf, lun, '</Icon>'
 
@@ -704 +713 @@
 ; 
 ; clean:
-; giss = LOAD_GISS_STATION_LIST()
-; ghcn = LOAD_GHCN_STATION_LIST()
-; data = load_station_data(giss,ghcn)
-; meta = make_meta_data(giss,ghcn,data)
+;; giss = LOAD_GISS_STATION_LIST()
+;; ghcn = LOAD_GHCN_STATION_LIST()
+;; data = load_station_data(giss,ghcn)
+;; meta = make_meta_data(giss,ghcn,data)
 
 ; quick:
@@ -717 +726 @@
 ;for i=0,100 do make_img, data[i]
 for i=0,n_elements(data)-1 do make_img, data[i]
-;makeKML, giss, ghcn, data, meta
+makeKML, giss, ghcn, data, meta
+
+; $ for v in *.ps; do convert $v $v.png; done
+; gzip the KML file
+; upload KML_net, images, and if changed, KML and CGI
 
 end