coskt(k)=dcos(k*theti2)
sinkt(k)=dsin(k*theti2)
enddo
- if (i.eq.3) then
- phii=0.0d0
- ityp1=nthetyp+1
- do k=1,nsingle
- cosph1(k)=0.0d0
- sinph1(k)=0.0d0
- enddo
- else
+cu if (i.eq.3) then
+cu phii=0.0d0
+cu ityp1=nthetyp+1
+cu do k=1,nsingle
+cu cosph1(k)=0.0d0
+cu sinph1(k)=0.0d0
+cu enddo
+cu else
if (i.gt.3 .and. itype(i-3).ne.ntyp1) then
#ifdef OSF
phii=phi(i)
sinph1(k)=0.0d0
enddo
endif
- endif
+cu endif
if (i.lt.nres .and. itype(i+1).ne.ntyp1) then
#ifdef OSF
phii1=phi(i+1)