@如何使用python处理nc数据制作Mike风场文件
NetCDF(network Common Data Form)网络通用数据格式是一种面向数组型并适于网络共享的数据的描述和编码标准。nc格式文件被广泛应用于气象、海洋、环境、地质等专业。目前可下载的大气再分析风场数据,如CCMP、ERA5、NCEP等均提供此格式的数据。本次我们以ERA5再分析数据为例演示如何通过python制作mike的dfs2风场文件。
1. ERA5风场数据得下载
ERA5数据下载的网站为https://cds.climate.copernicus.eu/cdsapp#!/dataset/reanalysis-era5-single-levels?tab=form,变量我们选择其中的10m u-component of wind、10m v-component of wind、Surface pressure下载,下载时间、空间范围根据实际需求自定义。
2. python读取nc数据
首先导入需要用到的模块,其中前面3个为制作Mike文件相关的模块,Dataset为读取nc文件的模块,datetime为时间处理模块。
from mikeio import Dfs2
from mikeio.eum import ItemInfo, EUMUnit
from mikeio.eum import EUMType
import numpy as np
from netCDF4 import Dataset
import datetime as dt
读取文件代码如下:
file = Dataset('.\\2015uvp.nc') #文件路径及名称
time =file['time'][: