Adding a pencil for soret diffusion

7 views
Skip to first unread message

Yuvraj Sinha

unread,
Sep 17, 2025, 6:06:00 PMSep 17
to pencil-code-discuss
Hello all,

I am trying to include the soret diffusion (thermal diffusion) term in the flux term of the species conservation equation. I want to introduce a pencil corresponding to the thermal diffusion coeffecient same as Mixture averaged diffusion coefficient. But the pencil 'Diff_penc_add'  is provided for the misture averaged diffusion coefficient in (chemistry.f90). Is it possible to provide a separate pencil 'Therm_Diff_penc_add'?

I checked cparam_pencils.inc including 200 pencils and realised it is generated by mkcparam. But couldn't figure out how to define additional pencil say using scripts/mkcparam.

Any help regarding this is appriciated.

Best,
Yuvraj

Kishore G

unread,
Sep 18, 2025, 12:57:20 AMSep 18
to pencil-code-discuss, Yuvraj Sinha
Hi Yuvraj.

You need to declare the pencil in the header of the relevant module (see
section 9.7 of the manual, entitled 'the pencil case').

--
Kishore G.


Yuvraj Sinha

unread,
Sep 18, 2025, 11:33:01 AMSep 18
to pencil-code-discuss
Hi Kishore,

Thanks for your help. I was able to declare a pencil.

Best,
Yuvraj 

Reply all
Reply to author
Forward
0 new messages