User: tbode Date: 2012/11/21 04:57 PM

Modified: /trunk/src/ GRHydro_PoloidalMagFieldM.F90

Log: Rename conflicting temporary Aphi

File Changes:

Directory: /trunk/src/ ======================

File [modified]: GRHydro_PoloidalMagFieldM.F90 Delta lines: +7 -7 =================================================================== --- trunk/src/GRHydro_PoloidalMagFieldM.F90 2012-11-09 01:59:28 UTC (rev 187) +++ trunk/src/GRHydro_PoloidalMagFieldM.F90 2012-11-21 22:57:45 UTC (rev 188) @@ -61,7 +61,7 @@ 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 @@ 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 @@ 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