/home/uryu/Cocal_svn_checkinout/Cocal/code/Subroutine/length_conversion.f90

Go to the documentation of this file.

subroutine length_conversion(i,inlen,outlen)
  use phys_constant, only : long,nmpt, g, c, solmas
  use def_matter_parameter_mpt
  implicit none
  real(long) :: inlen, outlen, R0
  integer :: i
  real(long) :: MM = solmas, LL = g*solmas/c**2, TT = g*solmas/c**3

  R0 = def_matter_param_real_(5,1)

! cocal units to km
  if(i.eq.1) outlen = inlen*LL*R0*1.0d-5

! km to cocal units
  if(i.eq.2) outlen = inlen/(LL*R0*1.0d-5)
!
!  write(6,'(a16,1p,e23.15)') 'Input length  = ', inlen
!  write(6,'(a16,1p,e23.15)') 'Output length = ', outlen
!
end subroutine length_conversion

---