-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathcurrentTS.m
65 lines (51 loc) · 1.79 KB
/
currentTS.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
%reads in all S and T, plots last 2 snapshots
close all
addpath('/nobackup1/gbrett/mStuff')
cd('/nobackup1/gbrett/spinupsteady')
load('geometrySpinupSteady.mat');
S=rdmds('S',NaN);
[nx,ny,nz,nt]=size(S);
figure; contourf(XC,YC,S(:,:,1,end),30:0.5:40); shading 'flat'; colorbar; %caxis([34 40])
hold all; plot(lonCoast,latCoast,'k')
fn1=strcat('SHour',num2str(nt));
save2pdf(fn1)
close all
figure; contourf(XC,YC,S(:,:,1,end-1),30:0.5:40); shading 'flat'; colorbar; %caxis([34 40])
hold all; plot(lonCoast,latCoast,'k')
fn1=strcat('SHour',num2str(nt-1));
save2pdf(fn1)
close all
figure; contourf(XC,YC,S(:,:,21,end),30:0.5:40); shading 'flat'; colorbar; %caxis([34 40])
hold all; plot(lonCoast,latCoast,'k')
fn7=strcat('S300mHour',num2str(nt));
save2pdf(fn7)
close all
figure; contourf(XC,YC,S(:,:,21,end-1),30:0.5:40); shading 'flat'; colorbar; %caxis([34 40])
hold all; plot(lonCoast,latCoast,'k')
fn7=strcat('S300mHour',num2str(nt-1));
save2pdf(fn7)
close all
clear S
load('geometrySpinupSteady.mat');
T=rdmds('T',NaN);
figure; contourf(XC,YC,T(:,:,1,end),30); shading 'flat'; colorbar; %caxis([34 40])
hold all; plot(lonCoast,latCoast,'k')
fn1=strcat('THour',num2str(nt));
save2pdf(fn1)
close all
figure; contourf(XC,YC,T(:,:,1,end-1),30); shading 'flat'; colorbar; %caxis([34 40])
hold all; plot(lonCoast,latCoast,'k')
fn1=strcat('THour',num2str(nt-1));
save2pdf(fn1)
close all
figure; contourf(XC,YC,T(:,:,21,end),30); shading 'flat'; colorbar; %caxis([34 40])
hold all; plot(lonCoast,latCoast,'k')
fn7=strcat('T300mHour',num2str(nt));
save2pdf(fn7)
close all
figure; contourf(XC,YC,T(:,:,21,end-1),30); shading 'flat'; colorbar; %caxis([34 40])
hold all; plot(lonCoast,latCoast,'k')
fn7=strcat('T300mHour',num2str(nt-1));
save2pdf(fn7)
close all
clear T