from __future__ import print_function import sys, os import PySCRIP as scrip from PySCRIP.config import PySCRIPConfig from netCDF4 import Dataset import numpy as np from matplotlib import pylab as plt a = PySCRIPConfig() # Mapping part ncfile = Dataset("testdata/cesmpifv1mts_cam_mapping_testdata.nc", "r") src_array = ncfile.variables["TS"][:, :] print(src_array.shape) dest = scrip.remap( src_array, "/Users/dchandan/Research/CESM/bc/cesmpifv1mts/cpl_s2/wrkdir/map_fv1-gx1_a_cesmpifv1_130513.nc", fformat="ncar-csm") print(dest.shape) print(src_array.min(), src_array.max()) print(dest.min(), dest.max()) plt.figure() plt.imshow(np.flipud(src_array[0, :])) plt.colorbar() plt.figure() plt.imshow(np.flipud(dest)) plt.colorbar() # plt.figure()
import PySCRIP as scrip from PySCRIP.config import PySCRIPConfig from netCDF4 import Dataset import numpy as np from matplotlib import pylab as plt a = PySCRIPConfig() # Mapping part ncfile = Dataset("testdata/cesmpifv1mts_pop_mapping_testdata.nc", "r") src_array = ncfile.variables["SSH"][:,:] print(src_array.shape) dest = scrip.remap(src_array, a.mapFile("cesmpifv1mts", "conservative", "gx1", "fv1")) print(dest.shape) print(src_array.min(), src_array.max()) print(dest.min(), dest.max()) plt.figure() plt.imshow(np.flipud(src_array[0,:])) plt.colorbar() # from IPython import embed # print type(dest) # embed() # plt.imshow(dest.mask) # plt.show() # print np.where(dest > src_array.max())
import PySCRIP as scrip from PySCRIP.config import PySCRIPConfig from netCDF4 import Dataset import numpy as np from matplotlib import pylab as plt a = PySCRIPConfig() # Mapping part ncfile = Dataset("testdata/cesmpifv1mts_pop_mapping_testdata.nc", "r") src_array = ncfile.variables["SSH"][:, :] print(src_array.shape) dest = scrip.remap(src_array, a.mapFile("cesmpifv1mts", "conservative", "gx1", "fv1")) print(dest.shape) print(src_array.min(), src_array.max()) print(dest.min(), dest.max()) plt.figure() plt.imshow(np.flipud(src_array[0, :])) plt.colorbar() # from IPython import embed # print type(dest) # embed() # plt.imshow(dest.mask) # plt.show() # print np.where(dest > src_array.max()) # print dest[-1,:]
# print(a.getmap(case, maptype, g1, g2).fname) # scrip.compute_remap_weights("/Users/dchandan/Research/PlioTopo/gencesmbc/prism4_N/cpl_s2/grids/fv0.9x1.25_070727.nc", # "/Users/dchandan/Development/PySCRIP/grids/ll1deg_grid.nc", # a.getmap(case, maptype, g1, g2).fname, # a.getmap(case, maptype, g2, g1).fname, # a.getmap(case, maptype, g1, g2).name, # a.getmap(case, maptype, g2, g1).name, # 2, # maptype, # normalize_opt="fracarea") scrip.compute_remap_weights( "/Users/dchandan/Research/PlioTopo/gencesmbc/prism4_N/cpl_s2/grids/prism4_N_gx1.nc", "/Users/dchandan/Development/PySCRIP/grids/ll1deg_grid.nc", a.getmap(case, maptype, g1, g2).fname, a.getmap(case, maptype, g2, g1).fname, a.getmap(case, maptype, g1, g2).name, a.getmap(case, maptype, g2, g1).name, 2, maptype, normalize_opt="fracarea", ) # # Testing part scrip.test_remap_weights(2, a.getmap(case, maptype, g1, g2).fname, "out1.nc") scrip.test_remap_weights(2, a.getmap(case, maptype, g2, g1).fname, "out2.nc")
from PySCRIP.config import PySCRIPConfig from netCDF4 import Dataset import numpy as np from matplotlib import pylab as plt a = PySCRIPConfig() # Mapping part ncfile = Dataset("testdata/cesmpifv1mts_cam_mapping_testdata.nc", "r") src_array = ncfile.variables["TS"][:, :] print(src_array.shape) dest = scrip.remap( src_array, "/Users/dchandan/Research/CESM/bc/cesmpifv1mts/cpl_s2/wrkdir/map_fv1-gx1_a_cesmpifv1_130513.nc", fformat="ncar-csm", ) print(dest.shape) print(src_array.min(), src_array.max()) print(dest.min(), dest.max()) plt.figure() plt.imshow(np.flipud(src_array[0, :])) plt.colorbar() plt.figure() plt.imshow(np.flipud(dest)) plt.colorbar()