Enter a species  (e.g. Al, Fe): Ar
Enter a hexagonal lattice type (e.g. hcp): hcp
Enter a model name: Morse_Shifted_Jelinek_1972_Ar__MO_831902330215_003
Enter hexagonal lattice constants 'a' and 'c' (m): 3.714955519697364e-10, 6.062858905024546e-10

Calculating isothermal bulk modulus

Optimization terminated successfully.
         Current function value: -0.185738
         Iterations: 72
         Function evaluations: 187

Finished calculation

C11 = 6.072974235767956 GPa
C12 = 1.609910917785639 GPa
  B = 3.0975986904464112 GPa