
namerun='erai51c' 
yd='1993' 
yf='1994' 
rep=[yd '0101_' yf '1231'] 
ydur=[yd '-' yf]

%namerun='erai51c';
%yd='1993'
%yf='1994'
%rep='19930101_19941231' 

champ='mldr10_1';

%root_dir='/L1/cholod/NEMO/PERU12/'
 root_dir=['/net/adonis/usr/adonis/phybiocar/cholod/peru12_nemo/']

 fs=[root_dir '/' namerun '/' rep '/peru12_5d_' rep '_grid_T_2D.' namerun '.nc']
%fs='/usr/temp/fclod/peru12_1d_1994_grid_T.erai3b.nc';
 
 fg=[root_dir 'mesh_mask_PERU12.nc']
 
nc=netcdf(fs,'r');
ng=netcdf(fg,'r');
lat=ng{'nav_lat'}(:);
lon=ng{'nav_lon'}(:);
 
%        float tos(time_counter, lat, lon) ;
champ2=nc{champ}(:,:,:);
mchamp=squeeze(mean(champ2,1));


opengl neverselect % pour pb cartes graphiques

figure;
pcolor(lon,lat,mchamp);
shading flat;
set(gca,'PlotBoxAspectRatio',[0.6 1 1]);
set(gca,'FontS',12);
caxis([13 31]);
title(['MLDR01-1 (m) - ' namerun ' (' ydur ')'], 'FontSize',16);
hcc=colorbar; set(hcc,'FontS',12);

print ('-dpng','figure1.png')

quit
