      SUBROUTINE GET_SABERO_CLIM(LUN,NHFILE,SHFILE,MZ,MD,NZ,ND,
     &ALT,NHO,SHO)
      IMPLICIT NONE
      INTEGER LUN,ID,IZ,MZ,MD,NZ,ND
      CHARACTER*(*) NHFILE,SHFILE
      REAL*4 ALT(MZ),NHO(MZ,MD),SHO(MZ,MD),DUM
      OPEN(LUN,FILE=NHFILE,FORM='FORMATTED',STATUS='OLD')
      READ(LUN,*) NZ,ND
      IF((NZ.GT.MZ).OR.(ND.GT.MD)) THEN
        WRITE(*,*) 'SABER NH CLIMATOLOGY FILE EXCEEDS DIMENSIONS'
        WRITE(*,*) 'NZ,ND= ',NZ,ND
        STOP 'ERROR IN DIMENSIONS OF SABER CLIMATOLOGY DATA'
      ENDIF
      DO ID=1,ND
        DO IZ=1,NZ
          READ(LUN,*) DUM,ALT(IZ),NHO(IZ,ID)
        ENDDO
      ENDDO
      CLOSE(LUN)
      OPEN(LUN,FILE=SHFILE,FORM='FORMATTED',STATUS='OLD')
      READ(LUN,*) NZ,ND
      IF((NZ.GT.MZ).OR.(ND.GT.MD)) THEN
        WRITE(*,*) 'SABER SH CLIMATOLOGY FILE EXCEEDS DIMENSIONS'
        WRITE(*,*) 'NZ,ND= ',NZ,ND
        STOP 'ERROR IN DIMENSIONS OF SABER CLIMATOLOGY DATA'
      ENDIF
      DO ID=1,ND
        DO IZ=1,NZ
          READ(LUN,*) DUM,ALT(IZ),SHO(IZ,ID)
        ENDDO
      ENDDO
      CLOSE(LUN)
      RETURN
      END
