!
!
!##################################################################
!##################################################################
!######                                                      ######
!######                SUBROUTINE INITADAS                   ######
!######                                                      ######
!######                     Developed by                     ######
!######     Center for Analysis and Prediction of Storms     ######
!######                University of Oklahoma                ######
!######                                                      ######
!##################################################################
!##################################################################
!


SUBROUTINE init3dvar(cntl_var_rh_out) 1
!
!-----------------------------------------------------------------------
!
!  PURPOSE:
!  Read in analysis variables in namelist format from standard input.
!
!  AUTHOR:
!  Jidong GAO, add 3dvar input parameter, 2001
!
!  MODIFICATION HISTORY:
!
!
!-----------------------------------------------------------------------
!
!  Variable Declarations:
!
!-----------------------------------------------------------------------
!
  IMPLICIT NONE
!
  INCLUDE 'varpara.inc'
  INCLUDE 'grid.inc'
!
  INTEGER :: nx,ny,nz          ! ARPS grid size
  INTEGER :: nstyps            ! Maximum number of soil types per grid point.
  INTEGER :: nt                ! Number of time levels of data
  INTEGER :: i                
!
!
!-----------------------------------------------------------------------
!
!  3DVAR namelists
!
!-----------------------------------------------------------------------
!
  NAMELIST /var_const/maxin

  NAMELIST /var_refil/ipass_filt,hradius,nradius_z

  NAMELIST /var_exprt/turn_chk,turn_3dda,cntl_var,cntl_var_rh

  NAMELIST /var_weight/turn_div,wei_div_h,wei_div_v

!
  INTEGER :: cntl_var_rh_out
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
!  Beginning of executable code...
!
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
!
!
!-----------------------------------------------------------------------
!
!  Assign default values to the ADAS input variables
!
!-----------------------------------------------------------------------
!
   maxin(1) = 50
   maxin(2) = 20
   maxin(3) = 30
   ipass_filt(1) = 1
   ipass_filt(2) = 1
   ipass_filt(3) = 1
   hradius(1) = 25.0
   hradius(2) = 25.0
   hradius(3) = 15.0
   nradius_z(1) = 2
   nradius_z(2) = 2
   nradius_z(3) = 2
    turn_chk = 0
    turn_3dda= 1
    cntl_var = 0
   turn_div = 0
   DO i=1,maxpass
     wei_div_h(i) = -1.0
     wei_div_v(i) = -1.0
   ENDDO
!
!-----------------------------------------------------------------------
!
!  read in ADAS namelists
!
!-----------------------------------------------------------------------
!
  READ(5, var_const,  END=350)
  WRITE(6,*) 'Namelist block 3dvar_const sucessfully read.'
  350 CONTINUE
  READ(5, var_refil,  END=360)
  WRITE(6,*) 'Namelist block 3dvar_refil sucessfully read.'
  360 CONTINUE
  READ(5, var_exprt,  END=370)
  WRITE(6,*) 'Namelist block 3dvar_exprt sucessfully read.'
  370 CONTINUE
  READ(5, var_weight,  END=390)
  WRITE(6,*) 'Namelist block 3dvar_weight sucessfully read.'
  390 CONTINUE

!
  DO i=1,maxpass
    hradius(i)= hradius(i)*1000.0/dx
    write(*,*) 'The influence radius:',hradius(i),'grid points'
  END DO

  cntl_var_rh_out=cntl_var_rh

  RETURN
END SUBROUTINE init3dvar