Cocal: /home/uryu/Cocal_svn_checkinout/Cocal/code/Main_utility/calc_crust

00001 include '../Module/phys_constant.f90'
00002 include '../EOS/Module/def_peos_parameter.f90'
00003 include '../EOS/Subroutine/peos_initialize.f90'
00004 include '../EOS/Subroutine/peos_lookup.f90'
00005 !
00006 program calc_crust_interface
00007 !
00008   use phys_constant                     !g,c,solmas,nnpeos
00009   use def_peos_parameter                !abc,abi,rhoi,qi,hi,nphase
00010   implicit none
00011 !
00012   real(8) :: rho_int, abccgs_crust, abi_crust
00013 !
00014   call peos_initialize
00015 !
00016   abi_crust = 1.35692d+0
00017   abccgs_crust = 3.99874d-08*c**2
00018   rho_int = (abccgs(1)/abccgs_crust)**(1.0d0/(abi_crust - abi(1)))
00019   write(6,'(a23,es13.5)') 'crust Gamma            ', abi_crust
00020   write(6,'(a23,es13.5)') 'crust adiabtic constant', abccgs_crust
00021   write(6,'(a23,es13.5)') 'crust interface density', rho_int
00022 !
00023 end program calc_crust_interface

/home/uryu/Cocal_svn_checkinout/Cocal/code/Main_utility/calc_crust_interface.f90