gkowal/amun-code
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Merge branch 'master' into reconnection

Browse Source
This commit is contained in:
Grzegorz Kowal 2017-03-29 10:44:06 -03:00
commit f1ef16be8a
1 changed files with 7 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
src/sources.F90
View File

@ -582,11 +582,17 @@ module sources
call curl(dh(:), pdata%q(ibx:ibz,1:im,1:jm,1:km) & call curl(dh(:), pdata%q(ibx:ibz,1:im,1:jm,1:km) &
, tmp(inz,inx:inz,1:im,1:jm,1:km)) , tmp(inz,inx:inz,1:im,1:jm,1:km))
! calculate J²
!
db(1:im,1:jm,1:km) = tmp(inz,inx,1:im,1:jm,1:km)**2 &
+ tmp(inz,iny,1:im,1:jm,1:km)**2 &
+ tmp(inz,inz,1:im,1:jm,1:km)**2
! add the second resistive source term to the energy equation, i.e. ! add the second resistive source term to the energy equation, i.e.
! d/dt E + ∇.F = η J² ! d/dt E + ∇.F = η J²
! !
du(ien,1:im,1:jm,1:km) = du(ien,1:im,1:jm,1:km) & du(ien,1:im,1:jm,1:km) = du(ien,1:im,1:jm,1:km) &
+ resistivity * sum(tmp(inz,inx:inz,1:im,1:jm,1:km)**2, 2) + resistivity * db(1:im,1:jm,1:km)
end if ! energy equation present end if ! energy equation present