# adding passive tracer to a buoyant plume

Hello,

I am trying to set up an idealized non-hydrostatic simulation of a buoyancy-driven plume rising from a heat source on the bottom. I would like to release a passive tracer at the bottom heat source but could not make it work. I specified the bottom heat and tracer flux (btflx(itemp),btflx(itpas)) in analytical.F. My CPP options are listed below. It is quite puzzling that the temperature profile seems okay but the passive tracer concentration quickly develops very large positive and negative values. I have attached two plots to demonstrate this.

   # define MPI
# define NBQ
# define SOLVE3D
# define NEW_S_COORD
# define SALINITY
# define PASSIVE_TRACER
# define ANA_GRID
# define ANA_INITIAL
# define ANA_SMFLUX
# define ANA_STFLUX
# define ANA_SRFLUX
# define ANA_BTFLUX
# define ANA_SSFLUX
# define ANA_BSFLUX
# define NO_FRCFILE
# define NONLIN_EOS
# define SPLIT_EOS
# define UV_VIS_SMAGO_3D


edit retag close merge delete

Sort by » oldest newest most voted

I am wondering if some of my CPP options are not appropriate when a passive tracer is added? I would appreciate it very much if anyone can provide some suggestions.

more