From 6fd6126abe84add517345d66d57e530da0866d3f Mon Sep 17 00:00:00 2001 From: jmlloren Date: Tue, 25 Jan 2005 18:11:15 +0000 Subject: [PATCH] Eliminado el bucle para el numero de QD, puesto que la energia no es acumulativa. --- pot_slat.f90 | 18 ++++-------------- 1 file changed, 4 insertions(+), 14 deletions(-) diff --git a/pot_slat.f90 b/pot_slat.f90 index f1d39f5..3f27a98 100755 --- a/pot_slat.f90 +++ b/pot_slat.f90 @@ -47,9 +47,9 @@ XD: DO I_X=1,XDim X_VEC=(/X,Y,Z/) -N3: DO I_N3=NMin_Z,NMax_Z -N2: DO I_N2=NMin_Y,NMax_Y -N1: DO I_N1=NMin_X,NMax_X + ENERGY=0.E0 + + I_N1=0; I_N2=0; I_N3=0 R_SL=REAL(I_N1)*A1_S+REAL(I_N2)*A2_S+REAL(I_N3)*A3_S @@ -90,10 +90,6 @@ N1: DO I_N1=NMin_X,NMax_X CALL POT_CALC_CYL_WZ(CHI,ENERGY) END IF - END DO N1 - END DO N2 - END DO N3 - EEL(I_X,I_Y,I_Z) = ENERGY(1) EHHUP(I_X,I_Y,I_Z) = ENERGY(2) EHHDW(I_X,I_Y,I_Z) = ENERGY(3) @@ -562,9 +558,7 @@ XD: DO I_X=1,XDim X_VEC=(/X,Y,Z/) -N3: DO I_N3=NMin_Z,NMax_Z -N2: DO I_N2=NMin_Y,NMax_Y -N1: DO I_N1=NMin_X,NMax_X + I_N1=0; I_N2=0; I_N3=0 R_SL=REAL(I_N1)*A1_S+REAL(I_N2)*A2_S+REAL(I_N3)*A3_S @@ -605,10 +599,6 @@ N1: DO I_N1=NMin_X,NMax_X CALL POT_CALC_CYL_ZB(CHI,ENERGY) END IF - END DO N1 - END DO N2 - END DO N3 - EEL(I_X,I_Y,I_Z) = ENERGY(1) EHHUP(I_X,I_Y,I_Z) = ENERGY(2) ELHUP(I_X,I_Y,I_Z) = ENERGY(3) -- 2.11.4.GIT