We recommend testing that TRACMASS was properly compiled by letting PROJECT and CASE be “Theoretical” in the Makefile (which is the default). In this case, TRACMASS will use a simple oscillating velocity field to trace trajectories.
-
Make sure in Makefile both PROJECT and CASE are set to Theoretical.
-
Recompile the code
make clean make
-
Run the code by typing
./runtracmass
-
If everything went without an error, you will get the following output in the terminal:
=============================================================================== TRACMASS Lagrangian off-line particle tracking =============================================================================== Start date : 2020-06-30 Start time : 20:57:55 ------------------------------------------------------------------------------- Model information: Project code : Theoretical Case name : Theoretical Namelist file : projects/Theoretical/namelist_Theoretical.in ------------------------------------------------------------------------------- Output information: Directory for output files : /Users/user/Tracmass_output Prefix for output files : test ------------------------------------------------------------------------------- Configuration options: - Two-dimensional trajectories, no change in depth ------------------------------------------------------------------------------- Start date in model-time : 2000-01-01 12:00 End date in model-time : 2000-01-01 21:00 Length of run in timesteps : 10 Number of seeding timesteps : 1 Steps between two GCM fields : 100 ------------------------------------------------------------------------------- t-step run out err tot model date ------------------------------------------------------------------------------- 1 | 1 | 0 | 0 | 1 | 2000-01-01 12:00:00 2 | 1 | 0 | 0 | 1 | 2000-01-01 13:00:00 3 | 1 | 0 | 0 | 1 | 2000-01-01 14:00:00 4 | 1 | 0 | 0 | 1 | 2000-01-01 15:00:00 5 | 1 | 0 | 0 | 1 | 2000-01-01 16:00:00 6 | 1 | 0 | 0 | 1 | 2000-01-01 17:00:00 7 | 1 | 0 | 0 | 1 | 2000-01-01 18:00:00 8 | 1 | 0 | 0 | 1 | 2000-01-01 19:00:00 9 | 1 | 0 | 0 | 1 | 2000-01-01 20:00:00 10 | 1 | 0 | 0 | 1 | 2000-01-01 21:00:00 =============================================================================== 1 particles calculated 0 particles exited the space and time domain 0 particles flagged with errors 1 particles in domain ------------------------------------------------------------------------------- =============================================================================== End date : 2020-06-30 End time : 20:57:55 ===============================================================================
-
You can rerun the code by modifying the namelist.in in the main directory.
How to configure a customised experiment
This is a basic guidance on how to implement a new setup from a GCM with a grid that does not exist in the present TRACMASS projects list. This requires good knowledge in the numerical discretisation of your GCM and of Fortran programming.
-
Make a new «project» which we call “mymodel” here.
Copy the directory from an existing project that has a similar grid (A,B or C-grid) or has a similar format (netcdf, grib, etc.). So that:
cd projects cp -r NEMO mymodel
or:
cp -r AVISO mymodel
or:
cp -r IFS mymodel
-
Rename the namelist file to namelist_mymodel.in or namelist_mycase.in if you are using a case within the project. Find more information about the namelist here.
-
Modify and adapt the files under /projects/mymodel, which are setup_grid.F90, read_field.F90, and kill_zones.F90. Find more information about the subroutines setup_grid.F90 and read_field.F90 here.
-
Adapt the Makefile.prj to adjust the compiling flags. The list of possible flags can be found here.
-
Adapt the Makefile in the main directory and set PROJECT = mymodel and CASE = mymodel.