Restore the original Suresh & Huynh limiter and add modification by
Ahn & Lee for compact schemes. Additionally, pass the variable
positivity flag and for positive variables modifie the lower limit
using the minimum value obtained from the cubic spline interpolation.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This is an improvement in the case a job restarted with a higher
number of processes. The reading of data blocks is spread among new
processes, so it better balances memory.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
By default this parameter is set to zero and the behavior is the same as
before. However, if this parameter is set to a value larger than zero,
it determines the driving time scale for the Ornstein–Uhlenbeck method,
and the injection power for the Alvelius method.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
The first line shows the location from which the message was called, and
the second one shows the actual message.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
New functions uppercase() and lowercase() change the input string to the
upper case and lower case, respectively.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
The possible failure could take place in the boundary update or
conversion to the conservative variables. In case of a failure, jump
directly to the end part of subroutine in which the %boundary flag is
cleared.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
After the boundary update, all respective blocks modified by the mesh or
boundary update are required to update their primitive variables. If a
data block contains a non-physical cell, it is corrected by
correct_unphysical_states(), if fix_unphysical_cells is switched on.
As the last steps, the shapes are updated if they are switched on.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This field should be set to .true. if all block cells are physical,
otherwise it should be set to .false..
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
If cons2prim() returns status indicating the existence of non-physical
cell, just update the subroutine status and continue with the variable
conversion.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
Do not quit immediately if a negative density or pressure is
encountered. Just mark the status flag non-zero and continue. The
non-physical cell can be corrected later.
Also, if the compilation flag DEBUG is present, print the conservative
variables for the non-physical cell.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
When the block has been (de)refined, its field %update is set in order
to update its primitive variables. However, the block which are
neighbors of the (de)refined blocks, since their boundaries are updated
too, were not selected for the primitive variable update.
Fix this by adding a new meta block field %boundary, which indicates
that a block boundary has been updated, and update the corresponding
data block primitive variables either if the %update or %boundary flags
are set.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
The specific boundaries (open, outflow, reflective, custom, etc.) work
on the primitive variables. Since the boundaries are updated using the
conservative variables, we have to update the primitive variables in the
block to which the specific boundary is applied. Then, the conservative
variables are update at the ghost zone where the specific boundary
condition was applied.
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
