Release the memory allocated in the mhd module

xcompact3d · Apr 22, 2024 · 380017c · 380017c
1 parent a1ec640
commit 380017c
Showing 1 changed file with 9 additions and 0 deletions.
diff --git a/src/mhd.f90 b/src/mhd.f90
@@ -87,8 +87,17 @@ subroutine mhd_fin()
 
      implicit none
 
+     ! Close the IO unit opened in the subroutine mhd_sta
      if (nrank==0) close(mhd_iounit)
 
+     ! Release memory
+     if (allocated(Bmean)) deallocate(Bmean)
+     if (allocated(Bm)) deallocate(Bm)
+     if (allocated(magelf)) deallocate(magelf)
+     if (allocated(Je)) deallocate(Je)
+     if (allocated(dBm)) deallocate(dBm)
+     if (allocated(elcpot)) deallocate(elcpot)
+
   end subroutine mhd_fin
   !
   subroutine int_time_magnet