读取RAMS输出文件(RAMS-ISAN文件)的基本程序

import numpy as np
from netCDF4 import Dataset
import xarray as xr

def read_rams_data(filename):
    """读取RAMS输出文件的基本函数"""
    try:
        # 使用xarray打开文件
        ds = xr.open_dataset(filename)
        
        # 提取基本变量
        temp = ds['THETA'].values  # 位温
        pressure = ds['PI'].values  # 压力
        u_wind = ds['UP'].values   # U风分量
        v_wind = ds['VP'].values   # V风分量
        rh = ds['RV'].values       # 相对湿度
        
        # 获取坐标信息
        lat = ds['GLAT'].values
        lon = ds['GLON'].values
        z = ds['TOPT'].values      # 地形高度
        
        return {
            'temperature': temp,
            'pressure': pressure,
            'u_wind': u_wind,
            'v_wind': v_wind,
            'relative_humidity': rh,
            'latitude': lat,
            'longitude': lon,
            'elevation': z
        }
        
    except Exception as e:
        print(f"Error reading file: {e}")
        return None

# 使用示例
data = read_rams_data("rams_output.nc")
if data:
    print("Data shapes:")
    for key, value in data.items():
        print(f"{key}: {value.shape}")

这个程序的主要功能是：

• 使用xarray库打开RAMS输出的NetCDF文件
• 提取基本气象要素：位温、气压、风场、相对湿度等
• 获取地理信息：经纬度、地形高度
• 将所有数据组织成字典形式返回
• 包含错误处理机制，确保文件读取安全