subroutine updown_gc(sl,zoncoef)
        use resol_def
        implicit none
        integer k,kd,ku
        real(kind=kind_evod) ysk(levs),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