Changeset 107 for trunk/src/flux_evaluation_tpr_timeseries.pro

Timestamp:

11/30/11 15:53:12 (12 years ago)

Author:

pinsard

Message:

first diag on pkb tools for timeseries

File:

• trunk/src/flux_evaluation_tpr_timeseries.pro (modified) (7 diffs)

trunk/src/flux_evaluation_tpr_timeseries.pro

@@ -1 @@
-;------------------------------------------------------------
+;+
+;
+; .. _flux_evaluation_tpr_timeseries.pro:
+;
+; ==================================
+; flux_evaluation_tpr_timeseries.pro
+; ==================================
+;
+; DESCRIPTION
+; ===========
+;
+; SEE ALSO
+; ========
+;
+; :ref:`project_profile.sh`
+; :ref:`project_init.pro`
+; :ref:`cm_project.pro`
+;
+; TODO
+; ====
+;
+; make it work : missing READ_TPR_NETFLUX
+;
+; EXAMPLES
+; ========
+;
+; ::
+;
+;   IDL> date1 = 20000101L
+;   IDL> date2 = 20091231L
+;   IDL> flux_evaluation_tpr_timeseries, date1, date2
+;
+; EVOLUTIONS
+; ==========
+;
+; $Id$
+;
+; $URL$
+;
+; - fplod 20111130T141341Z cratos (Linux)
+;
+;   * try to make it work on my account on cratos
+;   * remove x_site_location and  x_site_location because already exist
+;
+;-
 pro flux_evaluation_tpr_timeseries, $
 ;                    var,    $ ;;  flux variable (swr, lwr, lhf, shf) to calculate the statistics
@@ -6 +50 @@
 
 @common
-;------------------------------------------------------------
+@cm_project
 reinitplt, /z,/invert
 key_portrait = 0
-openps, FILENAME = 'idl.ps'
+openps, FILENAME = project_od_env+'flux_evaluation_tpr_timeseries.ps'
 ;------------------------------------------------------------
 ;; part to change
@@ -45 +89 @@
 close,/all
 
-fi='/Users/pkb/work/MY_SAXO/TropFlux_update/flux_stat.txt'
+fi=project_id_env + 'flux_stat.txt'
 openw,1,fi
 
@@ -52 +96 @@
 ;; First, this program reads the full TropFlux data and later extract it at specific TPR locations
 
-file="/Volumes/PAYASAM/TropFlux/TropFlux/shf_tropflux_1d_1989_2010.nc"
+file = project_id_env + "shf_tropflux_1d_1989_2010.nc"
 initncdf, file
 var=-1*read_ncdf("shf", date1, date2, file=file,/nostr)
@@ -74 +118 @@
     if (y ge 0. and y le 30.) then y=y+360.
     dx=0.5 & dy=0.5 & box=[y-dy, y+dy, x-dx, x+dx]
-     
+
     read_tpr_netflux, csite,date1,date2,nsmooth, $
             sw,lw,sh,lh
@@ -97 +141 @@
            cstat=string(cor, bias, std, rmsd, format='(f4.2,3x,f7.2,3x,f4.2,3x,f5.2)')
            print, cstat
- 
+
 ;;         PLOTTING THE TIME-SERIES
            array=[tpr, var_tpr] & mi=min(array,/nan) & ma=max(array,/nan) & int=(ma-mi)/3.
@@ -111 +155 @@
 
 close,/all
-;----------------------------------------------------------
 closeps
 
-fig='flux_evaluation_tpr_timeseries.ps'
-spawn, 'mv '+psdir+'idl.ps '+updatedir+fig
-spawn, 'gv '+updatedir+fig
-return
 end
-;--------------------------------------------------------------------------
-function x_site_location, site
-    n1=strpos(site, 's')
-if (n1 gt -1) then begin
-    ns=-1.
-    x=strmid(site, 0, n1)
-    x=float(x)*ns
-endif else begin
-    n1=strpos(site, 'n')
-    x=strmid(site, 0, n1)
-    ny=1.
-    x=float(x)*ny
-endelse
-return, float(x)
-end
-;--------------------------------------------------------------------------
-function y_site_location, site
-    n1=strpos(site, 'e')
-if (n1 gt -1) then begin
-    n=strpos(site, 's')
-    if (n gt -1) then begin
-        y=strmid(site, n+1, n1-n-1)
-    endif else begin
-        n=strpos(site, 'n')
-        y=strmid(site, n+1, n1-n-1)
-    endelse
-
-endif else begin
-    n1=strpos(site, 'w')
-    n=strpos(site, 's')
-    if (n gt -1) then begin
-        y=strmid(site, n+1, n1-n-1)
-        y=180+(180-float(y))
-    endif else begin
-        n=strpos(site, 'n')
-        y=strmid(site, n+1, n1-n-1)
-        y=180+(180-float(y))
-   endelse
-endelse
-return,float(y)
-end
-
-;--------------------------------------------------------------------------
-