Airborne doppler radars profiles: data format

Airborne doppler radars profiles data set description page

Airborne doppler radars profiles data set

Each file name davad_IOPx_By*.dat is made by appending the IOP number 
and the purl (circular trajectory) number to the method ident (davad). 
Each file contains a header of 2 lines followed by a data block of 
24 lines (altitude from 150m to 7050m, every 300m).

Header organization (line 1 to 2)
--------------------

Line No Parameter Meaning                                   Unit          Format
---- -- --------- ----------------------------------------- ------------  ------
 1    1 provider  Scientist and laboratory who provided the                 a15
                  data set 
 1    2 method                                                              a14
 1    3 version                                                             a11
 1    4 timestamp                                             DD/MM/YY      a8

 2    1 beg_date  Begin time for processing                    HHMISS       a6
 2    2 end_date  End time for processing                      HHMISS       a6
 2    3 lat       Latitude of reference for the vert. prof.   dec.degr.     f9.3
 2    4 lon       Longitude of reference for the vert. prof.  dec.degr.     f9.3

Data block organisation
-----------------------

 No Parameter Meaning                                  Unit        Format Miss.
 -- --------- ---------------------------------------- --------    ------ -----

  1 height    Altitude                                 m           f9.3   
  2 u0        West-east wind component                 m/s         f9.3    999
  3 v0        South-North wind component               m/s         f9.3    999
  4 ref       Reflectivity                             dBZ         f9.3    999
  5 div       Horizontal divergence                    10-4 s-1    f9.3    999
  6 vort      Vertical component of vorticity          10-4 s-1    f9.3    999
  7 stretch   Stretching deformation                   10-4 s-1    f9.3    999
  8 shear     Shearing deformation                     10-4 s-1    f9.3    999
  9 Vt        Terminal fall velocity of hydrometeors   m/s         f9.3    999
              (Convention : Positive downward)

Example to read data in FORTRAN format:

      program read_davad

      parameter (nplan = 24)                           ! there are 24 heights for all DAVAD files

      real lat,lon,height(nplan),u0(nplan),v0(nplan),ref(nplan),
     &  div(nplan),vort(nplan),stretch(nplan),shear(nplan),Vt(nplan)
      character provider*15, method*14, version*11, timestamp*8,
     &  beg_time*8,end_time*8

      open (1, file='davad_IOP12_B10',form='formatted')

      read(1, fmt='(A15, 3X, A14, 3X, A11, 3X, A8)' )
     &   provider, method, version, timestamp          ! timestamp is the time at which the data were made available

      print*,provider, method, version, timestamp

      read(1, fmt='(A8, 2X, A8, 2F9.3)')
     &   beg_time,                                     ! beginning time for processing
     &   end_time,                                     ! end time for processing
     &   lat,                                          ! latitude of reference for the vertical profile
     &   lon                                           ! longitude of reference for the vertical profile

      print*,beg_time,end_time,lat,lon

      do k=1,nplan
            read(1, fmt='( F9.3, 8(1X,F9.3) )')
     &        height(k),                               ! altitude (m)
     &        u0(k),                                   ! West-east wind component (m/s)
     &        v0(k),                                   ! South-North wind component (m/s)
     &        ref(k),                                  ! Reflectivity (dBZ)
     &        div(k),                                  ! Horizontal divergence (10-4 s-1)
     &        vort(k),                                 ! Vertical component of vorticity (10-4 s-1)
     &        stretch(k),                              ! Stretching deformation (10-4 s-1)
     &        shear(k),                                ! Shearing deformation (10-4 s-1)
     &        Vt(k)                                    ! Terminal fall velocity of hydrometeors (m/s)
! Convention : Positive downward
      enddo

      close(1)

      end