/[MITgcm]/MITgcm_contrib/verification_other/offline_cheapaml/input/gendata.m
ViewVC logotype

Diff of /MITgcm_contrib/verification_other/offline_cheapaml/input/gendata.m

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

revision 1.1 by jmc, Wed May 22 19:39:15 2013 UTC revision 1.2 by jmc, Tue Jun 11 01:23:52 2013 UTC
# Line 4  kwr=0; Line 4  kwr=0;
4  nx=80; ny=42; nr=3; nt=1;  nx=80; ny=42; nr=3; nt=1;
5    
6  xc=[1:nx]; xc=xc-mean(xc);  xc=[1:nx]; xc=xc-mean(xc);
7  yc=[1:ny]-.5;  yc=[1:ny]-.5; ymid=mean(yc);
8    yv=yc-0.5;
9    
10  %------------------------------------------------------  %------------------------------------------------------
11    
# Line 18  if kwr > 0, Line 19  if kwr > 0,
19   fid=fopen(namf,'w','b'); fwrite(fid,depth,'real*8'); fclose(fid);   fid=fopen(namf,'w','b'); fwrite(fid,depth,'real*8'); fclose(fid);
20  end  end
21    
22  namf='windx.bin';  namf=['windx_',int2str(windx),'ms.bin'];
23  wnd=windx*ones(nx,ny,nt);  uwind=windx*ones(nx,ny,nt);
24  if kwr > 0,  if kwr > 0,
25   fprintf('write to file: %s\n',namf);   fprintf('write to file: %s\n',namf);
26   fid=fopen(namf,'w','b'); fwrite(fid,wnd,'real*8'); fclose(fid);   fid=fopen(namf,'w','b'); fwrite(fid,uwind,'real*8'); fclose(fid);
27    end
28    
29    namf='windy_conv.bin';
30    dvdy=-1.e-6*5.e+3; %- uniform convergence: wWind = 10^-6 m/s
31    yy=yv-ymid; vwind=dvdy*yy; vwind(1)=0;
32    fld=ones(nx,1)*vwind;
33    if kwr > 0,
34     fprintf('write to file: %s\n',namf);
35     fid=fopen(namf,'w','b'); fwrite(fid,fld,'real*8'); fclose(fid);
36  end  end
37    
38  %- file name convention: "const_{xx}.bin" <-> uniform value = xx (in percent)  %- file name convention: "const_{xx}.bin" <-> uniform value = xx (in percent)
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22