This Package aims to identify discrete vortexs from the velocity field data, which is the result of DaVis PIV process.
Run these in terminal. Set INPUTDIR, OUTPUTDIR and CORE_NUM for your data.
#!/bin/bash#Set INPUTDIR and will create new directory in OUTPUTDIR to save resultINPUTDIR=~/Data/beta_2023_pivdata/cone25_4K/230527_0.4Hz_4K_4sOUTPUTDIR=.#Paralleled Thread Number CORE_NUMCORE_NUM=90make cleanmake DATADIR=$INPUTDIRmpirun -np $CORE_NUM ./vortex_mpi_cmd $INPUTDIR $OUTPUTDIR
To compute all vortexs, there are three steps. First, you have to set parametersinparameters.F90. Second, you have to compile the program. And the third is toprocessing the program.
Since for different original PIV data, different NX_DEF, NY_DEF (point number in x, y direction) should be set before compiling. The batch processingmethod is also provided in order to set NX_DEF and NY_DEF automatically, using shellscriptsetnxny.sh .
- X_MIN_DEF, X_MAX_DEF, Y_MIN_DEF, Y_MAX_DEF
These 4 parameters determine the range of field, a square area. The Unit is meter. For example, X_MIN_DEF -0.1
The points number between X_MIN_DEF and X_MAX_DEF, and also between Y_MIN_DEF and Y_MAX_DEF. The larger INTERPOLATIONPOINTS_DEF, the more accurate for vortexs position. For example, INTERPOLATIONPOINTS_DEF 2000
The points number in x and y direction of the original PIV data. The shell script *setnxny.sh can read NX_DEF and NY_DEF from the first line of data, and write intoparameters.F90.You can also set by yourself.
The vortexs out of these range will be deleted. e.g. EFF_RADIUS_DEF 0.1
Delete a ratio of vortexs. Larger to have fewer vortewx. e.g. FLITER_BETA_DEF 1
The program uses Makefile to compile.
To compile using parameters inparameters.F90, run
in terminal.
To set NX_DEF and NY_DEF forspecified data, run
make cleanmake DATADIR=/directory
in terminal, which can set NX_DEF and NY_DEF inparameters.F90 first and then compiling.
To runsingle processor version program usingmain_findvortex.f90, run
make cleanmake single_processor
in terminal.
Totest program usingmain_test.f90, run
in terminal.
The example is shown in./batch_processing_example.sh.
#!/bin/bash#Set INPUTDIR and will create new directory in OUTPUTDIR to save resultINPUTDIR=~/Data/beta_2023_pivdata/cone25_4K/230527_0.4Hz_4K_4sOUTPUTDIR=.#Paralleled Thread Number CORE_NUMCORE_NUM=90make cleanmake DATADIR=$INPUTDIRmpirun -np $CORE_NUM ./vortex_mpi_cmd $INPUTDIR $OUTPUTDIR
run in terminal.