Description of model and simulation 1: conduction
Start with this template file (f77) or this one (f90).Description of simulation 2: linear convection
This version of the f90 template program will run, store some results in an array, print the array using subroutine PRINT, and write the array to a file that can be read by MATLAB.
This MATLAB version of the f90 template program will run. Variables in the workspace.can be accessed.
To help you write your Fortran code, use this diagram and the diagrams in these slides of the arrangement of the model variables upon the staggered grid.
Learn about printing source code and the Fortran compiler..
Use subroutine PRINT to print the results of your simulations.
To help you write
your Fortran code, use this diagram
and the diagrams in these slides of the arrangement of the model variables upon the staggered grid.
Use the MATLAB program QCOM_movie.m to produce animations from 2D QCOM output.
Results for simulation 2: linear
convection
Simulation of Kelvin-Helmholtz instability using QCOM (256 x 256 grid points):
Evolution of potential tempeature field for Re = 900 and Re = 1400 (animation)