MPI shield

[unres4.git] / source / unres / control.F90

diff --git a/source/unres/control.F90 b/source/unres/control.F90
index f7cba30..4c413ec 100644 (file)
--- a/source/unres/control.F90
+++ b/source/unres/control.F90
@@ -512,7 +512,7 @@
             ind_eleint_old_nucl,ind_eleint_nucl,nele_int_tot_vdw_nucl,&
             my_ele_inds_vdw_nucl,my_ele_inde_vdw_nucl,ind_eleint_vdw_nucl,&
             ind_eleint_vdw_old_nucl,nscp_int_tot_nucl,my_scp_inds_nucl,&
-            my_scp_inde_nucl,ind_scpint_nucl,ind_scpint_old_nucl
+            my_scp_inde_nucl,ind_scpint_nucl,ind_scpint_old_nucl,impishi
 !      integer,dimension(5) :: nct_molec,nnt_molec
 !el      allocate(itask_cont_from(0:nfgtasks-1)) !(0:max_fg_procs-1)
 !el      allocate(itask_cont_to(0:nfgtasks-1)) !(0:max_fg_procs-1)
@@ -1383,6 +1383,13 @@
         call MPI_Type_contiguous(18,MPI_DOUBLE_PRECISION,MPI_UYZGRAD,&
           IERROR)
         call MPI_Type_commit(MPI_UYZGRAD,IERROR)
+        call MPI_Type_contiguous(maxcontsshi,MPI_INTEGER,MPI_I50,IERROR)
+        call MPI_Type_commit(MPI_I50,IERROR)
+         impishi=maxcontsshi*3
+        call MPI_Type_contiguous(impishi,MPI_DOUBLE_PRECISION, &
+        MPI_SHI,IERROR)
+        call MPI_Type_commit(MPI_SHI,IERROR)
+
         call MPI_Type_contiguous(2,MPI_DOUBLE_PRECISION,MPI_MU,IERROR)
         call MPI_Type_commit(MPI_MU,IERROR)
         call MPI_Type_contiguous(4,MPI_DOUBLE_PRECISION,MPI_MAT1,IERROR)