! if(.not.allocated(dc_norm2)) allocate(dc_norm2(3,0:nres2+2)) !(3,0:maxres2)
if(.not.allocated(dc_norm2)) then
allocate(dc_norm2(3,0:nres2+2)) !(3,0:maxres2)
- do i=0,nres2+2
- dc_norm2(1,i)=0.d0
- dc_norm2(2,i)=0.d0
- dc_norm2(3,i)=0.d0
- enddo
+ dc_norm2(:,:)=0.d0
endif
!
!el if(.not.allocated(dc_norm))
!elwrite(iout,*) "jestem w alloc geo 1"
if(.not.allocated(dc_norm)) then
allocate(dc_norm(3,0:nres2+2)) !(3,0:maxres2)
- do i=0,nres2+2
- dc_norm(1,i)=0.d0
- dc_norm(2,i)=0.d0
- dc_norm(3,i)=0.d0
- enddo
+ dc_norm(:,:)=0.d0
endif
!elwrite(iout,*) "jestem w alloc geo 1"
allocate(xloc(3,nres),xrot(3,nres))
!elwrite(iout,*) "jestem w alloc geo 1"
- do i=1,nres
- do j=1,3
- xloc(j,i)=0.0D0
- enddo
- enddo
+ xloc(:,:)=0.0D0
!elwrite(iout,*) "jestem w alloc geo 1"
allocate(dc_work(6*nres)) !(MAXRES6) maxres6=6*maxres
! common /rotmat/
#if defined(WHAM_RUN) || defined(CLUSTER)
allocate(vbld(2*nres))
- do i=1,2*nres
- vbld(i)=0.d0
- enddo
+ vbld(:)=0.d0
allocate(vbld_inv(2*nres))
- do i=1,2*nres
- vbld_inv(i)=0.d0
- enddo
+ vbld_inv(:)=0.d0
#endif
return