/[MITgcm]/MITgcm/verification/offline_exf_seaice/input/getdata.m
ViewVC logotype

Contents of /MITgcm/verification/offline_exf_seaice/input/getdata.m

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

Revision 1.1 - (show annotations) (download)
Sat Dec 8 00:36:34 2012 UTC (11 years, 4 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint64y, checkpoint64x, checkpoint64z, checkpoint64o, checkpoint64q, checkpoint64p, checkpoint64s, checkpoint64r, checkpoint64u, checkpoint64t, checkpoint64w, checkpoint64v, checkpoint66g, checkpoint66f, checkpoint66e, checkpoint66d, checkpoint66c, checkpoint66b, checkpoint66a, checkpoint64n, checkpoint66o, checkpoint66n, checkpoint66m, checkpoint66l, checkpoint66k, checkpoint66j, checkpoint66i, checkpoint66h, checkpoint64b, checkpoint64e, checkpoint64d, checkpoint64c, checkpoint64g, checkpoint64f, checkpoint65z, checkpoint65x, checkpoint65y, checkpoint65r, checkpoint65s, checkpoint65p, checkpoint65q, checkpoint65v, checkpoint65w, checkpoint65t, checkpoint65u, checkpoint65j, checkpoint65k, checkpoint65h, checkpoint65i, checkpoint65n, checkpoint65o, checkpoint65l, checkpoint65m, checkpoint65b, checkpoint65c, checkpoint65a, checkpoint65f, checkpoint65g, checkpoint65d, checkpoint65e, checkpoint64i, checkpoint64h, checkpoint64k, checkpoint65, checkpoint64j, checkpoint64m, checkpoint64l, HEAD 
change (for now, only fwd test) offline_exf_seaice experiment:
 re-entrant channel with idealized costline (trapezoidal shape) at southern
 boundary to explore seaice dynamics and thermodynamics with simple forcing.
1
2 it=720;
3
4 rDir='res_o16/'; sfx='_3c0.bin';
5 u=rdmds([rDir,'U'],it);
6 v=rdmds([rDir,'V'],it);
7 et=rdmds([rDir,'Eta'],it);
8 u1=u(:,:,1);
9 v1=v(:,:,1);
10
11 namf=['uVel',sfx]; var=u1;
12 fprintf('write to file: %s\n',namf);
13 fid=fopen(namf,'w','b'); fwrite(fid,var,'real*8'); fclose(fid);
14 namf=['vVel',sfx]; var=v1;
15 fprintf('write to file: %s\n',namf);
16 fid=fopen(namf,'w','b'); fwrite(fid,var,'real*8'); fclose(fid);
17 namf=['eta',sfx]; var=et;
18 fprintf('write to file: %s\n',namf);
19 fid=fopen(namf,'w','b'); fwrite(fid,var,'real*8'); fclose(fid);
20
21 nx=80; ny=42; nr=3;
22 xax=[1:nx]-.5 ; yax=[1:ny]-.5;
23 xtxt=1; ytxt=-2;
24
25 nf=1; ccB=[0 0];
26 figure(nf);clf;
27 subplot(211);
28 var=v1;
29 var=et;
30 imagesc(xax,yax,var'); set(gca,'YDir','normal');
31 if ccB(2) > ccB(1), caxis(ccB); end
32 change_colmap(-1);
33 colorbar;
34 grid
35 titv=['dir: ',strrep(rDir,'_','\_'),' ; sfx=',strrep(sfx,'_','\_')];
36 title(titv);
37
38 rDir='res_o06/'; sfx='_3c1.bin';
39 u=rdmds([rDir,'U'],it);
40 v=rdmds([rDir,'V'],it);
41 et=rdmds([rDir,'Eta'],it);
42 u1=u(:,:,1);
43 v1=v(:,:,1);
44
45 namf=['uVel',sfx]; var=u1;
46 fprintf('write to file: %s\n',namf);
47 fid=fopen(namf,'w','b'); fwrite(fid,var,'real*8'); fclose(fid);
48 namf=['vVel',sfx]; var=v1;
49 fprintf('write to file: %s\n',namf);
50 fid=fopen(namf,'w','b'); fwrite(fid,var,'real*8'); fclose(fid);
51 namf=['eta',sfx]; var=et;
52 fprintf('write to file: %s\n',namf);
53 fid=fopen(namf,'w','b'); fwrite(fid,var,'real*8'); fclose(fid);
54
55 subplot(212);
56 var=v1;
57 var=et;
58 imagesc(xax,yax,var'); set(gca,'YDir','normal');
59 if ccB(2) > ccB(1), caxis(ccB); end
60 change_colmap(-1);
61 colorbar;
62 grid
63 titv=['dir: ',strrep(rDir,'_','\_'),' ; sfx=',strrep(sfx,'_','\_')];
64 title(titv);
65
  ViewVC Help
Powered by ViewVC 1.1.22