Ignore:
Timestamp:
09/08/11 14:08:39 (14 years ago)
Author:
Mathieu Morlighem
Message:

Added geometry

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/m/model/basevert.m

    r9684 r9691  
    2020end
    2121
    22 hu=md.thickness.*md.initialization.vx;
    23 hv=md.thickness.*md.initialization.vy;
     22hu=md.geometry.thickness.*md.initialization.vx;
     23hv=md.geometry.thickness.*md.initialization.vy;
    2424
    2525summation=[1;1;1];
     
    3030velem=md.initialization.vy(md.elements)*summation/3;
    3131
    32 dbdx=(md.bed(md.elements).*alpha)*summation;
    33 dbdy=(md.bed(md.elements).*beta)*summation;
     32dbdx=(md.geometry.bed(md.elements).*alpha)*summation;
     33dbdy=(md.geometry.bed(md.elements).*beta)*summation;
    3434
    3535wb=-md.materials.rho_ice/md.materials.rho_water*(hux+hvy)+uelem.*dbdx+velem.*dbdy;
Note: See TracChangeset for help on using the changeset viewer.