Changeset 26863

Timestamp:

02/08/22 20:18:43 (3 years ago)

Author:

dmenemen

Message:

D test/MITgcm/input_4003/data_uncoupled
M test/NightlyRun/test4003.m

Location:

issm/trunk-jpl/test

Files:

• MITgcm/input_4003/data_uncoupled (deleted)
• NightlyRun/test4003.m (modified) (7 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:
@@ -41 +40 @@
     jEnd = Ny*3/4;  % where ice-shelf ends
     j2   = jEnd+1;
+    
+    % initial ocean conditions
+    T_sfc = -2;
+    T_bot = -1.9;
+    S_sfc = 34.2;
+    S_bot = 34.3;
 
     savedata(org, Nx, Ny, nPx, nPy, Nz, dLong, dLat, delZ, xgOrigin, ygOrigin, ...
-             rho_ice, rho_water, di, H, Hmin, Hmax, jEnd, j2, gravity, prec);
+             rho_ice, rho_water, di, H, Hmin, Hmax, jEnd, j2, gravity, prec, ...
+             T_sfc,T_bot,S_sfc,S_bot);
 end
 % }}}
@@ -82 +88 @@
         draft(i,:)=Hmin+dHdy*[-1:Ny-2]*dLat;
     end
+    draft(:,j2:Ny)=0;
 
     ice_mask=ones(Nx,Ny);
@@ -149 +156 @@
     fid=fopen(namF,'w','b'); fwrite(fid,dMdt,prec);fclose(fid);
 
-    T_sfc = -1.9;
-    T_bot = 2;
     del_T = (T_bot - T_sfc)/(59*delZ);
     tref=zeros(1,Nz);
@@ -161 +166 @@
     fid=fopen(namF,'w','b'); fwrite(fid,tref,prec);fclose(fid);
 
-    S_sfc = 34.2;
-    S_bot = 34.7;
     del_S = (S_bot - S_sfc)/(59*delZ);
     sref=zeros(1,Nz);
@@ -428 +431 @@
     fid=fopen(namF,'w','b'); fwrite(fid,dMdt,prec);fclose(fid);
 
-    T_sfc = -1.9;
-    T_bot = 2;
     del_T = (T_bot - T_sfc)/(59*delZ);
     tref=zeros(1,Nz);
@@ -440 +441 @@
     fid=fopen(namF,'w','b'); fwrite(fid,tref,prec);fclose(fid);
 
-    S_sfc = 34.2;
-    S_bot = 34.7;
     del_S = (S_bot - S_sfc)/(59*delZ);
     sref=zeros(1,Nz);