Update 11/2

  • Finished GRFP application on Friday. Waiting on two letters of recommendation (as of noon today).
  • WASP-12b: According to matlab script, no bow shock. Plot for ram pressures is

Initialized stellar wind out too far, so that planet was created inside:

​WASP-12b global simulation

  • MHD simulations of HD209458b: attempting to control refinement. On most recent attempt, segfaulted:

​Error message

Relevant code snippet: 

SUBROUTINE ProblemSetErrFlag(pos, Info)
    TYPE(InfoDef) :: Info
    REAL(KIND=qPREC), DIMENSION(:), INTENT(IN) :: pos
    REAL(KIND=qPREC) :: x(3), r
    x=pos-PlanetPos
    r=sqrt(sum(x**2))
!   Force lowest resolution outside of 5*Rp
    IF (r > 5*planetRadius) THEN
        Info%ErrFlag = 0
    END IF
  END SUBROUTINE ProblemSetErrFlag

Next thought is to write a separate subroutine that can be called in ProblemSetErrFlag, similar to the subroutines being called in the before step - but is this too simplistic?

  • Also time to register for next semester - high-energy astrophysics or compressible fluids?

Attachments (3)

Comments

No comments.