istp/gray
2
0
Fork 0
Code Issues Pull Requests Packages Projects Releases 1 Wiki Activity

src/gray_core.f90: remove unused variables in rhs

Browse Source
This commit is contained in:
Michele Guerini Rocco 2023-09-21 13:59:21 +02:00
parent f75cedbf2b
commit 2fed2ab5ec
Signed by: rnhmjoj
GPG Key ID: BFBAF4C975F76450
1 changed files with 8 additions and 13 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

21
src/gray_core.f90
View File

@ -979,38 +979,33 @@ contains
integer, intent(in) :: igrad,sox integer, intent(in) :: igrad,sox
real(wp_), dimension(6) :: yy,fk1,fk2,fk3,fk4 real(wp_), dimension(6) :: yy,fk1,fk2,fk3,fk4
real(wp_) :: gr2
real(wp_), dimension(3) :: dgr2
gr2 = dot_product(dgr, dgr)
dgr2 = 2 * matmul(ddgr, dgr)
fk1 = yp
fk1 = yp
yy = y + fk1*hh yy = y + fk1*hh
call rhs(sox,bres,xgcn,yy,gr2,dgr2,dgr,ddgr,fk2,igrad) call rhs(sox,bres,xgcn,yy,dgr,ddgr,fk2,igrad)
yy = y + fk2*hh yy = y + fk2*hh
call rhs(sox,bres,xgcn,yy,gr2,dgr2,dgr,ddgr,fk3,igrad) call rhs(sox,bres,xgcn,yy,dgr,ddgr,fk3,igrad)
yy = y + fk3*h yy = y + fk3*h
call rhs(sox,bres,xgcn,yy,gr2,dgr2,dgr,ddgr,fk4,igrad) call rhs(sox,bres,xgcn,yy,dgr,ddgr,fk4,igrad)
y = y + h6*(fk1 + 2*fk2 + 2*fk3 + fk4) y = y + h6*(fk1 + 2*fk2 + 2*fk3 + fk4)
end subroutine rkstep end subroutine rkstep
subroutine rhs(sox,bres,xgcn,y,gr2,dgr2,dgr,ddgr,dery,igrad) subroutine rhs(sox,bres,xgcn,y,dgr,ddgr,dery,igrad)
! Compute right-hand side terms of the ray equations (dery) ! Compute right-hand side terms of the ray equations (dery)
! used in R-K integrator ! used in R-K integrator
implicit none implicit none
! arguments ! arguments
real(wp_), dimension(6), intent(in) :: y real(wp_), dimension(6), intent(in) :: y
real(wp_), intent(in) :: bres,xgcn,gr2 real(wp_), intent(in) :: bres,xgcn
real(wp_), dimension(3), intent(in) :: dgr2,dgr real(wp_), dimension(3), intent(in) :: dgr
real(wp_), dimension(3,3), intent(in) :: ddgr real(wp_), dimension(3,3), intent(in) :: ddgr
real(wp_), dimension(6), intent(out) :: dery real(wp_), dimension(6), intent(out) :: dery
integer, intent(in) :: igrad,sox integer, intent(in) :: igrad,sox
! local variables ! local variables
real(wp_) :: psinv,dens,btot,xg,yg real(wp_) :: dens,btot,xg,yg
real(wp_), dimension(3) :: xv,anv,bv,derxg,deryg real(wp_), dimension(3) :: xv,anv,bv,derxg,deryg
real(wp_), dimension(3,3) :: derbv real(wp_), dimension(3,3) :: derbv