PRO plot_ssh_xy_obs_jrpaper8
@inithycom_woa025deg

    glamt[where(glamt LT 0.)]=glamt[where(glamt LT 0.)]+360.
    glamf[where(glamf LT 0.)]=glamf[where(glamf LT 0.)]+360.
;    glamu[where(glamu LT 0.)]=glamu[where(glamu LT 0.)]+360.
;    glamv[where(glamv LT 0.)]=glamv[where(glamv LT 0.)]+360.
	;; Read mean and variance
	   key_portrait = 1
	   PS = 1
	   io_nc = '/Net/gleam/abozec/HYCOM/JRA55-res_PAPER/'
	   ps_dir=io_nc+'PS_PAPER/'
	   model='Observations AVISO'
		;; low resolution
	   file='OBS/SDA_EKE_CMEMS_1993_2018_WOAgrid.nc'
	   eke=ncdf_lec(io_nc+file,var='sda')
		eke[where(eke LE 0)]=0.
		eke[1439,719,*]=0.
		 eke[where(eke GE 1e10)]=!values.f_nan
 

	bbox1=[0,360.,-78,89,0,10]
		  eke_surf=eke[*,*]^2
		  ;eke_surf[where(eke_surf LE 50.)]=0.
	   if (PS EQ 1 ) then openps,ps_dir+'sshvar_Obs_aviso.ps'
	   plt,eke_surf,box=bbox1,lct=64,small=[1,2,1],$
		 title=' Mean surface SSH variability AVISO Observations  (m2)',min=0, max=0.025,$
         /nocontour,window=1,format='(f7.3)'
	   if (PS EQ 1 ) then closeps

stop
end