Beispiel #1
0
# 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")