CMR_aijk Subroutine

public pure subroutine CMR_aijk(ai, daidt, daidt2, k, dkdt, dkdt2, a, dadt, dadt2)

`CMR_aijk`

Calculate the $a_{ijk}$ tensor of the CMR.

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=pr),	intent(in)			::	ai(:)	$\a_i$
real(kind=pr),	intent(in)			::	daidt(:)	$\frac{d\a_i}{dT}$
real(kind=pr),	intent(in)			::	daidt2(:)	$\frac{d^2\a_i}{dT^2}$
real(kind=pr),	intent(in)			::	k(:,:,:)	k_{ijk} matrix
real(kind=pr),	intent(in)			::	dkdt(:,:,:)	$\frac{d k_{ijk}}{dT}$
real(kind=pr),	intent(in)			::	dkdt2(:,:,:)	$\frac{d^2k_{ijk}}{dT^2}$
real(kind=pr),	intent(out)			::	a(:,:,:)	$a_{ijk}$ matrix
real(kind=pr),	intent(out)			::	dadt(:,:,:)	$\frac{da_{ijk}}{dT}$ matrix
real(kind=pr),	intent(out)			::	dadt2(:,:,:)	$\frac{d^2a_{ijk}{dT^2}$ matrix

Called by

Help

Variables

Type	Visibility	Attributes		Name		Initial
real(kind=pr),	public		::	aaa
real(kind=pr),	public		::	aijk_3
real(kind=pr),	public		::	daaadt
integer,	public		::	i
integer,	public		::	j
integer,	public		::	l
integer,	public		::	nc
real(kind=pr),	public,	parameter	::	third	=	1._pr/3._pr