source: issm/oecreview/Archive/21724-22754/ISSM-21789-21790.diff@ 22755

Last change on this file since 22755 was 22755, checked in by Mathieu Morlighem, 7 years ago

CHG: added 21724-22754
File size: 9.4 KB

  • ../trunk-jpl/test/NightlyRun/test441.m

     
     1%Test Name: MISMIP3D
     2md=triangle(model(),'../Exp/Square.exp',100000.);
     3md=setmask(md,'../Exp/SquareShelf.exp','');
     4md=parameterize(md,'../Par/SquareSheetShelf.par');
     5md.initialization.vx(:)=1.;
     6md.initialization.vy(:)=1.;
     7md.geometry.thickness(:)=500-md.mesh.x/10000;
     8md.geometry.bed =-100-md.mesh.x/1000;
     9md.geometry.base=-md.geometry.thickness*md.materials.rho_ice/md.materials.rho_water;
     10md.mask.groundedice_levelset=md.geometry.thickness+md.materials.rho_water/md.materials.rho_ice*md.geometry.bed;
     11pos=find(md.mask.groundedice_levelset>=0);
     12md.geometry.base(pos)=md.geometry.bed(pos);
     13md.geometry.surface=md.geometry.base+md.geometry.thickness;
     14md=setflowequation(md,'SSA','all');
     15
     16%Boundary conditions:
     17md.mask.ice_levelset=-ones(md.mesh.numberofvertices,1);
     18md.mask.ice_levelset(find(md.mesh.x==max(md.mesh.x)))=0;
     19md.stressbalance.spcvx(:)=NaN;
     20md.stressbalance.spcvy(:)=NaN;
     21md.stressbalance.spcvz(:)=NaN;
     22pos=find((md.mesh.y<1000000.1 & md.mesh.y>999999.9) | (md.mesh.y<0.1 & md.mesh.y>-0.1));
     23md.stressbalance.spcvy(pos)=0;
     24pos2=find(md.mesh.x<0.1 & md.mesh.x>-0.1);
     25md.stressbalance.spcvx(pos2)=0;
     26md.stressbalance.spcvy(pos2)=0;
     27
     28md.materials.rheology_B=1/((10^-25)^(1/3))*ones(md.mesh.numberofvertices,1);
     29md.materials.rheology_law='None';
     30md.friction.coefficient(:)=sqrt(10^7)*ones(md.mesh.numberofvertices,1);
     31md.friction.p=3*ones(md.mesh.numberofelements,1);
     32md.smb.mass_balance(:)=1;
     33md.basalforcings.groundedice_melting_rate(:)=0;
     34md.basalforcings.floatingice_melting_rate(:)=30;
     35md.transient.isthermal=0;
     36md.transient.isstressbalance=1;
     37md.transient.isgroundingline=1;
     38md.transient.ismasstransport=1;
     39md.transient.issmb=1;
     40md.transient.requested_outputs={'default','BasalforcingsFloatingiceMeltingRate'};
     41md.groundingline.migration='SubelementMigration2';
     42md.timestepping.final_time=30;
     43md.timestepping.time_step=10;
     44
     45md.cluster=generic('name',oshostname(),'np',3);
     46md=solve(md,'Transient');
     47
     48%Fields and tolerances to track changes
     49field_names     ={'Bed1','Surface1','Thickness1','Floatingice1','Vx1','Vy1','Pressure1','FloatingiceMeltingrate1',...
     50        'Bed2','Surface2','Thickness2','Floatingice2','Vx2','Vy2','Pressure2','FloatingiceMeltingrate2',...
     51        'Bed3','Surface3','Thickness3','Floatingice3','Vx3','Vy3','Pressure3','FloatingiceMeltingrate3'};
     52field_tolerances={2e-11,5e-12,2e-11,1e-11,5e-10,1e-08,1e-13,1e-13,...
     53        3e-11,3e-11,9e-10,7e-11,1e-09,5e-08,1e-10,1e-13,...
     54        1e-10,3e-11,1e-10,7e-11,1e-09,5e-08,1e-10,1e-13};
     55field_values={...
     56        (md.results.TransientSolution(1).Base),...
     57        (md.results.TransientSolution(1).Surface),...
     58        (md.results.TransientSolution(1).Thickness),...
     59        (md.results.TransientSolution(1).MaskGroundediceLevelset),...
     60        (md.results.TransientSolution(1).Vx),...
     61        (md.results.TransientSolution(1).Vy),...
     62        (md.results.TransientSolution(1).Pressure),...
     63        (md.results.TransientSolution(1).BasalforcingsFloatingiceMeltingRate),...
     64        (md.results.TransientSolution(2).Base),...
     65        (md.results.TransientSolution(2).Surface),...
     66        (md.results.TransientSolution(2).Thickness),...
     67        (md.results.TransientSolution(2).MaskGroundediceLevelset),...
     68        (md.results.TransientSolution(2).Vx),...
     69        (md.results.TransientSolution(2).Vy),...
     70        (md.results.TransientSolution(2).Pressure),...
     71        (md.results.TransientSolution(2).BasalforcingsFloatingiceMeltingRate),...
     72        (md.results.TransientSolution(3).Base),...
     73        (md.results.TransientSolution(3).Surface),...
     74        (md.results.TransientSolution(3).Thickness),...
     75        (md.results.TransientSolution(3).MaskGroundediceLevelset),...
     76        (md.results.TransientSolution(3).Vx),...
     77        (md.results.TransientSolution(3).Vy),...
     78        (md.results.TransientSolution(3).Pressure),...
     79        (md.results.TransientSolution(3).BasalforcingsFloatingiceMeltingRate),...
     80        };

  • ../trunk-jpl/test/NightlyRun/test442.m

     
     1%Test Name: MISMIP3DHO
     2md=triangle(model(),'../Exp/Square.exp',100000.);
     3md=setmask(md,'../Exp/SquareShelf.exp','');
     4md=parameterize(md,'../Par/SquareSheetShelf.par');
     5md.initialization.vx(:)=1.;
     6md.initialization.vy(:)=1.;
     7md.geometry.thickness(:)=500-md.mesh.x/10000;
     8md.geometry.bed =-100-md.mesh.x/1000;
     9md.geometry.base=-md.geometry.thickness*md.materials.rho_ice/md.materials.rho_water;
     10md.mask.groundedice_levelset=md.geometry.thickness+md.materials.rho_water/md.materials.rho_ice*md.geometry.bed;
     11pos=find(md.mask.groundedice_levelset>=0);
     12md.geometry.base(pos)=md.geometry.bed(pos);
     13md.geometry.surface=md.geometry.base+md.geometry.thickness;
     14md=extrude(md,4,1.);
     15md=setflowequation(md,'HO','all');
     16
     17%Boundary conditions:
     18md.mask.ice_levelset=-ones(md.mesh.numberofvertices,1);
     19md.mask.ice_levelset(find(md.mesh.x==max(md.mesh.x)))=0;
     20md.stressbalance.spcvx(:)=NaN;
     21md.stressbalance.spcvy(:)=NaN;
     22md.stressbalance.spcvz(:)=NaN;
     23pos=find((md.mesh.y<1000000.1 & md.mesh.y>999999.9) | (md.mesh.y<0.1 & md.mesh.y>-0.1));
     24md.stressbalance.spcvy(pos)=0;
     25pos2=find(md.mesh.x<0.1 & md.mesh.x>-0.1);
     26md.stressbalance.spcvx(pos2)=0;
     27md.stressbalance.spcvy(pos2)=0;
     28
     29md.materials.rheology_B=1/((10^-25)^(1/3))*ones(md.mesh.numberofvertices,1);
     30md.materials.rheology_law='None';
     31md.friction.coefficient(:)=sqrt(10^7)*ones(md.mesh.numberofvertices,1);
     32md.friction.p=3*ones(md.mesh.numberofelements,1);
     33md.smb.mass_balance(:)=1;
     34md.basalforcings.groundedice_melting_rate(:)=0;
     35md.basalforcings.floatingice_melting_rate(:)=30;
     36md.transient.isthermal=0;
     37md.transient.isstressbalance=1;
     38md.transient.isgroundingline=1;
     39md.transient.ismasstransport=1;
     40md.transient.issmb=1;
     41md.transient.requested_outputs={'default','BasalforcingsFloatingiceMeltingRate'};
     42md.groundingline.migration='SubelementMigration2';
     43md.timestepping.final_time=30;
     44md.timestepping.time_step=10;
     45
     46md.cluster=generic('name',oshostname(),'np',3);
     47md=solve(md,'Transient');
     48
     49%Fields and tolerances to track changes
     50field_names     ={'Bed1','Surface1','Thickness1','Floatingice1','Vx1','Vy1','Vz1','Pressure1','FloatingiceMeltingrate1',...
     51        'Bed2','Surface2','Thickness2','Floatingice2','Vx2','Vy2','Vz2','Pressure2','FloatingiceMeltingrate2',...
     52        'Bed3','Surface3','Thickness3','Floatingice3','Vx3','Vy3','Vz3','Pressure3','FloatingiceMeltingrate3',};
     53field_tolerances={2e-11,5e-12,2e-11,1e-11,5e-10,1e-08,6e-10,1e-13,1e-13,...
     54        3e-11,3e-11,9e-10,7e-11,1e-09,5e-08,1e-09,1e-10,1e-13,...
     55        1e-10,3e-11,1e-10,7e-11,1e-09,5e-08,1e-09,1e-10,1e-13};
     56field_values={...
     57        (md.results.TransientSolution(1).Base),...
     58        (md.results.TransientSolution(1).Surface),...
     59        (md.results.TransientSolution(1).Thickness),...
     60        (md.results.TransientSolution(1).MaskGroundediceLevelset),...
     61        (md.results.TransientSolution(1).Vx),...
     62        (md.results.TransientSolution(1).Vy),...
     63        (md.results.TransientSolution(1).Vz),...
     64        (md.results.TransientSolution(1).Pressure),...
     65        (md.results.TransientSolution(1).BasalforcingsFloatingiceMeltingRate),...
     66        (md.results.TransientSolution(2).Base),...
     67        (md.results.TransientSolution(2).Surface),...
     68        (md.results.TransientSolution(2).Thickness),...
     69        (md.results.TransientSolution(2).MaskGroundediceLevelset),...
     70        (md.results.TransientSolution(2).Vx),...
     71        (md.results.TransientSolution(2).Vy),...
     72        (md.results.TransientSolution(2).Vz),...
     73        (md.results.TransientSolution(2).Pressure),...
     74        (md.results.TransientSolution(2).BasalforcingsFloatingiceMeltingRate),...
     75        (md.results.TransientSolution(3).Base),...
     76        (md.results.TransientSolution(3).Surface),...
     77        (md.results.TransientSolution(3).Thickness),...
     78        (md.results.TransientSolution(3).MaskGroundediceLevelset),...
     79        (md.results.TransientSolution(3).Vx),...
     80        (md.results.TransientSolution(3).Vy),...
     81        (md.results.TransientSolution(3).Vz),...
     82        (md.results.TransientSolution(3).Pressure),...
     83        (md.results.TransientSolution(3).BasalforcingsFloatingiceMeltingRate),...
     84        };

  • ../trunk-jpl/test/Archives/Archive442.arch 

    Cannot display: file marked as a binary type.
svn:mime-type = application/octet-stream

  • ../trunk-jpl/test/Archives/Archive441.arch 

    Property changes on: ../trunk-jpl/test/Archives/Archive442.arch
___________________________________________________________________
Added: svn:mime-type
## -0,0 +1 ##
+application/octet-stream
\ No newline at end of property
Cannot display: file marked as a binary type.
svn:mime-type = application/octet-stream
Note: See TracBrowser for help on using the repository browser.