c---------------------------------------------------------------------
c include file for stump model
c
c*notes
c formulated for plane stress elements 
c   NDI = 2, NSHEAR = 1
c
c---------------------------------------------------------------------
      integer
     *    elm_no, elm_tot, int_points, maxcrd, ss_comps, 
     *    sim_num, elm_list, elm_list_prep, elm_material_check, 
     *    elm_remodel_check, maxdata_ST
      real*8
     *    rho, rho2, e_mod, poiss_ST, poiss_bone, 
     *    dr, phi, e_tissue, e_socket, e_skin,  beta, 
     *    eo, so, tot_time, d_t,  
     *    k_trabecular, k_cortical, k_soft_tissue, k_ho, 
     *    min_cortical_rho, max_cortical_rho, cutoff, stress_hydro,
     *    max_trabecular_rho, min_rho, cort_initial_R, trab_initial_R,
     *    ST_initial_R, non_remodel_R, SED_store, max_HO_rho, 
     *    skin_initial_R, sx_store, sy_store, coord_group,  
     *    inc_max, s_val_HO, s_val_femur, d_val, rho_recruit,
     *    el_area, el_mass, mass_HO, mass_bone, mass_femur, mass_all
      parameter
     *(   elm_no =  7446               ,elm_tot = 9748               ,
     *    int_points = 1               ,maxcrd = 2                   ,
     *    ss_comps = 3                 ,maxdata_ST = 16              ,
     *    d_t = 1.0d0                  ,             
     *    k_cortical = 0.0429d0        ,k_trabecular = 0.0125d0      ,
     *    k_ho = 0.0029d0              ,k_soft_tissue = 0.0029d0     ,
     *    poiss_ST = 0.47d0            ,poiss_bone = 0.3d0           ,
     *    e_socket = 1500000.0d0       ,s_val_HO = 0.0d0             ,
     *    d_val = 0.35d0               ,s_val_femur = 0.05d0         ,
     *    cutoff = -0.0d0              ,
     *    cort_initial_R = 1.7d0       ,trab_initial_R = 0.8d0       ,
     *    ST_initial_R = 0.95d0        ,skin_initial_R = 1.02d0      , 
     *    non_remodel_R = 0.0d0        ,
     *    rho_recruit = 1.2d0          ,
     *    min_cortical_rho = 1.2d0     ,max_cortical_rho = 2.0d0     ,
     *    max_trabecular_rho = 1.2d0   ,min_rho = 0.1d0              ,
     *    max_HO_rho = 2.2d0           ,
     *    inc_max= 1.0d0                                             )
      data
     *    eo                                                         /
     *    0.0d0      ,0.05d0     ,0.1d0      ,0.15d0     ,0.2d0      ,
     *    0.25d0     ,0.3d0      ,0.35d0     ,0.40d0     ,0.45d0     ,
     *    0.50d0     ,0.55d0     ,0.60d0     ,0.65d0     ,0.70d0     ,
     *    0.75d0                                                     /
      data
     *    so                                                         /
     *    0.0d0      ,0.39d0     ,1.08d0     ,1.88d0     ,3.00d0     ,
     *    4.66d0     ,7.06d0     ,10.42d0    ,14.96d0    ,20.88d0    ,
     *    28.39d0    ,37.71d0    ,49.05d0    ,62.62d0    ,78.64d0    ,
     *    97.31d0                                                    /
      dimension 
     *    elm_list_prep(elm_tot)       ,rho(elm_tot)                 ,
     *    rho2(elm_tot)                ,dr(elm_tot)                  ,
     *    phi(elm_tot)                 ,elm_list(elm_no)             ,
     *    e_mod(elm_tot)               ,elm_material_check(elm_tot)  ,
     *    elm_remodel_check(elm_tot)   ,beta(elm_tot)                ,
     *    eo(maxdata_ST)               ,so(maxdata_ST)               ,
     *    stress_hydro(elm_tot)        ,e_tissue(ss_comps)           ,
     *    SED_store(elm_tot)           ,sx_store(elm_tot)            ,
     *    sy_store(elm_tot)            ,el_area(elm_tot)             ,
     *    el_mass(elm_tot)             ,coord_group(int_points, maxcrd)
      common
     *   /user_nr/
     *   rho, rho2, dr, phi, elm_list, e_mod, SED_store, 
     *   elm_material_check, elm_remodel_check, tot_time, sim_num, 
     *   beta, stress_hydro, sx_store, sy_store, el_area

c---------------------------------------------------------------------
c*description
c d_t              - 
c k_trabecular     - 
c k_cortical       - 
c k_soft_tissue    - 
c k_ho             - 
c s_val            - 
c rho_recruit      - 
c                        used if want to use spacial influence, 
c                        might not be computationally worth it
c d_val            - 
c min_cortical_rho - 
c cutoff           -   
c                        material properties (hypela2)
c poiss            - poisson's ratio
c c1               - 
c e_socket         - youngs modulus for socket
c e_skin           - youngs modulus for skin
c---------------------------------------------------------------------
c