Skip to content

Commit

Permalink
Update oderhs.m
Browse files Browse the repository at this point in the history 
Adding the coding for external power source
  • Loading branch information
@umansky
umansky authored Oct 19, 2023
1 parent ec603c7 commit 2d04dec
Showing 1 changed file with 9 additions and 0 deletions.
9 changes: 9 additions & 0 deletions bbb/oderhs.m
View file
Original file line number Diff line number Diff line change
Expand Up @@ -4732,6 +4732,15 @@ cc elseif (ishosor .ne. 0)
if (istimingon .eq. 1) call timimpfj (tsimp, xc)
endif #loop for isimpon==2
ccc MVU - 18-oct-2023, adding external heating terms
do iy=0,ny+1
do ix=0,nx+1
resee(ix,iy) = resee(ix,iy) + pwre_use(ix,iy)*vol(ix,iy)
resei(ix,iy) = resei(ix,iy) + pwri_use(ix,iy)*vol(ix,iy)
enddo
enddo
ccc MVU
* -- joule heating --
Expand Down

2 comments on commit 2d04dec

@holm10
Copy link
Collaborator

@holm10 holm10 commented on 2d04dec Oct 19, 2023

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

@umansky just saw your email and had a look at this coding. I think this coding might be obsolete, and that columetric sources can be user-defined in the current version of UEDGE.

Try the following:

  • bbb.isvolsorext=1 # Tells UEDGE to use user-defined volumetric sources
  • bbb.pwrsore=... # Sets the volumetric electron energy source
  • bbb.pwrsori=... # Sets the volumetric ion energy source
  • bbb.volpsor=... # Sets the volumetric ion particle source
  • etc

This does not work if one needs to specify two sources separately, e.g. a Maxwellian heat source in addition to a user-specified source. However, in that case the net source could be calculated and used as user input.

@umansky
Copy link
Collaborator Author

@umansky umansky commented on 2d04dec Oct 19, 2023 via email

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Please sign in to comment.