istp/gray
2
0
Fork 0
Code Issues Pull Requests Packages Projects Releases 1 Wiki Activity

Fix label for Icd units in log output

Browse Source
This commit is contained in:
Lorenzo Figini 2023-10-04 16:56:05 +02:00
parent 526195ca94
commit cf8a37fec9
1 changed files with 3 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
src/gray_core.f90
View File

@ -553,7 +553,7 @@ contains
write(msg, '(3x,a,g0.3," MW")') 'absoption: P=', pabs_beam write(msg, '(3x,a,g0.3," MW")') 'absoption: P=', pabs_beam
call log_info(msg, mod='gray_core', proc='gray_main') call log_info(msg, mod='gray_core', proc='gray_main')
write(msg, '(3x,a,g0.3," MW")') 'current drive: I=', icd_beam*1.0e3_wp_ write(msg, '(3x,a,g0.3," kA")') 'current drive: I=', icd_beam*1.0e3_wp_
call log_info(msg, mod='gray_core', proc='gray_main') call log_info(msg, mod='gray_core', proc='gray_main')
if(ip < params%raytracing%ipass) then if(ip < params%raytracing%ipass) then
@ -1373,7 +1373,7 @@ contains
dbvdc(3,1) = dbvcdc(3,1) dbvdc(3,1) = dbvcdc(3,1)
dbvdc(1,2) = -bv(2) dbvdc(1,2) = -bv(2)
dbvdc(2,2) = bv(1) dbvdc(2,2) = bv(1)
dbvdc(3,2) = dbvcdc(3,2) dbvdc(3,2) = dbvcdc(3,2) ! = 0
dbvdc(1,3) = dbvcdc(1,3)*csphi - dbvcdc(2,3)*snphi dbvdc(1,3) = dbvcdc(1,3)*csphi - dbvcdc(2,3)*snphi
dbvdc(2,3) = dbvcdc(1,3)*snphi + dbvcdc(2,3)*csphi dbvdc(2,3) = dbvcdc(1,3)*snphi + dbvcdc(2,3)*csphi
dbvdc(3,3) = dbvcdc(3,3) dbvdc(3,3) = dbvcdc(3,3)
@ -1391,6 +1391,7 @@ contains
! B magnitude and derivatives ! B magnitude and derivatives
btot = norm2(bv) btot = norm2(bv)
! dbtot(i) = d |B| / dxvcart(i)
dbtot = matmul(bv, dbv)/btot dbtot = matmul(bv, dbv)/btot
yg = btot/Bres yg = btot/Bres