Changeset 204 for trunk/src/flux_evaluation_tpr_timeseries.pro

Timestamp:

03/26/14 11:58:56 (10 years ago)

Author:

pinsard

Message:

fix thanks to coding rules; typo

File:

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

trunk/src/flux_evaluation_tpr_timeseries.pro

@@ -41 +41 @@
 ;
 ;   * try to make it work on my account on cratos
-;   * remove x_site_location and  x_site_location because already exist
+;   * 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
-                    date1,  $ ;;  start date (in julian date. eg. 20000101)
-                    date2     ;;  end date (in julian date eg. 20091231)
+date1, $ ;;  start date (in Julian date. eg. 20000101)
+date2    ;;  end date (in Julian date eg. 20091231)
 
 @common
@@ -110 +110 @@
 for n=0, nn-1 do begin
 
-;; reading data from mooring
+    ;; reading data from mooring
 
     site=sitelist(n) & csite=site
@@ -122 +122 @@
             sw,lw,sh,lh
 
-;; select the appropriate variables for evaluation (trp = sw or lw or sh or lh)
+    ;; select the appropriate variables for evaluation (trp = sw or lw or sh or lh)
 
-   tpr=sh
+    tpr=sh
 
-   ind=where(finite(tpr)) & no_valid=n_elements(ind)
+    ind=where(finite(tpr)) & no_valid=n_elements(ind)
 
-   if (no_valid ge min_obs) then begin
+    if (no_valid ge min_obs) then begin
+        extract_tpr_location,var,box, $
+        var_tpr
+        var_tpr=reform(var_tpr)
 
+        stats_5d, tpr,var_tpr, $  ;; tpr=TPR observation and var_tpr=gridded product extracted at TPR location
+        cor, bias, std, rmsd
 
-       extract_tpr_location,var,box, $
-       var_tpr
-       var_tpr=reform(var_tpr)
+        printf, 1, x, y, cor, bias, std, rmsd, format='(f6.2, 3x, f6.2, 3x, f4.2, 3x, f7.2, 3x, f4.2, 3x, f5.2)'
+        cstat=string(cor, bias, std, rmsd, format='(f4.2,3x,f7.2,3x,f4.2,3x,f5.2)')
+        print, cstat
 
-       stats_5d, tpr,var_tpr, $  ;; tpr=TPR observation and var_tpr=gridded product extracted at TPR location
-       cor, bias, std, rmsd
-
-       printf, 1, x, y, cor, bias, std, rmsd, format='(f6.2, 3x, f6.2, 3x, f4.2, 3x, f7.2, 3x, f4.2, 3x, f5.2)'
-       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.
-       pltt, ts_smooth(tpr,5,/nan), "t",/rempl, small=[1,3,1], lct=65, $
-       title='Five day stats are shown below.  TPR (black) and Product (red) at'+csite+' ', charsize=1., $
-       subtitle=cstat
-       ind=where(finite(tpr,/nan)) & var_tpr(ind)=!Values.f_nan
-       pltt, ts_smooth(var_tpr,5,/nan), "t",/ov1d, color=250
-       erase
-       mooring=[mooring,tpr] & product=[product,var_tpr]
+        ;;     PLOTTING THE TIME-SERIES
+        array=[tpr, var_tpr] & mi=min(array,/nan) & ma=max(array,/nan) & int=(ma-mi)/3.
+        pltt, ts_smooth(tpr,5,/nan), "t",/rempl, small=[1,3,1], lct=65, $
+        title='Five day stats are shown below.  TPR (black) and Product (red) at'+csite+' ', charsize=1., $
+        subtitle=cstat
+        ind=where(finite(tpr,/nan)) & var_tpr(ind)=!Values.f_nan
+        pltt, ts_smooth(var_tpr,5,/nan), "t",/ov1d, color=250
+        erase
+        mooring=[mooring,tpr] & product=[product,var_tpr]
     endif
 endfor