From fba24a1753d874047237908fa4b71ab70293a970 Mon Sep 17 00:00:00 2001 From: tbode Date: Wed, 21 Nov 2012 22:57:45 +0000 Subject: Rename conflicting temporary Aphi git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/GRHydro_InitData/trunk@188 ac85fae7-cede-4708-beff-ae01c7fa1c26 --- src/GRHydro_PoloidalMagFieldM.F90 | 14 +++++++------- 1 file changed, 7 insertions(+), 7 deletions(-) diff --git a/src/GRHydro_PoloidalMagFieldM.F90 b/src/GRHydro_PoloidalMagFieldM.F90 index e4ac43f..52e9ca8 100644 --- a/src/GRHydro_PoloidalMagFieldM.F90 +++ b/src/GRHydro_PoloidalMagFieldM.F90 @@ -61,7 +61,7 @@ subroutine GRHydro_PoloidalMagFieldM(CCTK_ARGUMENTS) CCTK_REAL :: sdet CCTK_REAL :: dx,dy,dz CCTK_REAL :: rhofac, delPcut, maxP_Pcut - CCTK_REAL :: Aphi, Ax, Ay, Az + CCTK_REAL :: AphiL, Ax, Ay, Az CCTK_REAL :: rho_dx, rho_dy, rho_dz CCTK_REAL :: press_dx, press_dy, press_dz CCTK_REAL :: Aphi_dx, Aphi_dy, Aphi_dz @@ -84,14 +84,14 @@ subroutine GRHydro_PoloidalMagFieldM(CCTK_ARGUMENTS) rhofac = 1.0d0-rho(i,j,k)/poloidal_rho_max delPcut = press(i,j,k)-poloidal_P_cut maxP_Pcut = max(delPcut,0.0d0) - Aphi = poloidal_A_b*rhofac**poloidal_n_p*maxP_Pcut - Ax = -y(i,j,k)*Aphi - Ay = x(i,j,k)*Aphi + AphiL = poloidal_A_b*rhofac**poloidal_n_p*maxP_Pcut + Ax = -y(i,j,k)*AphiL + Ay = x(i,j,k)*AphiL Az = 0.0 !! write(*,*)'Before accessing rho(i,k,k)' !! write(*,*)'rho(',i,',',j,',',k,') = ', rho(i,j,k) -!! write(*,*)'Ax, Ay, Az, Aphi = ', Ax, Ay, Az,Aphi +!! write(*,*)'Ax, Ay, Az, Aphi = ', Ax, Ay, Az,AphiL !! write(*,*)'rhofac = ', rhofac !! write(*,*)'delPcut = ', delPcut !! write(*,*)'maxP_Pcut = ', maxP_Pcut @@ -110,9 +110,9 @@ subroutine GRHydro_PoloidalMagFieldM(CCTK_ARGUMENTS) Aphi_dz = poloidal_A_b*rhofac**(poloidal_n_p-1)*maxP_Pcut* & ( rhofac*press_dz/delPcut - poloidal_n_p*rho_dz/poloidal_rho_max ) - Ax_dy = -Aphi - y(i,j,k)*Aphi_dy + Ax_dy = -AphiL - y(i,j,k)*Aphi_dy Ax_dz = -y(i,j,k)*Aphi_dz - Ay_dx = Aphi + x(i,j,k)*Aphi_dx + Ay_dx = AphiL + x(i,j,k)*Aphi_dx Ay_dz = x(i,j,k)*Aphi_dz -- cgit v1.2.3