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

      subroutine xpyp_m (xp,yp,lonref,latref,ofi,ofj,con,ni,nj) 5
      implicit none
*     
      integer ni,nj,ofi,ofj
      real*8 xp(ni),yp(nj),lonref,latref,con
**
#include "consdyn_8.cdk"
*
      integer i,j
      real*8 c2,c3,d2r
*
*---------------------------------------------------------------------
*
      d2r = pi_8 / 180.
      c2  = 2. / d2r
      c3  = con * d2r  
*    
      do i=1,ni
         xp(i) = lonref + dble(i+ofi-1) * con
         if (xp(i).lt.0.) xp(i)=xp(i)+360.
      end do
*
      do j=1,nj
         yp(j) = c2*atan(tan((latref+90.)/c2)*exp(c3*dble(j+ofj-1)))-90.
      end do
*
*---------------------------------------------------------------------
*
      return
      end