Initial commit from WindTunnel_heat
 Replaced Grid_markRefineDerefine.F90 by FLASH4.2.1 version.
 Added Magnetic field.
  Compilable, but still problmatic (sound speed).
Compilable, runable. Boundary singularity problem.
Stable version. Added taper func. & init. mag. field
Stable HD version w/ initial larger cavity.
Nozzle solution is now mixed smoothly with Flash solution.
Replenish B field using E field.
New variables: sim_timeMHDon, sim_bPosZ.
Correctly replenish the toroidal field.
Poloidal field (w/ radial residuals).
Poloidal field w/o residuals, but wrong strength.
Move Heat_electricNozzle.F90 to hy_uhd_electricNozzle.F90
Modification of the electric field is now in hy_uhd_unsplit.F90
to exclude the nozzle region from FLASH calculation
Use jetPower to parametrize the jet.
Wind-driven bubble solution for the jet pressure.
Differential update of the nozzle solution in
hy_uhd_unsplitUpdate (instead of in Heat).
More outspread poloidal field structure.
Bigger electric nozzle.
Small radial outflow from the nozzle to stabilize the boundary.
Rename variables: bfeather_* -> rfeather_*.
Minor clean-ups.
Rotation and movement of the nozzle using
"Heat" routine.
- No boundary mixing
- Hard-coded angular veloctiy
Now use betaJet and helicityJet to parametrize the
magnetic field. Field strength changes according
to the pressure decrease.

Angular/linear velocity of the nozzle can be
configured in flash.par.

Renamed variables for consistency:
- bPosZ -> zTorInj
- *_old/*_inner/*_outer -> *Old/*In/*Out
- *feather* -> *Feather*
Added "sim_outflowRation" and "initGeometry".
Cylindrical or spherical initial condition.
Moved Etaper to hy_uhd_electricNozzle.F90
Added density profile.
New parameters:
 - sim_densityProfile
 - sim_densityBeta
 - sim_densityCoreR
Derefine according to jet heights.
New parameters: deref_z1, deref_z2