Enter a species  (e.g. Al, Fe): Fe
Enter a hexagonal lattice type (e.g. hcp): hcp
Enter a model name: Morse_Shifted_GirifalcoWeizer_1959LowCutoff_Fe__MO_331285495617_002
Enter hexagonal lattice constants 'a' and 'c' (m): 2.574566084107211e-10, 4.187412606010261e-10

Calculating isothermal bulk modulus

Optimization terminated successfully.
         Current function value: -7.806334
         Iterations: 71
         Function evaluations: 185

Finished calculation

C11 = 340.61099738151216 GPa
C12 = 84.77760972941411 GPa
  B = 170.0554056134468 GPa