import numpy as N

def u_star(u,v,w):
    '''
Compute the friction velocity, u_star, from the timeseries of the velocity \
components u, v, and w (an nD array)
    '''
    from metpy.bl.turb.fluxes import rs as R
    rs = R(u,v,w)
    uw = rs[3]
    vw = rs[4]

    us = N.power(N.power(uw,2)+N.power(vw,2),0.25)

    return us

def theta_star(u,v,w,T):
    '''
Compute the friction temperature, theta_star, from the timeseries of the velocity \
components u, v, and w, and temperature (an nD array)
    '''
    from metpy.bl.turb.fluxes import turb_covar as TC

    ts = -TC(w,T)/u_star(u,v,w)

    return ts

def obu_length(u,v,w,T):
    '''
Compute the Obukhov Length, L, using the timeseries of the velocity \
components u, v, and w, and temperature (an nD array)
    '''
    L = N.power(u_star(u,v,w),2)*N.average(T)/(0.4*9.81*theta_star(u,v,w,T))

    return L

import numpy as N

def u_star(u,v,w):
    '''
Compute the friction velocity, u_star, from the timeseries of the velocity \
components u, v, and w (an nD array)
    '''
    from metpy.bl.turb.fluxes import rs as R
    rs = R(u,v,w)
    uw = rs[3]
    vw = rs[4]

    us = N.power(N.power(uw,2)+N.power(vw,2),0.25)

    return us

def theta_star(u,v,w,T):
    '''
Compute the friction temperature, theta_star, from the timeseries of the velocity \
components u, v, and w, and temperature (an nD array)
    '''
    from metpy.bl.turb.fluxes import turb_covar as TC

    ts = -TC(w,T)/u_star(u,v,w)

    return ts

def obu_length(u,v,w,T):
    '''
Compute the Obukhov Length, L, using the timeseries of the velocity \
components u, v, and w, and temperature (an nD array)
    '''
    from metpy.constants import g
    L = N.power(u_star(u,v,w),2)*N.average(T)/(0.4*g*theta_star(u,v,w,T))

    return L

Change to import gravity constant from constants file.
Change to import gravity constant from constants file.

git-svn-id: acf0ef94bfce630b1a882387fc03ab8593ec6522@23 150532fb-1d5b-0410-a8ab-efec50f980d4

  import numpy as N
  
  def u_star(u,v,w):
      '''
  Compute the friction velocity, u_star, from the timeseries of the velocity \
  components u, v, and w (an nD array)
      '''
      from metpy.bl.turb.fluxes import rs as R
      rs = R(u,v,w)
      uw = rs[3]
      vw = rs[4]
  
      us = N.power(N.power(uw,2)+N.power(vw,2),0.25)
  
      return us
  
  def theta_star(u,v,w,T):
      '''
  Compute the friction temperature, theta_star, from the timeseries of the velocity \
  components u, v, and w, and temperature (an nD array)
      '''
      from metpy.bl.turb.fluxes import turb_covar as TC
  
      ts = -TC(w,T)/u_star(u,v,w)
  
      return ts
  
  def obu_length(u,v,w,T):
      '''
  Compute the Obukhov Length, L, using the timeseries of the velocity \
  components u, v, and w, and temperature (an nD array)
      '''
+     from metpy.constants import g
-     L = N.power(u_star(u,v,w),2)*N.average(T)/(0.4*9.81*theta_star(u,v,w,T))
+     L = N.power(u_star(u,v,w),2)*N.average(T)/(0.4*g*theta_star(u,v,w,T))
  
      return L
  
Change to import gravity constant from constants file.
## Code Before:
import numpy as N

def u_star(u,v,w):
    '''
Compute the friction velocity, u_star, from the timeseries of the velocity \
components u, v, and w (an nD array)
    '''
    from metpy.bl.turb.fluxes import rs as R
    rs = R(u,v,w)
    uw = rs[3]
    vw = rs[4]

    us = N.power(N.power(uw,2)+N.power(vw,2),0.25)

    return us

def theta_star(u,v,w,T):
    '''
Compute the friction temperature, theta_star, from the timeseries of the velocity \
components u, v, and w, and temperature (an nD array)
    '''
    from metpy.bl.turb.fluxes import turb_covar as TC

    ts = -TC(w,T)/u_star(u,v,w)

    return ts

def obu_length(u,v,w,T):
    '''
Compute the Obukhov Length, L, using the timeseries of the velocity \
components u, v, and w, and temperature (an nD array)
    '''
    L = N.power(u_star(u,v,w),2)*N.average(T)/(0.4*9.81*theta_star(u,v,w,T))

    return L

## Instruction:
Change to import gravity constant from constants file.
## Code After:
import numpy as N

def u_star(u,v,w):
    '''
Compute the friction velocity, u_star, from the timeseries of the velocity \
components u, v, and w (an nD array)
    '''
    from metpy.bl.turb.fluxes import rs as R
    rs = R(u,v,w)
    uw = rs[3]
    vw = rs[4]

    us = N.power(N.power(uw,2)+N.power(vw,2),0.25)

    return us

def theta_star(u,v,w,T):
    '''
Compute the friction temperature, theta_star, from the timeseries of the velocity \
components u, v, and w, and temperature (an nD array)
    '''
    from metpy.bl.turb.fluxes import turb_covar as TC

    ts = -TC(w,T)/u_star(u,v,w)

    return ts

def obu_length(u,v,w,T):
    '''
Compute the Obukhov Length, L, using the timeseries of the velocity \
components u, v, and w, and temperature (an nD array)
    '''
    from metpy.constants import g
    L = N.power(u_star(u,v,w),2)*N.average(T)/(0.4*g*theta_star(u,v,w,T))

    return L
