The code can be compiled as (for example):
clang++ vel_channel2D.cpp -std=c++11 -O3

It can then be run as
./a.out 0 0
where the numbers specify the location of the localized initial population. (0,0) gives a population at the center of the domain. 

The code reads in velocity files generated using an ocean model initialized with observational data. The complete filepath to these velocity files must be specified in the read_velocity function. Examples are included in the k46 folder.

The output of the model is concentration fields for species 0 and species 1 (file format c_[species]_[time]_[x_location]_[y_location].dat) as well as total concentration (file format ctots_[x_location]_[y_location]_[init time].dat). These files are stored in a data folder automatically. The model output can then be analyzed.