""" 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