----------------------------------------------
   Pochva test package
----------------------------------------------

This package contains the code of “Pochva” land-model (pochva.F90,
pochva_module.F90) used for testing the column model (model.F90). The
model reads data from 2 input data files: initial.dat2 (contains all
physical parameters of soil and vegetation for the soil column and
initial data for temperature and water content of soil) and
atm_forcing.dat2 (contains data for the parameters requested by Pochva
scheme to simulate the interaction between the surface and the
atmosphere: bottom atmospheric pressure, temperature, humidity,
turbulent exchange coefficient, solar radiation fluxes, precipitation
fluxes and time-varying vegetation parameters).

The model performs a time integration for 2,5 year calling at each
time step the Pochva scheme subroutine (time step values is fixed 60s
and determined by the input data prepared in atm_forcing.dat2).
Pochva subroutine writes output data in files: atm_01.dat,
surf_01.dat, tsoil_01.dat, qsoil_01.dat, frac_sice_01.dat,
flux_entropy_soil_01.dat, flux_water_soil_01.dat, lev_snow_01.dat,
tsnow_01.dat, fice_snow_01.dat, snow_age_01.dat, snow_melt_age_01.dat.

There are input data for 6 different geographical points:
1 - Leipzig (Germany), 2 - Moscow (European Russia),
3 - Rennes (France), 4 - Larisa (Greece),
5 - Diyrbakir (Turkey), 6 - Kutsho (Norway).
All data have been obtained by the Bolam NWP model simulation
described in the article.

The sample input data, compressed with gzip, can be downloaded at the
following address:
https://www.isac.cnr.it/dinamica/oxana/POCHVA/pochva_test_input_data/
(total download size about 3605MB).

Output data produced by the simulation *.dat may be visualized by
gnuplot tool using the script gnuplot_pochva_output.load.

You may perform the simulation, obtain the figures with output data
and compare the obtained figures with figures prepared by the author
(pochva_output_diagnostic_01_author.pdf, pochva_output_diagnostic_02_author.pdf,
pochva_output_diagnostic_03_author.pdf, pochva_output_diagnostic_04_author.pdf,
pochva_output_diagnostic_05_author.pdf, pochva_output_diagnostic_06_author.pdf)
available in the same input data address:
https://www.isac.cnr.it/dinamica/oxana/POCHVA/pochva_test_input_data/

You can set the compilation parameters in Makefile.

Compile with command make.

Execute test integration by using launch.script setting DIR_DATA
parameter which is the path to the directory with input data
(initial*.dat2, atm_forcing*.dat2).  Attention! time integration is
very fast, while elaboration of output data by gnuplot tool is very
slow, please be patient.