Tracers with Ariane

In addition to a particle’s trajectory (longitude and latitude) and depth, Ariane can also help us analyze tracers along the particle’s trajectory.

  • Temperature

  • Salinity

  • Density

We will be making changes in namelist.

Namelist: Add Sections

Ariane requires both salinity and temperature data as input. It also requires density data or an indication that density should instead be calculated using salinity and temperature.

Our model produces files with filenames that follow this format: SalishSea_t_yyyymmdd_yyyymmdd_grid_T.nc. These files contain salinity and temperature. Therefore, Ariane will be calculating density.

Add the following sections to namelist. Remember to modify the directory and filename of your input files.

Temperature

&TEMPERAT
    c_dir_te ='/ocean/nsoontie/MEOPAR/SalishSea/results/storm-surges/final/dec2006/all_forcing/30min/',
    c_prefix_te ='SalishSea_30m_20061214_20061215_grid_T.nc',
    ind0_te =-1,
    indn_te =-1,
    maxsize_te =-1,
    c_suffix_te ='NONE',
    nc_var_te ='votemper',
    nc_att_mask_te ='NONE',
/

Salinity

&SALINITY
    c_dir_sa ='/ocean/nsoontie/MEOPAR/SalishSea/results/storm-surges/final/dec2006/all_forcing/30min/',
    c_prefix_sa ='SalishSea_30m_20061214_20061215_grid_T.nc',
    ind0_sa =-1,
    indn_sa =-1,
    maxsize_sa =-1,
    c_suffix_sa ='NONE',
    nc_var_sa ='vosaline',
    nc_att_mask_sa ='NONE',
/

votemper and vosaline are the variable names for temperature and salinity in our input file.

Namelist: Add Parameters

Add the following parameters to the ARIANE section of the namelist:

  • key_computesigma: False if density is read into Ariane; True if Ariane calculates density.

  • zsigma: Reference depth for calculation of density.

And change:

  • key_alltracers: True if tracers are required.

Ariane

&ARIANE
        key_alltracers =.TRUE.,
        key_sequential =.FALSE.,
        key_ascii_outputs =.TRUE.,
        mode ='qualitative',
        forback ='forward',
        bin ='nobin',
        init_final ='init',
        nmax =5,
        tunit =3600.,
        ntfic =1,
        tcyc =0.,
        key_computesigma =.TRUE.,
        zsigma =1.,
/

Results

The variables names for the tracers are:

  • traj_temp

  • traj_salt

  • traj_dens

Plots

The results produced for the example above:

../_images/result_tracers_doc.png