Changeset 90 for trunk/src

Timestamp:

08/22/11 17:42:26 (13 years ago)

Author:

pinsard

Message:

progress on coare correction

Location:

trunk/src

Files:

• TropFlux_19890101_20091231.pro (modified) (7 diffs, 1 prop)
• TropFlux_NRT_ncdf.pro (modified) (4 diffs, 1 prop)
• TropFlux_swr_BLND_19890101_20091231.pro (modified) (1 diff, 1 prop)
• TropFlux_swr_DT_19890101_20071231.pro (modified) (1 diff, 1 prop)
• TropFlux_swr_NRT_19890101_20091231.pro (modified) (1 diff, 1 prop)
• cor30a.pro (modified) (2 diffs)
• interp_erai_msl_1989_2009.pro (modified) (1 diff, 1 prop)
• interp_olr_30n30s_1989_2009.pro (modified) (1 diff, 1 prop)
• lwr_correction_ncdf.pro (modified) (1 diff, 1 prop)
• q2m_correction_ncdf.pro (modified) (1 diff, 1 prop)
• sst_correction_ncdf.pro (modified) (1 diff, 1 prop)
• t2m_correction_ncdf.pro (modified) (8 diffs, 1 prop)
• ws_correction_ncdf.pro (modified) (1 diff, 1 prop)

trunk/src/TropFlux_19890101_20091231.pro

@@ -50 +50 @@
 ; :ref:`lwr_correction_ncdf.pro`.
 ;
-; :file:`${PROJECT_OD}/TropFlux_t2m_19890101_20091231.nc`
 ; containing
 ; t2m corrected on OAFLUX grid
@@ -110 +109 @@
 ;
 ; :func:`qsee`
-; :func:`calc_claud_vlat` 
+; :func:`lwnet_clark`
+; :func:`calc_claud_vlat`
 ; :func:`cor30a`
+;
+; .. note::
+;
+;    :func:`lwdown_clark` is not used here but available.
 ;
 ; EXAMPLES
@@ -124 +128 @@
 ; ====
 ;
-; missing lwnet_clark::
-;
-;    % Variable is undefined: LWNET_CLARK.
-;    % Execution halted at: TROPFLUX_19890101_20091231  434 /.autofs/home/fplod/incas/tropflux/tropflux_ws/src/TropFlux_19890101_20091231.pr
-;      o
-;    %                      $MAIN$          
-;    % Program caused arithmetic error: Floating overflow
-;    % Program caused arithmetic error: Floating illegal operand
-;
-; check if K or °C in input 
+; make it work ::
+;
+;     % Compiled module: EXTRAPOLATE.
+;     % Attempt to subscript LAND with OK is out of range.
+;     % Execution halted at: EXTRAPOLATE       176
+;       /usr/home/fplod/SAXO_DIR/SRC/Interpolation/extrapolate.pro
+;     %                      TROPFLUX_19890101_20091231  450
+;       /usr/home/fplod/incas/tropflux/tropflux_ws/src/TropFlux_19890101_20091231.p
+;       ro
+;     %                      $MAIN$
+;     % Program caused arithmetic error: Floating overflow
+;     % Program caused arithmetic error: Floating illegal operand
+;
+; the incriminated line is ::
+;
+;     tab=fltarr(jpi,jpj)+!values.f_nan & x=reform(y(*,1)) & tab(ocean)=x & m=fi    nite(tab) & tab=extrapolate(tab,m) & lwr(*,*,jt)=tab*msk+valmask*(1-msk)
+;
+; to make progress I DO NOT extrapolate anymore lwr, sen and lat 
+; !!!!! must be reactivated
+;
+; when I DO NOT extrapolate lwr, sen and lat , program end with ::
+;
+;     Written to 
+;     !/usr/work/incas/fplod/tropflux_d/TropFlux_19890101_20091231_coarev3.nc
+;     -------------------------
+;     % Program caused arithmetic error: Floating divide by 0
+;     % Program caused arithmetic error: Floating underflow
+;     % Program caused arithmetic error: Floating overflow
+;     % Program caused arithmetic error: Floating illegal operand
+;
+; explain why :func:`lwdown_clark` is not used
+;
+; check if K or °C in input
 ;
 ; get rid of::
@@ -143 +170 @@
 ;   da1=19880101
 ;
-; 1st date 19890101 no ? 
+; 1st date 19890101 no ?
 ;
 ; why da1-1 (with da1=19880101) when reading gustiness file ::
@@ -165 +192 @@
 ;
 ; $URL$
+;
+; - fplod 20110822T090838Z aedon.locean-ipsl.upmc.fr (Darwin)
+;
+;    * split some multiple lines statement to investigate
+;      pb of extropalation for lwr
 ;
 ; - fplod 20110819T144332Z aedon.locean-ipsl.upmc.fr (Darwin)
@@ -431 +463 @@
   wn=w(*,*,jt) & wn=wn(ocean)          ; wind speed (m/s)
   ts=sst(*,*,jt) & ts=ts(ocean)        ; Bulk sst (°C)
-  t=t2m(*,*,jt) & t=t(ocean)           ; 2m Air T (°C)
+  ; 2m Air T (°C)
+  t=t2m(*,*,jt) 
+  t=t(ocean)           
   qs=qsee(ts,P)                        ; Sea surface sat. spec. humidity (g/kg)
   q=q2m(*,*,jt) & q=q(ocean)           ; 2m AIr specific humidity  (g/kg)
@@ -449 +483 @@
 
 ; A few punctual missing values (coare does not converge): filled by spatial extrapolation
-  tab=fltarr(jpi,jpj)+!values.f_nan & x=reform(y(*,0)) & tab(ocean)=x & m=finite(tab) & tab=extrapolate(tab,m) & swr(*,*,jt)=tab*msk+valmask*(1-msk)
-  tab=fltarr(jpi,jpj)+!values.f_nan & x=reform(y(*,1)) & tab(ocean)=x & m=finite(tab) & tab=extrapolate(tab,m) & lwr(*,*,jt)=tab*msk+valmask*(1-msk)
-  tab=fltarr(jpi,jpj)+!values.f_nan & x=reform(y(*,2)) & tab(ocean)=x & m=finite(tab) & tab=extrapolate(tab,m) & lat(*,*,jt)=tab*msk+valmask*(1-msk)
-  tab=fltarr(jpi,jpj)+!values.f_nan & x=reform(y(*,3)) & tab(ocean)=x & m=finite(tab) & tab=extrapolate(tab,m) & sen(*,*,jt)=tab*msk+valmask*(1-msk)
-;  tab=fltarr(jpi,jpj)+!values.f_nan & x=reform(lw_clrk) & tab(ocean)=x & m=finite(tab) & tab=extrapolate(tab,m) & lwnet_clrk(*,*,jt)=tab*msk+valmask*(1-msk)
-;  tab=fltarr(jpi,jpj)+!values.f_nan & x=reform(y(*,6)) & tab(ocean)=x & m=finite(tab) & tab=extrapolate(tab,m) & Ch(*,*,jt)=tab*msk+valmask*(1-msk)
-;  tab=fltarr(jpi,jpj)+!values.f_nan & x=reform(y(*,7)) & tab(ocean)=x & m=finite(tab) & tab=extrapolate(tab,m) & Ce(*,*,jt)=tab*msk+valmask*(1-msk)
+  tab=fltarr(jpi,jpj)+!values.f_nan
+  x=reform(y(*,0))
+  tab(ocean)=x
+  m=finite(tab)
+  help, tab, m, swr
+  tab=extrapolate(tab,m)
+  swr(*,*,jt)=tab*msk+valmask*(1-msk)
+;
+  tab=fltarr(jpi,jpj)+!values.f_nan
+  x=reform(y(*,1))
+  tab(ocean)=x
+  m=finite(tab)
+  help, tab, m, lwr
+  ;!!!!!!!!!!!!+++++++++++++++  tab=extrapolate(tab,m)
+  print, 'www : extrapolation for lwr is not done'
+  tab=tab
+  lwr(*,*,jt)=tab*msk+valmask*(1-msk)
+;
+  tab=fltarr(jpi,jpj)+!values.f_nan
+  x=reform(y(*,2))
+  tab(ocean)=x
+  m=finite(tab)
+  help, tab, m, lat
+  ;!!!!!!!!!!!!+++++++++++++++  tab=extrapolate(tab,m)
+  print, 'www : extrapolation for lat is not done'
+  tab=tab
+  lat(*,*,jt)=tab*msk+valmask*(1-msk)
+  tab=fltarr(jpi,jpj)+!values.f_nan
+  x=reform(y(*,3))
+  tab(ocean)=x
+  m=finite(tab)
+  help, tab, m, sen
+  ;!!!!!!!!!!!!+++++++++++++++  tab=extrapolate(tab,m)
+  print, 'www : extrapolation for sen is not done'
+  tab=tab
+  sen(*,*,jt)=tab*msk+valmask*(1-msk)
+  ;  tab=fltarr(jpi,jpj)+!values.f_nan & x=reform(lw_clrk) & tab(ocean)=x & m=finite(tab) & tab=extrapolate(tab,m) & lwnet_clrk(*,*,jt)=tab*msk+valmask*(1-msk)
+  ;  tab=fltarr(jpi,jpj)+!values.f_nan & x=reform(y(*,6)) & tab(ocean)=x & m=finite(tab) & tab=extrapolate(tab,m) & Ch(*,*,jt)=tab*msk+valmask*(1-msk)
+  ;  tab=fltarr(jpi,jpj)+!values.f_nan & x=reform(y(*,7)) & tab(ocean)=x & m=finite(tab) & tab=extrapolate(tab,m) & Ce(*,*,jt)=tab*msk+valmask*(1-msk)
 endfor
 tt=time

trunk/src/TropFlux_NRT_ncdf.pro

@@ -62 +62 @@
 ; ::
 ;
+;  IDL> .compile TropFlux_NRT_ncdf
 ;  IDL> tropflux_nrt_ncdf
 ;
@@ -67 +68 @@
 ; ====
 ;
-; no yet tested because pb in TropFlux_19890101_20091231.pro
-;
 ; coding rules
+;
+; get rid of uppercase to avoid compile
 ;
 ; create 4 files instead of one now
@@ -76 +77 @@
 ; ==========
 ;
-; $Id$
+; $Id: TropFlux_NRT_ncdf.pro 88 2011-08-19 15:40:14Z pinsard $
 ; 
 ; $URL$
@@ -212 +213 @@
 lon=reform(glamt(0:jpi-1,0))
 
-ncfile='!/Users/pkb/data/TropFlux/TropFlux_19890101_20091231.nc'
+ncfile='!' + fullfilename_out
 lon_attr={units:'degrees_east',long_name:'Longitude'}
 lat_attr={units:'degrees_north',long_name:'Latitude'}

trunk/src/TropFlux_swr_BLND_19890101_20091231.pro

@@ -100 +100 @@
 ; ==========
 ;
-; $Id$
+; $Id: TropFlux_swr_BLND_19890101_20091231.pro 88 2011-08-19 15:40:14Z pinsard $
 ;
 ; $URL$

trunk/src/TropFlux_swr_DT_19890101_20071231.pro

@@ -87 +87 @@
 ; ==========
 ;
-; $Id$
+; $Id: TropFlux_swr_DT_19890101_20071231.pro 88 2011-08-19 15:40:14Z pinsard $
 ;
 ; $URL$

trunk/src/TropFlux_swr_NRT_19890101_20091231.pro

@@ -89 +89 @@
 ; ==========
 ;
-; $Id$
+; $Id: TropFlux_swr_NRT_19890101_20091231.pro 88 2011-08-19 15:40:14Z pinsard $
 ;
 ; $URL$

trunk/src/cor30a.pro

@@ -71 +71 @@
 ; used by :ref:`TropFlux_19890101_20091231.pro`
 ;
+; call :func:`psiu`, :func:`psit`.
+;
 ; TODO
 ; ====
@@ -86 +88 @@
 ;
 ; $URL$
+;
+; - fplod 20110822T085404Z aedon.locean-ipsl.upmc.fr (Darwin)
+;
+;   * add ref to psiu and psit called her and provided by pk today
 ;
 ; - fplod 20101214T093615Z aedon.locean-ipsl.upmc.fr (Darwin)

trunk/src/interp_erai_msl_1989_2009.pro

@@ -89 +89 @@
 ; ==========
 ;
-; $Id$
+; $Id: interp_erai_msl_1989_2009.pro 88 2011-08-19 15:40:14Z pinsard $
 ;
 ; $URL$

trunk/src/interp_olr_30n30s_1989_2009.pro

@@ -106 +106 @@
 ; ==========
 ;
-; $Id$
+; $Id: interp_olr_30n30s_1989_2009.pro 88 2011-08-19 15:40:14Z pinsard $
 ;
 ; $URL$

trunk/src/lwr_correction_ncdf.pro

@@ -82 +82 @@
 ; ==========
 ;
-; $Id$
+; $Id: lwr_correction_ncdf.pro 88 2011-08-19 15:40:14Z pinsard $
 ;
 ; $URL$

trunk/src/q2m_correction_ncdf.pro

@@ -73 +73 @@
 ; ==========
 ;
-; $Id$
+; $Id: q2m_correction_ncdf.pro 88 2011-08-19 15:40:14Z pinsard $
 ;
 ; $URL$

trunk/src/sst_correction_ncdf.pro

@@ -131 +131 @@
 ; ==========
 ;
-; $Id$
+; $Id: sst_correction_ncdf.pro 88 2011-08-19 15:40:14Z pinsard $
 ;
 ; $URL$

trunk/src/t2m_correction_ncdf.pro

@@ -10 +10 @@
 ; applied.
 ;
-; :file:`${PROJECT_ID}/erai_t2m_19890101_20091231_oafluxgrid.nc` have been produced by :ref:`interp_erai_t2m_1989_2009.pro`.
+; :file:`${PROJECT_ID}/erai_t2m_19890101_20091231_oafluxgrid.nc`
+; has been produced by
+; :ref:`interp_erai_t2m_1989_2009.pro`.
 ; It contains air temperature at 2 m height from ERA-I interpolated on OAFLUX grid.
 ;
-; Corrected air temperature at 2 m height is written in
+; Corrected air temperature at 2 m height
+; is written in
 ; :file:`${PROJECT_OD}/TropFlux_t2m_19890101_20091231.nc`
 ; if this file not already exists.
@@ -22 +25 @@
 ;
 ;        digraph t2m_correction_ncdf {
-;           graph [
-;           rankdir="LR",
-;           ]
 ;
 ;           file_in [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/erai_t2m_19890101_20091231_oafluxgrid.nc"];
@@ -67 +67 @@
 ; ====
 ;
+; note .5 in time value ... is this correct ? do others files
+; present the same feature ?::
+:
+;    $ ncks -v time -d time,0,1 /usr/work/incas/fplod/tropflux_d/TropFlux_t2m_19890101_20091231.nc
+;    time: type NC_DOUBLE, 1 dimension, 3 attributes, chunked? no, compressed? no, packed? no, ID = 3
+;    time RAM size is 7670*sizeof(NC_DOUBLE) = 7670*8 = 61360 bytes
+;    time dimension 0: time, size = 7670 NC_DOUBLE, dim. ID = 2 (CRD)(REC)
+;    time attribute 0: units, size = 30 NC_CHAR, value = days since 1950-01-01 00:00:00
+;    time attribute 1: long_name, size = 9 NC_CHAR, value = Time axis
+;    time attribute 2: time_origin, size = 20 NC_CHAR, value = 1950-JAN-01 00:00:00
+;
+;    time[0]=14244.5 days since 1950-01-01 00:00:00
+;    time[1]=14245.5 days since 1950-01-01 00:00:00
+;
 ; coding rules
 ;
@@ -94 +108 @@
 ; ==========
 ;
-; $Id$
+; $Id: t2m_correction_ncdf.pro 88 2011-08-19 15:40:14Z pinsard $
 ;
 ; $URL$
+;
+; - fplod 20110822T112525Z aedon.locean-ipsl.upmc.fr (Darwin)
+;
+;   * correction for time values in output file
 ;
 ; - fplod 20110808T130628Z aedon.locean-ipsl.upmc.fr (Darwin)
@@ -170 +188 @@
 ;
 ; build uncorrected t2 filename
-filename_t2_uncor='erai_t2m_'+ string(da1,format='(i8.8)')+'_'+ string(da2,format='(i8.8)')+'_oafluxgrid.nc'
+filename_t2_uncor = 'erai_t2m_' $
+    + string(da1,format='(i8.8)') $
+    + '_' $
+    + string(da2,format='(i8.8)') $
+    + '_oafluxgrid.nc'
 ;
 ; check if this file exists
@@ -181 +203 @@
 ;
 ; build output filename
-filename_out='TropFlux_t2m_'+ string(da1,format='(i8.8)')+'_'+ string(da2,format='(i8.8)')+'.nc'
+filename_out = 'TropFlux_t2m_' $
+    + string(da1,format='(i8.8)') $
+    + '_' $
+    + string(da2,format='(i8.8)') $
+    + '.nc'
+;
 fullfilename_out = iodirout + filename_out
 ; in order to avoid unexpected overwritten
@@ -256 +283 @@
 cda1=string(da2)
 
-time=time-julday(1,1,1950)
+% conversion for hours since 19570101 to days from 19500101
+time=julday(01, 01, 1957,timein)-julday(1,1,1950)
 jpt=n_elements(time)
 
@@ -275 +303 @@
 globattr={source:'Basic data obtained from ERAI.  Mean correction for air temperautre bias and correction for variability are applied',timerange:cda0+' - '+cda1}
 
-
 ncfields = 't2m[longitude,latitude,time]=t2m_new:t2m_attr; ' $
                       + 'longitude[]=lon:lon_attr; ' $

trunk/src/ws_correction_ncdf.pro

@@ -77 +77 @@
 ; ==========
 ;
-; $Id$
+; $Id: ws_correction_ncdf.pro 88 2011-08-19 15:40:14Z pinsard $
 ;
 ; $URL$