Solve the partial differential equation governing unsteady 3D homogeneous, isotropic molecular diffusion using 2nd-order-accurate central differences in space and 2nd-order Runge-Kutta time advancement.
Type | Attributes | Name | Initial | |||
---|---|---|---|---|---|---|
type(subdomain_t) | :: | phi | ||||
type(subdomain_t) | :: | phi_half | ||||
integer, | parameter | :: | steps | = | 1000 | |
integer | :: | step | ||||
real, | parameter | :: | D | = | 1. | |
real | :: | dt |