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

Go to the documentation of this file.

subroutine calc_surface(rsnew,emd)
  use phys_constant, only  : long
  use grid_parameter, only : nrf, ntf, npf
  use coordinate_grav_r, only : drg
  use def_matter, only : rs
!  use def_matter_parameter, only : emdc
  implicit none
  real(long), pointer :: emd(:,:,:)
  real(long), pointer :: rsnew(:,:)
  real(long) :: delta
  integer    :: itf, ipf
  do ipf = 0, npf
    do itf = 0, ntf
      delta = drg(nrf)*(0.0d0            - emd(nrf  ,itf,ipf)) &
    &                 /(emd(nrf,itf,ipf) - emd(nrf-1,itf,ipf))
      rsnew(itf,ipf) = rs(itf,ipf)*(1.0d0 + delta)
    end do
  end do
  rsnew(0,1:npf) = rsnew(0,0)
end subroutine calc_surface

---