c    this subroutine is used for retrieving all the data of a
c    specific station from wind_gen.wdb (master database) and 
c    putting them in a file

      subroutine write_wdb
     i                     (iu)
 
      integer iu
      include 'wind_db.fi'

      integer i,m

c     for the database header --- 2 lines
  100 format (a,1x,i5,1x,a,a,1x,a)
  110 format (i2,1x,i2,1x,a,1x,i3,1x,i2,1x,a,1x,i4,1x,a,1x,a,1x,a)

c     for dir(i,m) --- 17 lines
  200 format (12(1x,f4.1))

c     for scale(i,m) and shape(i,m) --- 32 lines 
  300 format (12(2x,f4.2))

c     for ratio(m) --- 1 line
  400 format (12(2x,f3.1))

c     for hr_max(m) --- 1 line
  500 format (12(3x,i2))

c     for the search_no
  800 format (i5)
  810 format (2x,i5)
  820 format (1x)

c     print the database header
         write (iu,100) sign, st_no, st_country, st_state, st_city
         write(iu,110) lati_deg, lati_min, lati_hem, long_deg, long_min,
     +                long_hem, elevation, s_date, e_date, obs

c        get dir(i,m) 
         do 40 i = 0, NDIR
            write(iu,200) (dir(i,m), m=0,NMONTHS-1)
   40    continue    
 
c        get scale(i,m) 
         do 72 i = 0, NDIR-1
            write (iu,300) (scale(i,m), m=0,NMONTHS-1)
   72    continue

c        get shape(i,m)
         do 73 i = 0, NDIR-1
            write (iu,300) (shape(i,m), m=0,NMONTHS-1)
   73    continue

c        get ratio(m) 
         write (iu,400) (ratio(m), m=0,NMONTHS-1)

c        get hr_max(m) 
         write (iu,500) (hr_max(m), m=0,NMONTHS-1)

       return
       end