      !-----------------------------------------------------------------
      !  ...input
      !-----------------------------------------------------------------
      character*72 dir, file_snd        ! directory and filename
      integer nt                        ! number of time levels
      character*15 var_snd(nvar_snd)    ! sounding variable names
      real dsnd(kd,nvar_snd)            ! sounding data
      !
      !  "z"            geopotential height (m)
      !  "u"            zonal wind speed (m/s)
      !  "v"            meridional wind speed (m/s)
      !  "th"           potential temperature (K)
      !  "qv"           water vapor mixing ratio (kg/kg)
      !  "p"            pressure (Pa)
      !
      real z(kd), u(kd), v(kd), th(kd), qv(kd), p(kd)
      !-----------------------------------------------------------------
      !  ...output
      !-----------------------------------------------------------------
      character*72 odir, ofile_snd      ! directory and filename
      integer nts, ntf                  ! starting and final time levels
      integer irec_snd                  ! sounding record
      character*15 varo_snd(nvaro_snd)  ! sounding variable names
      !
      !  "po"           pressure (Pa)
      !  "zo"           geopotential height (m)
      !  "uo"           zonal wind speed (m/s)
      !  "vo"           meridional wind speed (m/s)
      !  "tho"          potential temperature (K)
      !  "qvo"          water vapor mixing ratio (kg/kg)
      !
      real zo(kdo), uo(kdo), vo(kdo), tho(kdo), qvo(kdo), po(kdo)
      !
      !  "wspdo"        wind speed (m/s)
      !  "wdiro"        wind direction (degrees from N)
      !  "to"           air temperature (K)
      !  "tdo"          dewpoint air temperature (K)
      !
      real wspdo(kdo), wdiro(kdo), to(kdo), tdo(kdo)
