Changeset 46 for trunk/src/paper01/fig15

Timestamp:

04/20/11 16:45:53 (13 years ago)

Author:

pinsard

Message:

going on consolidation of paper01 materials

File:

• trunk/src/paper01/fig15/fig15_timeseries_mjo.pro (modified) (11 diffs)

trunk/src/paper01/fig15/fig15_timeseries_mjo.pro

@@ -6 +6 @@
 ; ========================
 ;
-;
 ; DESCRIPTION
 ; ===========
 ;
+; .. graphviz::
+;
+;    digraph fig15_timeseries_mjo {
+;       graph [
+;       rankdir="LR",
+;       ]
+;       tropflux_general [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/TropFlux_19890101_20091231.nc"];
+;       ncep2 [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/NCEP2_flux_19890101_20090729.nc"];
+;       erai [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/erai_fluxes_19930101_20090801_TROP_oafluxgrid.nc"];
+;       ncep1 [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/fluxe_ncep1_19890101_20091231.nc"];
+;       shf_oafluxgrid [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/shf_oafluxgrid_1985_2008.nc"];
+;       lhf_oafluxgrid [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/lhf_oafluxgrid_1985_2008.nc"];
+;       lwr_oafluxgrid [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/lwr_oafluxgrid_1985_2007.nc"];
+;       swr_oafluxgrid [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/swr_oafluxgrid_1985_2007.nc"];
+;       olr_oafluxgrid [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/olr_oafluxgrid_30n30s_19890101_20091231.nc"];
+;
+;       figure [shape=ellipse,fontname=Courier,label="${TROPFLUX_OD}/WHOTS_validation_net_2004_2007_all_v10.ps"];
+;
+;       fig15_timeseries_mjo [shape=box,
+;       fontname=Courier,
+;       color=blue,
+;       URL="http://forge.ipsl.jussieu.fr/tropflux/broswrer/trunk/src/paper01/fig15/fig15_timeseries_mjo.pro",
+;       label="${TROPFLUX}/src/paper01/fig15/fig15_timeseries_mjo.pro"];
+;
+;       {tropflux_general swr_oafluxgrid lwr_oafluxgrid lhf_oafluxgrid shf_oafluxgrid ncep2 ncep1} -> {fig15_timeseries_mjo} -> {figure}
+;    }
+;
 ; SEE ALSO
 ; ========
 ;
 ; :ref:`tropflux_profile.sh`
+; :ref:`tropflux_init.pro`
+; :ref:`cm_tropflux.pro`
 ;
 ; EXAMPLES
@@ -23 +51 @@
 ;  IDL> fig15_timeseries_mjo
 ;
+; TODO
+; ====
+;
+; make it work on cratos : missing data
+;
+; coding rules
+;
+; get rid of spwan
+;
+; complete description
+;
+; handle IO error
+;
 ; EVOLUTIONS
 ; ==========
@@ -28 +69 @@
 ; $Id$
 ;
+; - fplod 20110420T122720Z aedon.locean-ipsl.upmc.fr (Darwin)
+;
+;   * remove hard coding path
+;   * add graphviz
+;   * externalize functions
+;
 ; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin)
 ;
@@ -34 +81 @@
 ;-
 pro fig15_timeseries_mjo
-@common
-
+@cm_general
+@cm_tropflux
 reinitplt, /z,/invert
 key_portrait = 1
@@ -44 +91 @@
 box=[50,80,-12,-5]
 
-file_trop='/Users/pkb/data/TropFlux/TropFlux_19890101_20091231.nc'
+file_trop=tropflux_id_env+'TropFlux_19890101_20091231.nc'
 initncdf, file_trop
 shf_trop=read_ncdf("shf", date1, date2, box=box, file=file_trop,/nostr)
@@ -54 +101 @@
 trop=grossemoyenne(trop, "xy",/nan)
 
-file_ncep='/Volumes/Iomega_HDD/flux_reconstruction/ncep2/NCEP2_flux_19890101_20090729.nc'
+file_ncep=tropflux_id_env+'NCEP2_flux_19890101_20090729.nc'
 initncdf, file_ncep
 shf=read_ncdf("shf", date1, date2, box=box, file=file_ncep,/nostr)
@@ -65 +112 @@
 help, nce2
 
-file_erai='/Volumes/Iomega_HDD/work/flux_reconstruction/gridded_data/erai_fluxes_19930101_20090801_TROP_oafluxgrid.nc'
+file_erai=tropflux_id_env+'erai_fluxes_19930101_20090801_TROP_oafluxgrid.nc'
 initncdf, file_erai
 shf=read_ncdf("shf", date1, date2, box=box, file=file_erai,/nostr)
@@ -75 +122 @@
 help, era
 
-file='/Volumes/Iomega_HDD/flux_reconstruction/ncep/fluxe_ncep1_19890101_20091231.nc'
+file=tropflux_id_env+'fluxe_ncep1_19890101_20091231.nc'
 initncdf, file
 lwr=read_ncdf("lwr", date1, date2, box=box, file=file,/nostr)
@@ -86 +133 @@
 help, nce
 
-file_oaf='/Volumes/Iomega_HDD/work/flux_reconstruction/OAFLX_GRID/shf_oafluxgrid_1985_2008.nc'
+file_oaf=tropflux_id_env+'shf_oafluxgrid_1985_2008.nc'
 initncdf, file_oaf
 shf=read_ncdf("shf", date1, date2, box=box, file=file_oaf,/nostr)
-file_oaf='/Volumes/Iomega_HDD/work/flux_reconstruction/OAFLX_GRID/lhf_oafluxgrid_1985_2008.nc'
+file_oaf=tropflux_id_env+'lhf_oafluxgrid_1985_2008.nc'
 initncdf, file_oaf
 lhf=read_ncdf("lhf", date1, date2, box=box, file=file_oaf,/nostr)
-file_oaf='/Volumes/Iomega_HDD/work/flux_reconstruction/OAFLX_GRID/lwr_oafluxgrid_1985_2007.nc'
+file_oaf=tropflux_id_env+'lwr_oafluxgrid_1985_2007.nc'
 initncdf, file_oaf
 lwr=read_ncdf("lwr", date1, date2, box=box, file=file_oaf,/nostr)
-file_oaf='/Volumes/Iomega_HDD/work/flux_reconstruction/OAFLX_GRID/swr_oafluxgrid_1985_2007.nc'
+file_oaf=tropflux_id_env+'swr_oafluxgrid_1985_2007.nc'
 initncdf, file_oaf
 swr=read_ncdf("swr", date1, date2, box=box, file=file_oaf,/nostr) & swr_trop=swr
@@ -102 +149 @@
 oafl=grossemoyenne(oafl, "xy",/nan)
 
-file="/Users/pkb/data/heat_budget/NOAA_OLR/olr_oafluxgrid_30n30s_19890101_20091231.nc"
+file=tropflux_id_env+"olr_oafluxgrid_30n30s_19890101_20091231.nc"
 initncdf, file
 olr=read_ncdf("olr", date1-1, date2, box=box, file=file,/nostr)
@@ -161 +208 @@
 
 closeps
-fig='fig15_timeseries_mjo.ps'
-spawn, 'mv '+psdir+'idl.ps '+cpsdir+fig
-spawn, 'gv '+cpsdir+fig
-return
+fig=tropflux_od_env+'fig15_timeseries_mjo.ps'
+spawn, 'mv '+psdir+'idl.ps '+fig
+spawn, 'gv '+fig
 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