pro fig2_timeline_diagram
@common
;-----------------------------------------------
reinitplt, /z,/invert
key_portrait = 1
coefpalit=.9

openps, FILENAME = 'idl.ps'
;------------------------------------------------------------
; partie a changer
;------------------------------------------------------------
;; site locations

sitelist=['8s67e','12s55e', '8s55e', '8s80.5e', '1.5s80.5e', '0n80.5e', '1.5n80.5e', '1.5s90e', $
           '0n90e', '1.5n90e', '4n90e','8n90e','12n90e', '15n90e', '5s95e', $
           '8s165e', '8s180w',  '8s155w', '8s125w', '8s110w', '8s95w',  '5s156e', '5s165e', '5s180w', '5s170w', $
          '5s155w', '5s140w', '5s125w', '5s110w', '5s95w', '2s156e', '2s165e', '2s180w', '2s170w', '2s155w', '2s140w', $
          '2s125w', '2s110w', '2s95w', '0n147e', '0n156e', '0n165e', '0n180w', '0n170w', '0n155w', '0n140w', '0n125w', $
          '0n110w', '0n95w', '2n147e', '2n156e', '2n165e', '2n180w', '2n170w', '2n155w', '2n140w', '2n125w', '2n110w', $
          '2n95w', '5n147e', '5n156e', '5n165e', '5n170w', '5n155w', '5n140w', '5n125w', '5n110w', '5n95w', $
          '8n156e', '8n165e', '8n180w', '8n170w', '9n140w', '8n125w', '8n110w', '8n95w', $
          '0n0e', '0n10w', '0n23w', '0n35w', '10s10w', '12n23w', '12n38w', '14s32w', '15n38w', '19s34w', '20n38w', $
          '21n23w', '4n23w', '4n38w', '6s10w', '8n38w', '8s30w']


nn=n_elements(sitelist)
date1=19890101
date2=20091231

no_obs_lh=findgen(7670.)*0.
no_obs_lw=findgen(7670.)*0.
no_obs_sw=findgen(7670.)*0.

for n=0, nn-1 do begin

;; reading data from mooring

    site=sitelist(n) & csite=site
    print, csite
    x=x_site_location(site)
    y=y_site_location(site)
    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]

    nsmooth=1.
    read_lh, csite,date1,date2,nsmooth, $
                     lh

    ind=where(finite(lh)) & lh(ind)=1.
    ind=where(finite(lh,/nan)) & lh(ind)=0.
    no_obs_lh=no_obs_lh+lh

endfor

sitelist=['5n165e','8s67e','12s55e', '8s55e', '8s80.5e', '0n80.5e', '1.5n80.5e', '1.5s90e', $
           '0n90e', '1.5n90e', '4n90e','8n90e','12n90e', '15n90e', '5s95e', $
           '8s165e', '8s125w', '8s110w', '8s95w',  '5s156e', '5s165e', '5s180w', '5s170w', $
          '5s155w', '5s140w', '5s125w', '5s110w', '5s95w', '2s156e', '2s165e', '2s180w', '2s170w', '2s155w', '2s140w', $
          '2s125w', '2s110w', '2s95w', '0n147e', '0n156e', '0n165e', '0n180w', '0n170w', '0n155w', '0n140w', '0n125w', $
          '0n110w', '0n95w', '2n147e', '2n156e', '2n165e', '2n180w', '2n170w', '2n155w', '2n140w', '2n125w', '2n110w', $
          '2n95w', '5n147e', '5n156e', '5n170w', '5n155w', '5n140w', '5n125w', '5n110w', '5n95w', $
          '8n156e', '8n165e', '8n180w', '8n170w', '9n140w', '8n125w', '8n110w', '8n95w', $
          '0n0e', '0n10w', '0n23w', '0n35w', '10s10w', '12n23w', '12n38w', '14s32w', '15n38w', '19s34w', '20n38w', $
          '21n23w', '4n23w', '4n38w', '6s10w', '8n38w', '8s30w']

nn=n_elements(sitelist)
date1=date1
date2=date2

for n=0, nn-1 do begin

;; reading data from mooring

    site=sitelist(n) & csite=site
    print, csite
    x=x_site_location(site)
    y=y_site_location(site)
    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]

    nsmooth=1.
    read_sw, csite, date1, date2, box, $
                           sw
    swr=sw & ind=where(finite(swr)) & valid=n_elements(ind)
    if (valid ge 10) then begin
       ind=where(finite(sw)) & sw(ind)=1.
       ind=where(finite(sw,/nan)) & sw(ind)=0.
       no_obs_sw=no_obs_sw+sw
    endif
endfor
help, no_obs_sw

sitelist=['0n110w', '0n140w', '0n165e', '0n170w', '0n23w', '0n80.5e', '0n95w', '10s10w', '12n23w', '15n38w', $
          '15n90e', '2n95w', '2s95w', '5n95w', '5s95w', '8n95w', '8s67e', '8s95w']

nn=n_elements(sitelist)
date1=date1
date2=date2


for n=0, nn-1 do begin

;; reading data from mooring

    site=sitelist(n) & csite=site
    print, csite
    x=x_site_location(site)
    y=y_site_location(site)
    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]

    nsmooth=1.
    read_lw, csite, date1, date2, box, $
                           lw
    lwr=lw & ind=where(finite(lwr)) & valid=n_elements(ind)
    if (valid ge 10) then begin
    
       ind=where(finite(lw)) &  lw(ind)=1.
       ind=where(finite(lw,/nan)) & lw(ind)=0.
       no_obs_lw=no_obs_lw+lw
    endif

endfor

help,no_obs_lw
pltt, no_obs_lh, "t", small=[1,3,1],/rempl,0,80, $
     title='No of valid observations for turbulent flux validation', charsize=1., $
     ytitle="No. of observations",subtitle='black(turbulent), red(SWR), blue(LWR)'
pltt, no_obs_sw, "t",/ov1d, color=250
pltt, no_obs_lw, "t",/ov1d, color=50
pltt, no_obs_lh, "t",/ov1d
;------------------------------------------------------------
closeps
fig='fig2_timeline_diagram.ps'
spawn, 'mv '+psdir+'idl.ps '+cpsdir+fig
spawn, 'gv '+cpsdir+fig
;------------------------------------------------------------
return
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