Changeset 43 for trunk/src/paper01/fig7/ws_validation_scatter_2000_2009_basin.pro

Timestamp:

04/11/11 18:05:51 (13 years ago)

Author:

pinsard

Message:

add minimal header in paper01 IDL files

File:

• trunk/src/paper01/fig7/ws_validation_scatter_2000_2009_basin.pro (modified) (15 diffs)

trunk/src/paper01/fig7/ws_validation_scatter_2000_2009_basin.pro

@@ -1 @@
-;------------------------------------------------------------
+;+
+;
+; .. _ws_validation_scatter_2000_2009_basin.pro
+;
+; =========================================
+; ws_validation_scatter_2000_2009_basin.pro
+; =========================================
+;
+; DESCRIPTION
+; ===========
+;
+; SEE ALSO
+; ========
+;
+; :ref:`tropflux_profile.sh`
+;
+; EXAMPLES
+; ========
+;
+; ::
+;
+;  IDL> @tropflux_init
+;  IDL> date1=++
+;  IDL> date2=++
+;  IDL> ws_validation_scatter_2000_2009_basin, date1, date2
+;
+; EVOLUTIONS
+; ==========
+;
+; $Id$
+;
+; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin)
+;
+;   * minimal header
+;
+;-
 pro ws_validation_scatter_2000_2009_basin,date1,date2
 @common
-;------------------------------------------------------------
 
 reinitplt, /z,/invert
@@ -8 +42 @@
 
 openps, FILENAME = 'idl.ps'
-;------------------------------------------------------------
 ;; Before running this program, you have to compile the following subrutines
 ;;   1.  .r read_era_total
@@ -49 +82 @@
 sitelist=['8s125w', '8s110w', '8s95w','5s155w', '5s140w', '5s125w', '5s110w', '5s95w','2s155w', '2s140w', $
           '2s125w', '2s110w', '2s95w', '0n155w', '0n140w', '0n125w', '0n110w', '0n95w', '2n155w', '2n140w', $
-          '2n125w', '2n110w', '2n95w',  '5n155w', '5n140w', '5n125w', '5n110w', '5n95w'] 
+          '2n125w', '2n110w', '2n95w',  '5n155w', '5n140w', '5n125w', '5n110w', '5n95w']
 ocean='Pac_coldtongue'
 
@@ -65 +98 @@
 
 nsmooth=1
-;------------------------------------------------------------------------------------------------------------------------
 ;;   This program will create the following text files with statistics of respective variables
-;------------------------------------------------------------------------------------------------------------------------
 close,/all
 
@@ -91 +122 @@
 printf,6, 'x     y      cor    bias     std     rmsd      mean_tao'
 
-;------------------------------------------------------------------------------------------------------------------------
 ;; first reading the whole ERAI uncorrected and corrected data
 ;; date1=20000101 & date2=20090931
@@ -97 +127 @@
 file='/Volumes/Iomega_HDD/TropFlux/input_uncor/erai_ws_19890101_20091231_oafluxgrid.nc'
 initncdf, file
-u=read_ncdf('u10',date1,date2,file=file,/nostr) 
-;v=read_ncdf('v10',date1,date2,file=file,/nostr)     
-;unc=sqrt(u*u+v*v) 
+u=read_ncdf('u10',date1,date2,file=file,/nostr)
+;v=read_ncdf('v10',date1,date2,file=file,/nostr)
+;unc=sqrt(u*u+v*v)
 ;help, unc
 
 ;file='/Volumes/Iomega_HDD/TropFlux/input_cor/full_cor/TropFlux_ws_19890101_20091231_v20.nc'
 ;initncdf, file
-;cor=read_ncdf('ws',date1,date2,file=file,/nostr) 
+;cor=read_ncdf('ws',date1,date2,file=file,/nostr)
 ;help, cor
 
@@ -114 +144 @@
 ;fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/uwind_ncep2_oafluxgrid_19890101_20091231.nc'
 ;initncdf, fi
-;u=read_ncdf("u", date1-1, date2, file=fi,/nostr) 
+;u=read_ncdf("u", date1-1, date2, file=fi,/nostr)
 ;fi='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/vwind_ncep2_oafluxgrid_19890101_20091231.nc'
 ;initncdf, fi
-;v=read_ncdf("v", date1-1, date2, file=fi,/nostr)   
+;v=read_ncdf("v", date1-1, date2, file=fi,/nostr)
 ;nce=sqrt(u*u+v*v)
 ;help, nce
@@ -142 +172 @@
 help, unc, cor, oaf, nce, ws_tmi, nce1
 
-;------------------------------------------------------------------------------------------------------------------------
 nn=n_elements(sitelist)
 date1=date1
@@ -168 +197 @@
     uncr=tropflux
 
-     extract_flux_tropflux,cor,box, $  
+     extract_flux_tropflux,cor,box, $
         tropflux
      corr=tropflux
-    
+
      extract_flux_tropflux,oaf,box, $
         tropflux
@@ -188 +217 @@
      ncep1=tropflux
 
-    ind=where(finite(ws)) & ws=ws(ind) & uncr_ws=uncr(ind)  & corr_ws=corr(ind) 
+    ind=where(finite(ws)) & ws=ws(ind) & uncr_ws=uncr(ind)  & corr_ws=corr(ind)
     oafl=oafl(ind) & ncep=ncep(ind) & tmi=tmi(ind) & ncep1=ncep1(ind)
 
     mean_tao=total(ws)/n_elements(ws)
- 
+
     statistics_3var_v1, ws, uncr_ws, corr_ws, $
          cor1, cor2, bias1, bias2, std1, std2, rmsd1, rmsd2
- 
+
     printf, 1, x, y, cor1, bias1, std1, rmsd1, mean_tao, format='(f6.2, 3x, f6.2, 3x, f5.2,3x,f5.2,3x,f4.2,3x,f4.2,3x,f5.2)'
     printf, 2, x, y, cor2, bias2, std2, rmsd2, mean_tao, format='(f6.2, 3x, f6.2, 3x, f5.2,3x,f5.2,3x,f4.2,3x,f4.2,3x,f5.2)'
@@ -210 +239 @@
 
 endfor
-;-----------------------------------------------------------------------------------------
-close,/all  
+close,/all
 fi_ws_erai='/Users/pkb/work/MY_SAXO/flux_automat/ws_2000_2009_erai_'+ocean+'.txt'
 res=read_ascii(fi_ws_erai,data_start=1)
@@ -376 +404 @@
 oplot, mean_tao, yfit, color=250, thick=2
 
-;-----------------------------------------------------------
 closeps
 
@@ -384 +411 @@
 return
 end
-;--------------------------------------------------------------------------
 function x_site_location, site
     n1=strpos(site, 's')
@@ -399 +425 @@
 return, float(x)
 end
-;--------------------------------------------------------------------------
 function y_site_location, site
     n1=strpos(site, 'e')
@@ -425 +450 @@
 return,float(y)
 end
-
-;--------------------------------------------------------------------------
-