c print*,"ETOTAL Processor",fg_rank," absolute rank",myrank,
c & " nfgtasks",nfgtasks
call flush(iout)
-#ifndef DFA
- edfadis=0.0d0
- edfator=0.0d0
- edfanei=0.0d0
- edfabet=0.0d0
-#endif
if (nfgtasks.gt.1) then
#ifdef MPI
time00=MPI_Wtime()
c call int_from_cart1(.false.)
c endif
#endif
+#ifndef DFA
+ edfadis=0.0d0
+ edfator=0.0d0
+ edfanei=0.0d0
+ edfabet=0.0d0
+#endif
#ifdef TIMING
#ifdef MPI
time00=MPI_Wtime()
c sgtheta=-gtheta(k)*theta_diff(k)*sigma_theta(k,i)*waga_theta ! right functional form?
sum_sgtheta=sum_sgtheta+sgtheta ! cum variable
enddo
-c grad_theta3=sum_sgtheta/sum_gtheta 1/*theta(i)? s. line below
-c grad_theta3=sum_sgtheta/sum_gtheta
-c
c Final value of gradient using same var as in Econstr_back
gloc(nphi+i-2,icg)=gloc(nphi+i-2,icg)
& +sum_sgtheta/sum_gtheta*waga_theta