Changeset 43 for trunk/src/paper01/fig13

Timestamp:

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

Author:

pinsard

Message:

add minimal header in paper01 IDL files

File:

• trunk/src/paper01/fig13/net_flux_validation_scatter_2000_2007.pro (modified) (13 diffs)

trunk/src/paper01/fig13/net_flux_validation_scatter_2000_2007.pro

@@ -1 @@
-;------------------------------------------------------------
-pro net_flux_validation_scatter_2000_2007,date1,date2
+;+
+; .. _net_flux_validation_scatter_2000_2007.pro
+;
+; =========================================
+; net_flux_validation_scatter_2000_2007.pro
+; =========================================
+;
+; DESCRIPTION
+; ===========
+;
+; SEE ALSO
+; ========
+;
+; :ref:`tropflux_profile.sh`
+;
+; EXAMPLES
+; ========
+;
+; ::
+;
+;  IDL> @tropflux_init
+;  IDL> date1=++
+;  IDL> date2=++
+;  IDL> net_flux_validation_scatter_2000_2007, date1, date2
+;
+; EVOLUTIONS
+; ==========
+;
+; $Id$
+;
+; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin)
+;
+;   * minimal header
+;
+;-
+pro net_flux_validation_scatter_2000_2007, date1, date2
 @common
-;------------------------------------------------------------
 
 reinitplt, /z,/invert
@@ -8 +41 @@
 
 openps, FILENAME = 'idl.ps'
-;------------------------------------------------------------
 
 sitelist=['8s67e','12s55e', '8s55e', '8s80.5e', '1.5s80.5e', '0n80.5e', '1.5n80.5e', '1.5s90e', $
@@ -27 +59 @@
 
 ocean='net_flux'
-nsmooth=1. 
-;------------------------------------------------------------------------------------------------------------------------
+nsmooth=1.
 ;;   This program will create the following text files with statistics of respective variables
-;------------------------------------------------------------------------------------------------------------------------
 file_trop='/Users/pkb/data/TropFlux/TropFlux_19890101_20091231.nc'
 initncdf, file_trop
 shf_trop=read_ncdf("shf", date1, date2, file=file_trop,/nostr)
 lhf_trop=read_ncdf("lhf", date1, date2, file=file_trop,/nostr)
-lwr_trop=read_ncdf("lwr", date1, date2, file=file_trop,/nostr) 
+lwr_trop=read_ncdf("lwr", date1, date2, file=file_trop,/nostr)
 
 file_trop='/Users/pkb/data/TropFlux/TropFlux_swr_19890101_20071231_DT_v50.nc'
@@ -71 +101 @@
 file_oaf='/Volumes/Iomega_HDD/work/flux_reconstruction/OAFLX_GRID/shf_oafluxgrid_1985_2008.nc'
 initncdf, file_oaf
-shf=read_ncdf("shf", date1, date2, file=file_oaf,/nostr) 
+shf=read_ncdf("shf", date1, date2, file=file_oaf,/nostr)
 file_oaf='/Volumes/Iomega_HDD/work/flux_reconstruction/OAFLX_GRID/lhf_oafluxgrid_1985_2008.nc'
 initncdf, file_oaf
-lhf=read_ncdf("lhf", date1, date2, file=file_oaf,/nostr) 
+lhf=read_ncdf("lhf", date1, date2, file=file_oaf,/nostr)
 file_oaf='/Volumes/Iomega_HDD/work/flux_reconstruction/OAFLX_GRID/lwr_oafluxgrid_1985_2007.nc'
 initncdf, file_oaf
-lwr=read_ncdf("lwr", date1, date2, file=file_oaf,/nostr) 
+lwr=read_ncdf("lwr", date1, date2, file=file_oaf,/nostr)
 file_oaf='/Volumes/Iomega_HDD/work/flux_reconstruction/OAFLX_GRID/swr_oafluxgrid_1985_2007.nc'
 initncdf, file_oaf
@@ -94 +124 @@
 help, ncep1
 
-;------------------------------------------------------------------------------------------------------------------------
 close,/all
 
@@ -118 +147 @@
 printf,6, 'x     y      cor    bias     std     rmsd      mean_tao'
 
-;------------------------------------------------------------------------------------------------------------------------
 nn=n_elements(sitelist)
 
@@ -133 +161 @@
     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_variables_net_flux, csite,date1,date2,nsmooth, $
          net
     tao=net & ind=where(finite(tao)) & valid=n_elements(ind)
 
-;;    
+;;
     if (valid ge 180. ) then begin
 
@@ -155 +183 @@
            extract_flux_tropflux, trop, box, $
               tropflux
-           tropflux_new=tropflux 
+           tropflux_new=tropflux
 
            extract_flux_tropflux, ncep1, box, $
@@ -165 +193 @@
            tropdf=tropflux
 
-           ind=where(finite(tao)) & tao=tao(ind) & oaf_net=oaflux(ind) & nce_net=ncepflux(ind) & trop_net=tropflux_new(ind) 
+           ind=where(finite(tao)) & tao=tao(ind) & oaf_net=oaflux(ind) & nce_net=ncepflux(ind) & trop_net=tropflux_new(ind)
            erai_net=eraiflux(ind) & nce1_net=ncepflux1(ind) & tropd_net=tropdf(ind)
 
@@ -383 +411 @@
 oplot, [-100,200], [-100,200]
 
-;-----------------------------------------------------------
 closeps
 
@@ -391 +418 @@
 return
 end
-;--------------------------------------------------------------------------
 function x_site_location, site
     n1=strpos(site, 's')
@@ -406 +432 @@
 return, float(x)
 end
-;--------------------------------------------------------------------------
 function y_site_location, site
     n1=strpos(site, 'e')
@@ -432 +457 @@
 return,float(y)
 end
-
-;--------------------------------------------------------------------------
-