cws_read__readData_TPS.pro

;  reads TROPOS DATA
;

function _read_cop, file, var
  fileId = h5f_open(file)
  dataId = h5d_open(fileID,'ScienceData/'+var)
  data = h5d_read(dataID)
  return, rotate(data,7)          
  
end 



function cws_read::readData_TPS
  CASE strLowCase(self.product) OF 
          'cmb': begin
               self.infile = FILE_SEARCH(self->build_source_path() $
				      + '/TROPOS_MSI_CM_*'+ self->get_date(/string,/KMR) +'*.h5',c=c)
             
            fileId = h5f_open(self.infile)
		        dataId = h5d_open(fileID,'ScienceData/cloud_mask')
		        data = h5d_read(dataID)
            img = make_array(3712,3712,/byte,value=0.)
            img = (data eq 0 or data eq 1) + 100 * ( data eq 2 or data eq 3)
            img = rotate(img,7)
            
		        return,img   
          end 
          
          'ctp': begin
              self.infile = FILE_SEARCH(self->build_source_path() $
				      + '/TROPOS_MSI_COP_*'+ self->get_date(/string,/KMR) +'*.h5',c=c)
          
            img = _read_cop( self.infile, 'cloud_top_pressure')
            return,img
           
          end
          
          'cth': begin
              self.infile = FILE_SEARCH(self->build_source_path() $
				      + '/TROPOS_MSI_COP_*'+ self->get_date(/string,/KMR) +'*.h5',c=c)
          
            img = _read_cop( self.infile, 'cloud_top_height')
            return,img
           
          end
          
           'ctt': begin
              self.infile = FILE_SEARCH(self->build_source_path() $
				      + '/TROPOS_MSI_COP_*'+ self->get_date(/string,/KMR) +'*.h5',c=c)
          
            img = _read_cop( self.infile, 'cloud_top_temperature')
            return,img
           
          end
          
           'cod': begin
              self.infile = FILE_SEARCH(self->build_source_path() $
				      + '/TROPOS_MSI_COP_*'+ self->get_date(/string,/KMR) +'*.h5',c=c)
          
            img = _read_cop( self.infile, 'cloud_optical_thickness')
            return,img
           
           
           
           
          end
          
          'ref': begin
              self.infile = FILE_SEARCH(self->build_source_path() $
				      + '/TROPOS_MSI_COP_*'+ self->get_date(/string,/KMR) +'*.h5',c=c)
          
            img = _read_cop( self.infile, 'cloud_effective_radius')
            return,img
      
            end
             'cph': begin
               self.infile = FILE_SEARCH(self->build_source_path() $
				      + '/TROPOS_MSI_CM_*'+ self->get_date(/string,/KMR) +'*.h5',c=c)
             
            fileId = h5f_open(self.infile)
		        dataId = h5d_open(fileID,'ScienceData/cloud_phase')
		        data = h5d_read(dataID)
            img = make_array(3712,3712,/byte,value=0.)
            img = (data eq 0 or data eq 1) + 100 * ( data eq 2 or data eq 3)
            img = rotate(img,7)
            
		        return,img   
          end 
            else:
        endcase



end function