Changeset 203 for trunk/src/flux_evaluation_tpr_timeseries.pro
- Timestamp:
- 02/11/14 18:41:24 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/flux_evaluation_tpr_timeseries.pro
r180 r203 25 25 ; ======== 26 26 ; 27 ; ::27 ; .. code-block:: idl 28 28 ; 29 ; date1 = 20000101L30 ; date2 = 20091231L31 ; flux_evaluation_tpr_timeseries, date1, date229 ; date1 = 20000101L 30 ; date2 = 20091231L 31 ; flux_evaluation_tpr_timeseries, date1, date2 32 32 ; 33 33 ; EVOLUTIONS … … 131 131 132 132 133 134 135 133 extract_tpr_location,var,box, $ 134 var_tpr 135 var_tpr=reform(var_tpr) 136 136 137 138 137 stats_5d, tpr,var_tpr, $ ;; tpr=TPR observation and var_tpr=gridded product extracted at TPR location 138 cor, bias, std, rmsd 139 139 140 141 142 140 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)' 141 cstat=string(cor, bias, std, rmsd, format='(f4.2,3x,f7.2,3x,f4.2,3x,f5.2)') 142 print, cstat 143 143 144 ;; 145 146 147 148 149 150 151 152 144 ;; PLOTTING THE TIME-SERIES 145 array=[tpr, var_tpr] & mi=min(array,/nan) & ma=max(array,/nan) & int=(ma-mi)/3. 146 pltt, ts_smooth(tpr,5,/nan), "t",/rempl, small=[1,3,1], lct=65, $ 147 title='Five day stats are shown below. TPR (black) and Product (red) at'+csite+' ', charsize=1., $ 148 subtitle=cstat 149 ind=where(finite(tpr,/nan)) & var_tpr(ind)=!Values.f_nan 150 pltt, ts_smooth(var_tpr,5,/nan), "t",/ov1d, color=250 151 erase 152 mooring=[mooring,tpr] & product=[product,var_tpr] 153 153 endif 154 154 endfor
Note: See TracChangeset
for help on using the changeset viewer.