Skip to content

Commit

Permalink
Add Updated stuff fro Minicourse
Browse files Browse the repository at this point in the history
Updated vignette and related files.
  • Loading branch information
rmendels committed Aug 21, 2015
1 parent a8a4baf commit 8f79d65
Show file tree
Hide file tree
Showing 4 changed files with 12 additions and 6 deletions.
Binary file added Marlin38606.mat
Binary file not shown.
8 changes: 8 additions & 0 deletions erdap8601.m
Original file line number Diff line number Diff line change
@@ -0,0 +1,8 @@
function matTime=erdap8601(erddapTime)
timeSize=size(erddapTime);
matTime=zeros(timeSize(1),1);
for i=1:timeSize(1);
matTime(i,1)=datenum8601(erddapTime(i,1:end));
end;
matTime=squeeze(matTime);

4 changes: 3 additions & 1 deletion makeMap.m
Original file line number Diff line number Diff line change
@@ -1,14 +1,16 @@
function [] = makeMap( longitude, latitude, param )
%UNTITLED3 Summary of this function goes here
% Detailed explanation goes here
latitude=double(latitude);
longitude=double(longitude);
[Plg,Plt]=meshgrid(longitude,latitude);
lat=[min(latitude) max(latitude)];
lon=[min(longitude) max(longitude)];
m_proj('mercator','lon',lon, 'lat', lat);
m_pcolor(Plg,Plt,param);
shading flat;
m_gshhs_h('patch',[.7 .7 .7]);
m_grid;
m_gshhs_h('patch',[.7 .7 .7]);
colorbar;


Expand Down
6 changes: 1 addition & 5 deletions vignette.m
Original file line number Diff line number Diff line change
Expand Up @@ -15,12 +15,7 @@
%% gridded data examples
xpos = [235 240];
ypos = [36 39];
tpos = ['1998-01-01';'2014-11-30'];

% Seawifs chla this will cause an error
SeaWiFS= xtracto_3D(xpos, ypos, tpos, 'swchlamday');

%this fixes the error
tpos{1} = '1998-01-16';
tpos{2} = 'last';

Expand All @@ -32,6 +27,7 @@
tpos{1}='2003-01-16';
MODIS = xtracto_3D(xpos, ypos, tpos,'mhchlamday');
chla=double(squeeze(MODIS.chlorophyll(1,:,:)));
figure;
makeMap(MODIS.longitude-360, MODIS.latitude, log(chla));

% VIIS chla for same bounds
Expand Down

0 comments on commit 8f79d65

Please sign in to comment.