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

Diff of /MITgcm_contrib/shelfice_remeshing/MANUAL/input/Remeshing.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:38 2015 UTC revision 1.2 by dgoldberg, Tue Oct 13 15:56:34 2015 UTC
# Line 5  clc Line 5  clc
5  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
6  % Read in and reshape files  % Read in and reshape files
7  fid = fopen('pickup.OLD.data','r','b');  fid = fopen('pickup.OLD.data','r','b');
8  pickupold = fread(fid,inf,'real*8');  pickupold = fread(fid,inf,'float64');
9  fclose(fid);  fclose(fid);
10    
11  fid = fopen('pload.OLD.jmd95z','r','b');  fid = fopen('pload.OLD.jmd95z','r','b');
# Line 53  stpigold=reshape(stpigold,[nx ny]); Line 53  stpigold=reshape(stpigold,[nx ny]);
53  %Change eta, iceshelf topo and r0surfto reflect change in mesh  %Change eta, iceshelf topo and r0surfto reflect change in mesh
54    
55  splitthreshold=5;  splitthreshold=5;
56  mergethreshold=-6;  mergethreshold=-5.6;
57    
58    
59  [splitx,splity] =find(oldetan>splitthreshold & stroundold<0);  [splitx,splity] =find(oldetan>splitthreshold & stroundold<0);
# Line 231  for i=1:length(OMERGEX) Line 231  for i=1:length(OMERGEX)
231      newtheta(OMERGEX(i),OMERGEY(i),NMERGEZ(i))=((oldtheta(OMERGEX(i),OMERGEY(i),OMERGEZ(i))*(dz+oldetan(OMERGEX(i),OMERGEY(i))))...      newtheta(OMERGEX(i),OMERGEY(i),NMERGEZ(i))=((oldtheta(OMERGEX(i),OMERGEY(i),OMERGEZ(i))*(dz+oldetan(OMERGEX(i),OMERGEY(i))))...
232          +(dz*oldtheta(OMERGEX(i),OMERGEY(i),NMERGEZ(i))))/(2*dz+oldetan(OMERGEX(i),OMERGEY(i)));          +(dz*oldtheta(OMERGEX(i),OMERGEY(i),NMERGEZ(i))))/(2*dz+oldetan(OMERGEX(i),OMERGEY(i)));
233            
234         newsalt(OMERGEX(i),OMERGEY(i),OMERGEZ(i))=0;
235        newtheta(OMERGEX(i),OMERGEY(i),OMERGEZ(i))=0;
236        
237      newuvel(OMERGEX(i),OMERGEY(i),NMERGEZ(i))=((olduvel(OMERGEX(i),OMERGEY(i),OMERGEZ(i))*(dz+oldetan(OMERGEX(i),OMERGEY(i))))...      newuvel(OMERGEX(i),OMERGEY(i),NMERGEZ(i))=((olduvel(OMERGEX(i),OMERGEY(i),OMERGEZ(i))*(dz+oldetan(OMERGEX(i),OMERGEY(i))))...
238          +(dz*olduvel(OMERGEX(i),OMERGEY(i),NMERGEZ(i))))/(2*dz+oldetan(OMERGEX(i),OMERGEY(i)));          +(dz*olduvel(OMERGEX(i),OMERGEY(i),NMERGEZ(i))))/(2*dz+oldetan(OMERGEX(i),OMERGEY(i)));
239      %newuvel(OMERGEX(i)+1,OMERGEY(i),NMERGEZ(i))=((olduvel(OMERGEX(i)+1,OMERGEY(i),OMERGEZ(i))*(dz+oldetan(OMERGEX(i)+1,OMERGEY(i))))...      %newuvel(OMERGEX(i)+1,OMERGEY(i),NMERGEZ(i))=((olduvel(OMERGEX(i)+1,OMERGEY(i),OMERGEZ(i))*(dz+oldetan(OMERGEX(i)+1,OMERGEY(i))))...
# Line 268  stroundnew=reshape(stroundnew,[nx*ny 1]) Line 271  stroundnew=reshape(stroundnew,[nx*ny 1])
271  ploadnew=reshape(ploadnew,[nx*ny 1]);  ploadnew=reshape(ploadnew,[nx*ny 1]);
272    
273    
274  fid = fopen('pickup.0000012960.data','w','b'); fwrite(fid,pickupnew,'real*8'); fclose(fid);  fid = fopen('pickup.0000004320.data','w','b'); fwrite(fid,pickupnew,'float64'); fclose(fid);
275  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)
276  fid = fopen('shelftopo.round.bin','w','b'); fwrite(fid,stroundnew,'real*8'); fclose(fid)  fid = fopen('shelftopo.round.bin','w','b'); fwrite(fid,stroundnew,'real*8'); fclose(fid)
277    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22