c
c------------------------------------------------------------------
c Begin file /home/sting/foster/timegcm/mksrc/getglbm.f
c------------------------------------------------------------------
c
      subroutine getglbm(lu)
c
c Get global mean history volume and define fields:
c     common /glbmcom/ glbmfld(kmxglbm,nglbmfld)
c 7/91: kmxglbm=37, nglbmfld=96
c
      include 'glbm.h'
      character*80 errmsg
      character*8 locflnm
      data locflnm /'glbmvol '/
c
c Get global mean history from which to assign fields to timegcm:
c (use replace option to force read from mss even if local disk file
c  is there)
c
      write(6,"(' ')")
      write(6,"('getglbm calling msread for glbmvol=',a)") glbmvol
c     call msread(ier,locflnm,glbmvol,' ','REPLACE')
      call msread(ier,locflnm,glbmvol,' ',' ')
      if (ier.eq.3) then
        write(6,"('getglbm msread found file on disk -- change this',
     +   ' later -- use option REPLACE')")
      elseif (ier.ne.0) then
        call mserror(errmsg)
        write(6,"('getglbm: msread error ',i3,' errmsg=',a,
     +    /' glbmvol=',a)") ier,errmsg,glbmvol
        stop 'getglbm'
      else
        write(6,"('getglbm: successful msread of mss file ',a,' to ',
     +    /'  unicos file ',a)") glbmvol,locflnm
      endif
      call assgn(locflnm,lu)
c
c Read global mean fields:
c
      read(lu) glbmfld
c
c Set glbm (and timegcm) vertical grid:
c
      do k=1,kmxglbm
        glbmzp(k) = glbmzp1 + (k-1)*glbmdelz
      enddo

      return
      end