cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c     
c     subroutine firewx: compute fire weather danger probability 
c     Fireweather prob: 1hr apcp <= 0.01(in) and RH2m <= 15% and 10m Wind speed >= 20MPH and T2m >= 60.0 F    
c
c     Author: Binbin Zhou, May 17, 2011
c     Modification: 
c
cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
   	subroutine get_firewx (nv,rawdata,jf,iens,derv_pr,
     +               missing,wgt)

            include 'parm.inc'

C for variable table:
        Integer numvar
        Character*4 vname(maxvar)
        Integer k5(maxvar), k6(maxvar),k4(maxvar)
        Character*1 Msignal(maxvar), Psignal(maxvar)
        Integer Mlvl(maxvar), MeanLevel(maxvar,maxmlvl)
        Integer Plvl(maxvar), ProbLevel(maxvar,maxplvl)
        Integer Tlvl(maxvar)
        Character*1 op(maxvar)
        Real    Thrs(maxvar,maxtlvl)
                                            
c    for derived variables
        Character*4 dvname(maxvar)
        Integer dk5(maxvar), dk6(maxvar),dk4(maxvar)
        Character*1 dMsignal(maxvar), dPsignal(maxvar)
        Integer dMlvl(maxvar), dMeanLevel(maxvar,maxmlvl)
        Integer dPlvl(maxvar), dProbLevel(maxvar,maxplvl)
        Character*1 dop(maxvar)
        Integer dTlvl(maxvar)
        Real    dThrs(maxvar,maxtlvl)
        Integer MPairLevel(maxvar,maxmlvl,2)
        Integer PPairLevel(maxvar,maxplvl,2)
 
        common /tbl/numvar,
     +              vname,k4,k5,k6,Mlvl,Plvl,Tlvl,
     +              MeanLevel,ProbLevel,Thrs,
     +              Msignal,Psignal,op
                                                                                                                                                                                                                                                                        
        common /dtbl/nderiv,
     +              dvname,dk4,dk5,dk6,dMlvl,dPlvl,dTlvl,
     +              dMeanLevel,dProbLevel,dThrs,
     +              dMsignal,dPsignal,MPairLevel,PPairLevel,dop
                                                                                                                                                                
  
        INTEGER, intent(IN) :: nv, jf, iens
        REAL,dimension(jf,iens,maxmlvl),intent(IN) :: rawdata
        REAL,dimension(jf,maxplvl,maxtlvl),intent(INOUT) :: 
     +               derv_pr

        real FWapnt(iens),T2mapnt(iens),RH2mapnt(iens),
     +     APCPapnt(iens),U10apnt(iens),V10apnt(iens),W10apnt(iens)
        real aprob,wgt(30)

        integer, dimension(jf,iens),intent(IN) :: missing
        integer miss(iens)

        ID_T2m  = index_table(k5,k6,11,105,maxvar)   
        ID_RH2m = index_table(k5,k6,52,105,maxvar)  
        ID_APCP = index_table(k5,k6,61,2,maxvar)     
        ID_U10  = index_table(k5,k6,33,105,maxvar)   
        ID_V10  = index_table(k5,k6,34,105,maxvar)  

        write(*,*) 'In firewx -->'


        write(*,*) ID_T2m,ID_RH2m,ID_APCP,ID_U10,ID_V10


      if (ID_T2m .gt.0 .and. ID_RH2m .gt.0 .and.
     +    ID_APCP .gt.0 .and. ID_U10 .gt.0 .and.
     +    ID_V10  .gt.0 ) then

        do igrid = 1,jf

          T2mapnt=rawdata(igrid,:,ID_T2m,1)
          RH2mapnt=rawdata(igrid,:,ID_RH2m,1)
          APCPapnt=rawdata(igrid,:,ID_APCP,1)
          U10apnt=rawdata(igrid,:,ID_U10,1)
          V10apnt=rawdata(igrid,:,ID_V10,1)

          FWapnt = 0.

          do i = 1, iens

           W10apnt(i)=sqrt(U10apnt(i)*U10apnt(i) + 
     +      V10apnt(i)*V10apnt(i))
      
           TF=(T2mapnt(i)-273.15)*1.8 + 32.0        !from K to F 

            if(APCPapnt(i).le.0.10.and.RH2mapnt(i).le.15.0.and.
     +        TF.ge.60.0.and.W10apnt(i).ge.8.9) then
              FWapnt(i) = 1.0
            else
              FWapnt(i) = 0.0  
            end if

           end do

           do 30 lv=1,dPlvl(nv)
             do lt = 1, dTlvl(nv)
 
              thr1 = dThrs(nv,lt)
               miss=missing(igrid,:)

               call getprob(FWapnt,iens,
     +             thr1,thr2,dop(nv),aprob,miss,wgt)
          
               derv_pr(igrid,nv,lv,lt)=aprob 
        
               if(igrid.eq.738457) then
                write(*,*) 'Firewx@738457', dop(nv),
     +            dThrs(nv,lt),':',derv_pr(igrid,nv,lv,lt)
                write(*,'(15f8.1)') T2mapnt
                write(*,'(15f8.1)') RH2mapnt
                write(*,'(15f8.1)') APCPapnt
                write(*,'(15f8.1)') U10apnt
                write(*,'(15f8.1)') V10apnt
                write(*,'(15f8.1)') FWapnt
                write(*,'(15i8)') miss
               end if

             end do

30          continue  

         end do

        end if

        return
        end