/[MITgcm]/MITgcm_contrib/ocean_inversion_project/write_netCDF/lookat.m
ViewVC logotype

Contents of /MITgcm_contrib/ocean_inversion_project/write_netCDF/lookat.m

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph


Revision 1.1 - (show annotations) (download)
Wed Dec 17 18:14:27 2003 UTC (21 years, 7 months ago) by dimitri
Branch: MAIN
CVS Tags: HEAD
Added Files to ocean_inversion_project

1 % lookat at 4x4 netcdf output
2
3 f=netcdf('ECCO_MaskAreaBathy.nc');
4 for i=1:11
5 f{i}
6 minmax(f{i}(:))
7 pause
8 end
9
10 f=netcdf('ECCO_Stationary_1_phys.nc')
11 for i=1:10
12 f{i}
13 minmax(f{i}(:))
14 pause
15 end
16
17 f=netcdf('ECCO_Stationary_3001_phys.nc')
18 for i=1:10
19 f{i}
20 minmax(f{i}(:))
21 pause
22 end
23
24 f=netcdf('ECCO_Timedep_phys.nc')
25 for i=1:10
26 f{i}
27 minmax(f{i}(:))
28 pause
29 end
30
31 f=netcdf('ECCO_Stationary_BasisFNCTNS_3001.nc')
32 tmp=f{1}(:);
33 for i=1:30
34 clf
35 subplot(221),mypcolor(squeeze(tmp(i,1,:,:)));colorbar,title(i)
36 subplot(222),mypcolor(squeeze(tmp(i,3,:,:)));colorbar,title(i)
37 subplot(223),mypcolor(squeeze(tmp(i,6,:,:)));colorbar,title(i)
38 subplot(224),mypcolor(squeeze(tmp(i,9,:,:)));colorbar,title(i)
39 pause(1)
40 end
41
42 f=netcdf('ECCO_Timedep_BasisFNCTNS_1775.nc')
43 tmp=f{1}(:);
44 for i=1:30
45 clf
46 subplot(221),mypcolor(squeeze(tmp(i,1,:,:)));colorbar,title(i)
47 subplot(222),mypcolor(squeeze(tmp(i,3,:,:)));colorbar,title(i)
48 subplot(223),mypcolor(squeeze(tmp(i,6,:,:)));colorbar,title(i)
49 subplot(224),mypcolor(squeeze(tmp(i,9,:,:)));colorbar,title(i)
50 pause(1)
51 end
52
53 f=netcdf('ECCO_Timedep_BasisFNCTNS_2005.nc')
54 tmp=f{1}(:);
55 for i=1:30
56 clf
57 subplot(221),mypcolor(squeeze(tmp(i,1,:,:)));colorbar,title(i)
58 subplot(222),mypcolor(squeeze(tmp(i,3,:,:)));colorbar,title(i)
59 subplot(223),mypcolor(squeeze(tmp(i,6,:,:)));colorbar,title(i)
60 subplot(224),mypcolor(squeeze(tmp(i,9,:,:)));colorbar,title(i)
61 pause(1)
62 end
63
64 f=netcdf('ECCO_Stationary_diag_2D.nc')
65 tmp1=f{2}(:);
66 tmp2=f{3}(:);
67 for i=1:30
68 clf
69 subplot(211),mypcolor(squeeze(tmp1(i,:,:)));colorbar,title(i)
70 subplot(212),mypcolor(squeeze(tmp2(i,:,:)));colorbar,title(i)
71 pause(.1)
72 end
73
74 f=netcdf('ECCO_Timedep_diag_0D.nc')
75 tmp1=f{1}(:);
76 tmp2=f{2}(:);
77 tmp3=f{3}(:);
78 tmp4=f{4}(:);
79 for i=1:30
80 clf
81 subplot(311),plot(tmp1,tmp2(:,i))
82 subplot(312),plot(tmp1,tmp3(:,i))
83 subplot(313),plot(tmp1,tmp4(:,i))
84 title(i), pause(1)
85 end
86
87 f=netcdf('ECCO_Timedep_diag_2D.nc')
88 tmp1=f{2}(:);
89 tmp2=f{3}(:);
90 for i=1:30
91 clf
92 subplot(221),mypcolor(squeeze(tmp1(1,i,:,:)));colorbar,title(i)
93 subplot(222),mypcolor(squeeze(tmp1(55,i,:,:)));colorbar
94 subplot(223),mypcolor(squeeze(tmp2(1,i,:,:)));colorbar,title(i)
95 subplot(224),mypcolor(squeeze(tmp2(55,i,:,:)));colorbar
96 pause(.1)
97 end
98
99 f=netcdf('ECCO_Stationary_diag_0D.nc')
100 tmp1=f{1}(:);
101 tmp2=f{2}(:);
102 tmp3=f{3}(:);
103 tmp4=f{4}(:);
104 for i=1:30
105 clf
106 subplot(311),plot(tmp1,tmp2(:,i))
107 subplot(312),plot(tmp1,tmp3(:,i))
108 subplot(313),plot(tmp1,tmp4(:,i))
109 title(i), pause(1)
110 end

  ViewVC Help
Powered by ViewVC 1.1.22