; Reads data from the Scope experiment, stored in the file Scope2.new txt.

;  File Scope2b.new is tab delimited with one header line.  
;  The variables are:
;  Date day
;     Month
;  GMT   decimal hour
;  u    wind speed m/s
;  Tw   water temperature  Cent
;  Ta   air temperature Cent
;  Tw-Ta
;  Time Decimal day, september
;  qs   water humidity g/kg
;  q    air humidity g/kg
;  del q     qs-qa
;  Hsc  Covariance sensible heat W/m^2
;  Hsb  Bulk sensible heat W/m^2
;  Hsi  Inertial sensible heat W/m^2
; <Hs>      Average sensible heat
;  Hlc  Covariance latent heat W/m^2
;  Hlb  Bulk latent heat W/m^2
;  Hli  Inertial latent heat W/m^2
;  <Hl>      Average latent heat
;  Tuc  Covariance stress N/m^2
;  Tub  Bulk stress N/m^2
;  Tui  Inertial stress N/m^2
;  <Tu>      Average stress
;  U*   friction velocity m/s
;  Oph  quality number for humidity; good if less than 15 and greater than 0
;  Tilt flow tilt in deg
;  Zi   Inversion height m
;  -z/L Dimensionless monin obukhov similarity parameter
;  W*   Deardorff convective velocity scale
;  tj UT     same as Time
;  Gmin min value of rms waveheight during period cm
;  Gmax max value of rms waveheight during period cm
;  fmin min value of dominant wave frequency during period Hz
;  fmax max value of dominant wave frequency during period Hz
;  <G>  average of min and max
;  <f>  average of min and max
;  Co   wave phase speed for <f> m/s
;  Co/U*     wave age with u*
;  Co/Ua     wave age with U
;  H1/3 Significant wave height (4*rms) m
;  dir  True wind direction deg

  function read_SCOPE, spec_hum, U_star, wind_spd, sfc_spec_hum, wave_age, $
 Tair, Tsst, Hsig

; open data file 
  OPENR, SCOPE_data, "/usr/people/stewart/work/airsea/Scope2b.new", /GET_LUN
 
; skip first line (with headers)
  junk = ' '
  READF, SCOPE_data, junk

; define output arrays
  n_data = 305L
  wind_spd = FLTARR( n_data )
  spec_hum = wind_spd
  Tair = wind_spd
  Tsst = wind_spd
  sfc_spec_hum = wind_spd
  wave_age= wind_spd
  Hsig = wind_spd
  U_star = wind_spd
  other_wave_age= wind_spd
  z_over_L = wind_spd 

; read in all data (305 lines)
  For i = 1L, n_data do begin
    READF, SCOPE_data, date, gmt, ums, Tw, Ta, Tw_Ta, time, qs, q, delq, Hsc, Hsb, Hsi,$
 Hs, Hlc, Hlb, Hli, Hl, Tuc, Tub, Tui, Tu, Ustar, Oph, Tilt, Zi, zL, Wstar, tjUT,$
Gmin, Gmax, fmin, fmax, G, f, Co, CoUstar, CoUa, H13, dir

     spec_hum(i-1) = 0.001 * q    
     U_star(i-1) = Ustar 
     wind_spd(i-1) = ums 
     sfc_spec_hum(i-1) = 0.001 * qs
     wave_age(i-1) = CoUstar
     Tair(i-1) = Ta
     Tsst(i-1) = Tw
     Hsig(i-1) = H13
     other_wave_age(i-1)=CoUa
     z_over_L(i-1) = zL

  ENDFOR
; print,U_star, wave_age

; Find the points where QC conditions are met
 i_good = where( Oph LT 15.0 AND wave_age GT 0.0 AND ABS( z_over_L ) LT 0.18 )
;  i_good = where( sound_spd GT 332.0 AND sound_spd LE 338 AND $
;   stnd_dev_Hdg LT 14 )

; Save only the data that meet the QC conditions
  spec_hum = spec_hum(i_good)
  u_star = u_star( i_good )
  wind_spd = wind_spd(i_good)
  sfc_spec_hum = sfc_spec_hum(i_good)
  wave_age = wave_age( i_good )
  Tair = Tair( i_good )
  Tsst = Tsst( i_good )
  Hsig = Hsig( i_good )
;  shf = shf( i_good )
;  lhf = lhf( i_good )
  atmos_stab = z_over_L( i_good )

;  stop
  return, N_elements( wind_spd )
;print, n_data
END