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

corrected current sign assignment, more output in units 6,7

Browse Source
This commit is contained in:
Lorenzo Figini 2014-02-04 12:38:10 +00:00
parent 1a60867b5b
commit 0d971665aa
1 changed files with 7 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
src/gray.f
View File

@ -514,7 +514,7 @@ c
write(12,*) ' #i sst psi w1 w2' write(12,*) ' #i sst psi w1 w2'
write(7,*)'#Icd Pa Jphip dPdVp '// write(7,*)'#Icd Pa Jphip dPdVp '//
.'rhotj rhotjava rhotp rhotpav drhotjava drhotpav '// .'rhotj rhotjava rhotp rhotpav drhotjava drhotpav '//
.'ratjbmx stmx psipol chipol index_rt Jphimx dPdVmx' .'ratjbmx stmx psipol chipol index_rt Jphimx dPdVmx drhotj drhotp'
write(48,*) '#psi rhot Jphi Jcdb dPdV Icdins Pins P% index_rt' write(48,*) '#psi rhot Jphi Jcdb dPdV Icdins Pins P% index_rt'
else else
@ -1190,7 +1190,7 @@ c icocos mod 10 = 2,3,6,7: psi decreasing with CCW Ip
c c
c add check for Ip/psi and B0/Fpol sign consistency? c add check for Ip/psi and B0/Fpol sign consistency?
c c
current=sign(current,psiax-psiedge) current=sign(current,psiaxis-psiedge)
btrcen=sign(btrcen,fpol(nr)) btrcen=sign(btrcen,fpol(nr))
c c
@ -6371,16 +6371,16 @@ c dPdV [MW/m^3], Jcd [MA/m^2]
write(6,*)' ' write(6,*)' '
write(6,*)'#beta0 alpha0 Icd Pa Jphip dPdVp '// write(6,*)'#beta0 alpha0 Icd Pa Jphip dPdVp '//
.'rhotj rhotjava rhotp rhotpav drhotjava drhotpav '// .'rhotj rhotjava rhotp rhotpav drhotjava drhotpav '//
.'ratjbmx stmx psipol chipol index_rt Jphimx dPdVmx' .'ratjbmx stmx psipol chipol index_rt Jphimx dPdVmx drhotj drhotp'
write(6,99) beta0,alpha0,currtka,pabstot,ajphip,dpdvp, write(6,99) beta0,alpha0,currtka,pabstot,ajphip,dpdvp,
. rhotjfi,rhotjava,rhotp,rhotpav, . rhotjfi,rhotjava,rhotp,rhotpav,
. drhotjava,drhotpav, . drhotjava,drhotpav,stmx,psipol,chipol,real(index_rt),
. stmx,psipol,chipol,real(index_rt),ajmxfi,dpdvmx . ajmxfi,dpdvmx,drhotjfi,drhotp
write(7,99) currtka,pabstot,ajphip,dpdvp, write(7,99) currtka,pabstot,ajphip,dpdvp,
. rhotjfi,rhotjava,rhotp,rhotpav, . rhotjfi,rhotjava,rhotp,rhotpav,
. drhotjava,drhotpav,ratjbmx, . drhotjava,drhotpav,ratjbmx,stmx,psipol,chipol,real(index_rt),
. stmx,psipol,chipol,real(index_rt),ajmxfi,dpdvmx . ajmxfi,dpdvmx,drhotjfi,drhotp
do i=1,nd do i=1,nd
if (ipec.eq.0) then if (ipec.eq.0) then