Changeset 80 for trunk/src

Timestamp:

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

Author:

pinsard

Message:

progress on humidity processing

Location:

trunk/src

Files:

• d2m_to_q2m_erai.pro (modified) (9 diffs)
• interp_erai_msl_1989_2009.pro (modified) (4 diffs)
• lwr_correction_ncdf.pro (modified) (1 diff)
• q2m_correction_ncdf.pro (modified) (4 diffs)
• qsat.pro (modified) (2 diffs)
• rh_to_spechum.pro (modified) (2 diffs)

trunk/src/d2m_to_q2m_erai.pro

@@ -19 +19 @@
 ;
 ; q2m on OAFLUX grid is written in
-; :file:`${PROJECT_OD}/q2m_erai_19890101_20091231_oafluxgrid.nc`
+; :file:`${PROJECT_OD}/erai_q2m_19890101_20091231_oafluxgrid.nc`
 ; if this file not already exists.
 ;
@@ -34 +34 @@
 ;           file_t2m [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/erai_t2m_19890101_20091231_oafluxgrid.nc"];
 ;           file_msl [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/erai_msl_19890101_20091231_oafluxgrid.nc"];
-;           file_q2m [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/q2m_erai_19890101_20091231_oafluxgrid.nc"];
+;           file_q2m [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/erai_q2m_19890101_20091231_oafluxgrid.nc"];
 ;
 ;           d2m_to_q2m_erai [shape=box,
@@ -48 +48 @@
 ; ========
 ;
+; [Gill:AP:1982]_
+;
+; [AOMIP]_
+;
 ; :ref:`project_profile.sh`
 ;
@@ -74 +78 @@
 ; ====
 ;
-; es0 undefined
+; pb on cratos idl7 even if output file is written::
+;
+;   % Program caused arithmetic error: Floating overflow
+;   % Program caused arithmetic error: Floating illegal operand
+;
+; related to occurence of NaNf in ncdump output ? do no yet
+;
+; pk also see these error messages ...
+;
+; check for NaN value in P return by reading Netcdf msl file
+; solution is may be to apply *mskout to q2m like  in interp* tools
+; not yet done to keep close as possible to pk processing
 ;
 ; coding rules
@@ -84 +99 @@
 ; 
 ; $URL$
+;
+; - fplod 20110811T114932Z aedon.locean-ipsl.upmc.fr (Darwin)
+;
+;   * add es0 value
+;   * add reference
+;   * replace q2m_erai by erai_q2m in output filename
 ;
 ; - fplod 20110809T154249Z aedon.locean-ipsl.upmc.fr (Darwin)
@@ -186 +207 @@
 ;
 ; build output filename
-filename_out = 'q2m_erai_19890101_20091231_oafluxgrid.nc'
+filename_out = 'erai_q2m_19890101_20091231_oafluxgrid.nc'
 fullfilename_out = iodirout + filename_out
 ; in order to avoid unexpected overwritten
@@ -205 +226 @@
 help, t, td
 
-;;since the t2m and d2m are already in K, no need for conversion
-
+;since the t2m and d2m are already in K, no need for conversion
+es0=6.11
 e=es0*exp((25.e5/461.5)*((1/273.15)-(1/td)))
 es=es0*exp((25.e5/461.5)*((1/273.15)-(1/t)))
@@ -221 +242 @@
 P=read_ncdf("msl", d1, d2, file=fullfilename_msl,/nostr)
 help, P
+; ++ for debug ftp error 
+; ++print,P(0:1,0,0)
+; ++stop
 q2m=rh_to_spechum(rh,t,P)
 
@@ -239 +263 @@
 globattr={source:'Surface specific humidity at 2m calculated from dew point temperature',timerange:cda0+' - '+cda1}
 
-ncfields = 'q2m[longitude,latitude,time]=q2m:q2m_attr; ' $
+ncfields = i'q2m[longitude,latitude,time]=q2m:q2m_attr; ' $
                       + 'longitude[]=xlon:lon_attr; ' $
                       + 'latitude[]=ylat:lat_attr; ' $

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

trunk/src/lwr_correction_ncdf.pro

@@ -70 +70 @@
 ; 19890101 is not provided in ERA-Intermim file str : a possible workaround is
 ; to copy 19890102
-; dataset to 19890101 using nco before processing. see
-; ref:`compute_erai_daily_region_2d.sh`
+; dataset to 19890101 using nco before processing.
+;
+; here is the sequence used by pk wher erai_lwr_19890102_20091231_oafluxgrid.nc is the regridded (interpolated output)
+; 
+;   $ ncks -d time,0,0 erai_lwr_19890102_20091231_oafluxgrid.nc erai_lwr_step0.nc
+;
+; This will combine the two files ::
+;
+;   $ ncrcat erai_lwr_step0.nc erai_lwr_19890102_20091231_oafluxgrid.nc erai_lwr_19890101_20091231_oafluxgrid.nc.
+;
+; must be put in a shell script or find IDL equivalent
 ;
 ; here I (fp) replace hard coded date 19890102 by 19890101. must check

trunk/src/q2m_correction_ncdf.pro

@@ -21 +21 @@
 ;
 ;        digraph q2m_correction_ncdf {
-;           graph [
-;           rankdir="LR",
-;           ]
 ;
 ;           file_in [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/erai_q2m_19890101_20091231_oafluxgrid.nc"];
@@ -63 +60 @@
 ; ====
 ;
-; can be tested because d2m_to_q2m.pro pb
+; work on cratos idl7 even if NaNf values in erai_q2m_19890101_20091231_oafluxgrid.nc written by
+; d2m_to_q2m.pro.
+; 
+; No way ... NaNf also in this output !!
 ;
 ; coding rules
@@ -95 +95 @@
 @cm_4data
 @cm_project
-;
 ;
 ; test if ${PROJECT_OD} defined
@@ -133 +132 @@
 ;
 ; check if this file exists
+msg='iii : looking for ' + filename
+ras = report(msg)
 fullfilename = isafile(iodirout + filename, NEW=0, /MUST_EXIST)
 IF fullfilename[0] EQ '' THEN BEGIN

trunk/src/qsat.pro

@@ -1 @@
-;--------------------------------------------------------
-function qsat, T,P
+;+
 ;
-; Air humidity at saturation (g/kg)
-; T: SST (degC)
-; P : surface pressure (mb)
-;--------------------------------------------------------
+; ========
+; qsat.pro
+; ========
+;
+; .. function:: qsat(T,P)
+;
+;     :param T: SST (degC)
+;     :param P: surface pressure (mb)
+; 
+;     :returns: Air humidity at saturation (g/kg)
+;     
+; EXAMPLES
+; ========
+;
+; ::
+;
+;  IDL> T=+todo+
+;  IDL> P=+todo+
+;  IDL> result=qsat(T,P)
+;
+; TODO
+; ====
+;
+; coding rules
+;
+; complete example
+;
+; get rid of uppercase
+;
+; SEE ALSO
+; ========
+;
+; :func:`rh_to_spechum`
+;
+; EVOLUTIONS
+; ==========
+;
+; $Id$
+;
+; $URL$
+;
+; - fplod 20110811T131113Z aedon.locean-ipsl.upmc.fr (Darwin)
+;
+;   * add minimal header to pk release given today (!)
+;
+; - pbk 2008
+;
+;   * creation
+;
+;-
+function qsat, T, P
+;
 x=T
 es=6.112*exp(17.502*x/(x+241.))*(1.0007+3.46e-6*p)
@@ -13 +60 @@
 
 end
-;--------------------------------------------------------
-

trunk/src/rh_to_spechum.pro

@@ -1 +1 @@
 ;+
 ;
-; =============
-; rh_to_spechum
-; =============
+; =================
+; rh_to_spechum.pro
+; =================
 ;
 ; .. function:: rh_to_spechum(rh,ta,P)
@@ -26 +26 @@
 ; ====
 ;
-; hard coded directory - usage of ${PROJECT_ID}
+; coding rules
 ;
-; coding rules
+; complete example
+;
+; SEE ALSO
+; ========
+;
+; called by :ref:`d2m_to_q2m_erai.pro`
+;
+; use :func:`qsat`
 ;
 ; EVOLUTIONS
 ; ==========
+;
+; $Id$
+;
+; $URL$
+;
+; - fplod 20110811T130223Z aedon.locean-ipsl.upmc.fr (Darwin)
+;
+;   * add reference to qsat function
 ;
 ; - fplod 20101214T093615Z aedon.locean-ipsl.upmc.fr (Darwin)