/[MITgcm]/MITgcm_contrib/eh3/llc/ecco-godae/climatology/gen_bathy.m
ViewVC logotype

Diff of /MITgcm_contrib/eh3/llc/ecco-godae/climatology/gen_bathy.m

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

revision 1.3 by edhill, Thu Aug 24 15:30:24 2006 UTC revision 1.5 by heimbach, Wed Nov 15 01:27:17 2006 UTC
# Line 56  end Line 56  end
56    
57  %  Do a quick-and-dirty regrid  %  Do a quick-and-dirty regrid
58  disp('Performing the regrid...');  disp('Performing the regrid...');
59  s_lon = -37;  %ph(
60  s_lat =   0.0;  %%%s_lon = 0.0;
61    s_lon = -37.0;
62    %ph)
63    s_lat = 0.0;
64  for k = 1:5  for k = 1:5
65    disp(sprintf('  k = %d',k));    disp(sprintf('  k = %d',k));
66    %    %
# Line 142  for k = 1:5 Line 145  for k = 1:5
145  end  end
146  axis equal, view(2)  axis equal, view(2)
147    
148    %--- BREAK POINT HERE
149    %--- do e.g. view(80,45) to change view angle
150    
151    
152    
153  if do_print > 0  if do_print > 0
154    print -depsc llc_bathy.eps    print -depsc llc_bathy.eps
# Line 192  iti = [   1    1   1  ; ... Line 199  iti = [   1    1   1  ; ...
199  for k = 1:size(iti,1)  for k = 1:size(iti,1)
200    disp(sprintf('  k = %d',k));    disp(sprintf('  k = %d',k));
201    fname = sprintf('bathy.%03d.001.data',k);    fname = sprintf('bathy.%03d.001.data',k);
202      disp( fname )
203    fid = fopen(fname, 'w', 'ieee-be');    fid = fopen(fname, 'w', 'ieee-be');
   %  be careful of index orientation  
204    is = iti(k,3);  ie = is + 90 - 1;    is = iti(k,3);  ie = is + 90 - 1;
205    js = iti(k,2);  je = js + 90 - 1;    js = iti(k,2);  je = js + 90 - 1;
206    tmp = ginfo(iti(k,1)).b_out(is:ie,js:je)';    tmp = ginfo(iti(k,1)).b_out(is:ie,js:je);
   surf(tmp), view(2), axis equal, shading flat, pause(1)  
207    fwrite(fid,tmp,'real*8',0,'ieee-be');    fwrite(fid,tmp,'real*8',0,'ieee-be');
208    fclose(fid);    fclose(fid);
209      surf(tmp), view(2), axis equal, shading flat, pause(1)
210  end  end
211    
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22