Commit 67e7fba0 authored by Matthias Redies's avatar Matthias Redies

replace implicit allocation with explicit

parent 0db341f5
......@@ -206,6 +206,10 @@ CONTAINS
PotDen3%iter = PotDen1%iter
PotDen3%potdenType = PotDen1%potdenType
! implicit allocation would break the bounds staring at 0
if(.not. allocated(PotDen3%mt)) allocate(PotDen3%mt, mold=PotDen1%mt)
PotDen3%mt = PotDen1%mt + PotDen2%mt
PotDen3%pw = PotDen1%pw + PotDen2%pw
PotDen3%vacz = PotDen1%vacz + PotDen2%vacz
......@@ -224,6 +228,10 @@ CONTAINS
PotDen3%iter = PotDen1%iter
PotDen3%potdenType = PotDen1%potdenType
! implicit allocation would break the bounds staring at 0
if(.not. allocated(PotDen3%mt)) allocate(PotDen3%mt, mold=PotDen1%mt)
PotDen3%mt = PotDen1%mt - PotDen2%mt
PotDen3%pw = PotDen1%pw - PotDen2%pw
PotDen3%vacz = PotDen1%vacz - PotDen2%vacz
......@@ -242,6 +250,10 @@ CONTAINS
PotDenCopy%iter = PotDen%iter
PotDenCopy%potdenType = PotDen%potdenType
! implicit allocation would break the bounds staring at 0
if(.not. allocated(PotDenCopy%mt)) allocate(PotDenCopy%mt, mold=PotDen%mt)
PotDenCopy%mt = PotDen%mt
PotDenCopy%pw = PotDen%pw
PotDenCopy%vacz = PotDen%vacz
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment