1 |
|
$Header$ |
2 |
|
$Name$ |
3 |
|
|
4 |
27-Jul-2008, PH: |
27-Jul-2008, PH: |
5 |
--------------- |
--------------- |
6 |
|
|
53 |
ML: it's not a zonal section but a "meridional" (j-) average of grid |
ML: it's not a zonal section but a "meridional" (j-) average of grid |
54 |
points i=100:130, j=116:122, but maybe a section is a better idea? |
points i=100:130, j=116:122, but maybe a section is a better idea? |
55 |
I have now actually removed some of the "inlets" that I have |
I have now actually removed some of the "inlets" that I have |
56 |
carelessly averaged over. |
carelessly averaged over: |
57 |
|
|
58 |
|
i = 100:130; |
59 |
|
j = 116:121; |
60 |
|
|
61 |
|
rhoice = 0.92; |
62 |
|
rhosnw = 0.33; |
63 |
|
dxg=rdmds(fullfile('grid','DXG')); |
64 |
|
dyg=rdmds(fullfile('grid','DYG')); |
65 |
|
hf =rdmds(fullfile('grid','hFacC')); |
66 |
|
ar = dxg.*dyg.*hf(:,:,1); |
67 |
|
|
68 |
|
[nx,ny]=size(xc); |
69 |
|
|
70 |
|
arij = ar(i,j); |
71 |
|
% special masking to remove various inlets, etc. |
72 |
|
arij(14:end,6) = 0; |
73 |
|
arij(24:end,5) = 0; |
74 |
|
arij(1:15,1) = 0; |
75 |
|
arij(11,2:3) = 0; |