# Flood fill for phase field
# This is a FORTRAN code to identify region outlined by a phase-field variable