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 and exp2 output | ||
3 | cd ../../../verification/lab_sea/matlab | ||
4 | dimitri | 1.3 | exp1_path='../input/exp1/'; |
5 | exp2_path='../input/exp4/'; | ||
6 | heimbach | 1.2 | |
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 | % plot comparison figures | ||
24 | clf, orient tall, wysiwyg | ||
25 | subplot(421), mypcolor(lon,lat,area2'); colorbar | ||
26 | title('Ice concentration during first 10 hours') | ||
27 | subplot(422), mypcolor(lon,lat,area2'-area1'); colorbar | ||
28 | title('2cpu - 1cpu solution') | ||
29 | subplot(423), mypcolor(lon,lat,heff2'); colorbar | ||
30 | title('Effective ice thickness') | ||
31 | subplot(424), mypcolor(lon,lat,heff2'-heff1'); colorbar | ||
32 | title('2cpu - 1cpu solution') | ||
33 | subplot(425), mypcolor(lon,lat,uice2'); colorbar | ||
34 | title('Uice') | ||
35 | subplot(426), mypcolor(lon,lat,uice2'-uice1'); colorbar | ||
36 | title('2cpu - 1cpu solution') | ||
37 | dimitri | 1.3 | subplot(427), mypcolor(lon,lat,vice2'); colorbar |
38 | heimbach | 1.2 | title('Vice') |
39 | subplot(428), mypcolor(lon,lat,vice2'-vice1'); colorbar | ||
40 | title('2cpu - 1cpu solution') |
ViewVC Help | |
Powered by ViewVC 1.1.22 |