diff --git a/tsg_tools/climatology/test_m_map.m b/tsg_tools/climatology/test_m_map.m
new file mode 100644
index 0000000000000000000000000000000000000000..cc515d3c3f9e0d2e5c9673382ea67b34ff9742dc
--- /dev/null
+++ b/tsg_tools/climatology/test_m_map.m
@@ -0,0 +1,36 @@
+% 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
diff --git a/tsg_tools/climatology/test_opendap.m b/tsg_tools/climatology/test_opendap.m
new file mode 100644
index 0000000000000000000000000000000000000000..7dc0aba8aaa42587ea41c88e95eb66eff1554850
--- /dev/null
+++ b/tsg_tools/climatology/test_opendap.m
@@ -0,0 +1,22 @@
+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