copyright (C) 2001  MSC-RPN COMM  %%%MC2%%%
***s/r wrgeop4
*

      subroutine wrgeop4 ( topo_l,topo_h,geobus,geosize,ni,nj ) 1
      implicit none
*
      integer geosize,ni,nj
      real topo_l(ni,nj,2),topo_h(ni,nj,2),geobus(geosize)
*
#include "lesbus.cdk"
#include "geobus.cdk"
#include "lun.cdk"
#include "yomdyn1.cdk"
*
*MODULES
*
**
      integer i,j,ier,l,m,splice
      real wk1,wk2,busd
      pointer (pawk1 , wk1 (ni,nj)), (pawk2, wk2(ni,nj)),
     $        (pabusd, busd(*))
*
*----------------------------------------------------------------------
*
      write(un_geo) geotop,geospc
*
      write(un_geo) (geonm (i,1),i=1,geotop),(geonm (i,2),i=1,geotop)
      write(un_geo) (geopar(i,1),i=1,geotop),(geopar(i,2),i=1,geotop),
     $                (geopar(i,3),i=1,geotop)
*
      splice=ni-7
      m = geosize/splice
*
      do j=1,splice
         write(un_geo) (geobus(i),i=(j-1)*m+1,j*m)
      end do
      if (splice*m+1.le.geosize)
     $     write(un_geo) (geobus(i),i=splice*m+1,geosize)
*   
      write (un_geo) maxhh01_l,maxhh02_l,maxhh01_h,maxhh02_h
      write (un_geo) topo_l
      write (un_geo) topo_h 
*----------------------------------------------------------------------
      return
      end