Skip to content
Snippets Groups Projects
Commit 062bd2c4 authored by jacques.grelet_ird.fr's avatar jacques.grelet_ird.fr
Browse files

add test scripts for climato

parent 08813cfb
No related branches found
No related tags found
No related merge requests found
% test plot climato with m_map
%ncdisp('woa13_annual_surf.nc'
SSTP = ncread('woa13_annual_surf.nc', 'WOA_MEAN_SSTP');
SSPS = ncread('woa13_annual_surf.nc', 'WOA_MEAN_SSPS');
LATX = ncread('woa13_annual_surf.nc', 'WOA_LATX');
LONX = ncread('woa13_annual_surf.nc', 'WOA_LONX');
%latMin = -50; latMax = 65; lonMin = 40; lonMax = 280;
%latMin = -50; latMax = 65; lonMin = 180; lonMax = 360;
latMin = -50; latMax = 65; lonMin = 179; lonMax = 181;
m_proj('Mercator','lat',[latMin latMax],'long',[lonMin lonMax]);
climato = SSTP;
% permute matrix
climato = permute(climato, fliplr(1:4));
climato = squeeze(climato(1, 1,:,:));
indLat = find(LATX > latMin & LATX < latMax);
%LONX = wrapTo360_p(LONX);
%LONX = wrapTo180_p(LONX);
indn = find(LONX < 0);
indp = find(LONX >= 0);
LONX = LONX(indn)
%LONX = mod(LONX+180,360)-180;
indLon = find(LONX > lonMin & LONX < lonMax);
lon = LONX(indLon);
lat = LATX(indLat);
m_pcolor(lon, lat, climato(indLat, indLon));
shading flat;
m_coast('patch',[.7 .7 .7]);
m_grid('box', 'fancy', 'tickdir', 'in', 'Fontsize', 12);
disp('end')
\ No newline at end of file
url = 'https://data.nodc.noaa.gov/thredds/dodsC/woa/WOA13/DATAv2/salinity/netcdf/A5B2/0.25/woa13_A5B2_s00_04v2.nc'
meta = ncinfo(url);
disp({meta.Variables.Name}');
dtime = ncread(url,'time');
depth = ncread(url,'depth');
lat = ncread(url,'lat');
lon = ncread(url,'lon');
%s_an = ncread(url,'s_an')
s_an = ncread(url,'s_an',[1,1,2,1], [1440,720,2,1]);
sst = s_an(:,:,1,1)';
pcolor(lon,lat,sst)
shading interp
colorbar
figure
sst = s_an(:,:,2,1)';
pcolor(lon,lat,sst)
shading interp
colorbar
\ No newline at end of file
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment