binietoglou@5: Basic instructions
binietoglou@5: ==================
binietoglou@5: 
binietoglou@5: These classes can be used to handle lidar input data. They are still in a very initial stage. There are many features probably not working, but they work for some specific tasks. They work with Licel input files (and also with the Raymetrics modified format).
binietoglou@5: 
binietoglou@5: 
binietoglou@5: Set up
binietoglou@9: ------
binietoglou@5: 
binietoglou@5: 
binietoglou@5: Parameter file
binietoglou@9: ~~~~~~~~~~~~~~
binietoglou@9: 
binietoglou@5: Before using the classes you need to setup some channel parameters, that are used when converting the lidar data to Single Calculus Chain format. 
binietoglou@5: 
binietoglou@5: All the parameters are read from an external file stored in the same folder as the code. You can start by changing the file “cf_netcdf_parameters.py” that describe such parameters for the Clermont Ferrand  lidar. 
binietoglou@5: 
binietoglou@5: System class
binietoglou@9: ~~~~~~~~~~~~
binietoglou@5: The next thing you need to create a class that describes you system. This is very simple if your lidar data are in the Licel format, as you only need to specify the external file with the extra SCC parameters. You can use as an example the file “cf_raymetrics.py”::
binietoglou@9:    
binietoglou@9:    :::python
binietoglou@6:    from licel import LicelLidarMeasurement
binietoglou@6:    import cf_netcdf_parameters
binietoglou@5: 
binietoglou@6:    class CfLidarMeasurement(LicelLidarMeasurement):
binietoglou@6:       extra_netcdf_parameters = cf_netcdf_parameters
binietoglou@6: 
binietoglou@5: 
binietoglou@5: Using the class
binietoglou@9: ---------------
binietoglou@5: 
binietoglou@6: Once you have made the above setup you can start using it. The best way to understand how it works is through an interactive shell (I suggest ipython(http://ipython.org/)).  In the following example I use the cf_raymetrics setup::
binietoglou@8:    
binietoglou@9:    :::python
binietoglou@7:    import glob # This is needed to read a list of filenames
binietoglou@7:    from lidar import cf_raymetrics #If you have saved the files in a directrory called “lidar”
binietoglou@5: 
binietoglou@7:    # Go to the folder where you files are stored
binietoglou@7:    cd /path/to/lidar/files
binietoglou@5: 
binietoglou@7:    # Read the filenames
binietoglou@7:    files  = glob.glob('*') # The * reads all the files in the folder.
binietoglou@5: 
binietoglou@7:    #Read the files
binietoglou@7:    my_measurement = cf_raymetrics.CfLidarMeasurement(files)
binietoglou@5: 
binietoglou@7:    #Now the data have been read, and you have a measurement object to work with:
binietoglou@7:    # See what channels are present
binietoglou@7:    print my_measurement.channels
binietoglou@5: 
binietoglou@7:    #Quicklooks of all the channels
binietoglou@7:    my_measurements.plot()
binietoglou@5: 
binietoglou@8: 
binietoglou@5: Converting to SCC format
binietoglou@9: ------------------------
binietoglou@5: 
binietoglou@5: There are some extra info you need to put in before converting to SCC format, 'Measurement_ID', 'Temperature', 'Pressure'::
binietoglou@8:    
binietoglou@9:    :::python
binietoglou@7:    my_measurement.info['Measurement_ID'] = “20101229op00”
binietoglou@7:    my_measurement.info['Temperature'] = “14”
binietoglou@7:    my_measurement.info['Pressure'] = “1010”
binietoglou@5: 
binietoglou@8: 
binietoglou@6: You can use standard values of temperature and pressure by just calling::
binietoglou@8:    
binietoglou@9:    :::python
binietoglou@7:    my_measurement.get_PT() 
binietoglou@5: 
binietoglou@5: The standard values can be changed in generic.py. Search the  get_PT method and change of what is appropriate for your station. If you have an external source of temperature and pressure information (a meteorological station) you can automate this by overriding the get_PT method in your system's class (in our example in the cf_raymetrics.py file).
binietoglou@5: 
binietoglou@8: 
binietoglou@6: After you have used this extra input, you save the file using this command::
binietoglou@9: 
binietoglou@9:    :::python
binietoglou@7:    my_measurement.save_as_netcdf(“filename”)
binietoglou@5: 
binietoglou@5: where you change the filename to the filename you want to use.
binietoglou@5: