[Commits] [svn:einsteintoolkit] GRHydro_InitData/trunk/src/ (Rev. 188)

tanja.bode at physics.gatech.edu tanja.bode at physics.gatech.edu
Wed Nov 21 16:57:46 CST 2012


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
 
            



More information about the Commits mailing list