Changeset 26886

Timestamp:

02/13/22 14:07:27 (3 years ago)

Author:

dmenemen

Message:

first working check-in of test4003

Location:

issm/trunk-jpl/test

Files:

• MITgcm/code_4003/shelfice_step_icemass.F (modified) (1 diff)
• MITgcm/input_4003/data.shelfice (modified) (1 diff)
• MITgcm/tools/lookat_4003.m (modified) (2 diffs)
• NightlyRun/test4003.m (modified) (4 diffs)

issm/trunk-jpl/test/MITgcm/code_4003/shelfice_step_icemass.F

@@ -73 +73 @@
           DO j=1-OLy,sNy+OLy-1
            DO i=1-OLx+1,sNx+OLx-1
-             shelficeMass(i,j,bi,bj) = R_shelfIce(i,j,bi,bj) * rhoConst
+             shelficeMass(i,j,bi,bj) = -R_shelfIce(i,j,bi,bj) * rhoConst
            ENDDO
           ENDDO

issm/trunk-jpl/test/MITgcm/input_4003/data.shelfice

@@ -17 +17 @@
 #--
  SHELFICEtopoFile='shelficeTopo.Lin.bin',
- SHELFICEmassFile='shelficeMass.Lin.bin',
- SHELFICEMassDynTendFile='shelfice_dMdt.r02.bin',
-#SHELFICEwriteState = .TRUE.,
+ SHELFICEwriteState = .TRUE.,
  &

issm/trunk-jpl/test/MITgcm/tools/lookat_4003.m

@@ -1 +1 @@
 pn='~/mitgcm/slr/components/issm/trunk-jpl/test/NightlyRun/';
-p1=[pn 'RunUncoupledMITgcm/'];
-p2=[pn 'RunCoupledMITgcmISSM/'];
+p1=[pn 'RunUncoupled/'];
+p2=[pn 'run/'];
 
-fld='Eta';
-ts=1;
-v1=rdmds([p1 fld],ts);
-v2=rdmds([p2 fld],ts);
-clf
-subplot(311), mypcolor(v1); title(fld), colorbar
-subplot(312), mypcolor(v2); title('coupled'), colorbar
-subplot(313), mypcolor(v2-v1); title('coupled-uncoupled'), colorbar
-
-ts=8;
+ts=0;
 fld='R_shelfIce1_';
 v1=readbin([p2 fld myint2str(ts,10) '.data'],[3 200]);
 fld='R_shelfIce2_';
-v2=readbin([p2 fld myint2str(ts,10) '.data'],[3 200]);
-clf
-subplot(311), mypcolor(v1); title(fld), colorbar
-subplot(312), mypcolor(v2); title('coupled'), colorbar
-subplot(313), mypcolor(v2-v1); title('coupled-uncoupled'), colorbar
+for ts=0:8:184
+    v2=readbin([p2 fld myint2str(ts,10) '.data'],[3 200]);
+    clf
+    subplot(311), mypcolor(v1); title(ts-8), colorbar
+    subplot(312), mypcolor(v2); title(ts), colorbar
+    subplot(313), mypcolor(v2-v1); title('diff'), colorbar
+    pause
+    v1=v2;
+end
 
 fld='surfDiag';
@@ -36 +31 @@
     pause
 end
+
+fld='Eta';
+ts=1;
+v1=rdmds([p1 fld],ts);
+v2=rdmds([p2 fld],ts);
+clf
+subplot(311), mypcolor(v1); title(fld), colorbar
+subplot(312), mypcolor(v2); title('coupled'), colorbar
+subplot(313), mypcolor(v2-v1); title('coupled-uncoupled'), colorbar

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

@@ -144 +144 @@
     fid=fopen(namF,'w','b'); fwrite(fid,draft,prec);fclose(fid);
     
-    var=([1:Ny]-2)/(jEnd-2);
-    dMdt_fy=zeros(Nx,Ny);
-    for i=1:Nx
-        dMdt_fy(i,:)=-cos(pi*var);
-        dMdt_fy(i,1)=0; dMdt_fy(i,j2:Ny)=0;
-    end    
-    %- rate of change due to ice-stream dynamics
-    rateDyn=rho_water*0.1/3600;  sfx='r02';
-    dMdt=rateDyn*dMdt_fy;
-    namF=sprintf('%s.%s.%s','run/shelfice_dMdt',sfx,'bin');
-    fid=fopen(namF,'w','b'); fwrite(fid,dMdt,prec);fclose(fid);
-
     del_T = (T_bot - T_sfc)/(59*delZ);
     tref=zeros(1,Nz);
@@ -196 +184 @@
     namF='run/salt_ini.bin';
     fid=fopen(namF,'w','b'); fwrite(fid,var,prec);fclose(fid);
-
-    rhoAvr=rho_water-1.345;
-    mIce=-rhoAvr*draft;
-    namF='run/shelficeMass.Lin.bin';
-    fid=fopen(namF,'w','b'); fwrite(fid,mIce,prec);fclose(fid);
     % }}}
 
@@ -416 +399 @@
     fid=fopen(namF,'w','b'); fwrite(fid,bathymetry,prec);fclose(fid);
     
+    draft=md.results.TransientSolution(end).Base;
+    draft=reshape(draft,[Nx+1,Ny+1]);
+    pos=ones(Nx+1,Ny+1);
+    pos(find(md.mask.ice_levelset>0))=0;
+    draft=draft.*pos;
+    draft=draft(1:Nx,1:Ny)+draft(2:end,1:Ny)+draft(1:Nx,2:end)+draft(2:end,2:end);
+    pos=pos(1:Nx,1:Ny)+pos(2:end,1:Ny)+pos(1:Nx,2:end)+pos(2:end,2:end);
+    draft(find(pos))=draft(find(pos))./pos(find(pos));
     namF='run/shelficeTopo.Lin.bin';
     fid=fopen(namF,'w','b'); fwrite(fid,draft,prec);fclose(fid);
     
-    var=([1:Ny]-2)/(jEnd-2);
-    dMdt_fy=zeros(Nx,Ny);
-    for i=1:Nx
-        dMdt_fy(i,:)=-cos(pi*var);
-        dMdt_fy(i,1)=0; dMdt_fy(i,j2:Ny)=0;
-    end    
-    %- rate of change due to ice-stream dynamics
-    rateDyn=rho_water*0.1/3600;  sfx='r02';
-    dMdt=rateDyn*dMdt_fy;
-    namF=sprintf('%s.%s.%s','run/shelfice_dMdt',sfx,'bin');
-    fid=fopen(namF,'w','b'); fwrite(fid,dMdt,prec);fclose(fid);
-
     del_T = (T_bot - T_sfc)/(59*delZ);
     tref=zeros(1,Nz);
@@ -471 +450 @@
     namF='run/salt_ini.bin';
     fid=fopen(namF,'w','b'); fwrite(fid,var,prec);fclose(fid);
-
-    rhoAvr=rho_water-1.345;
-    mIce=-rhoAvr*draft;
-    namF='run/shelficeMass.Lin.bin';
-    fid=fopen(namF,'w','b'); fwrite(fid,mIce,prec);fclose(fid);
     % }}}
     % }}}