/[MITgcm]/MITgcm/verification/matrix_example/input/from_31x31_to_32x32.m
ViewVC logotype

Contents of /MITgcm/verification/matrix_example/input/from_31x31_to_32x32.m

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

Revision 1.1 - (show annotations) (download)
Thu Aug 25 22:36:38 2011 UTC (12 years, 8 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint64y, checkpoint64x, checkpoint64z, checkpoint64q, checkpoint64p, checkpoint64s, checkpoint64r, checkpoint64u, checkpoint64t, checkpoint64w, checkpoint64v, checkpoint64i, checkpoint64h, checkpoint64k, checkpoint64j, checkpoint64m, checkpoint64l, checkpoint64o, checkpoint64n, checkpoint64a, checkpoint64c, checkpoint64b, checkpoint64e, checkpoint64d, checkpoint64g, checkpoint64f, checkpoint63p, checkpoint63q, checkpoint63r, checkpoint63s, checkpoint63l, checkpoint63m, checkpoint63n, checkpoint63o, checkpoint63h, checkpoint63i, checkpoint63j, checkpoint63k, checkpoint63d, checkpoint63e, checkpoint63f, checkpoint63g, checkpoint63b, checkpoint63c, checkpoint64, checkpoint65, checkpoint66g, checkpoint66f, checkpoint66e, checkpoint66d, checkpoint66c, checkpoint66b, checkpoint66a, checkpoint66o, checkpoint66n, checkpoint66m, checkpoint66l, checkpoint66k, checkpoint66j, checkpoint66i, checkpoint66h, 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, HEAD 
add a matlab script to do the conversion 31x31 --> 32x32
(adding 1 row (on Southern side) + 1 column (on Western side) of land points)
1
2 %- To allow to split the domain in several tiles, change the domain size
3 % from 31x31 to 32x32 by adding 1 row (on Southern side) + 1 column (on
4 % Western side) of land points.
5
6 %- read old input files (31x31) and write new ones:
7 nx1=31; ny1=31;
8 nx2=1+nx1;
9 ny2=1+ny1;
10 prec='real*4';
11
12 fnam1='topog.box';
13 fnam2='topo_box.bin';
14 fprintf('converting %s to %s ...',fnam1,fnam2)
15 vv1=rdda(fnam1,[nx1 ny1],1,prec,'b');
16 vv2=zeros(nx2,ny2); vv2(2:nx2,2:ny2)=vv1;
17 fid=fopen(fnam2,'w','b'); fwrite(fid,vv2,prec); fclose(fid);
18 fprintf(' done\n');
19 %figure(1);clf; imagesc(vv2');set(gca,'YDir','normal');colorbar
20
21 fnam1='windx.cos_y';
22 fnam2='taux_cosY.bin';
23 fprintf('converting %s to %s ...',fnam1,fnam2)
24 vv1=rdda(fnam1,[nx1 ny1],1,prec,'b');
25 vv2=zeros(nx2,ny2); vv2(2:nx2,2:ny2)=vv1;
26 fid=fopen(fnam2,'w','b'); fwrite(fid,vv2,prec); fclose(fid);
27 fprintf(' done\n');
28
29 fnam1='trinit1';
30 fnam2='tr1_ini.bin';
31 fprintf('converting %s to %s ...',fnam1,fnam2)
32 vv1=rdda(fnam1,[nx1 ny1],1,prec,'b');
33 vv2=zeros(nx2,ny2); vv2(2:nx2,2:ny2)=vv1;
34 fid=fopen(fnam2,'w','b'); fwrite(fid,vv2,prec); fclose(fid);
35 fprintf(' done\n');
36
37 prec='real*8';
38 fnam1='pickup'; it=200000;
39 fnam2='pickup.32x32.data';
40 fprintf('converting %s to %s ...',fnam1,fnam2)
41 vv1=rdmds(fnam1,it); n3d=size(vv1,3);
42 vv2=zeros(nx2,ny2,n3d); vv2(2:nx2,2:ny2,:)=vv1;
43 fid=fopen(fnam2,'w','b'); fwrite(fid,vv2,prec); fclose(fid);
44 fprintf(' done\n');
45
  ViewVC Help
Powered by ViewVC 1.1.22