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

Calculating isothermal bulk modulus

Optimization terminated successfully.
         Current function value: -8.314815
         Iterations: 74
         Function evaluations: 197

Finished calculation

C11 = 348.36645980192884 GPa
C12 = 85.21055690619032 GPa
  B = 172.92919120476984 GPa