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

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

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

Revision 1.2 - (hide annotations) (download)
Thu Dec 5 08:43:03 2002 UTC (21 years, 5 months ago) by dimitri
Branch: MAIN
CVS Tags: checkpoint47j_post, checkpoint48d_pre, branch-exfmods-tag, checkpoint47e_post, checkpoint47i_post, checkpoint47f_post, checkpoint47c_post, checkpoint50e_post, checkpoint50c_post, checkpoint47d_post, checkpoint48e_post, checkpoint48d_post, checkpoint48f_post, checkpoint48c_post, checkpoint47d_pre, checkpoint48, checkpoint49, checkpoint47b_post, checkpoint48i_post, checkpoint48h_post, checkpoint50d_pre, checkpoint50, checkpoint47h_post, checkpoint50d_post, checkpoint47g_post, checkpoint48a_post, checkpoint50c_pre, checkpoint50b_pre, checkpoint48b_post, checkpoint50b_post, checkpoint50f_post, checkpoint50a_post, checkpoint50f_pre, checkpoint48c_pre, checkpoint50e_pre, checkpoint48g_post
Branch point for: branch-exfmods-curt, ecco-branch
Changes since 1.1: +58 -0 lines 
checkpoint47b_post
Merging from 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:
   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/global_ocean.90x40x15/input/eedata
   verification/lab_sea/README
   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:
   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:
   verification/lab_sea/code/KPP_OPTIONS.h
1 dimitri 1.2 % need to be in verification/lab_sea/matlab directory
2     % and to specify location of exp1 and exp2 output
3     cd ../../../verification/lab_sea/matlab
4     exp1_path='../input/exp6/';
5     exp2_path='../input/exp7/';
6    
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 'uVeltave.0000000010.data'],[20 16 23],1);
17     vice1=readbin([exp1_path 'vVeltave.0000000010.data'],[20 16 23],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 'uVeltave.0000000010.data'],[20 16 23],1);
21     vice2=readbin([exp2_path 'vVeltave.0000000010.data'],[20 16 23],1);
22     eta1=readbin([exp1_path 'ETAtave.0000000010.data'],[20 16 1],1);
23     eta2=readbin([exp2_path 'ETAtave.0000000010.data'],[20 16 1],1);
24     t1=readbin([exp1_path 'Ttave.0000000010.data'],[20 16 23],1);
25     t2=readbin([exp2_path 'Ttave.0000000010.data'],[20 16 23],1);
26     s1=readbin([exp1_path 'Stave.0000000010.data'],[20 16 23],1);
27     s2=readbin([exp2_path 'Stave.0000000010.data'],[20 16 23],1);
28    
29     % plot comparison figures
30     clf, orient tall, wysiwyg, k=1;
31     subplot(721), mypcolor(lon,lat,area2'); colorbar
32     title('Ice concentration during first 10 hours'),set(gca,'xticklabel',[])
33     subplot(722), mypcolor(lon,lat,area2'-area1'); colorbar
34     title('2cpu - 1cpu solution'),set(gca,'xticklabel',[])
35     subplot(723), mypcolor(lon,lat,heff2'); colorbar
36     title('Effective ice thickness'),set(gca,'xticklabel',[])
37     subplot(724), mypcolor(lon,lat,heff2'-heff1'); colorbar
38     title('2cpu - 1cpu solution'),set(gca,'xticklabel',[])
39     subplot(725), mypcolor(lon,lat,uice2(:,:,k)'); colorbar
40     title('uVel at level 1'),set(gca,'xticklabel',[])
41     subplot(726), mypcolor(lon,lat,uice2(:,:,k)'-uice1(:,:,k)'); colorbar
42     title('2cpu - 1cpu solution'),set(gca,'xticklabel',[])
43     subplot(727), mypcolor(lon,lat,vice2(:,:,k)'); colorbar
44     title('vVel at level 1'),set(gca,'xticklabel',[])
45     subplot(728), mypcolor(lon,lat,vice2(:,:,k)'-vice1(:,:,k)'); colorbar
46     title('2cpu - 1cpu solution'),set(gca,'xticklabel',[])
47     subplot(7,2,9), mypcolor(lon,lat,eta2'); colorbar
48     title('SSH'),set(gca,'xticklabel',[])
49     subplot(7,2,10), mypcolor(lon,lat,eta2'-eta1'); colorbar
50     title('2cpu - 1cpu solution'),set(gca,'xticklabel',[])
51     subplot(7,2,11), mypcolor(lon,lat,t2(:,:,k)'); colorbar
52     title('theta at level 1'),set(gca,'xticklabel',[])
53     subplot(7,2,12), mypcolor(lon,lat,t2(:,:,k)'-t1(:,:,k)'); colorbar
54     title('2cpu - 1cpu solution'),set(gca,'xticklabel',[])
55     subplot(7,2,13), mypcolor(lon,lat,s2(:,:,k)'); colorbar
56     title('salt at level 1')
57     subplot(7,2,14), mypcolor(lon,lat,s2(:,:,k)'-s1(:,:,k)'); colorbar
58     title('2cpu - 1cpu solution')
  ViewVC Help
Powered by ViewVC 1.1.22