README update

README.md

@@ -34,6 +34,40 @@ Option `COMPILER` enables set of compiler-specific options. Most common are `gnu
 ./sfx.exe
 ```
 
+## Running test cases
+We have prepared several dataseta to make the test calculations. Data is prepared on the basis of measurements.
+
+1) MOSAiC Shupe M. D. et al. Overview of the MOSAiC expedition: Atmosphere // Elem Sci Anth. – 2022. – Т. 10. – №. 1. – С. 00060.
+2) Irgason - Measurements on Mt. Elbrus
+3) Sheba https://data.eol.ucar.edu/project/SHEBA
+
+For MOSAiC put - 1, for Irgason - 2, for Sheba - 3.
+There are two input files for each test in folder ./data/. For example MOSAiC.txt and MOSAiC_zh.
+The first file contains data on:
+1. wind velocity  (m/s),  
+2. difference between potential temperature at constant flux layer height and at the surface  (deg. K)          
+3. semi-sum of potential temperature at constant flux layer height and at surface  (deg. K)          
+4. difference between humidity at constant flux layer height and a surface (kg/kg)    
+
+The number of rows corresponds to the number of time intervals.
+Files with endings _zh.txt contain two values that are time-independent: constant flux layer height (m) and roughness of surface (m)
+
+Then build the executable file using "makefile". If the build is successful, the file will appear sfx.exe. Run the file sfx.exe.
+The output values will appear in the folder ./output/. The name of the file will correspond to the name of the test.
+The output file contains data on:
+
+1. non-dimensional constant flux layer height (z/L)  L-Obukhov length                         
+2. Richardson number                                    
+3. Reynods number                                       
+4. ln(zu/zt)                                            
+5. dynamical roughness zu (m)                           
+6. thermal   roughness zt (m)                           
+7. critical Richardson number                           
+8. transfer coefficient for momentum                    
+9. transfer coefficient for heat                         
+10. eddy viscosity coefficient  (m**2/s)
+11. inverse turbulent Prandtl number
+
 ## Developing and working with issues
 Suppose you were assigned to an issue `introduce a feature` for one of the projects mentioned above and you need to make some changes in *nse-diurnal-les* project. Your typical workflow will be like this: