Ask Your Question

adding passive tracer to a buoyant plume

asked 2021-08-21 22:23:43 +0200


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 UV_ADV
   # define TS_HADV_WENO5
   # define TS_VADV_WENO5
   # define BIO_HADV_WENO5
   # define UV_HADV_WENO5
   # define UV_VADV_WENO5
   # define W_HADV_WENO5
   # define W_VADV_WENO5
   # 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

image description image description

edit retag flag offensive close merge delete

1 answer

Sort by ยป oldest newest most voted

answered 2021-08-22 18:12:37 +0200

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.

edit flag offensive delete link more
Login/Signup to Answer

Question Tools

1 follower


Asked: 2021-08-21 22:23:43 +0200

Seen: 10 times

Last updated: Aug 22