Ignore:
Timestamp:
10/28/13 11:27:31 (11 years ago)
Author:
Mathieu Morlighem
Message:

merged trunk and trunk-jpl

Location:
issm/trunk-jpl
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl

  • issm/trunk-jpl/examples/Jakobshavn/runme.m

    r11676 r16555  
    4242        md=loadmodel('JksPar.mdl');
    4343
    44         md=setflowequation(md,'macayeal','all');
     44        md=setflowequation(md,'SSA','all');
    4545
    4646        %Control general
     
    6464
    6565        %Additional parameters
    66         md.diagnostic.restol=0.01;
    67         md.diagnostic.reltol=0.1;
    68         md.diagnostic.abstol=NaN;
     66        md.stressbalance.restol=0.01;
     67        md.stressbalance.reltol=0.1;
     68        md.stressbalance.abstol=NaN;
    6969
    7070        %Go solve
    7171        md.cluster=generic('name',oshostname,'np',4);
    72         md.solver=addoptions(md.solver,NoneAnalysisEnum,asmoptions);
    73         md.solver=addoptions(md.solver,DiagnosticVertAnalysisEnum,jacobiasmoptions);
    74         md=solve(md,DiagnosticSolutionEnum);
     72        md=solve(md,StressbalanceSolutionEnum);
    7573       
    7674        save JksControl.mdl md
     
    8381                'FontSize#all',12,...
    8482                'data',md.inversion.vel_obs,'title','Observed velocity',...
    85                 'data',md.results.DiagnosticSolution.Vel,'title','Modeled Velocity',...
     83                'data',md.results.StressbalanceSolution.Vel,'title','Modeled Velocity',...
    8684                'colorbar#1','off','colorbar#2','on','colorbartitle#2','[m/yr]',...
    8785                'caxis#1-2',[0,7000],...
    8886                'data',md.geometry.bed,'title','Bed elevation',...
    89                 'data',md.results.DiagnosticSolution.FrictionCoefficient,...
     87                'data',md.results.StressbalanceSolution.FrictionCoefficient,...
    9088                'title','Friction Coefficient',...
    9189                'colorbar#3','on','colorbartitle#3','[m]','colorbar#4','on');
Note: See TracChangeset for help on using the changeset viewer.