/[MITgcm]/MITgcm/verification/lab_sea/matlab/lookat_exp2.m
ViewVC logotype

Annotation of /MITgcm/verification/lab_sea/matlab/lookat_exp2.m

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

Revision 1.1.2.2 - (hide annotations) (download)
Thu Dec 5 08:06:18 2002 UTC (21 years, 6 months ago) by dimitri
Branch: release1
CVS Tags: release1_p12, release1_p13, release1_p10, release1_p11, release1_p14, release1_p15, release1_p13_pre, release1_p12_pre, release1_p9
Branch point for: release1_50yr
Changes since 1.1.2.1: +14 -6 lines 
release1_p9
o pkg/seaice
  - removed GOTO's and added taf directives
  - double precision constants to reduce the g77 (Linux)
    to F77 (SGI) differences reported in release1_p8
o tools/genmake
  - added SGI options
o verification/testscript
  - updated to that of checkpoint47a_post
o verification/global_ocean.90x40x15/input/eedata
  - modified for SGI f77 compatibility
o verification/lab_sea
  - added description of sea-ice model
  - added missing matlab routines
  - added test of thermodynamics parallelization
Modified Files:
 Tag: release1
   doc/tag-index pkg/seaice/SEAICE_FFIELDS.h
   pkg/seaice/SEAICE_PARAMS.h pkg/seaice/adi.F
   pkg/seaice/advect.F pkg/seaice/budget.F pkg/seaice/diffus.F
   pkg/seaice/dynsolver.F pkg/seaice/groatb.F pkg/seaice/growth.F
   pkg/seaice/lsr.F pkg/seaice/ostres.F
   pkg/seaice/seaice_do_diags.F pkg/seaice/seaice_get_forcing.F
   pkg/seaice/seaice_init.F pkg/seaice/seaice_model.F
   pkg/seaice/seaice_readparms.F tools/genmake
   verification/testscript
   verification/global_ocean.90x40x15/input/eedata
   verification/lab_sea/README
   verification/lab_sea/code/CPP_EEOPTIONS.h
   verification/lab_sea/code/CPP_EEOPTIONS_MPI.h
   verification/lab_sea/code/CPP_OPTIONS.h
   verification/lab_sea/code/SIZE.h
   verification/lab_sea/code/SIZE_2x1.h
   verification/lab_sea/input/data
   verification/lab_sea/matlab/lookat_exp1.m
   verification/lab_sea/matlab/lookat_exp2.m
   verification/lab_sea/matlab/lookat_exp3.m
   verification/lab_sea/matlab/lookat_exp4.m
   verification/lab_sea/matlab/lookat_exp5.m
   verification/lab_sea/matlab/lookat_exp6.m
   verification/lab_sea/results/AREAtave.0000000010.data
   verification/lab_sea/results/HEFFtave.0000000010.data
   verification/lab_sea/results/UICEtave.0000000010.data
   verification/lab_sea/results/VICEtave.0000000010.data
   verification/lab_sea/results/output.txt
Added Files:
 Tag: release1
   verification/lab_sea/seaice.ps
   verification/lab_sea/matlab/lookat_exp7.m
   verification/lab_sea/matlab/mmax.m
   verification/lab_sea/matlab/mypcolor.m
   verification/lab_sea/matlab/myquiver.m
   verification/lab_sea/matlab/readbin.m
   verification/lab_sea/matlab/wysiwyg.m
Removed Files:
 Tag: release1
   verification/lab_sea/code/KPP_OPTIONS.h
1 heimbach 1.1.2.1 % need to be in verification/lab_sea/matlab directory
2     % and to specify location of exp1 output
3     cd ../../../verification/lab_sea/matlab
4 dimitri 1.1.2.2 exp1_path='../input/exp1/';
5     exp2_path='../input/exp2/';
6 heimbach 1.1.2.1
7     % load monthly-mean SMMR-SSM/I data
8     load SSMI
9    
10     % B-grid latitude for uice and vice
11     Blat=lat+1; Blon=lon+1;
12    
13     % load model output
14     area1=readbin([exp1_path 'AREAtave.0000000010.data'],[20 16 1],1);
15     heff1=readbin([exp1_path 'HEFFtave.0000000010.data'],[20 16 1],1);
16     uice1=readbin([exp1_path 'UICEtave.0000000010.data'],[20 16 1],1);
17     vice1=readbin([exp1_path 'VICEtave.0000000010.data'],[20 16 1],1);
18     area2=readbin([exp2_path 'AREAtave.0000000010.data'],[20 16 1],1);
19     heff2=readbin([exp2_path 'HEFFtave.0000000010.data'],[20 16 1],1);
20     uice2=readbin([exp2_path 'UICEtave.0000000010.data'],[20 16 1],1);
21     vice2=readbin([exp2_path 'VICEtave.0000000010.data'],[20 16 1],1);
22    
23     % compute ice speed and direction
24     udir1=0*uice1; vdir1=0*vice1;
25     udir2=0*uice2; vdir2=0*vice2;
26     icespeed1=100*sqrt(uice1.^2+vice1.^2);
27     in=find(icespeed1);
28     udir1(in)=uice1(in)./icespeed1(in);
29     vdir1(in)=vice1(in)./icespeed1(in);
30     icespeed2=100*sqrt(uice2.^2+vice2.^2);
31     in=find(icespeed2);
32     udir2(in)=uice2(in)./icespeed2(in);
33     vdir2(in)=vice2(in)./icespeed2(in);
34    
35     % plot comparison figures
36 dimitri 1.1.2.2 clf, subplot(321)
37 heimbach 1.1.2.1 mypcolor(lon,lat,area2'); caxis([0 1]), colorbar
38     title('Sea-ice concentration using LSR solver')
39 dimitri 1.1.2.2 set(gca,'xticklabel',[])
40    
41 heimbach 1.1.2.1 subplot(322)
42     mypcolor(lon,lat,area2'-area1'); colorbar
43     title('Difference with ADI solver')
44 dimitri 1.1.2.2 set(gca,'xticklabel',[])
45    
46 heimbach 1.1.2.1 subplot(323)
47     mypcolor(lon,lat,heff2'); caxis([0 .3]), colorbar
48     title('Effective sea-ice thickness (m)')
49 dimitri 1.1.2.2 set(gca,'xticklabel',[])
50    
51 heimbach 1.1.2.1 subplot(324)
52     mypcolor(lon,lat,heff2'-heff1'); colorbar
53     title('Difference with ADI solver')
54 dimitri 1.1.2.2 set(gca,'xticklabel',[])
55    
56 heimbach 1.1.2.1 subplot(325)
57     mypcolor(Blon,Blat,icespeed2'); caxis([0 10]), colorbar
58 dimitri 1.1.2.2 hold on, myquiver(Blon,Blat,udir2',vdir2','k');
59 heimbach 1.1.2.1 title('Sea-ice velocity (cm/s)')
60 dimitri 1.1.2.2
61 heimbach 1.1.2.1 subplot(326)
62     mypcolor(Blon,Blat,icespeed2'-icespeed1'); colorbar
63     if mmax(abs(udir2-udir1)) | mmax(abs(vdir2-vdir1))
64 dimitri 1.1.2.2 hold on, myquiver(Blon,Blat,udir2'-udir1',vdir2'-vdir1','k');
65 heimbach 1.1.2.1 end
66     title('Difference with ADI solver')
  ViewVC Help
Powered by ViewVC 1.1.22