Changeset 14124


Ignore:
Timestamp:
12/10/12 19:38:20 (12 years ago)
Author:
seroussi
Message:

Finishing up tutorial

Location:
issm/workshop/2012/Talks/07_Ice_flow_models
Files:
8 added
1 edited

Legend:

Unmodified
Added
Removed

  • issm/workshop/2012/Talks/07_Ice_flow_models/Code/runme.m

    r14072 r14124  
    1 >> md.flowequation
    2 
    3 ans =
    4 
    5    flow equation parameters:
    6                  ismacayealpattyn       : 0    -- is the macayeal or pattyn approximation used ?
    7                                 ishutter               : 0    -- is the shallow ice approximation used ?
    8                                 isstokes               : 0    -- are the Full-Stokes equations used ?
    9                                 vertex_equation        : N/A  -- flow equation for each vertex
    10                                 element_equation       : N/A  -- flow equation for each element
    11                                 bordermacayeal         : N/A  -- vertices on MacAyeal's border (for tiling)
    12                                 borderpattyn           : N/A  -- vertices on Pattyn's border (for tiling)
    13                                 borderstokes           : N/A  -- vertices on Stokes' border (for tiling)
    14 
    15 md=setflowequation(md,'hutter','all')
    16 md=setflowequation(md,'stokes','all')
    17 md=setflowequation(md,'macayeal','all')
    18 md=setflowequation(md,'pattyn','all')
    19 
    20 >> plotmodel(md,'data','elements_type')
    21 
    22 >> md.diagnostic
    23 
    24 ans =
    25 
    26    Diagnostic solution parameters:
    27 
    28      Convergence criteria:
    29              restol                 : 0.0001          -- mechanical equilibrium residue convergence criterion
    30                   reltol                 : 0.01            -- velocity relative convergence criterion, NaN -> not applied
    31                   abstol                 : 10              -- velocity absolute convergence criterion, NaN -> not applied
    32              maxiter                : 100             -- maximum number of nonlinear iterations
    33                   viscosity_overshoot    : 0               -- over-shooting constant new=new+C*(new-old)
    34 
    35      boundary conditions:
    36                   spcvx                  : N/A             -- x-axis velocity constraint (NaN means no constraint)
    37                   spcvy                  : N/A             -- y-axis velocity constraint (NaN means no constraint)
    38                   spcvz                  : N/A             -- z-axis velocity constraint (NaN means no constraint)
    39                   icefront               : N/A             -- segments on ice front list (last column 0-> Air, 1-> Water, 2->Ice
    40 
    41      Rift options:
    42                   rift_penalty_threshold : 0               -- threshold for instability of mechanical constraints
    43                   rift_penalty_lock      : 10              -- number of iterations before rift penalties are locked
    44 
    45      Penalty options:
    46                   penalty_factor         : 3               -- offset used by penalties: penalty = Kmax*10^offset
    47                   vertex_pairing         : N/A             -- pairs of vertices that are penalized
    48 
    49      Other:
    50                   shelf_dampening        : 0               -- use dampening for floating ice ?  Only for Stokes model
    51                   stokesreconditioning   : 10000000000000  -- multiplier for incompressibility equation.  Only for Stokes model
    52                   referential            : N/A             -- local referential
    53                   requested_outputs      : N/A             -- additional outputs requested
    54 
    55 >> md=SetIceSheetBC(md)
    56 >> md=SetIceShelfBC(md,'Front.exp')
    57 >> md=SetMarineIceSheefBC(md,'Front.exp')
    58 
    59 >> plotmodel(md,'data','BC')
    60 
    61 md=setflowequation(md,'pattyn',md.elementongroundedice,'fill','macayeal','coupling','penalties')
    62 md=setflowequation(md,'pattyn',md.elementongroundedice,'fill','macayeal','coupling','tiling')
    63 md=setflowequation(md,'stokes','Contour.exp','fill','pattyn')
    64 
    65 >> exptool('Contour.exp')
    66 
    67 >> plotmodel(md,'data','elements_type','edgecolor','k','expdisp','Contour.exp')
    68 
    69 md=setflowequation(md,'macayeal','FloatingIce.exp','fill','pattyn','coupling','penalties')
     1md=triangle(model(),'Square.exp',50000.);
     2md=setmask(md,'all','');
     3md=parameterize(md,'SquareShelf.par');
     4md=extrude(md,3,1);
     5%Set flow equation
     6...
     7%Solve
Note: See TracChangeset for help on using the changeset viewer.