
* INPUT :Time series 
* parameters marked with (E) not needed for pre-processing
* (E)  1 CFNR	record no			 
* (E)  2 CFST	station id			 
*      3 CFVU	year				 
*      4 CFKU   month				 
*      5 CFPV	day				 
*      6 CFKL 	time (HHMM)	note: Finnish time (UTC+2)		 
*      7 CFPR 	sea level pressure (mbar)		 
*      8 CFT2	temperature at 2m (K)			 
* (E)  9 CFUU	rel. humidity (%)   		 
*     10 CFSU	state of ground (0 -> 9)		 
*         0 dry 1 humid 2 wet 3 icy/dry 4 wet ice/ice
*         5 snow or melting snow <less than 50%  6  >50%
*         7 snow or melting snow 100 %
*         8 >50 % <100% covered by dry snow 9 100 % snowcover
*     11 CFNN  total cloudiness (0->1)                   
* (E) 12 CFTD	dew point temp (K)			 
* (E) 13 CFTW	wet-bulb temp (K)		 
*     14 CFRR   amount of rain (mm/h)                
* (E) 15 CFVV   visibility (m)                  	
*     16 CFWW	weather/synop code (0 - 99)              
*     17 CFW1 	weather of prev. hour (0 -> 9)		 
*     18 CFW2	weather of prev. 3 hours (0 -> 9)	 
*     19 CFNH	amount of low clouds (0 ->1)		 
* (E) 20 CFCL   type of low clouds (code) 	 
*     21 CFHH	height of low clouds (m) (60 - 6000)	 
* (E) 22 CFCM  type of middle clouds 	 
* (E) 23 CFCH  type of high clouds	 	 
*     24 CFDD    direction of flow (ARITM. degr)) 	 
*              NOTE!!!!! 
*              this is NOT the meteorological wind direction but 
*              the direction of the flow in x-y coordinates =(270 -met. angle)
*              e.g. Wind from north -> flow direction = 270 degrees
*     25 CFUA	windspeed at 10 meters (m/s)		 
********parameters ->(25) from basic synoptic measurements

*     26 CFAP	hourly amount of sunshine (rel) (0->1)   
*	! if missing counted from cloudiness data by regression
*****************************************************************************
*
OUTPUT (hourly values)
                    
*     27 CFAB	albedo (0->1)				 
*     28 CFAK	solar elevation (deg)			 
*     29 CFAS   Solar radiation (W/m**2)		 
*     30 CFAL   Moisture parameter (0->1)		 
*     31 CFPL	Inverse of M-O length (1/m)		 
*     32 CFTT	Temperature scale (K)			 
*     33 CFUT	Friction velocity (m/s)			 
*     34 CFH0	Turb. heat flux H (W/m**2)		 
*     35 CFQT	Net radiation Qt (W/m**2)		 
*     36 CFLE   Latent heat flux lE (W/m**2)	

+++ mixing height
INPUT: 2 soundings /day 

OUTPUT
*     37 CFH1	mixing height(m)
*     38 CFH2	height of wind shear layer (m)
*     39 CFWT	convective velocity scale (m/s)
*     40 CFGR 	potential temp. gradient ; Z=CFH2 (K/m)

