!--------1---------2---------3---------4---------5---------6---------7---------8
module parameters 
!--------1---------2---------3---------4---------5---------6---------7---------8
  use global_constants, only : rkind, itmax, dt, time0, time, year &
    & , set_mass_ratio_rock_H2O_init, imr, GO_NEXT, W_ALTERATION, D_REACTION &
    & , WV_PRESSURE, Rocky_Core, Ordinary_Chondrite
  use global_parameters, only : irb, iau, irbmax, iaumax &
    & , set_radius, set_temperature, set_radius_body, set_temperature_initial 
  use heat_radiogenic, only : heat_radiogenic_Al, heat_radiogenic_Fe &
    & , heat_radiogenic_K
  implicit none
  private
  public :: set_parameters
contains
  subroutine set_parameters
    if(.NOT.Ordinary_Chondrite) call set_mass_ratio_rock_H2O_init
    call set_radius_body
    call set_temperature_initial
    call heat_radiogenic_Al
    call heat_radiogenic_Fe
    call heat_radiogenic_K
  if(Ordinary_Chondrite) then
    Rocky_Core=.FALSE.
    W_ALTERATION=.FALSE.
    WV_PRESSURE=.FALSE.
  endif
    Rocky_Core=.FALSE.
    W_ALTERATION=.TRUE.
    if(W_ALTERATION) D_REACTION=.FALSE.
    if(.NOT.W_ALTERATION .AND. D_REACTION) stop
    WV_PRESSURE=.FALSE.
  end subroutine set_parameters
end module parameters