Changeset 26721

Timestamp:

12/08/21 19:01:45 (3 years ago)

Author:

dmenemen

Message:

moving ocean boundary to North edge of domain

File:

• issm/trunk-jpl/test/NightlyRun/test4003.m (modified) (4 diffs)

issm/trunk-jpl/test/NightlyRun/test4003.m

@@ -4 +4 @@
 %Script control parameters
 steps=1:11;
+steps=[1:3 6:11];
 
 %To download and recompile MITgcm from scratch:
-!rm -rf ${ISSM_DIR}/test/MITgcm/install
-!rm -rf ${ISSM_DIR}/test/MITgcm/build/*
-!rm -rf Models
+%!rm -rf ${ISSM_DIR}/test/MITgcm/install
+%!rm -rf ${ISSM_DIR}/test/MITgcm/build/*
+%!rm -rf Models
 
 %Organizer
@@ -61 +62 @@
     %create bathymetry:
     bathymetry = ones(Nx,Ny)*H;
-    bathymetry(:,1) = 0;
+    bathymetry(:,end) = 0;
 
     %save bathymetry file for MITgcm
@@ -80 +81 @@
     for i=1:Nx
         draft(i,:)=Hmin+dHdy*[-1:Ny-2]*dLat;
-        draft(i,1)=0; draft(i,j2:Ny)=0;
     end
 
@@ -86 +86 @@
     ice_mask(:,j2:Ny)=0;
     iceshelf_mask=ice_mask;
-    iceshelf_mask(:,1)=0;
     thickness=abs(draft)/di;
-    thickness(:,1)=thickness(:,2)+draft(:,2);
     
     savedata(org,ice_mask,iceshelf_mask,draft,thickness);
+    
+    close all, figure(2), clf
+    subplot(411), pcolorcen(bathymetry); colorbar, title('bathymetry')
+    subplot(412), pcolorcen(ice_mask); colorbar, title('ice and iceshelf mask')
+    subplot(413), pcolorcen(draft); colorbar, title('draft')
+    subplot(414), pcolorcen(thickness); colorbar, title('thickness')
+    
 end
 % }}}