"""
Test conservative regridding of strange grids
inmcm grid
$Id: testinmcmGrid.py 2472 2012-10-12 13:10:52Z pletzer $
"""
import cdms2
import numpy as P
import ESMP
import glob
import os
import unittest
ESMP.ESMP_LogSet(True)
class TestGrid(unittest.TestCase):
def setUp(self):
dir = "../cdat_data/"
filename = dir + "so_Omon_inmcm4_1pctCO2_r1i1p1_209001-209412_2timesteps.nc"
g = cdms2.open(filename)
self.so = g('so')[0, 0, ...]
gLat = cdms2.createGaussianAxis(64)
deltaLon = (360 / 128.)
gLon = cdms2.createUniformLongitudeAxis(0, 128, deltaLon)
self.gaussGrid = cdms2.grid.createGenericGrid(gLat[:], gLon[:],
gLat.getBounds(),
gLon.getBounds())
def test1(self):
soN = self.so.regrid(self.gaussGrid,
from tvariable import asVariable, createVariable, isVariable
from mvSphereMesh import SphereMesh
from mvBaseWriter import BaseWriter
from mvVsWriter import VsWriter
from mvVTKSGWriter import VTKSGWriter
from mvVTKUGWriter import VTKUGWriter
from mvCdmsRegrid import CdmsRegrid
# Gridspec is not installed by default so just pass on if it isn't installed
try:
from gsStaticVariable import StaticFileVariable
from gsTimeVariable import TimeFileVariable
except:
pass
from restApi import esgfConnection, esgfDataset, FacetConnection
MV = MV2
ESMP_HAS_BEEN_INITIALIZED = False
if not ESMP_HAS_BEEN_INITIALIZED:
try:
import ESMP
ESMP.ESMP_Initialize(ESMP.ESMP_LOGKIND_NONE)
# this turns off the PET file logs
ESMP.ESMP_LogSet(False)
ESMP_HAS_BEEN_INITIALIZED = True
except:
pass