distriREAC corrected

octavioroncero · Jun 28, 2021 · 3b35cca · 3b35cca
1 parent 1b016c1
commit 3b35cca
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 3 deletions.
diff --git a/SRC/AUXILIARy-CODES/distriREACwvp.f b/SRC/AUXILIARy-CODES/distriREACwvp.f
@@ -321,7 +321,7 @@ program distriREACwvp ! version v6 in progress
  do ir2=1,npunt
  r=rmis2+dble(ir2-1)*ahgauss
  arg=r*pini
- erot=hbr*hbr*pepe*0.5d0/(r*r*xmasa)
+ erot=hbr*hbr*pepe*0.5d0/(r2col*r2col*xmasa)
  if(ekinini.gt.erot)then
  CALL BESPH2(F,DF,G,DG,PEPE,ARG,KEY,0)
  zexpo=dcmplx(-g,f)
@@ -343,7 +343,8 @@ program distriREACwvp ! version v6 in progress
  else
  pfin=0.d0
  endif
- if(paqini.lt.1.d-30)then
+
+ if(paqini.lt.1.d-15)then
  S2prodfac(iv,j,iele)=0.d0
  else
  S2prodfac(iv,j,iele)=
@@ -407,7 +408,7 @@ program distriREACwvp ! version v6 in progress
  jjjj0=max0(iabs(iom),j00)
 
  do j=jjjj0,j1
- if(noBCstates(j,ielec).gt.0)then
+ if(noBCstates(j,ielec).ge.0)then
  zzz=zS2prod(iv,j,ielec,iom)/2.d0/pi
  else
  zzz=dcmplx(0.d0,0.d0)

diff --git a/SRC/photo/.cheby-spectra.f.swp b/SRC/photo/.cheby-spectra.f.swp