use mpi_data
! use mpi_
#endif
- use control, only:initialize
+ use control, only:initialize,hpb_partition
!el use io_config, only:parmread
!
#ifndef ISNAN
!#define WHAM_RUN
! call alloc_wham_arrays
!write(iout,*) "after alloc wham"
- write (*,"(A,/,A,/,A,/,A)") "Last Christmas", &
- "I gave you my heart","But the very next day", &
- "You gave it away"
#ifdef MPI
call MPI_Init( IERROR )
call MPI_Comm_rank( MPI_COMM_WORLD, me, IERROR )
!write(iout,*)"before proc_cont, define frag"
call proc_cont
call fragment_list
- if (constr_dist.gt.0) call read_dist_constr
+ endif
+ if (constr_dist.gt.0) then
+ call read_dist_constr
+ call hpb_partition
endif
write (iout,*) "Begin read_database"
call flush(iout)