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 | % 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/'; |
| 5 | |
| 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 |