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

Diff of /MITgcm/verification/offline_exf_seaice/input/gendata.m

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

revision 1.1 by jmc, Sat Dec 8 00:36:34 2012 UTC revision 1.2 by gforget, Thu Dec 27 14:10:03 2012 UTC
# Line 12  H0=-100.; Line 12  H0=-100.;
12    
13  namf='channel.bin';  namf='channel.bin';
14  depth=H0*ones(nx,ny); depth(:,1)=0.;  depth=H0*ones(nx,ny); depth(:,1)=0.;
15  %if kwr > 0,  if kwr > 0,
16  % fprintf('write to file: %s\n',namf);   fprintf('write to file: %s\n',namf);
17  % fid=fopen(namf,'w','b'); fwrite(fid,depth,'real*8'); fclose(fid);   fid=fopen(namf,'w','b'); fwrite(fid,depth,'real*8'); fclose(fid);
18  %end  end
19    
20  namf='bathy_3c.bin';  namf='bathy_3c.bin';
21  msk=abs(xc)'*ones(1,ny)+ones(nx,1)*yc;  msk=abs(xc)'*ones(1,ny)+ones(nx,1)*yc;
# Line 74  end Line 74  end
74  %end  %end
75  %------------------------------------------------------  %------------------------------------------------------
76    
77    namf='edge100.bin'; w0=1.;
78    var=w0*ones(nx,ny);
79    var(:,2)=0*var(:,end);
80    var(:,3)=0.01*var(:,end-1);
81    var(:,4)=0.1*var(:,end-2);
82    var(:,end)=0*var(:,end);
83    var(:,end-1)=0.01*var(:,end-1);
84    var(:,end-2)=0.1*var(:,end-2);
85    if kwr >0,
86     fprintf('write to file: %s\n',namf);
87     fid=fopen(namf,'w','b'); fwrite(fid,var,'real*8'); fclose(fid);
88    end
89    
90    namf='edge+20.bin'; w0=0.2;
91    var=w0*ones(nx,ny);
92    var(:,2)=0*var(:,end);
93    var(:,3)=0.01*var(:,end-1);
94    var(:,4)=0.1*var(:,end-2);
95    var(:,end)=0*var(:,end);
96    var(:,end-1)=0.01*var(:,end-1);
97    var(:,end-2)=0.1*var(:,end-2);
98    if kwr >0,
99     fprintf('write to file: %s\n',namf);
100     fid=fopen(namf,'w','b'); fwrite(fid,var,'real*8'); fclose(fid);
101    end
102    
103    %------------------------------------------------------
104    
105    
106  dsw0=100;  dsw0=100;
107  namf=['dsw_',int2str(dsw0),'.bin'];  namf=['dsw_',int2str(dsw0),'.bin'];
108  fld=dsw0*ones(nx,ny,nt);  fld=dsw0*ones(nx,ny,nt);
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22