Parent Directory | Revision Log | Revision Graph
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 | heimbach | 1.2 | % 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.3 | exp1_path='../input/'; |
5 | heimbach | 1.2 | |
6 | % load monthly-mean SMMR-SSM/I data | ||
7 | load SSMI | ||
8 | |||
9 | % B-grid latitude for uice and vice | ||
10 | Blat=lat+1; Blon=lon+1; | ||
11 | |||
12 | % load model output | ||
13 | area=readbin([exp1_path 'AREAtave.0000010950.data'],[20 16 1],1); | ||
14 | heff=readbin([exp1_path 'HEFFtave.0000010950.data'],[20 16 1],1); | ||
15 | uice=readbin([exp1_path 'UICEtave.0000010950.data'],[20 16 1],1); | ||
16 | vice=readbin([exp1_path 'VICEtave.0000010950.data'],[20 16 1],1); | ||
17 | |||
18 | % compute ice speed and direction | ||
19 | udir=0*uice; vdir=0*vice; | ||
20 | icespeed=100*sqrt(uice.^2+vice.^2); | ||
21 | in=find(icespeed); | ||
22 | udir(in)=uice(in)./icespeed(in); | ||
23 | vdir(in)=vice(in)./icespeed(in); | ||
24 | |||
25 | % plot comparison figures | ||
26 | clf | ||
27 | subplot(221) | ||
28 | mypcolor(lon,lat,SSMI(:,:,4)'/100); caxis([0 1]), colorbar | ||
29 | title('Observed March sea-ice concentration') | ||
30 | subplot(222), mypcolor(lon,lat,area'); caxis([0 1]), colorbar | ||
31 | title('Modeled March sea-ice concentration') | ||
32 | subplot(223), mypcolor(lon,lat,heff'); caxis([0 1]), colorbar | ||
33 | title('Effective thickness in m') | ||
34 | subplot(224) | ||
35 | mypcolor(Blon,Blat,icespeed'); caxis([0 10]), colorbar | ||
36 | hold on,myquiver(Blon,Blat,udir',vdir','k'); | ||
37 | title('Sea-ice velocity (cm/s)') |
ViewVC Help | |
Powered by ViewVC 1.1.22 |