Changeset 26744 for issm/trunk/test/NightlyRun/test2001.m

Timestamp:

12/22/21 10:39:44 (3 years ago)

Author:

Mathieu Morlighem

Message:

merged trunk-jpl and trunk for revision 26742

Location:

issm/trunk

Files:

• . (modified) (1 prop)
• test (modified) (1 prop)
• test/NightlyRun/test2001.m (modified) (1 diff)

issm/trunk/test/NightlyRun/test2001.m

@@ -1 +1 @@
 %Test Name: SquareSheetConstrainedGia2d
-%GIA test, inspired on test101. Running default GIA Ivins class.
+%GIA test, based off of test101. Running default GIA Ivins class.
 md=triangle(model(),'../Exp/Square.exp',100000.);
 md=setmask(md,'','');
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
 
-%GIA: 
-md.gia=giaivins();
-md.gia.lithosphere_thickness=100.*ones(md.mesh.numberofvertices,1); % in km
-md.gia.mantle_viscosity=1.0*10^21*ones(md.mesh.numberofvertices,1); % in Pa.s
-md.materials.lithosphere_shear_modulus=6.7*10^10;                   % in Pa
-md.materials.lithosphere_density=3.32;                              % in g/cm^-3
-md.materials.mantle_shear_modulus=1.45*10^11;                       % in Pa
-md.materials.mantle_density=3.34;                                   % in g/cm^-3
+%GIA Ivins, 2 layer model.
+md.solidearth.settings.grdmodel=2;
+md.solidearth.settings.isgrd=1;
 
-%% indicate what you want to compute 
-md.gia.cross_section_shape=1;    % for square-edged x-section 
+md.materials=materials('litho','ice');
+md.materials.numlayers=2;
+md.materials.radius =  [10 6271e3 6371e3];
+md.materials.density=  [3.34e3 3.32e3];
+md.materials.lame_mu=  [1.45e11         6.7e10];
+md.materials.viscosity=[1e21            0];
+md.initialization.sealevel=zeros(md.mesh.numberofvertices,1);
+md.solidearth.settings.cross_section_shape=1;    % for square-edged x-section 
+md.solidearth.settings.grdocean=0;  %do not compute sea level, only deformation
+md.solidearth.settings.sealevelloading=0;  %do not compute sea level, only deformation
+md.solidearth.requested_outputs={'Sealevel','BedGRD'};
 
-%% define loading history 
-md.timestepping.start_time=2400000; %2,400 kyr :: EVALUATION TIME
+%Loading history 
+md.timestepping.start_time=-2400000; %4,800 kyr :: EVALUATION TIME
+md.timestepping.time_step= 2400000; %2,400 kyr :: EVALUATION TIME
 % to get rid of default final_time: make sure final_time>start_time
-md.timestepping.final_time=2500000; %2,500 kyr
-md.geometry.thickness=[...
-        [md.geometry.thickness; 0.0],...
-        [md.geometry.thickness; md.timestepping.start_time],...
+md.timestepping.final_time=2400000; %2,500 kyr
+md.masstransport.spcthickness=[...
+        [md.geometry.thickness; 0],...
+        [md.geometry.thickness; 2400000]...
         ];
 
-%% solve for GIA deflection 
+%geometry at 0 initially: 
+md.geometry.thickness=zeros(md.mesh.numberofvertices,1);
+md.geometry.surface=zeros(md.mesh.numberofvertices,1);
+md.geometry.base=zeros(md.mesh.numberofvertices,1);
+
+%Physics: 
+md.transient.issmb=0; 
+md.transient.isstressbalance=0;
+md.transient.isthermal=0;
+md.transient.ismasstransport=1;
+md.transient.isslc=1;
+
+% Solve for GIA deflection 
 md.cluster=generic('name',oshostname(),'np',3);
-md.verbose=verbose('1111111');
-md=solve(md,'Gia');
+md.verbose.solver=0;
+md=solve(md,'Transient');
 
 %Fields and tolerances to track changes
-field_names     ={'UGia','UGiaRate'};
-field_tolerances={1e-13,1e-13};
-field_values={...
-        (md.results.GiaSolution.UGia),...
-        (md.results.GiaSolution.UGiaRate),...
-        };
+field_names     ={'UGrd'};
+field_tolerances={1e-13};
+field_values={md.results.TransientSolution(2).BedGRD};