Source codes


FLake is freely available under the terms of the MIT license.
Download tar.gz file or zip file. These files contain
(i) routines of the lake model FLake,
(ii) routines of the atmospheric surface-layer parameterisation scheme SfcFlx,
(iii) an example of FLake Interface for single-column applications.
A brief description of the source code is given in the synopsis of FLake routines (click Documentation and download PDF file).
Test runs provide examples of the use of FLake in a stand-alone mode.

If you download the FLake code, please send the following information to Arkady Terzhevik:
(i) your name,
(ii) your address,
(iii) your email address,
(iv) a few words to give an idea of how you are going to use FLake.
Arkady will then add you to the list of FLake users (this does not assume any obligations).



As a lake parameterisation module, FLake is implemented into a number of NWP and climate models (see menu item "Applications"). The source codes of those models are not necessarily freely available. If you are interested in details of the implementation of FLake into one or the other NWP or climate model, please contact the code managers of the respective model.


  • Download the precompiled Windows version of the model

    The windows-version of the FLake is compiled using Intel Fortran 8.1 for Windows under WinXP and should run on virtually any 32-bit Windows system. The main program in source files was slightly changed in comparison to the original Unix-version in order to provide a quick model reconfiguration without editing the source files.
  • Configuring the model

    To start the model on Your lake, apart from the model executable file you will need two additional text files describing the problem:
    1. 'namelist' file with the model configuration. Here all necessary information about your lake should be given. The file is in FORTRAN95 namelist format (see FORTRAN help for more information). It includes 4 groups of parameters, among them simulation period, time step, lake depth, typical fetch, filename for the input meteorology, output filename etc. Three example files Mueggelsee80-96.nml, Heiligensee80-96.nml and Stechlin94-98.nml can be downloaded, for three model cases described in Test Runs. The files have self-explaining comments and can give an overall idea for configuring the model without detailed knowledge of the namelist format. The files can also be used as a starting point for further editing.
    2. The file with the external forcing (meteorology input). The file should have the same name as that given in the configuration file under the METEOFILE variable (see examples in .nml files).

      The file consists of 6 columns divided by space characters with following values (See example files Potsdam80-96.dat, and Stechlin94-98.dat from Test Runs):

      1 2 3 4 5 6
      Sequential number Solar Radiation (W/m2) Air Temperature (oC) Air Humidity (mb) Wind Speed (m/s) Cloudiness (0-1)
      NOTE: The number of rows should be equal to or greater than the value of TIME_STEP_NUMBER in the namelist file.
  • Model run

    Put all 3 files in the same directory, open the directory in the Windows Explorer and type in the command line (press "Windows key"+R to call the command window):

    flake < YOURFILE >, 
    where < YOURFILE > is the name of Your namelist configuration file. The file with modeling results will be created after a successful model run.
  • Model output

    The output file name is that you have given in the namelist file as the value for the OUTPUTFILE variable. Variables are in columns, variable names are given in the first row. These are:
    No - sequential number
    time - time in days
    Ts - surface temperature (oC)
    Tm - mean water temperature (oC)
    Tb - bottom temperature (oC)
    ufr_a - friction velocity in air (m/s)
    ufr_w - friction velocity in surface water (m/s)
    Wconv - convective velocity scale (m/s)
    Q_se - Sensible surface heat flux (W/m2)
    Q_la - latent surface heat flux (W/m2)
    I_w - short-wave radiation (W/m2)
    Q_lwa - long-wave radiative flux from the atmosphere (W/m2)
    Q_lww - long-wave radiative flux from the water (W/m2)
    h_ML - mixed layer depth (m)
    C_T - stratification "shape factor" (dimensionless)
    H_B1 - depth of the thermal wave penetration in sediments (m)
    T_B1 - temperature at the crest of the thermal wave in sediments (oC)
    Q_bot - heat flux across the water-sediments boundary (W/m2)
    H_ice - ice thickness (m)
    H_snow - snow thickness (m)
    T_ice - ice temperature (oC)
    T_snow - snow temperature (oC)

    If you use MATLAB for postprocessing and analysis of results the matlab-script ReadFlake.m can help importing data directly to MATLAB. Just run it from the matlab prompt.

  • Please send any feedback concerning the stand-alone windows-version of FLake to Georgiy Kirillin