# 
# crust model file
# PREM + mmc crust setup
#
#
#   format :
#   upper_radius[m] density[g/cm^3] label[string]
#    - no space before first column
#    - delimiter must be A whitespace
#    - data must include at least 4 lines with labels
#      "moho_boundary", "rockice_boundary",
#      "earth_radius", and "atmo_radius"
#    - "earth_radius" must be the surface of ice
#    - lines start from '#' or whitespace are treated
#      as comments
#    - density of "moho_boundary" should be same as
#      the density of next layer
#

6346600 2.900 moho_boundary
6356000 2.900 inner_crust           # 
6371324 2.650 rockice_boundary      # surface of bed rock
6373934 0.917 light_ice_boundary    # 200m below ice surface
6374134 0.757 earth_radius          # surface of ice, 0.825 x 0.917
6478000 0.000811 atmo_radius        # 0.673atm x 1.205e-3(g/cm3 for 1atm)