Merge pull request #934 from SALMON-TDDFT/bugfix

Bugfix to release version 2.
SALMON-TDDFT · Jul 21, 2020 · 08c1ea9 · 08c1ea9
2 parents 864f391 + b79cee3
commit 08c1ea9
Show file tree

Hide file tree

Showing 2 changed files with 5 additions and 5 deletions.
diff --git a/src/common/structures.f90 b/src/common/structures.f90
@@ -539,7 +539,7 @@ subroutine allocate_vector(rg,field)
     do iz=rg%is(3),rg%ie(3)
     do iy=rg%is(2),rg%ie(2)
     do ix=rg%is(1),rg%ie(1)
-      field%v = 0d0
+      field%v(:,ix,iy,iz) = 0d0
     end do
     end do
     end do
@@ -555,7 +555,7 @@ subroutine allocate_vector_with_ovlp(rg,field)
     do iz=rg%is_array(3),rg%ie_array(3)
     do iy=rg%is_array(2),rg%ie_array(2)
     do ix=rg%is_array(1),rg%ie_array(1)
-      field%v = 0d0
+      field%v(:,ix,iy,iz) = 0d0
     end do
     end do
     end do

diff --git a/src/ms/main_ms.f90 b/src/ms/main_ms.f90
@@ -265,8 +265,6 @@ subroutine initialization_ms()
     end do
 
     ! incident field
-    allocate(Ac_inc(1:3, -1:itotNtime+2))
-
     call Weyl_init(fs, fw)
 
     allocate(fs%imedia(fs%mg%is_array(1):fs%mg%ie_array(1), &
@@ -348,7 +346,9 @@ subroutine initialization_ms()
         call comm_sync_all()
     end do
 
-    itt = mit; call incident()
+    itt = mit
+    allocate(Ac_inc(1:3, -1:itotNtime+2))
+    call incident()
 
     ! Experimental implementation
     if (comm_is_root(ms%id_ms_world)) then