# 
# crust model file
# modified PREM + mmc crust setup
# total column depth of crust is estimated 
# 38km x 2.7g/cm3 
# (Adrienne E.Block et.al, Earth and Planetary Science Letters 288, 2009, 194-203)
# To adjust earth's size, moho-boundary radius is reduced.
# (which is not correct, but may not affect for our purpose...)
#
# You must use mediadef_pon_38km too!
#
#   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
#

6334024 2.900 moho_boundary         # inner crust: 2.9g/cm3, 9.4km
6343424 2.900 inner_crust           # outer crust: 2.7g/cm3, 27.9km
6371324 1.890 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)