program calcmaxbg

!***************************************************************************************
! abstract: use previous URMA temperature first guess analysis and guess
!           files to compute maximum temperatuer according to NDFD definition
!           (7A-7P).  Account for local time zone using gridded time zone file (fixed file).
!
!! paramaters: names of grids via grids_input
!
! inputs: *_ges_valid${hchar}.bin - last 24 guess files for temperature.  hchar = 2 character hour
!         *_anl_valid${hchar}.bin - last 24 analysis files for temperature.  hchar = 2 character hour
!
! outputs: urma2p5.${PDY}.maxT.bin, urma2p5.${PDY}.maxT.grb2
!          akurma.${PDY}.maxT.bin, akurma.${PDY}.maxT.grb2
!          hiurma.${PDY}.maxT.bin, hiurma.${PDY}.maxT.grb2
!          prurma.${PDY}.maxT.bin, prurma.${PDY}.maxT.grb2
!
! program history log:
!    2017-04-24   levine
!
!***************************************************************************************

  implicit none

  interface
     subroutine readmaxb(nx,ny,hours,threehours,dname,fnum,bnum)
       implicit none
       include 'param.incl'
       !input variables
       integer(4),intent(in)::nx,ny,fnum,bnum!,nhours,max !,singledim calc singledim on the fly
       integer(4),intent(in)::hours(nhours),threehours(8)
       character(len=7),intent(in)::dname
       
       !internal variables
       integer(4)::i,j,k,l,lun,ptnum,singledim !counters
       real(4),allocatable::tzsec(:,:),tzhour(:,:)
       real(4),allocatable::temps(:,:)
       real(4),allocatable::hourlytempbg(:,:,:),hourlytempanl(:,:,:)
       real(4),allocatable::threehour(:,:,:)
       real(4),allocatable::sixhour(:,:,:)
       real(4),allocatable::allhrlybg(:,:),allhrlyanl(:,:),allsix(:,:)
       real(4),allocatable::finalmax(:,:)
       integer(4),allocatable::alltz(:)
       real(4),allocatable::maxtonebg(:),maxttwobg(:,:),maxtwoanl(:,:),maxoneanl(:)
       character(len=20)::desc(max)
       character(len=8)::names(max)
       logical::fexist
       character(len=2)::hchar

     end subroutine readmaxb
  end interface

  include 'param.incl'

  integer(4)::g,i,j,k,l,nx,ny,lun,fnum!,singledim,ptnum
  real(4)::ds
  integer(4)::hours(nhours),threehours(8)
  integer(4)::bnum
  integer(4),parameter::grids_total=10  !process up to this number of grids
  character(len=40)::cgrid
  character(len=2)::hchar
  character(len=7)::dname
  character(len=10)::gridnames(grids_total)
  logical::fexist

  integer(4) n,grids

  namelist/blendinput/bnum
  namelist/gridsinfo/grids,gridnames

  open(15,file="blend_input",form="formatted")
  read(15,blendinput)
  close(15)

  if (bnum.ne.1.and.bnum.ne.2.and.bnum.ne.3) then
     print*, "FATAL ERROR: Invalid bnum option.  Must be 1,2, or 3.  Current value is:",bnum
     stop
  endif

  inquire(file='gridsinfo_input',exist=fexist)
  if(fexist) then
     open (15,file='gridsinfo_input',form='formatted')
     read (15,gridsinfo)
     close(15)
     if (grids.gt.10) then
        print*, "FATAL ERROR: Maximum of 10 grids (grid names) allowed in dominfo_input file.  Current value is:",grids
        stop
     endif
  else
     print*, "WARNING: no dominfo_input file present.  Defaulting to cohreswexp, akhres, hawaii, and prico grids."
     print*, "Check ex-script to make sure dominfo_input file is being written properly."
     grids=4
     gridnames(1)="cohreswexp"
     gridnames(2)="akhres"
     gridnames(3)="hawaii"
     gridnames(4)="prico"
  endif

  print*,'in calcmaxbg: grids=',grids

  do n=1, grids
     print*,'in calcmaxbg: n,gridnames(n)=',n,gridnames(n)
  enddo

  do g=1,grids
     cgrid=trim(gridnames(g))
     
     call domain_dims(cgrid,nx,ny,ds)
     print*, "In create_firstguess: nx,ny,ds=",nx,ny,ds

     if (cgrid.eq."cohresext") then
        dname="urma2p5"
     elseif (cgrid.eq."cohreswexp") then
        dname="urma2p5"
     elseif (cgrid.eq."cohres") then
        dname="urma2p5"
     elseif (cgrid.eq."akhres") then
        dname="akurma"
     elseif (cgrid.eq."prico") then
        dname="prurma"
     elseif (cgrid.eq."hawaii") then
        dname="hiurma"
     elseif (cgrid.eq."guam") then
        dname="guurma"
     endif

     if (cgrid.ne."guam") then

        hours=(/6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,0,1,2,3,4,5,6/)
        threehours=(/9,12,15,18,21,0,3,6/)

        if (dname.eq."urma2p5") then
           fnum=71
           call readmaxb(nx,ny,hours,threehours,dname,fnum,bnum)
        elseif (dname.eq."akurma") then
           fnum=73
           call readmaxb(nx,ny,hours,threehours,dname,fnum,bnum)
        elseif (dname.eq."prurma") then
           fnum=75
           call readmaxb(nx,ny,hours,threehours,dname,fnum,bnum)
        elseif (dname.eq."hiurma") then
           fnum=77
           call readmaxb(nx,ny,hours,threehours,dname,fnum,bnum)
        else
           print*, "FATAL ERROR: Invalid domain name:",dname
           stop
        endif
     else
        print*, "Skipping Guam for now..."
     endif
  enddo

end program calcmaxbg