!copyright (C) 2001  MSC-RPN COMM  %%%RPNPHY%%%
***s/r raddel

      subroutine raddel( del, s, sh, n, nk, nkp) 2
#include "impnone.cdk"

      integer n, nk, nkp
      real del(n,nk), s(n,nkp), sh(n,nk)
*
*Author:
*          Marc Gagnon
*
*Arguments:
*          - output -
* del      Thickness between levels.
* s        Flux sigma levels.
*
*          - input -
* sh       Sigma levels.
* n        Number of points.
* nk       Number of levels.
* nkp      Number of levels including ground.
*
*Object:
*          To compute the thickness between levels and flux levels.
*
**
      integer k,i

      do i=1,n
         s(i,1)=2.*sh(i,1)-((sh(i,1)+sh(i,2))*0.5)

c        s(i,1)=amax1(s(i,1),0.0003)

         s(i,1)=amax1(s(i,1),sh(i,1)/2.)
         s(i,nkp)=1.0
      enddo

      do k=2,nk
         do i=1,n
            s(i,k)=(sh(i,k)+sh(i,k-1))*0.5
            del(i,k-1)=s(i,k)-s(i,k-1)
         enddo
      enddo

      do i=1,n
         del(i,nk)=1.-s(i,nk)
      enddo

      end