文章目錄
- 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}'################################