+ integer, dimension(:,:),allocatable :: icont
+ real(kind=4),dimension(:),allocatable :: dcont
+ if (.not.allocated(icont)) then
+ allocate(icont(2,100*nres_molec(1)+1))
+ endif
+ if (.not.allocated(dcont)) then
+ allocate(dcont(100*nres_molec(1)+1))
+ endif