Parent Directory
| 
Revision Log
| 
Revision Graph
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 | % need to be in verification/lab_sea/matlab directory |
| 2 | % and to specify location of exp1 output |
| 3 | cd ../../../verification/lab_sea/matlab |
| 4 | exp1_path='../input/exp1/'; |
| 5 | exp2_path='../input/exp2/'; |
| 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 '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 | clf, subplot(321) |
| 37 | mypcolor(lon,lat,area2'); caxis([0 1]), colorbar |
| 38 | title('Sea-ice concentration using LSR solver') |
| 39 | set(gca,'xticklabel',[]) |
| 40 | |
| 41 | subplot(322) |
| 42 | mypcolor(lon,lat,area2'-area1'); colorbar |
| 43 | title('Difference with ADI solver') |
| 44 | set(gca,'xticklabel',[]) |
| 45 | |
| 46 | subplot(323) |
| 47 | mypcolor(lon,lat,heff2'); caxis([0 .3]), colorbar |
| 48 | title('Effective sea-ice thickness (m)') |
| 49 | set(gca,'xticklabel',[]) |
| 50 | |
| 51 | subplot(324) |
| 52 | mypcolor(lon,lat,heff2'-heff1'); colorbar |
| 53 | title('Difference with ADI solver') |
| 54 | set(gca,'xticklabel',[]) |
| 55 | |
| 56 | subplot(325) |
| 57 | mypcolor(Blon,Blat,icespeed2'); caxis([0 10]), colorbar |
| 58 | hold on, myquiver(Blon,Blat,udir2',vdir2','k'); |
| 59 | title('Sea-ice velocity (cm/s)') |
| 60 | |
| 61 | subplot(326) |
| 62 | mypcolor(Blon,Blat,icespeed2'-icespeed1'); colorbar |
| 63 | if mmax(abs(udir2-udir1)) | mmax(abs(vdir2-vdir1)) |
| 64 | hold on, myquiver(Blon,Blat,udir2'-udir1',vdir2'-vdir1','k'); |
| 65 | end |
| 66 | title('Difference with ADI solver') |
| ViewVC Help | |
| Powered by ViewVC 1.1.22 |