/[MITgcm]/MITgcm_contrib/shelfice_remeshing/MANUAL/input/rwpu.m
ViewVC logotype

Diff of /MITgcm_contrib/shelfice_remeshing/MANUAL/input/rwpu.m

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

revision 1.1 by dgoldberg, Thu Sep 10 14:56:39 2015 UTC revision 1.2 by dgoldberg, Tue Oct 13 15:56:35 2015 UTC
# Line 59  threshold=5; Line 59  threshold=5;
59  %etax=etax(3);  %etax=etax(3);
60  %etay=etay(3);  %etay=etay(3);
61    
62  for i = 1:length(etax)  for i = 1:size(etax,1)
63      etaz(i)=((shelftopoold(etax(i),etay(i)))/10*-1)+1;      etaz(i)=((shelftopoold(etax(i),etay(i)))/10*-1)+1;
64  end  end
65            
# Line 187  newgvnm1=oldgvnm1; Line 187  newgvnm1=oldgvnm1;
187    
188    
189    
190  for i=1:length(OLDX)  for i=1:size(OLDX)
191      newsalt(OLDX(i),OLDY(i),NEWZ(i))=oldsalt(OLDX(i),OLDY(i),OLDZ(i));      newsalt(OLDX(i),OLDY(i),NEWZ(i))=oldsalt(OLDX(i),OLDY(i),OLDZ(i));
192      newtheta(OLDX(i),OLDY(i),NEWZ(i))=oldtheta(OLDX(i),OLDY(i),OLDZ(i));      newtheta(OLDX(i),OLDY(i),NEWZ(i))=oldtheta(OLDX(i),OLDY(i),OLDZ(i));
193        
# Line 224  end Line 224  end
224     ploadnew=reshape(ploadnew,[nx*ny 1]);     ploadnew=reshape(ploadnew,[nx*ny 1]);
225    
226    
227   fid = fopen('pickup.0000012960.data','w','b'); fwrite(fid,pickupnew,'real*8'); fclose(fid);   fid = fopen('pickup.0000025920.data','w','b'); fwrite(fid,pickupnew,'real*8'); fclose(fid);
228   fid = fopen('pload.pig.jmd95z','w','b'); fwrite(fid,ploadnew,'real*8'); fclose(fid)   fid = fopen('pload.pig.jmd95z','w','b'); fwrite(fid,ploadnew,'real*8'); fclose(fid)
229   fid = fopen('shelftopo.round.bin','w','b'); fwrite(fid,shelftoponew,'real*8'); fclose(fid)   fid = fopen('shelftopo.round.bin','w','b'); fwrite(fid,shelftoponew,'real*8'); fclose(fid)
230    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22