sgdecker

sgdecker

 GDATTIM  = f0
 GLEVEL   = 500
 GVCORD   = PRES
 GFUNC    = vor(wnd)
 GDFILE   = $MODEL/gfs/18072300_gfs003.gem
 CINT     = 4
 LINE     = 1//2
 MAP      = 8
 MSCALE   = 0
 TITLE    = 1

sgdecker

from datetime import datetime
import cartopy.crs as ccrs
import cartopy.feature as cfeature
import matplotlib.pyplot as plt
import numpy as np
import xarray as xr
from metpy.units import units
import metpy.calc as mpcalc

# Grab a recent GFS run on Grid 3

sgdecker

import numpy as np
import matplotlib.pyplot as plt
from metpy.units import units
from metpy.calc import saturation_vapor_pressure
import timeit

def goffgratch46(Tin):
    """Computes saturation vapor pressure (mb) from temperature using the
    Goff and Gratch (1946) formulation as quoted by Alduchov and Eskridge
    (1996)"""