MODULE module_wind_generic

  IMPLICIT NONE

  TYPE windturbine_specs
     INTEGER id             
     REAL    lat, lon       
     REAL    i, j           
     REAL    hubheight      
     REAL    diameter       
     REAL    stdthrcoef     
     REAL    power          
     REAL    cutinspeed     
     REAL    cutoutspeed    
  END TYPE windturbine_specs

  TYPE(windturbine_specs), TARGET, ALLOCATABLE, DIMENSION(:) :: windturbines
  INTEGER :: nwindturbines

  INTEGER, PARAMETER :: WIND_TURBINES_OFF      = 0
  INTEGER, PARAMETER :: WIND_TURBINES_IDEAL    = 1
  INTEGER, PARAMETER :: WIND_TURBINES_FROMLIST = 2

  INTEGER windspec

  LOGICAL , EXTERNAL      :: wrf_dm_on_monitor

CONTAINS

  SUBROUTINE read_windturbines_in








    IMPLICIT NONE

    CHARACTER*256  fname, message
    CHARACTER*512 inline
    INTEGER i,istat
    INTEGER id
    INTEGER n,lineno,ig,jg
    REAL  lat,lon,hubheight,diameter,stdthrcoef,power,cutinspeed,cutoutspeed

    CALL nl_get_windturbines_spec( 1, fname )
    windspec = WIND_TURBINES_OFF
    IF ( TRIM(fname) .EQ. "none" ) THEN
      RETURN
    ELSE IF ( TRIM(fname) .EQ. "ideal" ) THEN
     
     
      windspec = WIND_TURBINES_IDEAL
    ELSE
      
      
      
      IF ( wrf_dm_on_monitor() ) THEN
        OPEN(file=TRIM(fname),unit=19,FORM='FORMATTED',STATUS='OLD',IOSTAT=istat)
        IF ( istat .EQ. 0 ) THEN
          
          n = 0
          DO WHILE (.true.)
            READ(19,'(A256)',END=30)inline
            IF ( index(inline,'!') .EQ. 0 ) n = n + 1
          ENDDO
 30       CONTINUE
          nwindturbines = n
          IF ( .NOT. ALLOCATED(windturbines) ) ALLOCATE(windturbines(nwindturbines))
          REWIND(19)
          i = 1 
          lineno = 0
          DO WHILE (.true.)
            lineno = lineno + 1
            READ(19,'(A256)',END=120)inline
            IF ( i .LE. nwindturbines .AND. index(inline,'!') .EQ. 0 ) THEN
              READ(inline,*,ERR=130)id,lat,lon,hubheight,diameter,stdthrcoef,power,cutinspeed,cutoutspeed
              windturbines(i)%id = id
              windturbines(i)%lat = lat
              windturbines(i)%lon = lon
              windturbines(i)%i = -999   
              windturbines(i)%j = -999   
              windturbines(i)%hubheight = hubheight
              windturbines(i)%diameter = diameter
              windturbines(i)%stdthrcoef = stdthrcoef
              windturbines(i)%power = power
              windturbines(i)%cutinspeed = cutinspeed
              windturbines(i)%cutoutspeed = cutoutspeed
              i = i + 1
            ENDIF
          ENDDO
 120      CONTINUE
          CLOSE(19)
          GOTO 150
 130      CONTINUE
          CLOSE(19)   
          istat = 150150
          GOTO 150
        ENDIF
      ENDIF
 150  CONTINUE
      CALL wrf_dm_bcast_integer(istat,1) 
      IF ( istat .NE. 0 ) THEN
        WRITE(message,*)'Unable to open or read ',TRIM(fname),'. Proceeding without wind-turbine parameterization.'
        CALL wrf_message(message)
        IF ( istat .EQ. 150150 ) THEN
          WRITE(message,*)'Perhaps bad syntax at line ',lineno,' of ',TRIM(fname)
          CALL wrf_message(message)
        ENDIF
        IF ( ALLOCATED(windturbines) ) DEALLOCATE(windturbines)
        RETURN
      ENDIF
      CALL wrf_dm_bcast_integer(nwindturbines,1) 
      IF ( .NOT. wrf_dm_on_monitor() ) THEN
        IF ( .NOT. ALLOCATED(windturbines) ) ALLOCATE(windturbines(nwindturbines))
      ENDIF
      DO i = 1, nwindturbines 
        CALL wrf_dm_bcast_integer(windturbines(i)%id,1) 
        CALL wrf_dm_bcast_real(windturbines(i)%lat,1) 
        CALL wrf_dm_bcast_real(windturbines(i)%lon,1) 
        CALL wrf_dm_bcast_real(windturbines(i)%hubheight,1) 
        CALL wrf_dm_bcast_real(windturbines(i)%diameter,1) 
        CALL wrf_dm_bcast_real(windturbines(i)%stdthrcoef,1) 
        CALL wrf_dm_bcast_real(windturbines(i)%power,1) 
        CALL wrf_dm_bcast_real(windturbines(i)%cutinspeed,1) 
        CALL wrf_dm_bcast_real(windturbines(i)%cutoutspeed,1) 
      ENDDO
      windspec = WIND_TURBINES_FROMLIST
      RETURN
    ENDIF
  END SUBROUTINE read_windturbines_in

  SUBROUTINE init_module_wind_generic
    IMPLICIT NONE
    CALL read_windturbines_in 
  END SUBROUTINE init_module_wind_generic

END MODULE module_wind_generic