Reads geometric information from heart.iod for Auckland's full heart mesh (geometry fibres and sheets).
Heart H27 nodal parameters. (Note: this is H27 geometry with sheets and H02 fibres) There are 99 nodes which are used in 60 elements using the prolate spheroidal coordinate system (lambda, mu, theta). See Nielsen, LeGrice, Smaill and Hunter "Mathematical model of geometry and fibrous structure of the heart" AJP 260 H1365-H1378, 1991. x-axis is main axis of heart from base to apex y-axis goes through middle of RV z-axis completes right-handed orthogonal system XP(1,1,np) is lambda coord at node np XP(2,1,np) is derivative of lambda wrt s1 at node np XP(3,1,np) is derivative of lambda wrt s2 at node np XP(4,1,np) is cross deriv of lambda wrt s1 & s2 at node np XP(1,2,np) is mu coord (radians) at node np XP(1,3,np) is theta coord (radians) at node np XP(1,4,np) is fibre angle (radians) at node np XP(2,4,np) is deriv of fibre angle wrt s3 (radians) at node np XP(1,5,np) is imbrication angle (not used here) XP(1,6,np) is sheet angle XP(2,6,np) is derivative of sheet angle Basis nb=1 is bicubic-linear (used for lambda) Basis nb=2 is trilinear (used for mu & theta) Basis nb=3 is bilinear-cubic (used for fibre field) Basis nb=8 is linear-bicubic (used for sheet field) Nodes on the epicardium are, in layers from base to apex and in clockwise order (viewed from base towards apex) starting at z=0: 1st layer: 1,44, 2,56,47, 3,50,59, 4,53 2nd layer: 5,43, 6,55,46, 7,49,58, 8,52 3rd layer: 9,42,10,54,45,11,48,57,12,51 apex: 13(10 versions for theta) Nodes on the LV endocardium are, in layers from base to apex and in clockwise order (viewed from base towards apex) starting at z=0: 1st layer: 29,62,30,74,65,31,68,77,32,71 2nd layer: 33,61,34,73,64,35,67,76,36,70 3rd layer: 37,60,38,72,63,39,66,75,40,69 apex: 41(10 versions for theta) Nodes on the LV midwall are, in layers from base to apex and in clockwise order (viewed from base towards apex): 1st layer: 15,90,93,16,96,99,17 (continues on to RV endo) 2nd layer: 20,89,92,21,95,96,22 (continues on to RV endo) 3rd layer: 24,78,25,88,91,26,94,97,27,81 apex: 28(10 versions for theta) Nodes on the RV free wall endo are, in layers from base to apex and in clockwise order (viewed from base towards apex) starting at z=0: 1st layer: 14,80,15(2 versions for fibres),17(2 versions for fibres),83 2nd layer: 19,79,20(2 versions for fibres),22(2 versions for fibres),82 3rd layer: 24(2 versions for fibres),78(2 versions for fibres), 25(2 versions for fibres),27(2 versions for fibres), 81(2 versions for fibres) Nodes on the RV septal wall endo are, in layers from base to apex and in clockwise order (viewed from base towards apex) starting at z=0: 1st layer: 18,85,15(2 versions for fibres),17(2 versions for fibres),87 2nd layer: 23,84,20(2 versions for fibres),22(2 versions for fibres),86 3rd layer: 24(2 versions for fibres),78(2 versions for fibres), 25(2 versions for fibres),27(2 versions for fibres), 81(2 versions for fibres) The RV free wall has elements: 1, 2, 9,10 11,12, 19,20 The septal wall has elements: 31,32, 39,40 41,42, 49,50, The LV free wall has elements: (subepi) (subendo) 3.. 8 33..38 13..18 43..48 21,22,23..28,29,30 51,52,53..58,59,60 NPE(nn,nb,ne) is the global node# (np) associated with local node nn of element ne for basis nb. These (see below) give the connectivity of the mesh. SE(ns,nb,ne) is the scaling factor associated with element variable ns for basis nb of element ne where, for nb=1, ns is 1..4 for the 4 bicubic derivs in the Xi_1 & 2 dir.s at node nn=1 ns is 5..8 " " " " at node nn=2 . . . . ns is 29..32 " " " " at node nn=8 and for nb=3, ns is 1,2 for the 2 cubic derivs in the Xi_3 direction at node nn=1 . . ns is 15,16 " " " " at node nn=8 The way in which the scaling factors are used is explained in the AJP paper.