Changeset 80 for trunk/src/interp_erai_msl_1989_2009.pro

Timestamp:

08/11/11 17:29:34 (13 years ago)

Author:

pinsard

Message:

progress on humidity processing

File:

• trunk/src/interp_erai_msl_1989_2009.pro (modified) (4 diffs)

trunk/src/interp_erai_msl_1989_2009.pro

@@ -73 +73 @@
 ; ====
 ;
-; no input file available
-;
-; hard coded scale factor ? mslin=mslin/100
-;
 ; msl_attr={units:'milibars' .. millibar at least, hPa will be probably
 ; more SI
@@ -82 +78 @@
 ; why 19880101,20100930 as dates
 ;
+; first step time 14245.5 (days since 1950-01-01 00:00:00) : correct ?
+;
 ; KNOWN ISSUES
 ; ============
@@ -94 +92 @@
 ;
 ; $URL$
+;
+; - fplod 20110811T113646Z aedon.locean-ipsl.upmc.fr (Darwin)
+;
+;   * check for units in original file before conversion
 ;
 ; - fplod 20110809T145539Z aedon.locean-ipsl.upmc.fr (Darwin)
@@ -217 +219 @@
 print, 'lon grid ',min(lonout),max(lonout),lonout(1)-lonout(0)
 mskout=read_ncdf("msk", file=fullfilename_msk,/nostr)
-mslin=mslin/100
+
+; conversion to hPa
+ncdf_getatt, fullfilename, 'msl', units=units
+CASE units OF
+    'Pascal': BEGIN
+        mslin=mslin/100
+    END
+    ELSE: BEGIN
+        msg = 'eee : ' + units + ' unknown'
+        ras = report(msg)
+        STOP
+    END
+ENDCASE
 
 help, mslin,lonin,latin,mskin,lonout,latout,mskout