clear all;
close all;

ww3dir='/home1/datawork/xcouvela/WW3_GLOBAL_ST4/ECMWF_ERAIM/';
chafile='WW3_JFM.nc';
%chafile='WW3_JAS.nc';

lat=double(ncread([ww3dir,chafile],'latitude'));
lon=double(ncread([ww3dir,chafile],'longitude'));

%used for m_proj
minLat=min(lat(:));
maxLat=max(lat(:));
minLon=min(lon(:));
maxLon=max(lon(:));

CHA=squeeze(ncread([ww3dir,chafile],'cha'));

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


m_proj('Gall-Peters','lon',[minLon maxLon],'lat',[minLat maxLat])

%
figure
m_pcolor(lon,lat,CHA'-0.018);shading interp;%colorbar('h')
m_coast('patch',0.3*[.6 .6 .6],'edgecolor','k');
m_grid('box','fancy','tickdir','in','FontSize',15);
colormap(b2r(-.007,.007))
title('JFM avegraged Charnock-0.018')
set(gca,'fontweight','bold','fontsize',13)
%print -dpng cha.png
%!cp cha.png $HOME/Dropbox
