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

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

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

revision 1.7 by jmc, Mon May 2 19:35:38 2011 UTC revision 1.8 by jmc, Thu May 19 00:03:06 2011 UTC
# Line 81  uEast = cat(3,(u0-du)*ones(ny,nz),(u0+du Line 81  uEast = cat(3,(u0-du)*ones(ny,nz),(u0+du
81  w1=[0:nz-1]*pi/nz; w1=-w0*sin(w1);  w1=[0:nz-1]*pi/nz; w1=-w0*sin(w1);
82  wZonal = cat (3, ones(ny,1)*w1, zeros(ny,nz));  wZonal = cat (3, ones(ny,1)*w1, zeros(ny,nz));
83    
84    % to test prescribing Eta in NonLin-FreeSurf formulation:
85    et1=0.1;
86    etWest = cat(2,+et1*ones(ny,1),-et1*ones(ny,1));
87    etEast = cat(2,-et1*ones(ny,1),+et1*ones(ny,1));
88    
89   fid=fopen('OBmeridU.bin','w',ieee); fwrite(fid,uMerid,prec); fclose(fid);   fid=fopen('OBmeridU.bin','w',ieee); fwrite(fid,uMerid,prec); fclose(fid);
90  %fid=fopen('OBzonalU.bin','w',ieee); fwrite(fid,uZonal,prec); fclose(fid);  %fid=fopen('OBzonalU.bin','w',ieee); fwrite(fid,uZonal,prec); fclose(fid);
91   fid=fopen('OB_WestU.bin','w',ieee); fwrite(fid,uWest ,prec); fclose(fid);   fid=fopen('OB_WestU.bin','w',ieee); fwrite(fid,uWest ,prec); fclose(fid);
92   fid=fopen('OB_EastU.bin','w',ieee); fwrite(fid,uEast ,prec); fclose(fid);   fid=fopen('OB_EastU.bin','w',ieee); fwrite(fid,uEast ,prec); fclose(fid);
93   fid=fopen('OBzonalS.bin','w',ieee); fwrite(fid,sZonal,prec); fclose(fid);   fid=fopen('OBzonalS.bin','w',ieee); fwrite(fid,sZonal,prec); fclose(fid);
94   fid=fopen('OBzonalW.bin','w',ieee); fwrite(fid,wZonal,prec); fclose(fid);   fid=fopen('OBzonalW.bin','w',ieee); fwrite(fid,wZonal,prec); fclose(fid);
95     fid=fopen('OB_WestH.bin','w',ieee); fwrite(fid,etWest,prec); fclose(fid);
96     fid=fopen('OB_EastH.bin','w',ieee); fwrite(fid,etEast,prec); fclose(fid);
97    
98  %- rbcs mask & restauring tracer field:  %- rbcs mask & restauring tracer field:
99  msk=ones(nx,ny,nz);  msk=ones(nx,ny,nz);
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22