copyright (C) 2001  MSC-RPN COMM  %%%MC2%%%

      subroutine ref_h2 ( hgeow,hgeot,hgeom,hh0,ogl,ng,nk ) 3
      implicit none
*
      integer ng,nk
      real hgeow (ng,*), hgeot (ng,*), hgeom(ng,*), hh0(ng,1,2), ogl(*)
*
      integer i,j,k,iscal(2)
      real ozt(nk),ozm(nk+1),htop
#include "h_geop.cdk"
*
*---------------------------------------------------------------------
*
      htop    = ogl(nk+1)
      iscal(1)= ogl(nk+2)
      iscal(2)= ogl(nk+3)
*     
      do k=2,nk
         ozm(k) = (ogl(k)+ogl(k-1))*0.5
         ozt(k) = ogl(k)
      end do
*
      ozm (1   ) =  ogl(1)
      ozm (nk+1) =  ogl(nk)
      ozt (1   ) =  ozm(2)*0.5
      ozt (nk  ) = (ozm(nk+1)+ozm(nk))*0.5
*
      do k=1,nk
      do i=1,ng
         hgeow(i,k) = h_geop (ogl(k),i,1)
         hgeot(i,k) = h_geop (ozt(k),i,1)
         hgeom(i,k) = h_geop (ozm(k),i,1)
      end do
      end do
*
      do i=1,ng
         hgeom(i,nk+1) = h_geop (ozm(nk+1),i,1)
      end do
*
*---------------------------------------------------------------------
      return
      end