Re: [FLASH-USERS] which program/data input file to specifying y coordinates via the nyb parameter?

Hi Michal and flash users,
 
I managed to run Michal's setup for the jet nozzle and built a movie out of it.
Michal, pls take a look and see if it is correct.
 
I've also tried adapting the mhd modules to the jet with different require statements set in the Config file (according to the rt_mhd example) but found out that even with refinemax=1, the program aborted itself on the third checkpoint writeout.
Any suggestion?
 
Also, in the rt_mhd or kh_mhd examples the magnetic field specified in flash.par are Mx and My.
If I want to do Mtheta (azimuthal direction for the nozzle), should I do solndata (Mtheta)=Mx^2 + My^2 or stick with solndata(My) = function of x or Mx.
 
Many Thanks!
 
Bst,
Victor

Michal Hanasz <Michal.Hanasz@astri.uni.torun.pl> wrote:
Dear Victor,

Victor Shum wrote:
> I saw that you used formulations (in terms of y, rj, and rho_a) for the
> idens, ipres variables in the init_block.F90 and set the rj, rho_a
> params in the flash.par as constants.
> However, I still have some questions.
> For the y variable you used in init_block.F90, did you specify the
> boundary values wrt x in one of your setup files (I don't see that
> anywhere in your setups) or somewhere else (like in a dBase.F90 program)?

The inflow boundary conditions are specified in user_bnd.F90.
I did not specify boundary values with respect to x, but it shouldn't be a
problem to do it by modifying init_block.F90, parts following the comment line

"! -Y user-defined boundary"

in a way analogous to the part of the code following the comment line

"! -X user-defined boundary"

>
> If I have a set of values for the y boundaries at the top and bottom,
> can I then use that as an array of values in the flash.par (like
> replacing your rj=3 with rj = [ 1,2, 6, 5] ) ?

"rj" is the jet radius - see description of parameters in the init_block.F90
Then, I don't understand your question.

Best wishes,

Michal

>
> Many Thanks!
>
> Bst,
> Victor
>