Change of Con2Prim and Prim2Con interface.

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/GRHydro_InitData/trunk@126 ac85fae7-cede-4708-beff-ae01c7fa1c26

3 files changed, 40 insertions, 15 deletions

diff --git a/src/GRHydro_C2P2CM.F90 b/src/GRHydro_C2P2CM.F90
index 6c4cd43..dcea38f 100644
--- a/src/GRHydro_C2P2CM.F90
+++ b/src/GRHydro_C2P2CM.F90
@@ -41,6 +41,7 @@ subroutine c2p2cM(CCTK_ARGUMENTS)
   CCTK_REAL :: uxx,uxy,uxz,uyy,uyz,uzz
   CCTK_REAL :: gxx_send,gxy_send,gxz_send,gyy_send,gyz_send,gzz_send
   CCTK_REAL :: dens_send,sx_send,sy_send,sz_send,tau_send
+  CCTK_REAL :: bconsx_send, bconsy_send, bconsz_send
   CCTK_REAL :: rho_send,velx_send,vely_send,velz_send,eps_send
   CCTK_REAL :: press_send,w_lorentz_send,x_send,y_send,z_send,r_send
   CCTK_REAL :: bvcx_send, bvcy_send, bvcz_send, b2_send
@@ -93,6 +94,9 @@ subroutine c2p2cM(CCTK_ARGUMENTS)
   bvcx_send   = Bx_init
   bvcy_send   = By_init
   bvcz_send   = Bz_init
+  bconsx_send   = Bx_init
+  bconsy_send   = By_init
+  bconsz_send   = Bz_init
 
   eps_send = 1.0d-6
   press_send = 6.666666666666667d-7
@@ -124,6 +128,9 @@ subroutine c2p2cM(CCTK_ARGUMENTS)
   write(*,*) '   sy  : ',sy_send
   write(*,*) '   sz  : ',sz_send
   write(*,*) '   tau : ',tau_send
+  write(*,*) '   Bconsx  : ',bconsx_send
+  write(*,*) '   Bconsy  : ',bconsy_send
+  write(*,*) '   Bconsz  : ',bconsz_send
   write(*,*) '   eps : ',eps_send
   write(*,*) '   W   : ',w_lorentz_send
   write(*,*) '   Bvecx  : ',bvcx_send
@@ -132,11 +139,10 @@ subroutine c2p2cM(CCTK_ARGUMENTS)
   
   write(*,*) 'C2P2CM test: getting the associated primitive variables.'
   call Con2PrimGenM(GRHydro_eos_handle,local_gam(1),dens_send,sx_send,sy_send,sz_send, &
-       tau_send,rho_send,velx_send,vely_send,velz_send, &
-       eps_send,press_send,w_lorentz_send, &
+       tau_send,bconsx_send,bconsy_send,bconsz_send,rho_send,velx_send,vely_send,velz_send, &
+       eps_send,press_send,bvcx_send,bvcy_send,bvcz_send,b2_send,w_lorentz_send, &
        gxx_send,gxy_send,gxz_send,gyy_send,gyz_send,gzz_send,&
        uxx,uxy,uxz,uyy,uyz,uzz,det,&
-       bvcx_send,bvcy_send,bvcz_send,b2_send,&
        epsnegative,C2P_failed)
   
   write(*,*) 'C2P2CM test: the primitive variables are'
@@ -155,8 +161,8 @@ subroutine c2p2cM(CCTK_ARGUMENTS)
   
   write(*,*) 'C2P2CM test: converting back to conserved variables.'
   call Prim2ConGenM(GRHydro_eos_handle,gxx_send, gxy_send, gxz_send, gyy_send, gyz_send, gzz_send, det, &
-       dens_send, sx_send, sy_send, sz_send, tau_send, bvcx_send, bvcy_send, bvcz_send, rho_send, &
-       velx_send, vely_send, velz_send, eps_send, press_send, w_lorentz_send) 
+       dens_send, sx_send, sy_send, sz_send, tau_send, bconsx_send, bconsy_send, bconsz_send, rho_send, &
+       velx_send, vely_send, velz_send, eps_send, press_send, bvcx_send, bvcy_send, bvcz_send, w_lorentz_send) 
   
   write(*,*) 'C2P2CM test: the conserved variables are'
   write(*,*) '   conservative variables: '
@@ -165,6 +171,9 @@ subroutine c2p2cM(CCTK_ARGUMENTS)
   write(*,*) '   sy  : ',sy_send
   write(*,*) '   sz  : ',sz_send
   write(*,*) '   tau : ',tau_send
+  write(*,*) '   Bconsx  : ',bconsx_send
+  write(*,*) '   Bconsy  : ',bconsy_send
+  write(*,*) '   Bconsz  : ',bconsz_send
   write(*,*) '   eps : ',eps_send
   write(*,*) '   W   : ',w_lorentz_send
   write(*,*) '   Bvecx  : ',bvcx_send
diff --git a/src/GRHydro_P2C2PM.F90 b/src/GRHydro_P2C2PM.F90
index 62c998b..15cacf0 100644
--- a/src/GRHydro_P2C2PM.F90
+++ b/src/GRHydro_P2C2PM.F90
@@ -40,6 +40,7 @@ subroutine p2c2pm(CCTK_ARGUMENTS)
   CCTK_REAL :: uxx,uxy,uxz,uyy,uyz,uzz
   CCTK_REAL :: gxx_send,gxy_send,gxz_send,gyy_send,gyz_send,gzz_send
   CCTK_REAL :: dens_send,sx_send,sy_send,sz_send,tau_send
+  CCTK_REAL :: bconsx_send,bconsy_send,bconsz_send
   CCTK_REAL :: rho_send(1),velx_send,vely_send,velz_send,eps_send(1)
   CCTK_REAL :: press_send(1),w_lorentz_send,x_send,y_send,z_send,r_send
   CCTK_REAL :: bvcx_send,bvcy_send,bvcz_send,b2_send
@@ -86,6 +87,9 @@ subroutine p2c2pm(CCTK_ARGUMENTS)
   bvcx_send   = Bx_init
   bvcy_send   = By_init
   bvcz_send   = Bz_init
+  bconsx_send   = Bx_init
+  bconsy_send   = By_init
+  bconsz_send   = Bz_init
 
   w_lorentz_send = 1.d0/sqrt(1.0d0-velx_send*velx_send-vely_send*vely_send-velz_send*velz_send)
 
@@ -127,8 +131,9 @@ subroutine p2c2pm(CCTK_ARGUMENTS)
   call Prim2ConGenM(GRHydro_eos_handle,&
        gxx_send,gxy_send,gxz_send,gyy_send,gyz_send,gzz_send,det, &
        dens_send,sx_send,sy_send,sz_send,tau_send, &
-       bvcx_send,bvcy_send,bvcz_send,rho_send(1), &
-       velx_send,vely_send,velz_send,eps_send(1),press_send(1),w_lorentz_send) 
+       bconsx_send,bcibsy_send,bconsz_send, &
+       rho_send(1),velx_send,vely_send,velz_send,eps_send(1),press_send(1), &
+       bvcx_send,bvcy_send,bvcz_send,w_lorentz_send) 
 
   write(*,*) 'P2C2PM test: initial values.'
   write(*,*) '   conservative variables: '
@@ -137,6 +142,9 @@ subroutine p2c2pm(CCTK_ARGUMENTS)
   write(*,*) '   sy  : ',sy_send
   write(*,*) '   sz  : ',sz_send
   write(*,*) '   tau : ',tau_send
+  write(*,*) '   Bconsx  : ',bconsx_send
+  write(*,*) '   Bconsy  : ',bconsy_send
+  write(*,*) '   Bconsz  : ',bconsz_send
   write(*,*) '   eps : ',eps_send
   write(*,*) '   W   : ',w_lorentz_send
   write(*,*) '   Bvecx  : ',bvcx_send
@@ -146,11 +154,13 @@ subroutine p2c2pm(CCTK_ARGUMENTS)
   write(*,*) 'P2C2PM test: getting the associated primitive variables.'
   
   call Con2PrimGenM(GRHydro_eos_handle,local_gam(1),dens_send,sx_send,sy_send,sz_send, &
-       tau_send,rho_send(1),velx_send,vely_send,velz_send, &
-       eps_send(1),press_send(1),w_lorentz_send, &
+       tau_send,bconsx_send,bconsy_send,bconsz_send,&
+       rho_send(1),velx_send,vely_send,velz_send, &
+       eps_send(1),press_send(1), &      
+       bvcx_send,bvcy_send,bvcz_send,b2_send,&
+       w_lorentz_send, &
        gxx_send,gxy_send,gxz_send,gyy_send,gyz_send,gzz_send,&
        uxx,uxy,uxz,uyy,uyz,uzz,det,&
-       bvcx_send,bvcy_send,bvcz_send,b2_send,&
        epsnegative,C2P_failed)
   
   write(*,*) 'P2C2PM test: the primitive variables are'
diff --git a/src/GRHydro_P2C2PM_polytype.F90 b/src/GRHydro_P2C2PM_polytype.F90
index a8d0e77..eec4d13 100644
--- a/src/GRHydro_P2C2PM_polytype.F90
+++ b/src/GRHydro_P2C2PM_polytype.F90
@@ -39,6 +39,7 @@ subroutine p2c2pm_polytype(CCTK_ARGUMENTS)
   CCTK_REAL uxx,uxy,uxz,uyy,uyz,uzz
   CCTK_REAL gxx_send,gxy_send,gxz_send,gyy_send,gyz_send,gzz_send
   CCTK_REAL dens_send,sx_send,sy_send,sz_send,tau_send
+  CCTK_REAL bconsx_send,bconsy_send,bconsz_send
   CCTK_REAL rho_send,velx_send,vely_send,velz_send,eps_send
   CCTK_REAL press_send,w_lorentz_send,x_send,y_send,z_send,r_send
   CCTK_REAL bvcx_send,bvcy_send,bvcz_send,b2_send
@@ -93,6 +94,9 @@ subroutine p2c2pm_polytype(CCTK_ARGUMENTS)
   bvcx_send   = Bx_init
   bvcy_send   = By_init
   bvcz_send   = Bz_init
+  bconsx_send   = Bx_init
+  bconsy_send   = By_init
+  bconsz_send   = Bz_init
 
   w_lorentz_send = 1.d0/sqrt(1.0d0-velx_send*velx_send-vely_send*vely_send-velz_send*velz_send)
 
@@ -133,8 +137,8 @@ subroutine p2c2pm_polytype(CCTK_ARGUMENTS)
   
   write(*,*) 'P2C2PM test: converting back to conserved variables.'
   call prim2conpolytypeM(GRHydro_polytrope_handle,gxx_send, gxy_send, gxz_send, gyy_send, gyz_send, gzz_send, det, &
-       dens_send, sx_send, sy_send, sz_send, tau_send, bvcx_send, bvcy_send, bvcz_send, rho_send, &
-       velx_send, vely_send, velz_send, eps_send, press_send, w_lorentz_send) 
+       dens_send, sx_send, sy_send, sz_send, tau_send, bconsx_send, bconsy_send, bconsz_send, rho_send, &
+       velx_send, vely_send, velz_send, eps_send, press_send, bvcx_send, bvcy_send, bvcz_send, w_lorentz_send) 
 
   write(*,*) 'P2C2PM test: initial values.'
   write(*,*) '   conservative variables: '
@@ -143,6 +147,9 @@ subroutine p2c2pm_polytype(CCTK_ARGUMENTS)
   write(*,*) '   sy  : ',sy_send
   write(*,*) '   sz  : ',sz_send
   write(*,*) '   tau : ',tau_send
+  write(*,*) '   Bconsx  : ',bconsx_send
+  write(*,*) '   Bconsy  : ',bconsy_send
+  write(*,*) '   Bconsz  : ',bconsz_send
   write(*,*) '   eps : ',eps_send
   write(*,*) '   W   : ',w_lorentz_send
   write(*,*) '   Bvecx  : ',bvcx_send
@@ -154,11 +161,10 @@ subroutine p2c2pm_polytype(CCTK_ARGUMENTS)
 
   write(*,*) 'P2C2PM test: getting the associated primitive variables.'
   call GRHydro_Con2PrimM_Polytype_pt(GRHydro_polytrope_handle,local_gam,dens_send,sx_send,sy_send,sz_send, &
-       sc_send,rho_send,velx_send,vely_send,velz_send, &
-       eps_send,press_send,w_lorentz_send, &
+       sc_send,bconsx_send,bconsy_send,bconsz_send,rho_send,velx_send,vely_send,velz_send, &
+       eps_send,press_send,bvcx_send,bvcy_send,bvcz_send,b2_send,w_lorentz_send, &
        gxx_send,gxy_send,gxz_send,gyy_send,gyz_send,gzz_send,&
        uxx,uxy,uxz,uyy,uyz,uzz,det,&
-       bvcx_send,bvcy_send,bvcz_send,b2_send,&
        epsnegative,C2P_failed)
   
   write(*,*) 'P2C2PM test: the primitive variables are'