--- MITgcm_contrib/articles/ceaice/comments_todo.txt 2008/07/28 12:34:43 1.2 +++ MITgcm_contrib/articles/ceaice/comments_todo.txt 2008/07/29 08:33:23 1.3 @@ -1,3 +1,6 @@ +$Header: /home/ubuntu/mnt/e9_copy/MITgcm_contrib/articles/ceaice/comments_todo.txt,v 1.3 2008/07/29 08:33:23 mlosch Exp $ +$Name: $ + 27-Jul-2008, PH: --------------- @@ -50,4 +53,23 @@ ML: it's not a zonal section but a "meridional" (j-) average of grid points i=100:130, j=116:122, but maybe a section is a better idea? I have now actually removed some of the "inlets" that I have -carelessly averaged over. +carelessly averaged over: + +i = 100:130; +j = 116:121; + +rhoice = 0.92; +rhosnw = 0.33; +dxg=rdmds(fullfile('grid','DXG')); +dyg=rdmds(fullfile('grid','DYG')); +hf =rdmds(fullfile('grid','hFacC')); +ar = dxg.*dyg.*hf(:,:,1); + +[nx,ny]=size(xc); + +arij = ar(i,j); +% special masking to remove various inlets, etc. +arij(14:end,6) = 0; +arij(24:end,5) = 0; +arij(1:15,1) = 0; +arij(11,2:3) = 0;