diff options
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 |