Reads geometric information from ipfiles for Auckland's full heart mesh (geometry fibres and sheets) and exports node and element files to cmgui.
 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.