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 |