;+
; .. _fig1_no_obs_map.pro:
;
; ===================
; fig1_no_obs_map.pro
; ===================
;
; DESCRIPTION
; ===========
;
; :file:`${PROJECT_OD}/mask_oaflux_30N30S.nc`
; containing OAFLUX grid
; has been produced by
; :func:`oaflux_mask_30N30S`.
;
; Figure of
; ++
; is saved in
; :file:`${PROJECT_OD}/fig1_no_obs_map.ps`.
;
; .. graphviz::
;
;    digraph fig1_no_obs_map {
;
;       mask [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/mask_oaflux_30N30S.nc"];
;       lhf_no_obs [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/lhf_no_obs.txt"];
;       sw_no_obs [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/sw_no_obs.txt"];
;       lw_no_obs [shape=ellipse,fontname=Courier,label="${PROJECT_ID}/lw_no_obs.txt"];
;       figure [shape=ellipse,fontname=Courier,label="${PROJECT_OD}/fig1_no_obs_map.ps"];
;
;       fig1_no_obs_map [shape=box,
;       fontname=Courier,
;       color=blue,
;       URL="http://forge.ipsl.jussieu.fr/tropflux/browser/trunk/src/paper01/fig1/fig1_no_obs_map.pro",
;       label="${TROPFLUX}/src/paper01/fig1/fig1_no_obs_map.pro"];
;
;       {mask lhf_no_obs sw_no_obs lw_no_obs} -> {fig1_no_obs_map} -> {figure}
;
;    }
;
; SEE ALSO
; ========
;
; :ref:`data_in_oaflux`
;
; :ref:`project_profile.sh`
; :ref:`project_init.pro`
;
; :func:`x_site_location`
; :func:`y_site_location`
;
; :func:`oaflux_mask_30N30S`
;
; EXAMPLES
; ========
;
; ::
;
;  fig1_no_obs_map
;
; TODO
; ====
;
; how text files have been produced ?
;
; coding rules
;
; complete description
;
; handle IO error
;
; EVOLUTIONS
; ==========
;
; $Id$
;
; - fplod 20110817T104358Z aedon.locean-ipsl.upmc.fr (Darwin)
;
;   * get rid of spawn
;
; - pinsard 2011-04-20T15:08:28Z loholt1.ipsl.polytechnique.fr (Linux)
;
;   * make it work on loholt1
;
; - fplod 20110420T134821Z aedon.locean-ipsl.upmc.fr (Darwin)
;
;   * externalize functions
;
; - fplod 20110412T073818Z aedon.locean-ipsl.upmc.fr (Darwin)
;
;   * make it work on cratos
;   * remove hard coding path
;   * add graphviz
;
; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin)
;
;   * minimal header
;
;-
pro fig1_no_obs_map
;
@cm_4cal
@cm_4data
@cm_4mesh
@cm_4data
@cm_project
;
reinitplt, /z,/invert
key_portrait = 1
;
openps, FILENAME = project_od_env + 'fig1_no_obs_map.ps'
close,/all
;
yr=365.
lh_mi=365/yr
lh_ma=3285/yr
lh_int=0.25
sw_mi=365/yr
sw_ma=3285/yr
sw_int=0.25
lw_mi=365/yr
lw_ma=1460/yr
lw_int=0.125
lct=65
;
file=project_id_env + 'mask_oaflux_30N30S.nc'
initncdf, file
domdef, 30,390,-30,30
msk=ncdf_lec(file,var='msk')
;
fi=project_id_env + 'lhf_no_obs.txt'
res=read_ascii(fi,data_start=1)
ff=res.field1
lat=reform(ff(0,*))
lon=reform(ff(1,*))
no=reform(ff(2,*))
no=no/yr
;ind=where(no ge lh_ma)
; no(ind)=no-0.5
;
marge1=[0,0,-5,0]
plt, .4+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,box=box,lct=lct, lh_mi, lh_ma, int=lh_int, $
   title='a) Turbulent Flux validation sites', subtitle='', small=[1,4,1],/rempl, marge=marge1, format='(i5)'
NN=n_elements(lat)
;
usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill
;
for n=0,NN-1 do begin
  x=lon[n]
  y=lat[n]
  c=no[n]
  cmi=lh_mi
  cma=lh_ma
  dc=cma-cmi
  col=((10+244*(c-cmi)/dc) > 10) < 254
  plots, x,y,psym=8,symsize=1.5,color=col
endfor
;
usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5]
;
for n=0,NN-1 do begin
  x=lon[n]
  y=lat[n]
  c=no[n]
;
  plots, x,y,psym=8,symsize=1.5,color=0
endfor
;
fi=project_id_env + 'sw_no_obs.txt'
res=read_ascii(fi,data_start=1)
ff=res.field1
lat=reform(ff(0,*))
lon=reform(ff(1,*))
no=reform(ff(2,*))
no=no/yr
;ind=where(no ge sw_ma)
;no(ind)=sw_ma-0.5
marge1=[0,0,-5,0]
plt, .4+msk*1.e-5,realcont=2,/nocont,/nofill,xminor=1,yminor=1,box=box,lct=lct, sw_mi, sw_ma, int=sw_int, $
   title='b) SWR validation sites', subtitle='', small=[1,4,2],/noer,/rempl, marge=marge1, format='(i5)'
NN=n_elements(lat)
;
usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill
;
for n=0,NN-1 do begin
  x=lon[n]
  y=lat[n]
  c=no[n]
  cmi=sw_mi
  cma=sw_ma
  dc=cma-cmi
  col=((10+244*(c-cmi)/dc) > 10) < 254
  plots, x,y,psym=8,symsize=1.5,color=col
endfor
;
usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5]
;
for n=0,NN-1 do begin
  x=lon[n]
  y=lat[n]
  c=no[n]
;
  plots, x,y,psym=8,symsize=1.5,color=0
endfor
;
fi=project_id_env + 'lw_no_obs.txt'
res=read_ascii(fi,data_start=1)
ff=res.field1
lat=reform(ff(0,*))
lon=reform(ff(1,*))
no=reform(ff(2,*))
no=no/yr
;ind=where(no ge lw_ma)
;no(ind)=lw_ma-0.25
;
marge1=[0,0,-5,0]
plt, .4+msk*1.e-5,realcont=2,/nocont,/nofill,/noer,xminor=1,yminor=1,box=box,lct=lct,lw_mi, lw_ma, int=lw_int, $
   title='c) LWR validation sites', subtitle='', small=[1,4,3],/rempl, marge=marge1, format='(i5)'
NN=n_elements(lat)
;
usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5],/fill
;
for n=0,NN-1 do begin
  x=lon[n]
  y=lat[n]
  c=no[n]
  cmi=lw_mi
  cma=lw_ma
  dc=cma-cmi
  col=((10+244*(c-cmi)/dc) > 10) < 254
  plots, x,y,psym=8,symsize=1.5,color=col
endfor
;
usersym, [-.5,.5,.5,-.5,-.5],[-.5,-.5,.5,.5,-.5]
;
for n=0,NN-1 do begin
  x=lon[n]
  y=lat[n]
  c=no[n]
;
  plots, x,y,psym=8,symsize=1.5,color=0
endfor
;
; independent validation sites
;
plots, 61.5, 15.5,psym=2      ; as location
plots, 331, 25.5,psym=2       ; subduction (central)
plots, 326, 18,psym=2         ; subduction (sw)
plots, 338, 18,psym=2         ; subduction (se)
plots, 156, -1.7 , psym=2     ; coare
plots, 274.2, -21, psym=2     ; stratus
plots, 202.46, 22.45, psym=2  ; WHOTS
;
closeps
;
end