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

BTW, the setups I used for mhd jet is attached.
See if you can see how I can get across this problem I
described (early abort and azimuthal magnetic field)

Thanks!

Bst,
Victor
--- Victor Shum <vshum12@pacbell.net> wrote:

> Date: Tue, 1 Feb 2005 15:03:16 -0800 (PST)
> From: Victor Shum <vshum12@pacbell.net>
> Subject: Re: [FLASH-USERS] which program/data input
> file to specifying y coordinates via the nyb
> parameter?
> To: Michal Hanasz
> <Michal.Hanasz@astri.uni.torun.pl>,
> flash <flash-users@flash.uchicago.edu>
> CC: Michal Hanasz <Michal.Hanasz@astri.uni.torun.pl>
>
> 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
> >
>
>

> ATTACHMENT part 2 video/avi name=jetnozzle.avi