Update 5/24

MHD proposal figs

http://www.pas.rochester.edu/~adebrech/HD209458b/MHD/MHD_proposal_fig_0_nopress.png

http://www.pas.rochester.edu/~adebrech/HD209458b/MHD/MHD_proposal_fig_0_press.png

http://www.pas.rochester.edu/~adebrech/HD209458b/MHD/MHD_proposal_fig_90.png

Also, a test of the streamlines. I like the LICs better, since the streamlines have a habit of not ending.

http://www.pas.rochester.edu/~adebrech/HD209458b/MHD/streamline_test.png

​Radiation Pressure Paper

Found issue with postprocessing. Reprocessing synthetic observations. Here's a first look at the difference. The difference is significant, but doesn't change our conclusions.

http://www.pas.rochester.edu/~adebrech/HD209458b/HD209458b_obs_comparison_three_zoom_diff.png

Charge Exchange

Current frame

http://www.pas.rochester.edu/~adebrech/ChargeExchange/HD209458b_charge_exchange_10003.png

We also need to decide on the high-med-low for the stellar wind. Here are the current parameters (taken from John's simulation, except the mass loss rate, which is the solar mass loss rate - which corresponds to a 4x higher wind density than their strongest wind):

Parameter Value
cm
K
cm/s
g/s

AMR line transfer

Timing run is next up (for comparison with non-AMR and new optimizations), along with work on optimizations.

Also working on a way to output propagation of rays across grid for a step (just for a nice visualization):

Need: ray position, synchronized time

Ray position is easy.

Synchronized time pseudocode: 

call MPI_barrier     ! To synchronize time - unless MPI_wtime_is_global, in which case all processors are synchronized already and this becomes much easier
initTime = MPI_wtime
call PropagateRay

subroutine PropagateRay

time = MPI_wtime - initTime
MPI_file_write(time,position)    ! Or might be easier to save an array of time, position and write it out after all rays are processed, in processor order (so we don't have to calculate the buffer)

Comments

No comments.