subroutine read_cmaq_netcdf(mype)
  use kinds, only : i_kind,r_kind
  use gridmod, only : nlon,nlat,grd_a
  use guess_grids, only : ges_z,ges_ps,ges_u,ges_v,ges_tv,ges_q,ges_oz,nfldsig,ifilesig,ntguessig
!  use guess_grids, only : ges_no2
  use guess_grids, only : ifilesfc
  use guess_grids, only : hrdifsig,hrdifsfc,hrdifsig_all,hrdifsfc_all
  use constants, only : zero
  use cmaqmod, only : cmaqpsol,ncidcmaqd
  use chemmod, only : laeroana_gocart
  implicit none
  integer(i_kind)mype,iret,i,k
  character(24)filename
  integer(i_kind) :: it,ncidcmaqdum
!  ntguessig=1
!  nfldsig=1
!  write(6,*)'readcmaq set nfldsig',nfldsig
!  call flush(6)
!  call create_gesfinfo
  hrdifsfc=zero
!  hrdifsig(1)=0.
!  hrdifsig_all(1)=hrdifsig(1)
  hrdifsfc_all(1)=hrdifsfc(1)
!  ifilesig(1)=0
!  write(6,*)mype,'nfldsig',nfldsig,'sz',shape(ges_z),'sp',shape(ges_ps)
!  call flush(6)
!  ifilesig=0
  ifilesfc=0
  ncidcmaqum=0
!  if(mype==0) write(6,*)'find ajl read_cmaq ',(hrdifsig(i),i=1,nfldsig)
  
  do it=1,nfldsig
    write(filename,'("sigf",i2.2)')ifilesig(it)
!    if(mype.eq.0)write(6,*)'ifilesig',it,ifilesig(it)
    call get_cmaq_netcdf(grd_a,filename,mype,ges_z(1,1,it),ges_ps(1,1,it),ges_u(1,1,1,it), &
    ges_v(1,1,1,it),ges_tv(1,1,1,it),ges_q(1,1,1,it),ges_oz(1,1,1,it),it,iret)
!    ges_v(1,1,1,it),ges_tv(1,1,1,it),ges_q(1,1,1,it),ges_oz(1,1,1,it),ges_no2(1,1,1,it),it,iret)
  end do
  return
end subroutine read_cmaq_netcdf
