Jeff Whitaker jswhit

## test_localsolver_multiscale.py
"""non-cycled 1d test of LETKF with multiscale localization"""
import numpy as np
from scipy.linalg import eigh, lapack, solve
from scipy.fft import rfft, irfft, rfftfreq
from argparse import ArgumentParser

# set random seed for reproducibility
#np.random.seed(42)

def syminv(C):

## test_enkfsolvers.py
"""non-cycled 2d test of LETKF solver with B or R localization"""
import numpy as np
from scipy.linalg import eigh, inv, pinvh
from scipy.special import gamma, kv
from argparse import ArgumentParser

# function definitions.

def cartdist(x1,y1,x2,y2,xmax,ymax):
    """cartesian distance on doubly periodic plane"""

## Spharmt.py
import shtns
import numpy as np
class Spharmt(object):
    """
    wrapper class for commonly used spectral transform operations in
    atmospheric models.
    """
    def __init__(self,nlons,nlats,ntrunc,rsphere=6.3712e6,gridtype='gaussian'):
        self._shtns = shtns.sht(ntrunc, ntrunc, 1, \
                                norm=shtns.sht_orthonormal+shtns.SHT_NO_CS_PHASE)

## parallel_nn_sphere.py
from __future__ import print_function # Python 3 compatible print function
"""find nearest neighbors for points on a sphere when locations to be compared
 are distributed across MPI tasks"""
# to run: mpirun -np 4 python parallel_nn_sphere.py
# requires mpi4py ('conda install mpi4py' in anaconda python)
from mpi4py import MPI
import numpy as np
import sys

# generate lat,lon values for random points on the surface of a unit sphere

## writing_netCDF.ipynb

      
              1 file
            
          
              0 forks
            
          
              0 comments
            
          
              0 stars
            
          
                jswhit
                / writing_netCDF.ipynb
            
            
              Created
              October 18, 2014 17:34
            
              
                writing netcdf ipython notebook
              
          
      Sorry, something went wrong. Reload?
      Sorry, we cannot display this file.
      Sorry, this file is invalid so it cannot be displayed.
      
          Viewer requires iframe.
      
    
## reading_netCDF.ipynb

      
              1 file
            
          
              4 forks
            
          
              2 comments
            
          
              5 stars
            
          
                jswhit
                / reading_netCDF.ipynb
            
            
              Created
              October 18, 2014 17:32
            
              
                reading netcdf ipython notebook
              
          
      Sorry, something went wrong. Reload?
      Sorry, we cannot display this file.
      Sorry, this file is invalid so it cannot be displayed.
      
          Viewer requires iframe.
      
    
## reading_netCDF.slides.html

<!DOCTYPE html>
<html>
<head>

<meta charset="utf-8" />
<meta http-equiv="X-UA-Compatible" content="chrome=1" />

<meta name="apple-mobile-web-app-capable" content="yes" />
<meta name="apple-mobile-web-app-status-bar-style" content="black-translucent" />

## basemap_pygrib.ipynb

      
              1 file
            
          
              3 forks
            
          
              0 comments
            
          
              2 stars
            
          
                jswhit
                / basemap_pygrib.ipynb
            
            
              Last active
              August 30, 2023 17:22
            
          
      Sorry, something went wrong. Reload?
      Sorry, we cannot display this file.
      Sorry, this file is invalid so it cannot be displayed.
      
          Viewer requires iframe.
      
    
## basemap_netcdf4-python.ipynb

      
              1 file
            
          
              0 forks
            
          
              0 comments
            
          
              0 stars
            
          
                jswhit
                / basemap_netcdf4-python.ipynb
            
            
              Last active
              January 13, 2019 17:44
            
              
                ipython notebook: matplotlib + netcdf4-python
              
          
      Sorry, something went wrong. Reload?
      Sorry, we cannot display this file.
      Sorry, this file is invalid so it cannot be displayed.
      
          Viewer requires iframe.
      
    
## swm_test.py
import shtns
import numpy as np

class Spharmt(object):
    """
    wrapper class for commonly used spectral transform operations in
    atmospheric models.  Provides an interface to shtns compatible
    with pyspharm (pyspharm.googlecode.com).
    """
    def __init__(self,nlons,nlats,ntrunc,rsphere,gridtype='gaussian'):
	"""non-cycled 1d test of LETKF with multiscale localization"""
	import numpy as np
	from scipy.linalg import eigh, lapack, solve
	from scipy.fft import rfft, irfft, rfftfreq
	from argparse import ArgumentParser

	# set random seed for reproducibility
	#np.random.seed(42)

	def syminv(C):
	"""non-cycled 2d test of LETKF solver with B or R localization"""
	import numpy as np
	from scipy.linalg import eigh, inv, pinvh
	from scipy.special import gamma, kv
	from argparse import ArgumentParser

	# function definitions.

	def cartdist(x1,y1,x2,y2,xmax,ymax):
	"""cartesian distance on doubly periodic plane"""
	import shtns
	import numpy as np
	class Spharmt(object):
	"""
	wrapper class for commonly used spectral transform operations in
	atmospheric models.
	"""
	def __init__(self,nlons,nlats,ntrunc,rsphere=6.3712e6,gridtype='gaussian'):
	self._shtns = shtns.sht(ntrunc, ntrunc, 1, \
	norm=shtns.sht_orthonormal+shtns.SHT_NO_CS_PHASE)
	from __future__ import print_function # Python 3 compatible print function
	"""find nearest neighbors for points on a sphere when locations to be compared
	are distributed across MPI tasks"""
	# to run: mpirun -np 4 python parallel_nn_sphere.py
	# requires mpi4py ('conda install mpi4py' in anaconda python)
	from mpi4py import MPI
	import numpy as np
	import sys

	# generate lat,lon values for random points on the surface of a unit sphere

	<!DOCTYPE html>
	<html>
	<head>

	<meta charset="utf-8" />
	<meta http-equiv="X-UA-Compatible" content="chrome=1" />

	<meta name="apple-mobile-web-app-capable" content="yes" />
	<meta name="apple-mobile-web-app-status-bar-style" content="black-translucent" />