In [9]:
import os
os.environ['PATH'] = os.environ['PATH'] + ':' + '/home/fedor/Checkouts/esmf/DEFAULTINSTALLDIR/bin/binO/Linux.gfortran.64.mpiuni.default'
In [10]:
# This should show a missing source
!ESMF_RegridWeightGen
In [58]:
import netCDF4
src = netCDF4.Dataset('src.nc', mode='w', clobber=True)
src.createDimension('nodes', 10)
src.createDimension('elems', 13)
src.createDimension('three', 3)
mesh = src.createVariable('mesh', 'int32')
mesh.standard_name = 'mesh_topology'
mesh.dimension = 2
mesh.node_coordinates = "lon lat"
mesh.face_node_connectivity = "nv"
nv = src.createVariable("nv", "int32", dimensions=("elems", "three"))
nv.standard_name = "face_node_connectivity"
nv.start_index = 1
lon = src.createVariable("lon", "float", dimensions=("nodes",))
lon.standard_name = "longitude"
lon.units = "degrees_east"
lat = src.createVariable("lat", "float", dimensions=("nodes",))
lat.standard_name = "latitude"
lat.units = "degrees_north"
src.close()
dst = netCDF4.Dataset('dst.nc', mode='w', clobber=True)
dst.createDimension('nodes', 10)
dst.createDimension('elems', 13)
dst.createDimension('three', 3)
mesh = dst.createVariable('mesh', 'int32')
mesh.standard_name = 'mesh_topology'
mesh.dimension = 2
mesh.node_coordinates = "lon lat"
mesh.face_node_connectivity = "nv"
nv = dst.createVariable("nv", "int32", dimensions=("elems", "three"))
nv.standard_name = "face_node_connectivity"
nv.start_index = 1
lon = dst.createVariable("lon", "float", dimensions=("nodes",))
lon.standard_name = "longitude"
lon.units = "degrees_east"
lat = dst.createVariable("lat", "float", dimensions=("nodes",))
lat.standard_name = "latitude"
lat.units = "degrees_north"
dst.close()
In [60]:
%%bash
ncdump -h src.nc
rm PET*.RegridWeightGen.Log
ESMF_RegridWeightGen -s src.nc -d dst.nc --src_type UGRID --src_meshname mesh -m conserve --dst_type UGRID --dst_meshname mesh -w w.nc
cat PET*.RegridWeightGen.Log
In [ ]: