#!/bin/sh
cd ${0%/*} || exit 1    # Run from this directory

# Source tutorial run functions
. $WM_PROJECT_DIR/bin/tools/RunFunctions

wmake calcRadiusField
wclean calcRadiusField

runApplication blockMesh

##### Procedure for special refinement over Z

# We need the 0 folder to exist for these steps
mkdir 0

# Refine over Z, in 6 passes
for index in 1 2 3 4 5 6; do

  runApplication -s tier$index  calcRadiusField

  runApplication -s tier$index \
      topoSet -dict system/topoSetDict.tier$index

  ## foamToVTK -cellSet tier$index

  runApplication -s tier$index \
      refineMesh -dict system/refineMeshDict.tier$index -overwrite

  rm -r 0/*

done

# Refine over cylindrical coordinates, in 3 passes
for index in 1 2 3; do

  runApplication -s range$index  calcRadiusField -calcDirections

  runApplication -s range$index \
      topoSet -dict system/topoSetDict.range$index

  ## foamToVTK -cellSet tier$index

  runApplication -s range$index \
      refineMesh -dict system/refineMeshDict.range$index -overwrite

  rm -r 0/*

done

#------------------------------------------------------------------------------
