subroutine updown_gc(sl,zoncoef)
        use resol_def
        implicit none
        integer k,kd,ku
        real(kind=kind_evod) sl(levs),zoncoef(levs)
        real(kind=kind_evod) xd(levs),xu(levs),sf(levs)
 

!       print *,' enter getysk_gc_fd '
 
        do k=1,levs
          kd=max(k-1,1)
          ku=min(k+1,levs)
          sf(k)=1./(sl(ku)-sl(kd))/sqrt(2.)
        enddo
 
 
          xu(levs)=zoncoef(levs)
        do k=1,levs-1
          xu(k)=zoncoef(k+1)
        enddo
        xd(1)=zoncoef(1)
        do k=2,levs
          xd(k)=zoncoef(k-1)
        enddo
 
      do k=1,levs
        zoncoef(k)=(xu(k)-xd(k))*sf(k)/101.316
      enddo
 
!       print *,' leave getysk_gc_fd '

      return
      end