Changeset 13663
- Timestamp:
- 10/15/12 11:31:48 (12 years ago)
- Location:
- issm/trunk-jpl/test/NightlyRun
- Files:
-
- 115 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/NightlyRun/ad.m
r13508 r13663 23 23 md.autodiff.driver='fos_reverse'; 24 24 25 md=solve(md,TransientSolutionEnum );25 md=solve(md,TransientSolutionEnum()); 26 26 27 27 %plot -
issm/trunk-jpl/test/NightlyRun/ad_maxalloc.m
r13535 r13663 17 17 md.autodiff.driver='fov_forward'; 18 18 19 md=solve(md,TransientSolutionEnum );19 md=solve(md,TransientSolutionEnum()); -
issm/trunk-jpl/test/NightlyRun/test1101.m
r11186 r13663 12 12 nx=20; %numberof nodes in x direction 13 13 ny=20; 14 md=model ;14 md=model(); 15 15 md=squaremesh(md,L,L,nx,ny); 16 16 md=setmask(md,'',''); %ice sheet test … … 40 40 %Compute the diagnostic 41 41 md.cluster=generic('name',oshostname(),'np',8); 42 md=solve(md,DiagnosticSolutionEnum );42 md=solve(md,DiagnosticSolutionEnum()); 43 43 44 44 %Plot the results and save them -
issm/trunk-jpl/test/NightlyRun/test1102.m
r10976 r13663 12 12 nx=20; %numberof nodes in x direction 13 13 ny=20; 14 md=model ;14 md=model(); 15 15 md=squaremesh(md,L,L,nx,ny); 16 16 md=setmask(md,'',''); %ice sheet test … … 49 49 md.diagnostic.restol=1; 50 50 md.cluster=generic('name',oshostname(),'np',8); 51 md=solve(md,DiagnosticSolutionEnum );51 md=solve(md,DiagnosticSolutionEnum()); 52 52 53 53 %Plot the results and save them -
issm/trunk-jpl/test/NightlyRun/test1103.m
r11193 r13663 10 10 nx=20; %numberof nodes in x direction 11 11 ny=20; 12 md=model ;12 md=model(); 13 13 md=squaremesh(md,L,L,nx,ny); 14 14 md=setmask(md,'',''); %ice sheet test … … 37 37 %Compute the diagnostic 38 38 md.cluster=generic('name',oshostname(),'np',8); 39 md=solve(md,DiagnosticSolutionEnum );39 md=solve(md,DiagnosticSolutionEnum()); 40 40 41 41 %Plot the results and save them -
issm/trunk-jpl/test/NightlyRun/test1104.m
r11550 r13663 9 9 nx=20; %numberof nodes in x direction 10 10 ny=20; 11 md=model ;11 md=model(); 12 12 md=squaremesh(md,L,L,nx,ny); 13 13 md=setmask(md,'',''); %ice sheet test … … 37 37 md.diagnostic.abstol=NaN; 38 38 md.cluster=generic('name',oshostname(),'np',8); 39 md=solve(md,DiagnosticSolutionEnum );39 md=solve(md,DiagnosticSolutionEnum()); 40 40 pos=find(md.mesh.x==0 | md.mesh.y==0 | md.mesh.x==max(md.mesh.x) | md.mesh.y==max(md.mesh.y)); 41 41 md.diagnostic.spcvx(pos)=md.results.DiagnosticSolution.Vx(pos); … … 43 43 md.diagnostic.vertex_pairing=[]; 44 44 md=setflowequation(md,'stokes','all'); 45 md=solve(md,DiagnosticSolutionEnum );45 md=solve(md,DiagnosticSolutionEnum()); 46 46 47 47 %Plot the results and save them -
issm/trunk-jpl/test/NightlyRun/test1105.m
r11550 r13663 12 12 nx=30; %number of nodes in x direction 13 13 ny=30; 14 md=model ;14 md=model(); 15 15 md=squaremesh(md,L,L,nx,ny); 16 16 md=setmask(md,'',''); %ice sheet test … … 63 63 %Compute the diagnostic 64 64 md.cluster=generic('name',oshostname(),'np',8); 65 md=solve(md,DiagnosticSolutionEnum );65 md=solve(md,DiagnosticSolutionEnum()); 66 66 67 67 %Plot the results and save them -
issm/trunk-jpl/test/NightlyRun/test1106.m
r11323 r13663 41 41 %Compute the diagnostic 42 42 md.cluster=generic('name',oshostname(),'np',8); 43 md=solve(md,DiagnosticSolutionEnum );43 md=solve(md,DiagnosticSolutionEnum()); 44 44 45 45 %Plot the results and save them -
issm/trunk-jpl/test/NightlyRun/test1107.m
r12811 r13663 12 12 nx=30; %numberof nodes in x direction 13 13 ny=30; 14 md=model ;14 md=model(); 15 15 md=squaremesh(md,L,L,nx,ny); 16 16 md=setmask(md,'',''); %ice sheet test … … 58 58 %Compute the diagnostic 59 59 md.cluster=generic('name',oshostname(),'np',8); 60 md=solve(md,DiagnosticSolutionEnum );60 md=solve(md,DiagnosticSolutionEnum()); 61 61 62 62 %Plot the results and save them -
issm/trunk-jpl/test/NightlyRun/test1108.m
r12811 r13663 9 9 nx=30; %numberof nodes in x direction 10 10 ny=30; 11 md=model ;11 md=model(); 12 12 md=squaremesh(md,L,L,nx,ny); 13 13 md=setmask(md,'',''); %ice sheet test … … 39 39 md.cluster=generic('name',oshostname(),'np',8); 40 40 md.verbose=verbose('convergence',true); 41 md=solve(md,DiagnosticSolutionEnum );41 md=solve(md,DiagnosticSolutionEnum()); 42 42 md.diagnostic.reltol=NaN; 43 43 md.diagnostic.abstol=NaN; … … 51 51 md.diagnostic.spcvy(pos)=md.results.DiagnosticSolution.Vy(pos); 52 52 md=setflowequation(md,'stokes','all'); 53 md=solve(md,DiagnosticSolutionEnum );53 md=solve(md,DiagnosticSolutionEnum()); 54 54 55 55 %Plot the results and save them -
issm/trunk-jpl/test/NightlyRun/test1109.m
r11550 r13663 13 13 nx=3; %number of nodes in x direction 14 14 ny=51; 15 md=model ;15 md=model(); 16 16 md=squaremesh(md,Lx,Ly,nx,ny); 17 17 md=setmask(md,'',''); %ice sheet test … … 49 49 %Compute the diagnostic 50 50 md.cluster=generic('name',oshostname(),'np',8); 51 md=solve(md,DiagnosticSolutionEnum );51 md=solve(md,DiagnosticSolutionEnum()); 52 52 53 53 vx=(md.results.DiagnosticSolution.Vx); -
issm/trunk-jpl/test/NightlyRun/test1110.m
r11550 r13663 7 7 nx=30; %numberof nodes in x direction 8 8 ny=30; 9 md=model ;9 md=model(); 10 10 md=squaremesh(md,L,L,nx,ny); 11 11 %md=triangle(md,'../Exp/SquareISMIP.exp',5500); … … 57 57 md.cluster=generic('name',oshostname(),'np',8); 58 58 md.verbose=verbose('convergence',true,'solution',true); 59 md=solve(md,TransientSolutionEnum );59 md=solve(md,TransientSolutionEnum()); 60 60 61 61 %save the results -
issm/trunk-jpl/test/NightlyRun/test1201.m
r11194 r13663 38 38 md.settings.output_frequency=500; %keep only last step 39 39 md.verbose=verbose(); 40 md=solve(md,TransientSolutionEnum );40 md=solve(md,TransientSolutionEnum()); 41 41 results{stabilization}=(md.results.TransientSolution(end).Thickness); 42 42 end -
issm/trunk-jpl/test/NightlyRun/test1202.m
r11014 r13663 3 3 4 4 %tests 3 and 4: using Glen's flow law 5 md=model ;5 md=model(); 6 6 md=triangle(md,'../Exp/SquareEISMINT.exp',3550); 7 7 md=setmask(md,'all',''); … … 11 11 %Compute solution for MacAyeal's model 12 12 md.cluster=generic('name',oshostname(),'np',8); 13 md=solve(md,DiagnosticSolutionEnum );13 md=solve(md,DiagnosticSolutionEnum()); 14 14 15 15 %plot results -
issm/trunk-jpl/test/NightlyRun/test1203.m
r11014 r13663 3 3 4 4 %test 5 and 6 : 5 md=model ;5 md=model(); 6 6 md=triangle(md,'../Exp/SquareEISMINT.exp',5100); %test3 7 7 md=setmask(md,'all',''); … … 15 15 %Compute solution for MacAyeal's model 16 16 md.cluster=generic('name',oshostname(),'np',8); 17 md=solve(md,DiagnosticSolutionEnum );17 md=solve(md,DiagnosticSolutionEnum()); 18 18 19 19 vx=(md.results.DiagnosticSolution.Vx); -
issm/trunk-jpl/test/NightlyRun/test1204.m
r11014 r13663 3 3 4 4 %tests 3 and 4: using Glen's flow law 5 md=model ;5 md=model(); 6 6 md=triangle(md,'../Exp/SquareEISMINT.exp',3550); 7 7 md=setmask(md,'all',''); … … 15 15 %Compute solution for MacAyeal's model 16 16 md.cluster=generic('name',oshostname(),'np',8); 17 md=solve(md,DiagnosticSolutionEnum );17 md=solve(md,DiagnosticSolutionEnum()); 18 18 19 19 %plot results … … 24 24 md.timestepping.final_time=5000; 25 25 md.prognostic.stabilization=1; 26 md=solve(md,TransientSolutionEnum );26 md=solve(md,TransientSolutionEnum()); 27 27 28 28 plotmodel(md,'data',(md.results.TransientSolution(end).Vx)) -
issm/trunk-jpl/test/NightlyRun/test1205.m
r10976 r13663 6 6 7 7 %To begin with the numerical model 8 md=model ;8 md=model(); 9 9 md=roundmesh(md,750000,resolution); 10 10 md=setmask(md,'',''); %We can not test iceshelves nor ice rises with this analytical solution … … 29 29 %Now we can solve the problem 30 30 md.cluster=generic('name',oshostname(),'np',8); 31 md=solve(md,DiagnosticSolutionEnum );31 md=solve(md,DiagnosticSolutionEnum()); 32 32 33 33 %Calculate the depth averaged velocity field (2d): -
issm/trunk-jpl/test/NightlyRun/test1206.m
r11193 r13663 6 6 7 7 %To begin with the numerical model 8 md=model ;8 md=model(); 9 9 md=roundmesh(md,750000,resolution); 10 10 md=setmask(md,'',''); %We can not test iceshelves nor ice rises with this analytical solution … … 29 29 %Now we can solve the problem 30 30 md.cluster=generic('name',oshostname(),'np',8); 31 md=solve(md,DiagnosticSolutionEnum );31 md=solve(md,DiagnosticSolutionEnum()); 32 32 33 33 %Calculate the depth averaged velocity field (2d): -
issm/trunk-jpl/test/NightlyRun/test1207.m
r10976 r13663 6 6 7 7 %To begin with the numerical model 8 md=model ;8 md=model(); 9 9 md=roundmesh(md,750000,resolution); 10 10 md=setmask(md,'',''); %We can not test iceshelves nor ice rises with this analytical solution … … 29 29 %Now we can solve the problem 30 30 md.cluster=generic('name',oshostname(),'np',8); 31 md=solve(md,DiagnosticSolutionEnum );31 md=solve(md,DiagnosticSolutionEnum()); 32 32 33 33 %Calculate the depth averaged velocity field (2d): -
issm/trunk-jpl/test/NightlyRun/test1208.m
r11204 r13663 27 27 %Now we can solve the problem 28 28 md.cluster=generic('name',oshostname(),'np',8); 29 md=solve(md,TransientSolutionEnum );29 md=solve(md,TransientSolutionEnum()); 30 30 31 31 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test1301.m
r11193 r13663 5 5 printingflag=false; 6 6 7 md=model ;7 md=model(); 8 8 md=triangle(md,'../Exp/Square.exp',100000); 9 9 md=setmask(md,'',''); … … 26 26 %modeled results 27 27 md.cluster=generic('name',oshostname(),'np',2); 28 md=solve(md,ThermalSolutionEnum );28 md=solve(md,ThermalSolutionEnum()); 29 29 30 30 %plot results -
issm/trunk-jpl/test/NightlyRun/test1302.m
r11172 r13663 4 4 printingflag=false; 5 5 6 md=model ;6 md=model(); 7 7 md=triangle(md,'../Exp/Square.exp',100000); 8 8 md=setmask(md,'',''); … … 28 28 %modeled results 29 29 md.cluster=generic('name',oshostname(),'np',2); 30 md=solve(md,ThermalSolutionEnum );30 md=solve(md,ThermalSolutionEnum()); 31 31 32 32 %plot results -
issm/trunk-jpl/test/NightlyRun/test1303.m
r11024 r13663 5 5 printingflag=false; 6 6 7 md=model ;7 md=model(); 8 8 md=triangle(md,'../Exp/Square.exp',100000); 9 9 md=setmask(md,'all',''); … … 22 22 %modeled results 23 23 md.cluster=generic('name',oshostname(),'np',2); 24 md=solve(md,ThermalSolutionEnum );24 md=solve(md,ThermalSolutionEnum()); 25 25 26 26 %plot results -
issm/trunk-jpl/test/NightlyRun/test1304.m
r11024 r13663 5 5 printingflag=false; 6 6 7 md=model ;7 md=model(); 8 8 md=triangle(md,'../Exp/Square.exp',100000); 9 9 md=setmask(md,'',''); … … 23 23 %modeled results 24 24 md.cluster=generic('name',oshostname(),'np',2); 25 md=solve(md,ThermalSolutionEnum );25 md=solve(md,ThermalSolutionEnum()); 26 26 27 27 %plot results -
issm/trunk-jpl/test/NightlyRun/test1401.m
r9850 r13663 7 7 nx=70; %numberof nodes in x direction 8 8 ny=70; 9 md=model ;9 md=model(); 10 10 11 11 %mesh adaptation loop YAMS -
issm/trunk-jpl/test/NightlyRun/test1402.m
r9850 r13663 6 6 nx=30; %numberof nodes in x direction 7 7 ny=30; 8 md=model ;8 md=model(); 9 9 10 10 %mesh adaptation loop YAMS -
issm/trunk-jpl/test/NightlyRun/test1501.m
r11014 r13663 17 17 md.basalforcings.melting_rate= smb; 18 18 19 md=solve(md,PrognosticSolutionEnum );19 md=solve(md,PrognosticSolutionEnum()); 20 20 21 21 for i=1:10 22 md=solve(md,PrognosticSolutionEnum );22 md=solve(md,PrognosticSolutionEnum()); 23 23 md.surfaceforcings.mass_balance= md.surfaceforcings.mass_balance - ((md.results.PrognosticSolution.Thickness)-md.geometry.thickness); 24 24 end … … 33 33 md.surfaceforcings.mass_balance(end+1,:)=[1:2000]; 34 34 35 md=solve(md,TransientSolutionEnum );35 md=solve(md,TransientSolutionEnum()); 36 36 37 37 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test1502.m
r11254 r13663 18 18 md.basalforcings.melting_rate= smb; 19 19 20 md=solve(md,PrognosticSolutionEnum );20 md=solve(md,PrognosticSolutionEnum()); 21 21 22 22 for i=1:10 23 md=solve(md,PrognosticSolutionEnum );23 md=solve(md,PrognosticSolutionEnum()); 24 24 md.surfaceforcings.mass_balance= md.surfaceforcings.mass_balance - ((md.results.PrognosticSolution.Thickness)-md.geometry.thickness); 25 25 end … … 34 34 md.surfaceforcings.mass_balance(end+1,:)=[1:2000]; 35 35 36 md=solve(md,TransientSolutionEnum );36 md=solve(md,TransientSolutionEnum()); 37 37 38 38 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test1601.m
r11024 r13663 9 9 10 10 md.cluster=generic('name',oshostname,'np',2); 11 md=solve(md,DiagnosticSolutionEnum );11 md=solve(md,DiagnosticSolutionEnum()); 12 12 vel0=md.results.DiagnosticSolution.Vel; 13 13 … … 20 20 md.diagnostic.referential(:,1:3)=repmat([cos(theta),sin(theta),0],md.mesh.numberofvertices,1); 21 21 md.diagnostic.referential(:,4:6)=repmat([0,0,1],md.mesh.numberofvertices,1); 22 md=solve(md,DiagnosticSolutionEnum );22 md=solve(md,DiagnosticSolutionEnum()); 23 23 vel1=md.results.DiagnosticSolution.Vel; 24 24 … … 31 31 md.diagnostic.referential(pos,1:3)=repmat([cos(theta),sin(theta),0],size(pos,1),1); 32 32 md.diagnostic.referential(pos,4:6)=repmat([0,0,1],size(pos,1),1); 33 md=solve(md,DiagnosticSolutionEnum );33 md=solve(md,DiagnosticSolutionEnum()); 34 34 vel2=md.results.DiagnosticSolution.Vel; 35 35 -
issm/trunk-jpl/test/NightlyRun/test1602.m
r11024 r13663 10 10 11 11 md.cluster=generic('name',oshostname,'np',3); 12 md=solve(md,DiagnosticSolutionEnum );12 md=solve(md,DiagnosticSolutionEnum()); 13 13 vel0=md.results.DiagnosticSolution.Vel; 14 14 … … 21 21 md.diagnostic.referential(:,1:3)=repmat([cos(theta),sin(theta),0],md.mesh.numberofvertices,1); 22 22 md.diagnostic.referential(:,4:6)=repmat([0,0,1],md.mesh.numberofvertices,1); 23 md=solve(md,DiagnosticSolutionEnum );23 md=solve(md,DiagnosticSolutionEnum()); 24 24 vel1=md.results.DiagnosticSolution.Vel; 25 25 -
issm/trunk-jpl/test/NightlyRun/test3001.m
r13433 r13663 4 4 md=setflowequation(md,'macayeal','all'); 5 5 md.cluster=generic('name',oshostname(),'np',3); 6 md.diagnostic.requested_outputs=StressTensorEnum ;6 md.diagnostic.requested_outputs=StressTensorEnum(); 7 7 md.autodiff.isautodiff=true; 8 8 md.verbose=verbose('autodiff',true); 9 md=solve(md,DiagnosticSolutionEnum );9 md=solve(md,DiagnosticSolutionEnum()); 10 10 11 11 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test3002.m
r13433 r13663 7 7 md.autodiff.isautodiff=true; 8 8 md.verbose=verbose('autodiff',true); 9 md=solve(md,DiagnosticSolutionEnum );9 md=solve(md,DiagnosticSolutionEnum()); 10 10 11 11 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test3003.m
r13433 r13663 5 5 md=setflowequation(md,'pattyn','all'); 6 6 md.cluster=generic('name',oshostname(),'np',3); 7 md.diagnostic.requested_outputs=StressTensorEnum ;7 md.diagnostic.requested_outputs=StressTensorEnum(); 8 8 md.autodiff.isautodiff=true; 9 9 md.verbose=verbose('autodiff',true); 10 md=solve(md,DiagnosticSolutionEnum );10 md=solve(md,DiagnosticSolutionEnum()); 11 11 12 12 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test3004.m
r13433 r13663 7 7 md.autodiff.isautodiff=true; 8 8 md.verbose=verbose('autodiff',true); 9 md=solve(md,DiagnosticSolutionEnum );9 md=solve(md,DiagnosticSolutionEnum()); 10 10 11 11 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test3005.m
r13433 r13663 6 6 md.autodiff.isautodiff=true; 7 7 md.verbose=verbose('autodiff',true); 8 md=solve(md,PrognosticSolutionEnum );8 md=solve(md,PrognosticSolutionEnum()); 9 9 10 10 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test3006.m
r13433 r13663 9 9 md.autodiff.isautodiff=true; 10 10 md.verbose=verbose('autodiff',true); 11 md=solve(md,PrognosticSolutionEnum );11 md=solve(md,PrognosticSolutionEnum()); 12 12 13 13 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test3007.m
r13433 r13663 7 7 md.autodiff.isautodiff=true; 8 8 md.verbose=verbose('autodiff',true); 9 md=solve(md,PrognosticSolutionEnum );9 md=solve(md,PrognosticSolutionEnum()); 10 10 11 11 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test3008.m
r13433 r13663 8 8 md.autodiff.isautodiff=true; 9 9 md.verbose=verbose('autodiff',true); 10 md=solve(md,ThermalSolutionEnum );10 md=solve(md,ThermalSolutionEnum()); 11 11 12 12 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test3009.m
r13433 r13663 11 11 md.autodiff.isautodiff=true; 12 12 md.verbose=verbose('autodiff',true); 13 md=solve(md,TransientSolutionEnum );13 md=solve(md,TransientSolutionEnum()); 14 14 15 15 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test3010.m
r13433 r13663 8 8 md.autodiff.isautodiff=true; 9 9 md.verbose=verbose('autodiff',true); 10 md=solve(md,TransientSolutionEnum );10 md=solve(md,TransientSolutionEnum()); 11 11 12 12 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test3015.m
r13552 r13663 8 8 md=setflowequation(md,'macayeal','all'); 9 9 md.cluster=generic('name',oshostname(),'np',3); 10 md.prognostic.requested_outputs=IceVolumeEnum ;10 md.prognostic.requested_outputs=IceVolumeEnum(); 11 11 md.verbose=verbose('autodiff',true); 12 12 … … 41 41 md=SetIceShelfBC(md); 42 42 43 md=solve(md,PrognosticSolutionEnum );43 md=solve(md,PrognosticSolutionEnum()); 44 44 V0=md.results.PrognosticSolution.IceVolume; 45 45 … … 52 52 md=SetIceShelfBC(md); 53 53 54 md=solve(md,PrognosticSolutionEnum );54 md=solve(md,PrognosticSolutionEnum()); 55 55 V2=md.results.PrognosticSolution.IceVolume; 56 56 … … 66 66 md=SetIceShelfBC(md); 67 67 68 md=solve(md,PrognosticSolutionEnum );68 md=solve(md,PrognosticSolutionEnum()); 69 69 %retrieve directly 70 70 dVdh_ad=md.results.PrognosticSolution.AutodiffJacobian; -
issm/trunk-jpl/test/NightlyRun/test3019.m
r13516 r13663 14 14 md.autodiff.driver='fos_reverse'; 15 15 16 md=solve(md,TransientSolutionEnum );16 md=solve(md,TransientSolutionEnum()); 17 17 18 18 %recover jacobian: … … 25 25 md.autodiff.driver='fov_forward'; 26 26 27 md=solve(md,TransientSolutionEnum );27 md=solve(md,TransientSolutionEnum()); 28 28 29 29 %recover jacobian: -
issm/trunk-jpl/test/NightlyRun/test3020.m
r13516 r13663 8 8 md=setflowequation(md,'macayeal','all'); 9 9 md.cluster=generic('name',oshostname(),'np',3); 10 md.transient.requested_outputs=[IceVolumeEnum ;MaxVelEnum];10 md.transient.requested_outputs=[IceVolumeEnum();MaxVelEnum]; 11 11 md.verbose=verbose('autodiff',true); 12 12 md.diagnostic.restol=.000001; … … 43 43 md=SetIceShelfBC(md); 44 44 45 md=solve(md,TransientSolutionEnum );45 md=solve(md,TransientSolutionEnum()); 46 46 V0=md.results.TransientSolution(end).IceVolume; 47 47 MaxV0=md.results.TransientSolution(end).MaxVel; … … 55 55 md=SetIceShelfBC(md); 56 56 57 md=solve(md,TransientSolutionEnum );57 md=solve(md,TransientSolutionEnum()); 58 58 V2=md.results.TransientSolution(end).IceVolume; 59 59 MaxV2=md.results.TransientSolution(end).MaxVel; … … 72 72 md=SetIceShelfBC(md); 73 73 74 md=solve(md,TransientSolutionEnum );74 md=solve(md,TransientSolutionEnum()); 75 75 %retrieve directly 76 76 dVdh_ad=md.results.TransientSolution(1).AutodiffJacobian(1); -
issm/trunk-jpl/test/NightlyRun/test328.m
r13529 r13663 22 22 md.surfaceforcings.smbref=smbref'; 23 23 md.cluster=generic('name',oshostname(),'np',3); 24 md=solve(md,TransientSolutionEnum );24 md=solve(md,TransientSolutionEnum()); 25 25 26 26 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test329.m
r13529 r13663 23 23 md.transient.requested_outputs=TotalSmbEnum(); 24 24 md.cluster=generic('name',oshostname(),'np',3); 25 md=solve(md,TransientSolutionEnum );25 md=solve(md,TransientSolutionEnum()); 26 26 27 27 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test401.m
r12804 r13663 4 4 md=setflowequation(md,'macayeal','all'); 5 5 md.cluster=generic('name',oshostname(),'np',3); 6 md=solve(md,DiagnosticSolutionEnum );6 md=solve(md,DiagnosticSolutionEnum()); 7 7 8 8 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test402.m
r12804 r13663 5 5 md=setflowequation(md,'macayeal','all'); 6 6 md.cluster=generic('name',oshostname(),'np',3); 7 md=solve(md,DiagnosticSolutionEnum );7 md=solve(md,DiagnosticSolutionEnum()); 8 8 9 9 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test403.m
r12804 r13663 5 5 md=setflowequation(md,'pattyn','all'); 6 6 md.cluster=generic('name',oshostname(),'np',3); 7 md=solve(md,DiagnosticSolutionEnum );7 md=solve(md,DiagnosticSolutionEnum()); 8 8 9 9 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test404.m
r12804 r13663 5 5 md=setflowequation(md,'stokes','all'); 6 6 md.cluster=generic('name',oshostname(),'np',3); 7 md=solve(md,DiagnosticSolutionEnum );7 md=solve(md,DiagnosticSolutionEnum()); 8 8 9 9 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test405.m
r12804 r13663 5 5 md=setflowequation(md,'macayeal','../Exp/SquareHalfRight.exp','fill','pattyn','coupling','penalties'); 6 6 md.cluster=generic('name',oshostname(),'np',3); 7 md=solve(md,DiagnosticSolutionEnum );7 md=solve(md,DiagnosticSolutionEnum()); 8 8 9 9 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test406.m
r12804 r13663 6 6 md.cluster=generic('name',oshostname(),'np',3); 7 7 md.timestepping.time_step=0; 8 md=solve(md,ThermalSolutionEnum );8 md=solve(md,ThermalSolutionEnum()); 9 9 10 10 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test407.m
r12804 r13663 9 9 md.transient.isthermal=1; 10 10 md.transient.isgroundingline=0; 11 md=solve(md,TransientSolutionEnum );11 md=solve(md,TransientSolutionEnum()); 12 12 13 13 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test408.m
r12804 r13663 4 4 md=setflowequation(md,'macayeal','all'); 5 5 md.cluster=generic('name',oshostname(),'np',3); 6 md=solve(md,TransientSolutionEnum );6 md=solve(md,TransientSolutionEnum()); 7 7 8 8 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test409.m
r12804 r13663 5 5 md=setflowequation(md,'macayeal','../Exp/SquareHalfRight.exp','fill','pattyn','coupling','penalties'); 6 6 md.cluster=generic('name',oshostname(),'np',3); 7 md=solve(md,TransientSolutionEnum );7 md=solve(md,TransientSolutionEnum()); 8 8 9 9 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test410.m
r12804 r13663 6 6 md.cluster=generic('name',oshostname(),'np',3); 7 7 md.timestepping.time_step=0; 8 md=solve(md,SteadystateSolutionEnum );8 md=solve(md,SteadystateSolutionEnum()); 9 9 10 10 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test411.m
r12804 r13663 6 6 md.cluster=generic('name',oshostname(),'np',3); 7 7 md.timestepping.time_step=0; 8 md=solve(md,SteadystateSolutionEnum );8 md=solve(md,SteadystateSolutionEnum()); 9 9 10 10 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test412.m
r13021 r13663 33 33 34 34 %solve 35 md=solve(md,DiagnosticSolutionEnum ,'overwrite','y');35 md=solve(md,DiagnosticSolutionEnum(),'overwrite','y'); 36 36 37 37 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test413.m
r13021 r13663 34 34 35 35 %solve 36 md=solve(md,DiagnosticSolutionEnum ,'overwrite','y');36 md=solve(md,DiagnosticSolutionEnum(),'overwrite','y'); 37 37 38 38 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test414.m
r13021 r13663 48 48 49 49 %solve 50 md=solve(md,DiagnosticSolutionEnum ,'overwrite','y');50 md=solve(md,DiagnosticSolutionEnum(),'overwrite','y'); 51 51 md.qmu.results=md.results.dakota; 52 52 -
issm/trunk-jpl/test/NightlyRun/test415.m
r12804 r13663 20 20 21 21 md.cluster=generic('name',oshostname(),'np',3); 22 md=solve(md,SteadystateSolutionEnum );22 md=solve(md,SteadystateSolutionEnum()); 23 23 24 24 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test416.m
r12804 r13663 20 20 21 21 md.cluster=generic('name',oshostname(),'np',3); 22 md=solve(md,SteadystateSolutionEnum );22 md=solve(md,SteadystateSolutionEnum()); 23 23 24 24 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test417.m
r13021 r13663 50 50 md.diagnostic.reltol=10^-5; %tighten for qmu analyses 51 51 52 md=solve(md,DiagnosticSolutionEnum ,'overwrite','y');52 md=solve(md,DiagnosticSolutionEnum(),'overwrite','y'); 53 53 54 54 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test419.m
r12804 r13663 5 5 md=setflowequation(md,'macayeal','../Exp/SquareHalfRight.exp','fill','pattyn'); 6 6 md.cluster=generic('name',oshostname(),'np',3); 7 md=solve(md,DiagnosticSolutionEnum );7 md=solve(md,DiagnosticSolutionEnum()); 8 8 9 9 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test420.m
r13021 r13663 32 32 33 33 %solve 34 md=solve(md,DiagnosticSolutionEnum ,'overwrite','y');34 md=solve(md,DiagnosticSolutionEnum(),'overwrite','y'); 35 35 md.qmu.results=md.results.dakota; 36 36 -
issm/trunk-jpl/test/NightlyRun/test421.m
r13340 r13663 5 5 md=setflowequation(md,'stokes','../Exp/SquareHalfRight.exp','fill','pattyn'); 6 6 md.cluster=generic('name',oshostname(),'np',3); 7 md=solve(md,DiagnosticSolutionEnum );7 md=solve(md,DiagnosticSolutionEnum()); 8 8 9 9 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test422.m
r12804 r13663 6 6 md.cluster=generic('name',oshostname(),'np',3); 7 7 md.diagnostic.reltol=0.4; 8 md=solve(md,DiagnosticSolutionEnum );8 md=solve(md,DiagnosticSolutionEnum()); 9 9 10 10 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test423.m
r12804 r13663 25 25 %test different grounding line dynamics. 26 26 md.groundingline.migration='AgressiveMigration'; 27 md=solve(md,TransientSolutionEnum );27 md=solve(md,TransientSolutionEnum()); 28 28 element_on_iceshelf_agressive=(md.results.TransientSolution.MaskElementonfloatingice); 29 29 30 30 md.groundingline.migration='SoftMigration'; 31 md=solve(md,TransientSolutionEnum );31 md=solve(md,TransientSolutionEnum()); 32 32 element_on_iceshelf_soft=(md.results.TransientSolution.MaskElementonfloatingice); 33 33 -
issm/trunk-jpl/test/NightlyRun/test424.m
r12804 r13663 16 16 17 17 md.cluster=generic('name',oshostname(),'np',3); 18 md=solve(md,TransientSolutionEnum );18 md=solve(md,TransientSolutionEnum()); 19 19 20 20 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test425.m
r12804 r13663 16 16 17 17 md.cluster=generic('name',oshostname(),'np',3); 18 md=solve(md,TransientSolutionEnum );18 md=solve(md,TransientSolutionEnum()); 19 19 20 20 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test426.m
r12804 r13663 18 18 md.cluster=generic('name',oshostname(),'np',3); 19 19 20 md=solve(md,TransientSolutionEnum );20 md=solve(md,TransientSolutionEnum()); 21 21 22 22 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test427.m
r12804 r13663 17 17 md.groundingline.migration='SoftMigration'; 18 18 md.cluster=generic('name',oshostname(),'np',3); 19 md=solve(md,TransientSolutionEnum );19 md=solve(md,TransientSolutionEnum()); 20 20 21 21 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test428.m
r12804 r13663 6 6 md.diagnostic.restol=0.0001; 7 7 md.cluster=generic('name',oshostname(),'np',3); 8 md=solve(md,DiagnosticSolutionEnum );8 md=solve(md,DiagnosticSolutionEnum()); 9 9 10 10 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test429.m
r12804 r13663 7 7 md.diagnostic.restol=0.0001; 8 8 md.cluster=generic('name',oshostname(),'np',3); 9 md=solve(md,DiagnosticSolutionEnum );9 md=solve(md,DiagnosticSolutionEnum()); 10 10 11 11 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test430.m
r12804 r13663 7 7 md.diagnostic.restol=0.0001; 8 8 md.cluster=generic('name',oshostname(),'np',3); 9 md=solve(md,DiagnosticSolutionEnum );9 md=solve(md,DiagnosticSolutionEnum()); 10 10 11 11 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test431.m
r12804 r13663 8 8 md.thermal.isenthalpy=1; 9 9 md.initialization.waterfraction=zeros(md.mesh.numberofvertices,1); 10 md=solve(md,SteadystateSolutionEnum );10 md=solve(md,SteadystateSolutionEnum()); 11 11 12 12 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test432.m
r12804 r13663 8 8 md.thermal.isenthalpy=1; 9 9 md.initialization.waterfraction=zeros(md.mesh.numberofvertices,1); 10 md=solve(md,SteadystateSolutionEnum );10 md=solve(md,SteadystateSolutionEnum()); 11 11 12 12 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test501.m
r13657 r13663 4 4 md=setflowequation(md,'macayeal','all'); 5 5 md.cluster=generic('name',oshostname(),'np',3); 6 md=solve(md,DiagnosticSolutionEnum );6 md=solve(md,DiagnosticSolutionEnum()); 7 7 8 8 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test502.m
r12805 r13663 1 md=triangle(model,'../Exp/Pig.exp',20000) ;1 md=triangle(model,'../Exp/Pig.exp',20000).; 2 2 md=setmask(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp'); 3 3 md=parameterize(md,'../Par/Pig.par'); … … 5 5 md=setflowequation(md,'pattyn','all'); 6 6 md.cluster=generic('name',oshostname(),'np',3); 7 md=solve(md,DiagnosticSolutionEnum );7 md=solve(md,DiagnosticSolutionEnum()); 8 8 9 9 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test502.py
r13657 r13663 8 8 from MatlabFuncs import * 9 9 10 md=triangle(model(),'../Exp/Pig.exp',20000 )10 md=triangle(model(),'../Exp/Pig.exp',20000.) 11 11 md=setmask(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp') 12 12 md=parameterize(md,'../Par/Pig.py') -
issm/trunk-jpl/test/NightlyRun/test503.m
r12805 r13663 1 md=triangle(model,'../Exp/Pig.exp',20000 );1 md=triangle(model,'../Exp/Pig.exp',20000.); 2 2 md=setmask(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp'); 3 3 md=parameterize(md,'../Par/Pig.par'); … … 5 5 md=setflowequation(md,'stokes','all'); 6 6 md.cluster=generic('name',oshostname(),'np',3); 7 md=solve(md,DiagnosticSolutionEnum );7 md=solve(md,DiagnosticSolutionEnum()); 8 8 9 9 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test503.py
r13657 r13663 8 8 from MatlabFuncs import * 9 9 10 md=triangle(model(),'../Exp/Pig.exp',20000 )10 md=triangle(model(),'../Exp/Pig.exp',20000.) 11 11 md=setmask(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp') 12 12 md=parameterize(md,'../Par/Pig.py') -
issm/trunk-jpl/test/NightlyRun/test504.m
r12805 r13663 1 md=triangle(model,'../Exp/Pig.exp',20000 );1 md=triangle(model,'../Exp/Pig.exp',20000.); 2 2 md=setmask(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp'); 3 3 md=parameterize(md,'../Par/Pig.par'); 4 4 md=setflowequation(md,'macayeal','all'); 5 5 md.cluster=generic('name',oshostname(),'np',3); 6 md=solve(md,TransientSolutionEnum );6 md=solve(md,TransientSolutionEnum()); 7 7 8 8 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test504.py
r13657 r13663 8 8 from MatlabFuncs import * 9 9 10 md=triangle(model(),'../Exp/Pig.exp',20000 )10 md=triangle(model(),'../Exp/Pig.exp',20000.) 11 11 md=setmask(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp') 12 12 md=parameterize(md,'../Par/Pig.py') -
issm/trunk-jpl/test/NightlyRun/test505.m
r12805 r13663 1 md=triangle(model,'../Exp/Pig.exp',30000 );1 md=triangle(model,'../Exp/Pig.exp',30000.); 2 2 md=setmask(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp'); 3 3 md=parameterize(md,'../Par/Pig.par'); 4 md=extrude(md,3,1 );4 md=extrude(md,3,1.); 5 5 md=setflowequation(md,'macayeal','all'); 6 6 md.cluster=generic('name',oshostname(),'np',3); 7 md=solve(md,TransientSolutionEnum );7 md=solve(md,TransientSolutionEnum()); 8 8 9 9 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test505.py
r13657 r13663 8 8 from MatlabFuncs import * 9 9 10 md=triangle(model(),'../Exp/Pig.exp',30000 )10 md=triangle(model(),'../Exp/Pig.exp',30000.) 11 11 md=setmask(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp') 12 12 md=parameterize(md,'../Par/Pig.py') 13 md.extrude(3,1 )13 md.extrude(3,1.) 14 14 md=setflowequation(md,'macayeal','all') 15 15 md.cluster=generic('name',oshostname(),'np',3) -
issm/trunk-jpl/test/NightlyRun/test506.m
r12805 r13663 1 md=triangle(model,'../Exp/Pig.exp',30000 );1 md=triangle(model,'../Exp/Pig.exp',30000.); 2 2 md=setmask(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp'); 3 3 md=parameterize(md,'../Par/Pig.par'); 4 md=extrude(md,2,1 );4 md=extrude(md,2,1.); 5 5 md=setflowequation(md,'pattyn','all'); 6 6 md.cluster=generic('name',oshostname(),'np',3); 7 md=solve(md,TransientSolutionEnum );7 md=solve(md,TransientSolutionEnum()); 8 8 9 9 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test506.py
r13658 r13663 8 8 from MatlabFuncs import * 9 9 10 md=triangle(model(),'../Exp/Pig.exp',30000 )10 md=triangle(model(),'../Exp/Pig.exp',30000.) 11 11 md=setmask(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp') 12 12 md=parameterize(md,'../Par/Pig.py') 13 md.extrude(2,1 )13 md.extrude(2,1.) 14 14 md=setflowequation(md,'pattyn','all') 15 15 md.cluster=generic('name',oshostname(),'np',3) -
issm/trunk-jpl/test/NightlyRun/test507.m
r13658 r13663 1 md=triangle(model,'../Exp/Pig.exp',30000 );1 md=triangle(model,'../Exp/Pig.exp',30000.); 2 2 md=setmask(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp'); 3 3 md=parameterize(md,'../Par/Pig.par'); 4 md=extrude(md,2,1 );4 md=extrude(md,2,1.); 5 5 md=setflowequation(md,'stokes','all'); 6 6 md.cluster=generic('name',oshostname(),'np',3); 7 md=solve(md,TransientSolutionEnum );7 md=solve(md,TransientSolutionEnum()); 8 8 9 9 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test507.py
r13658 r13663 8 8 from MatlabFuncs import * 9 9 10 md=triangle(model(),'../Exp/Pig.exp',30000 )10 md=triangle(model(),'../Exp/Pig.exp',30000.) 11 11 md=setmask(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp') 12 12 md=parameterize(md,'../Par/Pig.py') 13 md.extrude(2,1 )13 md.extrude(2,1.) 14 14 md=setflowequation(md,'stokes','all') 15 15 md.cluster=generic('name',oshostname(),'np',3) -
issm/trunk-jpl/test/NightlyRun/test508.m
r13658 r13663 1 md=triangle(model,'../Exp/Pig.exp',30000 );1 md=triangle(model,'../Exp/Pig.exp',30000.); 2 2 md=setmask(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp'); 3 3 md=parameterize(md,'../Par/Pig.par'); … … 5 5 md=setflowequation(md,'macayeal','all'); 6 6 md.cluster=generic('name',oshostname(),'np',3); 7 md.timestepping.time_step=0 ;8 md=solve(md,SteadystateSolutionEnum );7 md.timestepping.time_step=0.; 8 md=solve(md,SteadystateSolutionEnum()); 9 9 10 10 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test508.py
r13658 r13663 8 8 from MatlabFuncs import * 9 9 10 md=triangle(model(),'../Exp/Pig.exp',30000 )10 md=triangle(model(),'../Exp/Pig.exp',30000.) 11 11 md=setmask(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp') 12 12 md=parameterize(md,'../Par/Pig.py') … … 14 14 md=setflowequation(md,'macayeal','all') 15 15 md.cluster=generic('name',oshostname(),'np',3) 16 md.timestepping.time_step=0 16 md.timestepping.time_step=0. 17 17 md=solve(md,SteadystateSolutionEnum()) 18 18 -
issm/trunk-jpl/test/NightlyRun/test509.m
r12805 r13663 1 md=triangle(model,'../Exp/Pig.exp',30000 );1 md=triangle(model,'../Exp/Pig.exp',30000.); 2 2 md=setmask(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp'); 3 3 md=parameterize(md,'../Par/Pig.par'); 4 md=extrude(md,3,1 );4 md=extrude(md,3,1.); 5 5 md=setflowequation(md,'pattyn','all'); 6 6 md.cluster=generic('name',oshostname(),'np',3); 7 md.timestepping.time_step=0 ;8 md=solve(md,SteadystateSolutionEnum );7 md.timestepping.time_step=0.; 8 md=solve(md,SteadystateSolutionEnum()); 9 9 10 10 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test509.py
r13658 r13663 8 8 from MatlabFuncs import * 9 9 10 md=triangle(model(),'../Exp/Pig.exp',30000 )10 md=triangle(model(),'../Exp/Pig.exp',30000.) 11 11 md=setmask(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp') 12 12 md=parameterize(md,'../Par/Pig.py') 13 md.extrude(3,1 )13 md.extrude(3,1.) 14 14 md=setflowequation(md,'pattyn','all') 15 15 md.cluster=generic('name',oshostname(),'np',3) 16 md.timestepping.time_step=0 16 md.timestepping.time_step=0. 17 17 md=solve(md,SteadystateSolutionEnum()) 18 18 -
issm/trunk-jpl/test/NightlyRun/test510.m
r13662 r13663 1 md=triangle(model,'../Exp/Pig.exp',20000 );1 md=triangle(model,'../Exp/Pig.exp',20000.); 2 2 md=setmask(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp'); 3 3 md=parameterize(md,'../Par/Pig.par'); 4 md=extrude(md,2,1 );4 md=extrude(md,2,1.); 5 5 md=setflowequation(md,'stokes','all'); 6 6 md.cluster=generic('name',oshostname(),'np',3); 7 md.timestepping.time_step=0 ;8 md=solve(md,SteadystateSolutionEnum );7 md.timestepping.time_step=0.; 8 md=solve(md,SteadystateSolutionEnum()); 9 9 10 10 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test510.py
r13662 r13663 8 8 from MatlabFuncs import * 9 9 10 md=triangle(model(),'../Exp/Pig.exp',20000 )10 md=triangle(model(),'../Exp/Pig.exp',20000.) 11 11 md=setmask(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp') 12 12 md=parameterize(md,'../Par/Pig.py') 13 md.extrude(2,1 )13 md.extrude(2,1.) 14 14 md=setflowequation(md,'stokes','all') 15 15 md.cluster=generic('name',oshostname(),'np',3) 16 md.timestepping.time_step=0 16 md.timestepping.time_step=0. 17 17 md=solve(md,SteadystateSolutionEnum()) 18 18 -
issm/trunk-jpl/test/NightlyRun/test511.m
r13021 r13663 1 md=triangle(model,'../Exp/Pig.exp',10000 );1 md=triangle(model,'../Exp/Pig.exp',10000.); 2 2 md=setmask(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp'); 3 3 md=parameterize(md,'../Par/Pig.par'); … … 6 6 md.geometry.bed=-md.materials.rho_ice/md.materials.rho_water*md.geometry.thickness; 7 7 md.geometry.surface=md.geometry.bed+md.geometry.thickness; 8 md=extrude(md,3,1 );8 md=extrude(md,3,1.); 9 9 md=setflowequation(md,'stokes','all'); 10 10 md=extract(md,md.mask.elementonfloatingice); … … 13 13 md.inversion.iscontrol=1; 14 14 md.inversion.control_parameters={'MaterialsRheologyBbar'}; 15 md.inversion.min_parameters=10 ^6*ones(md.mesh.numberofvertices,1);16 md.inversion.max_parameters=2 *10^9*ones(md.mesh.numberofvertices,1);15 md.inversion.min_parameters=10.^6*ones(md.mesh.numberofvertices,1); 16 md.inversion.max_parameters=2.*10^9*ones(md.mesh.numberofvertices,1); 17 17 md.inversion.nsteps=2; 18 18 md.inversion.cost_functions=101*ones(md.inversion.nsteps,1); 19 19 md.inversion.cost_functions_coefficients=ones(md.mesh.numberofvertices,1); 20 md.inversion.gradient_scaling=10 ^8*ones(md.inversion.nsteps,1);21 md.inversion.maxiter_per_step=2 *ones(md.inversion.nsteps,1);20 md.inversion.gradient_scaling=10.^8*ones(md.inversion.nsteps,1); 21 md.inversion.maxiter_per_step=2.*ones(md.inversion.nsteps,1); 22 22 md.inversion.step_threshold=0.99*ones(md.inversion.nsteps,1); 23 23 md.inversion.vx_obs=md.initialization.vx; md.inversion.vy_obs=md.initialization.vy; 24 24 25 25 md.cluster=generic('name',oshostname(),'np',3); 26 md=solve(md,DiagnosticSolutionEnum );26 md=solve(md,DiagnosticSolutionEnum()); 27 27 28 28 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test512.m
r13662 r13663 1 md=triangle(model,'../Exp/Pig.exp',20000 );1 md=triangle(model,'../Exp/Pig.exp',20000.); 2 2 md=setmask(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp'); 3 3 md=parameterize(md,'../Par/Pig.par'); 4 md=extrude(md,3,1 );4 md=extrude(md,3,1.); 5 5 md=setflowequation(md,'pattyn','all'); 6 6 … … 12 12 md.inversion.nsteps=2; 13 13 md.inversion.cost_functions=[103*ones(md.inversion.nsteps,1) 501*ones(md.inversion.nsteps,1)]; 14 md.inversion.cost_functions_coefficients=ones(md.mesh.numberofvertices,2); md.inversion.cost_functions_coefficients(:,2)=2 *10^-7;15 md.inversion.gradient_scaling=3 *ones(md.inversion.nsteps,1);14 md.inversion.cost_functions_coefficients=ones(md.mesh.numberofvertices,2); md.inversion.cost_functions_coefficients(:,2)=2.*10^-7; 15 md.inversion.gradient_scaling=3.*ones(md.inversion.nsteps,1); 16 16 md.inversion.maxiter_per_step=2*ones(md.inversion.nsteps,1); 17 17 md.inversion.step_threshold=0.99*ones(md.inversion.nsteps,1); … … 19 19 20 20 md.cluster=generic('name',oshostname(),'np',3); 21 md=solve(md,DiagnosticSolutionEnum );21 md=solve(md,DiagnosticSolutionEnum()); 22 22 23 23 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test513.m
r13662 r13663 1 md=triangle(model,'../Exp/Pig.exp',30000 );1 md=triangle(model,'../Exp/Pig.exp',30000.); 2 2 md=setmask(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp'); 3 3 md=parameterize(md,'../Par/Pig.par'); 4 md=extrude(md,3,1 );4 md=extrude(md,3,1.); 5 5 md=setflowequation(md,'macayeal','all'); 6 6 … … 8 8 md.inversion.iscontrol=1; 9 9 md.inversion.control_parameters={'FrictionCoefficient'}; 10 md.inversion.min_parameters=1 *ones(md.mesh.numberofvertices,1);11 md.inversion.max_parameters=200 *ones(md.mesh.numberofvertices,1);10 md.inversion.min_parameters=1.*ones(md.mesh.numberofvertices,1); 11 md.inversion.max_parameters=200.*ones(md.mesh.numberofvertices,1); 12 12 md.inversion.nsteps=2; 13 13 md.inversion.cost_functions=[103*ones(md.inversion.nsteps,1) 501*ones(md.inversion.nsteps,1)]; 14 md.inversion.cost_functions_coefficients=ones(md.mesh.numberofvertices,2); md.inversion.cost_functions_coefficients(:,2)=2 *10^-7;15 md.inversion.gradient_scaling=3 *ones(md.inversion.nsteps,1);16 md.inversion.maxiter_per_step=2 *ones(md.inversion.nsteps,1);14 md.inversion.cost_functions_coefficients=ones(md.mesh.numberofvertices,2); md.inversion.cost_functions_coefficients(:,2)=2.*10^-7; 15 md.inversion.gradient_scaling=3.*ones(md.inversion.nsteps,1); 16 md.inversion.maxiter_per_step=2.*ones(md.inversion.nsteps,1); 17 17 md.inversion.step_threshold=0.99*ones(md.inversion.nsteps,1); 18 18 md.inversion.vx_obs=md.initialization.vx; md.inversion.vy_obs=md.initialization.vy; 19 md.timestepping.time_step=0 ;19 md.timestepping.time_step=0.; 20 20 21 21 md.thermal.penalty_lock=5; 22 22 md.cluster=generic('name',oshostname(),'np',3); 23 md=solve(md,SteadystateSolutionEnum );23 md=solve(md,SteadystateSolutionEnum()); 24 24 25 25 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test513.py
r13662 r13663 28 28 md.inversion.vx_obs=md.initialization.vx 29 29 md.inversion.vy_obs=md.initialization.vy 30 md.timestepping.time_step=0 30 md.timestepping.time_step=0. 31 31 32 32 md.thermal.penalty_lock=5 -
issm/trunk-jpl/test/NightlyRun/test515.m
r12805 r13663 1 md=triangle(model,'../Exp/Pig.exp',30000 );1 md=triangle(model,'../Exp/Pig.exp',30000.); 2 2 md=setmask(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp'); 3 3 md=parameterize(md,'../Par/Pig.par'); 4 md=extrude(md,3,1 );4 md=extrude(md,3,1.); 5 5 md=setflowequation(md,'pattyn','all'); 6 6 md.thermal.stabilization=2; … … 10 10 md.transient.isthermal=1; 11 11 md.transient.isgroundingline=0; 12 md=solve(md,TransientSolutionEnum );12 md=solve(md,TransientSolutionEnum()); 13 13 14 14 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test515.py
r13662 r13663 8 8 from MatlabFuncs import * 9 9 10 md=triangle(model(),'../Exp/Pig.exp',30000 )10 md=triangle(model(),'../Exp/Pig.exp',30000.) 11 11 md=setmask(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp') 12 12 md=parameterize(md,'../Par/Pig.py') -
issm/trunk-jpl/test/NightlyRun/test516.m
r12805 r13663 1 md=triangle(model,'../Exp/Pig.exp',30000 );1 md=triangle(model,'../Exp/Pig.exp',30000.); 2 2 md=setmask(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp'); 3 3 md=parameterize(md,'../Par/Pig.par'); 4 md=extrude(md,3,1 );4 md=extrude(md,3,1.); 5 5 md=setflowequation(md,'pattyn','all'); 6 6 md.thermal.stabilization=2; … … 8 8 md.timestepping.time_step=0; 9 9 md.thermal.penalty_threshold=40; 10 md=solve(md,ThermalSolutionEnum );10 md=solve(md,ThermalSolutionEnum()); 11 11 12 12 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test601.m
r12806 r13663 4 4 md=setflowequation(md,'macayeal','all'); 5 5 md.cluster=generic('name',oshostname(),'np',3); 6 md=solve(md,PrognosticSolutionEnum );6 md=solve(md,PrognosticSolutionEnum()); 7 7 8 8 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test602.m
r12806 r13663 7 7 md.prognostic.spcthickness=md.geometry.thickness; 8 8 md.cluster=generic('name',oshostname(),'np',3); 9 md=solve(md,PrognosticSolutionEnum );9 md=solve(md,PrognosticSolutionEnum()); 10 10 11 11 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test603.m
r12806 r13663 5 5 md=setflowequation(md,'macayeal','all'); 6 6 md.cluster=generic('name',oshostname(),'np',3); 7 md=solve(md,PrognosticSolutionEnum );7 md=solve(md,PrognosticSolutionEnum()); 8 8 9 9 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test604.m
r12806 r13663 4 4 md=setflowequation(md,'macayeal','all'); 5 5 md.cluster=generic('name',oshostname(),'np',3); 6 md=solve(md,SurfaceSlopeSolutionEnum );6 md=solve(md,SurfaceSlopeSolutionEnum()); 7 7 8 8 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test605.m
r12806 r13663 5 5 md=setflowequation(md,'pattyn','all'); 6 6 md.cluster=generic('name',oshostname(),'np',3); 7 md=solve(md,SurfaceSlopeSolutionEnum );7 md=solve(md,SurfaceSlopeSolutionEnum()); 8 8 9 9 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test606.m
r12806 r13663 4 4 md=setflowequation(md,'macayeal','all'); 5 5 md.cluster=generic('name',oshostname(),'np',3); 6 md=solve(md,BedSlopeSolutionEnum );6 md=solve(md,BedSlopeSolutionEnum()); 7 7 8 8 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test607.m
r12806 r13663 5 5 md=setflowequation(md,'macayeal','all'); 6 6 md.cluster=generic('name',oshostname(),'np',3); 7 md=solve(md,BedSlopeSolutionEnum );7 md=solve(md,BedSlopeSolutionEnum()); 8 8 9 9 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test608.m
r12806 r13663 4 4 md=setflowequation(md,'macayeal','all'); 5 5 md.cluster=generic('name',oshostname(),'np',3); 6 md=solve(md,BalancethicknessSolutionEnum );6 md=solve(md,BalancethicknessSolutionEnum()); 7 7 8 8 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test609.m
r12806 r13663 6 6 md.balancethickness.stabilization=3; 7 7 md.cluster=generic('name',oshostname(),'np',3); 8 md=solve(md,BalancethicknessSolutionEnum );8 md=solve(md,BalancethicknessSolutionEnum()); 9 9 10 10 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test610.m
r12806 r13663 5 5 md=setflowequation(md,'macayeal','all'); 6 6 md.cluster=generic('name',oshostname(),'np',3); 7 md=solve(md,BalancethicknessSolutionEnum );7 md=solve(md,BalancethicknessSolutionEnum()); 8 8 9 9 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test611.m
r12806 r13663 19 19 20 20 md.cluster=generic('name',oshostname(),'np',3); 21 md=solve(md,BalancethicknessSolutionEnum );21 md=solve(md,BalancethicknessSolutionEnum()); 22 22 23 23 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test612.m
r12806 r13663 20 20 21 21 md.cluster=generic('name',oshostname(),'np',3); 22 md=solve(md,BalancethicknessSolutionEnum );22 md=solve(md,BalancethicknessSolutionEnum()); 23 23 24 24 %Fields and tolerances to track changes -
issm/trunk-jpl/test/NightlyRun/test613.m
r13021 r13663 27 27 md.verbose.control=1; 28 28 md.cluster=generic('name',oshostname(),'np',3); 29 md=solve(md,BalancethicknessSolutionEnum );29 md=solve(md,BalancethicknessSolutionEnum()); 30 30 31 31 %Fields and tolerances to track changes
Note:
See TracChangeset
for help on using the changeset viewer.