4 use input, only : ifmhd, ifaxis, ifcvode
5 use tstep, only : ifield
8 logical :: ifturb,if_conv_std
11 if (ifmhd .AND. ifaxis) if_conv_std = .false.
13 if(ifcvode .AND. ifield == 2) call
setprop
18 if (ifturb) call maketq
21 if ( .NOT. ifturb .AND. if_conv_std) call
makeuq
subroutine makeuq
Fill up user defined forcing function and collocate will the mass matrix on the Gauss-Lobatto mesh...
subroutine setprop
Set variable property arrays.