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

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

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

Revision 1.2 - (hide annotations) (download)
Tue Feb 18 05:33:56 2003 UTC (21 years, 4 months ago) by dimitri
Branch: MAIN
CVS Tags: checkpoint48i_post, checkpoint50b_post, checkpoint48g_post, checkpoint50c_pre, checkpoint50, checkpoint50f_post, checkpoint50a_post, checkpoint50f_pre, checkpoint50d_pre, checkpoint50d_post, checkpoint50c_post, checkpoint50e_post, checkpoint50b_pre, checkpoint48h_post, checkpoint50e_pre, checkpoint49, checkpoint48f_post
Changes since 1.1: +66 -0 lines 
Merging from release1_p12:
o Modifications for using pkg/exf with pkg/seaice
  - improved description of the various forcing configurations
  - added basic radiation bulk formulae to pkg/exf
  - units/sign fix for evap computation in exf_getffields.F
  - updated verification/global_with_exf/results/output.txt
o Added pkg/sbo for computing IERS Special Bureau for the Oceans
  (SBO) core products, including oceanic mass, center-of-mass,
  angular, and bottom pressure (see pkg/sbo/README.sbo).
o Lower bound for viscosity/diffusivity in pkg/kpp/kpp_routines.F
  to avoid negative values in shallow regions.
  - updated verification/natl_box/results/output.txt
  - updated verification/lab_sea/results/output.txt
o MPI gather, scatter: eesupp/src/gather_2d.F and scatter_2d.F
o Added useSingleCpuIO option (see PARAMS.h).
o Updated useSingleCpuIO option in mdsio_writefield.F to
  work with multi-field files, e.g., for single-file pickup.
o pkg/seaice:
  - bug fix in growth.F: QNET for no shortwave case
  - added HeffFile for specifying initial sea-ice thickness
  - changed SEAICE_EXTERNAL_FLUXES wind stress implementation
o Added missing /* */ to CPP comments in pkg/seaice, pkg/exf,
  kpp_transport_t.F, forward_step.F, and the_main_loop.F
o pkg/seaice:
  - adjoint-friendly modifications
  - added a SEAICE_WRITE_PICKUP at end of the_model_main.F
1 dimitri 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     exp8_path='../input/exp8/';
5     exp10_path='../input/exp10/';
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([exp8_path 'AREAtave.0000000010.data'],[20 16 1],1);
15     heff1=readbin([exp8_path 'HEFFtave.0000000010.data'],[20 16 1],1);
16     uice1=readbin([exp8_path 'UICEtave.0000000010.data'],[20 16 1],1);
17     vice1=readbin([exp8_path 'VICEtave.0000000010.data'],[20 16 1],1);
18     area2=readbin([exp10_path 'AREAtave.0000000010.data'],[20 16 1],1);
19     heff2=readbin([exp10_path 'HEFFtave.0000000010.data'],[20 16 1],1);
20     uice2=readbin([exp10_path 'UICEtave.0000000010.data'],[20 16 1],1);
21     vice2=readbin([exp10_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 with pkg/exf evaporation')
39     set(gca,'xticklabel',[])
40    
41     subplot(322)
42     mypcolor(lon,lat,area2'-area1'); colorbar
43     title('Difference with Experiment 8')
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 Experiment 8')
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 Experiment 8')
  ViewVC Help
Powered by ViewVC 1.1.22