diff options
author | goodale <goodale@5301f0c2-dbc4-4cee-b2f5-8d7afba4d129> | 2004-05-17 11:51:13 +0000 |
---|---|---|
committer | goodale <goodale@5301f0c2-dbc4-4cee-b2f5-8d7afba4d129> | 2004-05-17 11:51:13 +0000 |
commit | 58d8c0e9ba4217b011ed5fc572973dca842ee5f2 (patch) | |
tree | 8821b09629b3320f2ea5755ee44f7dcc000dbe5f /src/D2_extract.F | |
parent | d18babc7883e3673b4c109c292ce8f9e38ce4f7e (diff) |
Untabified.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/Extract/trunk@96 5301f0c2-dbc4-4cee-b2f5-8d7afba4d129
Diffstat (limited to 'src/D2_extract.F')
-rw-r--r-- | src/D2_extract.F | 24 |
1 files changed, 12 insertions, 12 deletions
diff --git a/src/D2_extract.F b/src/D2_extract.F index 69f65f1..4a21398 100644 --- a/src/D2_extract.F +++ b/src/D2_extract.F @@ -119,8 +119,8 @@ c ------------------ IF (do_ADMmass(1) == 1) THEN c print*,'1' -c Calculate ADM mass using standard formula - DO j = 2, Np+1 +c Calculate ADM mass using standard formula + DO j = 2, Np+1 DO i = 2, Nt+1 temp(i,j) = ADMmass_int1(i-1,j-1)*SIN(theta(i-1)) ENDDO @@ -132,7 +132,7 @@ c Calculate ADM mass using standard formula c print*,'2' c Calculate ADM mass using conformal formula - DO j = 2, Np+1 + DO j = 2, Np+1 DO i = 2, Nt+1 temp(i,j) = ADMmass_int2(i-1,j-1)*SIN(theta(i-1)) ENDDO @@ -146,7 +146,7 @@ c ------------------ IF (do_momentum == 1) THEN c print*,'3' - DO j = 2, Np+1 + DO j = 2, Np+1 DO i = 2, Nt+1 temp(i,j) = momentum_int1(i-1,j-1)*SIN(theta(i-1)) ENDDO @@ -154,7 +154,7 @@ c print*,'3' momentum(1) = sphere_int(temp,igrid,Nt+1,Np+1,Dt,Dp) c print*,'4' - DO j = 2, Np+1 + DO j = 2, Np+1 DO i = 2, Nt+1 temp(i,j) = momentum_int2(i-1,j-1)*SIN(theta(i-1)) ENDDO @@ -162,7 +162,7 @@ c print*,'4' momentum(2) = sphere_int(temp,igrid,Nt+1,Np+1,Dt,Dp) c print*,'5' - DO j = 2, Np+1 + DO j = 2, Np+1 DO i = 2, Nt+1 temp(i,j) = momentum_int3(i-1,j-1)*SIN(theta(i-1)) ENDDO @@ -176,7 +176,7 @@ c -------------- IF (do_spin == 1) THEN c print*,'6' - DO j = 2, Np+1 + DO j = 2, Np+1 DO i = 2, Nt+1 temp(i,j) = spin_int1(i-1,j-1)*SIN(theta(i-1)) ENDDO @@ -184,7 +184,7 @@ c print*,'6' spin(1) = sphere_int(temp,igrid,Nt+1,Np+1,Dt,Dp) c print*,'7' - DO j = 2, Np+1 + DO j = 2, Np+1 DO i = 2, Nt+1 temp(i,j) = spin_int2(i-1,j-1)*SIN(theta(i-1)) ENDDO @@ -192,7 +192,7 @@ c print*,'7' spin(2) = sphere_int(temp,igrid,Nt+1,Np+1,Dt,Dp) c print*,'8' - DO j = 2, Np+1 + DO j = 2, Np+1 DO i = 2, Nt+1 temp(i,j) = spin_int3(i-1,j-1)*SIN(theta(i-1)) ENDDO @@ -557,13 +557,13 @@ c ------------------------------------------------------------------ WRITE(104,*) eta,K WRITE(105,*) eta,c1 WRITE(106,*) eta,c2 - WRITE(107,*) eta,dG + WRITE(107,*) eta,dG WRITE(108,*) eta,dK WRITE(109,*) eta,dc2 WRITE(201,*) eta,mass - WRITE(202,*) eta,rsch + WRITE(202,*) eta,rsch WRITE(203,*) eta,drschdri - WRITE(301,*) eta,Qeven(:,lmax,mmax) + WRITE(301,*) eta,Qeven(:,lmax,mmax) WRITE(302,*) eta,Qodd(:,lmax,mmax) ENDIF |