/[MITgcm]/MITgcm/verification/global_ocean.90x40x15/diags_matlab/mit_oceanmasks.m
ViewVC logotype

Contents of /MITgcm/verification/global_ocean.90x40x15/diags_matlab/mit_oceanmasks.m

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

Revision 1.1 - (show annotations) (download)
Tue Oct 22 13:30:40 2002 UTC (21 years, 8 months ago) by mlosch
Branch: MAIN
CVS Tags: checkpoint57t_post, checkpoint58l_post, checkpoint48f_post, checkpoint46k_post, checkpoint51k_post, checkpoint53f_post, checkpoint47j_post, checkpoint54a_pre, checkpoint55c_post, checkpoint53b_pre, checkpoint48d_pre, checkpoint51l_post, checkpoint51j_post, branch-exfmods-tag, checkpoint47e_post, checkpoint57m_post, checkpoint47i_post, checkpoint52l_pre, checkpoint48i_post, checkpoint52e_pre, checkpoint58e_post, checkpoint57v_post, checkpoint57g_pre, checkpoint52j_post, checkpoint47f_post, checkpoint48d_post, checkpoint51o_pre, checkpoint57f_post, checkpoint50e_post, checkpoint52e_post, checkpoint50c_post, checkpoint57s_post, checkpoint51n_pre, checkpoint47d_post, checkpoint57j_post, checkpoint58b_post, checkpoint58m_post, checkpoint47a_post, checkpoint57b_post, checkpoint46l_pre, checkpoint53c_post, checkpoint53d_post, checkpoint57f_pre, checkpoint48a_post, checkpoint55d_pre, checkpoint46n_post, checkpoint51f_pre, checkpoint57g_post, checkpoint48e_post, checkpoint57a_post, checkpoint48h_post, checkpoint55j_post, checkpoint56b_post, checkpoint50c_pre, checkpoint57h_pre, checkpoint57y_post, branchpoint-genmake2, checkpoint46l_post, checkpoint58g_post, checkpoint57x_post, checkpoint52j_pre, checkpoint54a_post, branch-netcdf, checkpoint50d_pre, checkpoint55h_post, checkpoint51r_post, checkpoint52b_pre, checkpoint52n_post, checkpoint54b_post, checkpoint58h_post, checkpoint51i_post, checkpoint57e_post, checkpoint54d_post, checkpoint47h_post, checkpoint48c_post, checkpoint56c_post, checkpoint54e_post, checkpoint58j_post, checkpoint55b_post, checkpoint51e_post, checkpoint51b_post, checkpoint57h_post, checkpoint51l_pre, checkpoint52m_post, checkpoint51c_post, checkpoint57y_pre, checkpoint55, checkpoint53a_post, checkpoint55a_post, checkpoint57c_pre, checkpoint48, checkpoint49, checkpoint47b_post, checkpoint53b_post, checkpoint57o_post, checkpoint55g_post, checkpoint57r_post, checkpoint51o_post, checkpoint48g_post, checkpoint57k_post, checkpoint57d_post, checkpoint55f_post, checkpoint57i_post, checkpoint51q_post, checkpoint52l_post, checkpoint52k_post, checkpoint58, checkpoint57a_pre, checkpoint54, checkpoint57, checkpoint56, checkpoint51, checkpoint50, checkpoint53, checkpoint52, checkpoint50d_post, checkpoint52d_post, checkpoint46m_post, checkpoint51b_pre, checkpoint52a_post, checkpoint57h_done, checkpoint47g_post, checkpoint58f_post, checkpoint52b_post, checkpoint53g_post, checkpoint52f_post, checkpoint57n_post, checkpoint52c_post, checkpoint58d_post, checkpoint57w_post, checkpoint57p_post, checkpint57u_post, checkpoint51h_pre, checkpoint58a_post, checkpoint50g_post, checkpoint58i_post, checkpoint57q_post, checkpoint51g_post, ecco_c52_e35, checkpoint57z_post, checkpoint54f_post, checkpoint51f_post, checkpoint48b_post, checkpoint50b_post, eckpoint57e_pre, checkpoint58c_post, checkpoint58k_post, checkpoint57c_post, checkpoint50f_post, checkpoint50a_post, checkpoint50f_pre, checkpoint52a_pre, checkpoint47d_pre, checkpoint51d_post, checkpoint48c_pre, checkpoint51m_post, checkpoint51t_post, checkpoint53d_pre, checkpoint47, checkpoint55e_post, checkpoint54c_post, checkpoint50h_post, checkpoint52i_post, checkpoint51a_post, checkpoint50e_pre, checkpoint50i_post, checkpoint51p_post, checkpoint51n_post, checkpoint55i_post, checkpoint51i_pre, checkpoint57l_post, checkpoint52i_pre, checkpoint51u_post, checkpoint52h_pre, checkpoint52f_pre, checkpoint50b_pre, checkpoint56a_post, checkpoint51s_post, checkpoint55d_post
Branch point for: netcdf-sm0, branch-genmake2, branch-nonh, tg2-branch, checkpoint51n_branch, release1, branch-exfmods-curt 
 fixed the verification/global_ocean.90x40x15 experiment:
 - new bathymetry (the world according to A., JMC, and M.)
 - new initial fields and forcing fields (*.bin files)
 - new POLY3.COEFFS (for the next release one should switch to a full
   equation of state: JMD95P or MDJWF)
 - fixed several errors and redundancies in the data file
 - experiment uses looped cells
 - added matlab directory with diagnostic scripts for plotting of output
1 function grid = mit_oceanmasks(grid,plotit)
2 % function grid = mit_oceanmasks(grid)
3 % create masks for the 4x4 degree run for the individual oceans
4 % hfacc/s/w have to be available for this
5
6 if grid.nx ~= 90 | grid.ny ~= 40
7 error(sprintf( ...
8 '%s\n%s\n%s\n%s', ...
9 'Your horizontal grid has dimensions different from 90x40.', ...
10 'mit_oceanmasks creates masks for the individual ocean basins,', ...
11 'therefore it needs to be customized if you change the', ...
12 'horizontal dimensions.'))
13 end
14
15 % masks for the zonal velocities:
16 % hfacw
17
18 % Atlantic mask
19 grid.atlantic_hfacw = grid.hfacw;
20 ixsow = 6:73; % southern ocean and most of pacific
21 iysow = 1:24;
22 grid.atlantic_hfacw(ixsow,iysow,:) = NaN;
23 ixnpw = 15:67; % northern pacific
24 iynpw = 25:37;
25 grid.atlantic_hfacw(ixnpw,iynpw,:) = NaN;
26 % Pacific mask
27 grid.pacific_hfacw = grid.hfacw;
28 ixaw = 74:90; % most of the atlantic
29 iyaw = 1:40;
30 grid.pacific_hfacw(ixaw,iyaw,:) = NaN;
31 ixnaw = 68:73; % north atlantic and gulf of Mexico
32 iynaw = 27:40;
33 grid.pacific_hfacw(ixnaw,iynaw,:) = NaN;
34 ixcaw = 70:73; % carribean
35 iycaw = 24:26;
36 grid.pacific_hfacw(ixcaw,iycaw,:) = NaN;
37 ixiw = 1:27; % most of the indic
38 iyiw = 1:40;
39 grid.pacific_hfacw(ixiw,iyiw,:) = NaN;
40 ixauw = 28:35; % water south and west of Australia
41 iyauw = 1:20;
42 grid.pacific_hfacw(ixauw,iyauw,:) = NaN;
43
44
45 % masks for the meridional overturning stream functions:
46 % hfacs
47
48 % Atlantic mask
49 grid.atlantic_hfacs = grid.hfacs;
50 ixso = 5:73; % southern ocean
51 iyso = 1:24;
52 grid.atlantic_hfacs(ixso,iyso,:) = NaN;
53 ixnp = 14:65; % northern pacific
54 iynp = 25:36;
55 grid.atlantic_hfacs(ixnp,iynp,:) = NaN;
56 ixsoc = 1:90; % complete southern ocean
57 iysoc = 1:12; % south of latg(12) = -34 degN;
58 grid.atlantic_hfacs(ixsoc,iysoc,:) = NaN;
59
60 % Pacific mask
61 grid.pacific_hfacs = grid.hfacs;
62 ixatl = [1:26 74:90]; % most of altantic and indian ocean
63 iyatl = 1:grid.ny;
64 grid.pacific_hfacs(ixatl,iyatl,:) = NaN;
65 ixatlf = 73; % small fix
66 iyatlf = 31;
67 grid.pacific_hfacs(ixatlf,iyatlf,:) = NaN;
68 ixgm = 66:73; % Gulf of Mexico
69 iygm = 25:30;
70 grid.pacific_hfacs(ixgm,iygm,:) = NaN;
71 ixhb = 68:73; % Hudson Bay and Baffin Bay
72 iyhb = 36:40;
73 grid.pacific_hfacs(ixhb,iyhb,:) = NaN;
74 ixsop = 1:grid.nx; % southern ocean
75 iysop = 1:12;
76 grid.pacific_hfacs(ixsop,iysop,:) = NaN;
77 ixio = 27:32; % rest of indian ocean
78 iyio = 13:19;
79 grid.pacific_hfacs(ixio,iyio,:) = NaN;
80 %pacific_hfacs(33:34,12,:) = NaN; % two points south of Australia
81
82 % Indian ocean
83 grid.indic_hfacs = grid.hfacs;
84 ixsoi = 1:90; % southern ocean
85 iysoi = 1:12;
86 grid.indic_hfacs(ixsoi,iysoi,:) = NaN;
87 ixap = [1:4 33:90]; % atlantic and pacific
88 iyap = 1:40;
89 grid.indic_hfacs(ixap,iyap,:) = NaN;
90 iynps = 35:40; % north polar sea
91 grid.indic_hfacs(:,iynps,:) = NaN;
92 ixmeds = 1:10; % mediterraenan
93 iymeds = 29:31;
94 grid.indic_hfacs(ixmeds,iymeds,:) = NaN;
95 ixscs = 27:32; % south china sea
96 iyscs = 21:29;
97 grid.indic_hfacs(ixscs,iyscs,:) = NaN;
98 grid.indic_hfacs(5,39,:) = NaN; % singular point near Murmansk
99
100 % Southern Ocean
101 grid.so_hfacs = change( ...
102 change(grid.atlantic_hfacs,'==',NaN,-1) ...
103 + change(grid.pacific_hfacs,'==',NaN,-1) ...
104 + change(grid.indic_hfacs,'==',NaN,-1) ...
105 + grid.hfacs-grid.hfacs,'==',-1,NaN);
106 iso = find(~isnan(grid.so_hfacs));
107 grid.so_hfacs(iso)= grid.hfacs(iso);
108
109 if exist('plotit','var')
110 figure;
111 spy(isnan(grid.hfacs(:,:,1)'),'kx');
112 hold on;
113 spy(~isnan(grid.atlantic_hfacs(:,:,1)'));
114 spy(~isnan(grid.indic_hfacs(:,:,1)'),'g.');
115 spy(~isnan(grid.so_hfacs(:,:,1)'),'r.');
116 spy(~isnan(grid.pacific_hfacs(:,:,1)'),'y.');
117 title('hfacs')
118 axis xy
119 end
120
121 % masks for zonal averages of C-properties:
122 % hfacc
123
124 % Atlantic mask
125 grid.atlantic_hfacc = grid.hfacc;
126 ixe = 6:65; % indian ocean
127 iye = 1:28;
128 grid.atlantic_hfacc(ixe,iye,:) = NaN;
129 ixep = 26:65; % eastern pacific
130 iyep = 1:grid.ny;
131 grid.atlantic_hfacc(ixep,iyep,:) = NaN;
132 ixw = 6:73; % western pacific
133 iyw = 1:23;
134 grid.atlantic_hfacc(ixw,iyw,:) = NaN;
135 ixgl = 66:68; % galapagos
136 iygl = 24;
137 grid.atlantic_hfacc(ixgl,iygl,:) = NaN;
138 ixch = 73; % coast of Chile
139 iych = 13:15;
140 grid.atlantic_hfacc(ixch,iych,:) = NaN;
141
142 % Pacific mask
143 grid.pacific_hfacc = grid.hfacc;
144 ixatl = [1:25 74:90]; % most of altantic and indian ocean
145 iyatl = 1:grid.ny;
146 grid.pacific_hfacc(ixatl,iyatl,:) = NaN;
147 ixatlf = 73; % small fix
148 iyatlf = 31;
149 grid.pacific_hfacc(ixatlf,iyatlf,:) = NaN;
150 ixgme = 69:73; % Gulf of Mexico
151 iygme = 24:30;
152 ixgmw = 66:69;
153 iygmw = 26:28;
154 grid.pacific_hfacc(ixgme,iygme,:) = NaN;
155 grid.pacific_hfacc(ixgmw,iygmw,:) = NaN;
156 ixhbb = 68:73; % Hudson and Baffin Bay
157 iyhbb = 35:40;
158 grid.pacific_hfacc(ixhbb,iyhbb,:) = NaN;
159 ixio = 26:34; % rest of indian ocean
160 iyio = 1:19;
161 grid.pacific_hfacc(ixio,iyio,:) = NaN;
162 % grid.pacific_hfacc(33:34,12,:) = NaN; % two points south of Australia
163
164 % Indian ocean
165 grid.indic_hfacc = ...
166 change(change(grid.atlantic_hfacc,'==',NaN,0),'>',0,NaN) ...
167 + change(change(grid.pacific_hfacc,'==',NaN,0),'>',0,NaN) ...
168 + grid.hfacc;
169
170 if exist('plotit','var')
171 figure;
172 spy(isnan(grid.hfacc(:,:,1)'),'kx'); axis xy
173 hold on;
174 spy(~isnan(grid.atlantic_hfacc(:,:,1)'));
175 spy(~isnan(grid.pacific_hfacc(:,:,1)'),'r.');
176 spy(~isnan(grid.indic_hfacc(:,:,1)'),'g.');
177 axis xy
178 title('hfacc');
179
180 end
181
182 return
  ViewVC Help
Powered by ViewVC 1.1.22