!#define VERBOSE
subroutine gsisub(mype,init_pass,last_pass)
!$$$ subprogram documentation block
! . . . .
! subprogram: gsisub high level driver for gridpoint
! statistical interpolation
! prgmmr: parrish org: np22 date: 1992-09-16
!
! abstract: This subroutine embodies the main driver for the gsi analysis
! code. This routine
! a) get analysis time and vertical coordinate information
! b) generate grid information
! c) initialize arrays containing information regarding use
! of radiance, precipitation rate, and ozone observations
! d) decomponses the domain into subdomains
! e) perform various initializations (subdomain variables,
! mpi communications, random number for pcp assimilation)
! f) if (.not.regional) iniitialize spectral <--> grid transforms
! g) call outer/inner loop driver
! h) deallocate arrays
!
! program history log:
! 1991-xx-xx parrish/derber
! 1991-12-10 parrish/derber fixed coding error in near sfc anal
! 1992-09-14 derber improved version of global analysis
! 1998-05-15 weiyu yang mpp version
! 1999-08-24 derber, j., treadon, r., yang, w., first frozen mpp version
! 2003-10-31 kleist, d. add capability to run hybrid or sigma
! 2003-12-22 derber, j. comments
! 2004-06-21 treadon - update documentation
! 2004-07-28 treadon - add only to module use, add intent in/out
! 2004-12-29 treadon - replace wrf_nmm_regional and wrf_mass_regional logicals
! with regional
! 2005-02-16 xu li - add sub destroy_sst_an
! 2005-03-07 dee - support gmao model interface
! 2005-04-18 treadon - remove destroy_sst_an
! 2005-05-27 pondeca - bypass radinfo_ and pcpinfo_write when twodvar_regional=.t.!
! 2005-05-25 guo - added interfaces to handle GMAO first guess gridded fields
! 2005-07-25 treadon - remove redundant call to gengrid_vars
! 2005-09-08 derber - modify to use input group time window and simplify data set handling
! 2005-10-18 treadon - remove obs_load and dload, move deter_subdomain before read_obs
! 2005-11-04 guo - moved call to gengrid_vars to ncep initialization part
! 2005-11-28 derber - move read_obs to glbsoi
! 2005-12-09 guo - minor comment change
! 2006-01-10 treadon - consolidate query guess file code in gesinfo
! 2006-04-20 kistler - moved conv_read from read_obs here to parallel other *info modules
! 2006-04-21 parrish - changes for new processing of level 2 radar wind data
! 2007-03-15 todling - merged in da Silva/Cruz ESMF changes
! 2007-10-03 todling - add observer call
! 2009-01-28 todling - update observer calling procedure
! 2009-08-19 guo - #ifdef out destroy_gesfinfo() call for multi-pass observer.
!
! input argument list:
! mype - mpi task id
!
! comments:
!
! attributes:
! language: f90
! machine: ibm RS/6000 SP
!
!$$$
use kinds, only: i_kind
use obsmod, only: iadate,lobserver
use observermod, only: observer_init,observer_run,observer_finalize
use gridmod, only: twodvar_regional,regional,&
create_grid_vars,create_mapping,init_subdomain_vars,&
destroy_mapping,destroy_grid_vars
use gridmod, only: wrf_mass_regional,wrf_nmm_regional,nems_nmmb_regional
use mpimod, only: npe,mpi_comm_world,ierror
use radinfo, only: radinfo_read
use pcpinfo, only: pcpinfo_read,create_pcp_random,&
destroy_pcp_random
use convinfo, only: convinfo_read
use ozinfo, only: ozinfo_read
use read_l2bufr_mod, only: radar_bufr_read_all
use oneobmod, only: oneobtest,oneobmakebufr
#ifndef HAVE_ESMF
use guess_grids, only: destroy_gesfinfo
#endif
use mpeu_util, only: die,tell
implicit none
! Declare passed variables
integer(i_kind),intent(in ) :: mype
logical ,intent(in) :: init_pass
logical ,intent(in) :: last_pass
#ifdef VERBOSE
call tell('gsisub','entered ..')
call tell('gsisub','init_pass =',init_pass)
call tell('gsisub','last_pass =',last_pass)
call tell('gsisub','iadate(1)=',iadate(1))
call tell('gsisub','iadate(2)=',iadate(2))
call tell('gsisub','iadate(3)=',iadate(3))
call tell('gsisub','iadate(4)=',iadate(4))
call tell('gsisub','iadate(5)=',iadate(5))
#endif
#ifndef HAVE_ESMF
! Allocate grid arrays.
call create_grid_vars
! Get date, grid, and other information from model guess files
call gesinfo(mype)
! If single ob test, create prep.bufr file with single ob in it
if (oneobtest) then
if(mype==0)call oneobmakebufr
call mpi_barrier(mpi_comm_world,ierror)
end if
! Create analysis subdomains and initialize subdomain variables
call create_mapping(npe)
call deter_subdomain(mype)
call init_subdomain_vars
#endif /* !HAVE_ESMF */
! Process any level 2 bufr format land doppler radar winds and create radar wind superob file
if(wrf_nmm_regional.or.wrf_mass_regional.or.nems_nmmb_regional) call radar_bufr_read_all(npe,mype)
! Read info files for assimilation of various obs
if (init_pass) then
if (.not.twodvar_regional) then
call radinfo_read
call ozinfo_read(mype)
call pcpinfo_read(mype)
endif
call convinfo_read(mype)
#ifdef VERBOSE
call tell('gsisub','returned from convinfo_read()')
#endif
endif
#ifndef HAVE_ESMF
! Set communicators between subdomain and global/horizontal slabs
call init_commvars(mype)
#endif /* HAVE_ESMF */
! Compute random number for precipitation forward model.
if(init_pass) then
call create_pcp_random(iadate,mype)
#ifdef VERBOSE
call tell('gsisub','returned from create_pcp_random()')
#endif
endif
! Complete setup and execute external and internal minimization loops
#ifdef VERBOSE
call tell('gsisub','lobserver=',lobserver)
#endif
if (lobserver) then
if(init_pass) call observer_init()
#ifdef VERBOSE
call tell('gsisub','calling observer_run()')
#endif
call observer_run(init_pass=init_pass,last_pass=last_pass)
#ifdef VERBOSE
call tell('gsisub','returned from observer_run()')
#endif
if(last_pass) call observer_finalize()
#ifndef HAVE_ESMF
call destroy_gesfinfo() ! paired with gesinfo()
#endif
else
call glbsoi(mype)
endif
if(last_pass) then
! Deallocate arrays
call destroy_pcp_random
#ifndef HAVE_ESMF
call destroy_mapping
call destroy_grid_vars
#endif /* HAVE_ESMF */
endif
#ifdef VERBOSE
call tell('gsisub','exiting ..')
#endif
! End of gsi driver routine
return
end subroutine gsisub