      SUBROUTINE RR_vuy_Qiao_init(id)
      IMPLICIT NONE
      INCLUDE 'ferret_cmn/EF_Util.cmn'
      INTEGER id

      CALL ef_set_desc             (id,'v du/dy from Qiao (1997)')
      CALL ef_set_num_args         (id, 3)
      CALL ef_set_axis_inheritance (id, IMPLIED_BY_ARGS, 
     . IMPLIED_BY_ARGS,IMPLIED_BY_ARGS,IMPLIED_BY_ARGS)

      CALL ef_set_arg_name         (id, ARG1, 'v4') 
      CALL ef_set_arg_desc         (id, ARG1, 'v4')
      CALL ef_set_arg_unit         (id, ARG1, 'cm/sec')
      CALL ef_set_arg_type         (id, ARG1, FLOAT_ARG)
      CALL ef_set_axis_influence   (id, ARG1, YES, YES, YES, YES)

      CALL ef_set_arg_name         (id, ARG2, 'u2')
      CALL ef_set_arg_desc         (id, ARG2, 'u2')
      CALL ef_set_arg_unit         (id, ARG2, 'cm/sec')
      CALL ef_set_arg_type         (id, ARG2, FLOAT_ARG)
      CALL ef_set_axis_influence   (id, ARG2, YES,  NO, YES, YES)

      CALL ef_set_arg_name         (id, ARG3, 'u5')
      CALL ef_set_arg_desc         (id, ARG3, 'u5')
      CALL ef_set_arg_unit         (id, ARG3, 'cm/sec')
      CALL ef_set_arg_type         (id, ARG3, FLOAT_ARG)
      CALL ef_set_axis_influence   (id, ARG3, YES,  NO, YES, YES)

      END


      SUBROUTINE RR_vuy_Qiao_compute(id,arg_1,arg_2,arg_3,result)
      IMPLICIT NONE
      INCLUDE 'ferret_cmn/EF_Util.cmn'
      INCLUDE 'ferret_cmn/EF_mem_subsc.cmn'

      INTEGER        id
      REAL           arg_1 (mem1lox:mem1hix, mem1loy:mem1hiy, 
     .                      mem1loz:mem1hiz, mem1lot:mem1hit)
      REAL           arg_2 (mem2lox:mem2hix, mem2loy:mem2hiy, 
     .                      mem2loz:mem2hiz, mem2lot:mem2hit)
      REAL           arg_3 (mem3lox:mem3hix, mem3loy:mem3hiy, 
     .                      mem3loz:mem3hiz, mem3lot:mem3hit)
      REAL           result (memreslox:memreshix,memresloy:memreshiy,
     .                       memresloz:memreshiz,memreslot:memreshit)
      INTEGER        res_lo_ss (4)
      INTEGER        res_hi_ss (4)
      INTEGER        res_incr  (4)
      INTEGER        arg_lo_ss (4,EF_MAX_ARGS)
      INTEGER        arg_hi_ss (4,EF_MAX_ARGS)
      INTEGER        arg_incr  (4,EF_MAX_ARGS)

      INTEGER        i,  j,  k,  l
      INTEGER        i1, j1, k1, l1
      INTEGER        i3, j3, k3, l3
      INTEGER        i2, j2, k2, l2

      REAL           dy

      CALL ef_get_res_subscripts (id, res_lo_ss, res_hi_ss, res_incr)
      CALL ef_get_arg_subscripts (id, arg_lo_ss, arg_hi_ss, arg_incr)

      dy = 212.7e5  ! 2 degrees of latitude in cm

      i  = res_lo_ss(X_AXIS)
      i1 = arg_lo_ss(X_AXIS,ARG1) 
      i2 = arg_lo_ss(X_AXIS,ARG2) 
      i3 = arg_lo_ss(X_AXIS,ARG3) 

      j  = res_lo_ss(Y_AXIS)
      j1 = arg_lo_ss(Y_AXIS,ARG1) 
      j2 = arg_lo_ss(Y_AXIS,ARG2) 
      j3 = arg_lo_ss(Y_AXIS,ARG3) 

      l1 = arg_lo_ss(T_AXIS,ARG1) 
      l2 = arg_lo_ss(T_AXIS,ARG2) 
      l3 = arg_lo_ss(T_AXIS,ARG3) 
      DO l = res_lo_ss(T_AXIS), res_hi_ss(T_AXIS)

         k1 = arg_lo_ss(Z_AXIS,ARG1)
         k2 = arg_lo_ss(Z_AXIS,ARG2)
         k3 = arg_lo_ss(Z_AXIS,ARG3)

         DO k = res_lo_ss(Z_AXIS), res_hi_ss(Z_AXIS)

           !  ARG_1 = v4 
           !  ARG_2 = u2
           !  ARG_3 = u5

           result(i,j,k,l) = ARG_1(i1,j1,k1,l1) *
     .        (ARG_3(i3,j3,k3,l3)-ARG_2(i2,j2,k2,l2))/dy

           k1 = k1 + arg_incr(Z_AXIS,ARG1)
           k2 = k2 + arg_incr(Z_AXIS,ARG2)
           k3 = k3 + arg_incr(Z_AXIS,ARG3)
         ENDDO 

         l1 = l1 + arg_incr(T_AXIS,ARG1)
         l2 = l2 + arg_incr(T_AXIS,ARG2)
         l3 = l3 + arg_incr(T_AXIS,ARG3)
      ENDDO               
      END				






