/[MITgcm]/MITgcm/verification/seaice_obcs/input/compare_lab_sea.m
ViewVC logotype

Diff of /MITgcm/verification/seaice_obcs/input/compare_lab_sea.m

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

revision 1.1 by dimitri, Thu Jan 17 20:48:08 2008 UTC revision 1.2 by dimitri, Wed Dec 19 02:44:21 2012 UTC
# Line 2  Line 2 
2  % results to ../lab_sea/input.salt_plume  % results to ../lab_sea/input.salt_plume
3    
4  % preamble  % preamble
5  nx=20; ny=16; nz=18; nt=6; ix=8:17; iy=4:11;  %cd MITgcm/verification/seaice_obcs/input
6    nx=20; ny=16; nz=18; nt=10; ix=8:17; iy=4:11;
7  p1='../../lab_sea/tr_run.salt_plume/';  p1='../../lab_sea/tr_run.salt_plume/';
8  p2='../run/';  p2='../run/';
9    
# Line 20  for f={'T','S','U','V'} Line 21  for f={'T','S','U','V'}
21          title(['seaice obcs ' f{1} ', level ' int2str(k)]);          title(['seaice obcs ' f{1} ', level ' int2str(k)]);
22          tmp3=tmp2(:,:,k)-tmp1(:,:,k);          tmp3=tmp2(:,:,k)-tmp1(:,:,k);
23          subplot(313), mypcolor(tmp3'); thincolorbar          subplot(313), mypcolor(tmp3'); thincolorbar
24          title('difference'), pause(.1)          title('difference'), pause(.2)
25          J=J+sum(sum(tmp3(2:(end-1),2:(end-1)).^2));          J=J+sum(sum(tmp3(2:(end-1),2:(end-1)).^2));
26      end      end
27      disp([f{1} ': ' num2str(J)])      disp([f{1} ': ' num2str(J)])
# Line 40  for f={'UICE','VICE','AREA','HSNOW','HSA Line 41  for f={'UICE','VICE','AREA','HSNOW','HSA
41      title(['seaice obcs ' f{1} ', time step ' int2str(t)]);      title(['seaice obcs ' f{1} ', time step ' int2str(t)]);
42      tmp3=tmp2-tmp1;      tmp3=tmp2-tmp1;
43      subplot(313), mypcolor(tmp3'); thincolorbar      subplot(313), mypcolor(tmp3'); thincolorbar
44      title('difference'), pause(.1)      title('difference'), pause(.2)
45      J=J+sum(sum(tmp3(2:(end-1),2:(end-1)).^2));      J=J+sum(sum(tmp3(2:(end-1),2:(end-1)).^2));
46    end    end
47    disp([f{1} ': ' num2str(J)])    disp([f{1} ': ' num2str(J)])
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22