source: trunk/src/paper01/fig11/swr_isccp_tropflux_new_v1.pro @ 43

;+
;
; .. _swr_isccp_tropflux_new_v1.pro:
;
; =============================
; swr_isccp_tropflux_new_v1.pro
; =============================
;
; DESCRIPTION
; ===========
;
; SEE ALSO
; ========
;
; :ref:`tropflux_profile.sh`
;
; EXAMPLES
; ========
;
; ::
;
;  IDL> @tropflux_init
;  IDL> swr_isccp_tropflux_new_v1
;
; EVOLUTIONS
; ==========
;
; $Id$
;
; - fplod 20110411T142955Z aedon.locean-ipsl.upmc.fr (Darwin)
;
;   * minimal header
;
;-
pro swr_isccp_tropflux_new_v1
@common
reinitplt, /z,/invert
key_portrait = 1

openps, FILENAME = 'idl.ps'
st=19890101 & en=20071231

file='/Volumes/Iomega_HDD/work/flux_reconstruction/OAFLX_GRID/swr_oafluxgrid_1985_2007.nc'
initncdf, file
sw_isccp=read_ncdf("swr", st, en, file=file,/nostr)

file='/Volumes/Iomega_HDD/flux_reconstruction/basic_data/sw_reconstred_1989_2009.nc'
initncdf, file
sw_trop=read_ncdf("sw", st, en+1, file=file,/nostr)

file='/Volumes/Iomega_HDD/flux_reconstruction/basic_data/olr_oafluxgrid_30n30s.nc'
initncdf, file
olr=read_ncdf("olr", st, en, file=file,/nostr)

help, sw_isccp, sw_trop , olr

;;creating climatology

ny=(2007-1989)+1
nt=365*ny

olri=total(reform(olr(*,*,0:nt-1),nxt,nyt,365,ny),4)/ny
swi=total(reform(sw_isccp(*,*,0:nt-1),nxt,nyt,365,ny),4)/ny
help, olri, swi

;; applying smoothing
nsmooth=45

swi=[[[swi]],[[swi]],[[swi]]] & olri=[[[olri]],[[olri]],[[olri]]]
si=size(sw) & n1=si(1) & n2=si(2)

for j=0,n2-1 do begin
  for i=0,n1-1 do begin
    x=reform(swi(i,j,*))
    y=reform(olri(i,j,*))
    swi(i,j,*)=smooth(x,nsmooth)
    olri(i,j,*)=smooth(y,nsmooth)
  endfor
endfor

swi=swi[*,*,365:365+364] & olri=olri[*,*,365:365+364]
help, swi, olri

caldat, time,mon,day,yea
swm=olr*0. & olrm=olr*0.
for jt=0,jpt-1 do begin
  jtt=(time(jt)-julday(1,1,yea(jt))) < 364
  swm(*,*,jt)=swi(*,*,jtt)
  olrm(*,*,jt)=olri(*,*,jtt)
  print, jt
endfor
swa=sw_isccp-swm        ; swr anomaly
olra=olr-olrm     ; olr anomaly
help, swa, olra
m=reform(sw_trop(*,*,0))
tmask=finite(m)
cor_olr=fltarr(jpi,jpj)+!values.f_nan & rms_olr=fltarr(jpi,jpj)+!values.f_nan

for jj=0,jpj-1 do begin
  for ji=0,jpi-1 do begin
    if (tmask(ji,jj)) then begin
      print, 'over',jj, '/', jpj-1

      swr_ano=reform(swa(ji,jj,*))
      olr_ano=reform(olra(ji,jj,*))

      cor_olr(ji,jj)=correlate(swr_ano,olr_ano)
      rms_olr(ji,jj)=(total((swr_ano-olr_ano)^2)/(jpt-1l))^.5

    endif
  endfor
endfor

cor=fltarr(jpi,jpj)+!values.f_nan & rms=fltarr(jpi,jpj)+!values.f_nan
for jj=1,jpj-2 do begin
  for ji=1,jpi-2 do begin
    if (tmask(ji,jj)) then begin
       cor(ji,jj)=(cor_olr(ji-1,jj)+cor_olr(ji+1,jj)+cor_olr(ji,jj-1)+cor_olr(ji,jj+1))*tmask(ji,jj)/(tmask(ji-1,jj)+tmask(ji+1,jj)+tmask(ji,jj-1)+tmask(ji,jj+1))
    endif
  endfor
endfor

save, cor, cor_olr, filename='correlation_isccp_olr.idl'

lct=64 & cs=0.9
marge1=[0,0,-5,0]
help, cor

plt, cor_olr,0.5,1, lct=lct, $
     small=[1,4,1],/rempl,/noer,/nocont, $
     charsize=cs,/nocolorb,marge=marge, $
     subt='',title='a) Correlation (OLRA, SWRA)',xminor=1,yminor=1

file="/Volumes/Iomega_HDD/TropFlux/input_cor/full_cor/TropFlux_sst_19890101_20091231_v20.nc"
initncdf, file
sst=read_ncdf("sst", 19890101, 20091231, file=file,/nostr)
sst=grossemoyenne(sst, "t",/nan)

plt, cor_olr,0.5,1, lct=lct, $
     small=[1,4,3],/rempl,/noer,/nocont, $
     charsize=cs,marge=marge, $
     subt='',title='a) Correlation (OLRA, SWRA)',xminor=1,yminor=1, $
     contour=sst, contmin=26.5, contmax=28.5, contint=1.

plt, rms_olr, lct=lct, $
     small=[1,4,4],/rempl,/noer,/nocont, $
     charsize=cs, marge=marge,/nocolorb, $
     subt='',title='b) RMS Difference (OLRA, SWRA)',xminor=1,yminor=1

closeps
fig='correlation_rmsd_olra_swra_smooth.ps'
spawn, 'mv '+psdir+'idl.ps '+psdir+fig
spawn, 'gv '+psdir+fig

return

end