文章目录
- 1. 按年输出数据
- 2. 按月输出数据
将PIOMAS标量的二进制数据格式转化成nc格式。
1. 按年输出数据
# 按年输出数据
import numpy as np
import pandas as pd
import struct
import xarray as xr
import matplotlib.pyplot as plt
# from cartoplot import cartoplotgrids = {}dims = (120,360)for i in ['lon','lat']:grid = np.array(pd.read_csv(f'/home/mw/project/grid_data/{i}grid.dat',header=None,delim_whitespace=True))flat_grid = grid.ravel()# if i == 'lon':shaped_grid = flat_grid.reshape(dims)# else:# shaped_grid = flat_grid.reshape((360,120))grids[i] = shaped_griddef process_piomas(year):binary_dir = f'/home/mw/project/bin_data/iceprod.H{year}'################################