解决ERA5资料驱动WRF real 时报错smooth_cg_topo=T但缺少SOILHGT的问题

离域

运行real.exe时报错
--- ERROR: NML option smooth_cg_topo=T
            But found no soil elevation / terrain / topography data in metgrid files
            The field SOILHGT is required when smoothing the CG topography on d01

在met_em中缺少了SOILHGT变量，基本又是ERA5 single layer下载的变量少了导致的问题。

参考气象家园：进行UNGRIB操作所需的物理场
https://bbs.06climate.com/forum.php?mod=viewthread&tid=30940&fromuid=73002

三维数据（例如，气压层的数据）
温度
风场u,v分量
位势高度
相对湿度（/比湿）


二维数据
地面气压
平均海平面气压
地表温度
2米温度
2米相对湿度或比湿
10风场u,v分量

如果需要用到masked，建议提供LANDSEA场。
如果在进行物理过程参数化时需要用到Noah LSM陆面模式，需要提供两层Soil Temperature和Soil Moisture。
如果模式运行时需要更新SST，需要提供SST。
建议提供Soil Height（可选），这样就可以调整土壤温度达到与WRF模式地形高度向适应。
建议提供SNOW场。
建议对高纬度冬季个例，提供SEAICE场。



在ERA5 single layer的数据尤其是soil数据中，没有找到soil height，询问chatgpt：

下面的回答是确认METGRID.TBL中土壤高度的相关设置，默认的就有，不截取了。

在ERA5 single layer资料中添加 Other-Geopotential 后重新下载数据并ungrib和metgrid，新的met_em*.nc中产生了SOILHGT变量，后续步骤也顺利进行。

ps：现在驱动WRF用的ERA5 single layer变量，由于不跑冬季个例没选SEAICE场
        'variable': [
            '10m_u_component_of_wind', '10m_v_component_of_wind', '2m_dewpoint_temperature',
            '2m_temperature', 'geopotential', 'land_sea_mask',
            'mean_sea_level_pressure', 'sea_surface_temperature', 'skin_temperature',
            'snow_albedo', 'snow_depth', 'snowfall',
            'soil_temperature_level_1', 'soil_temperature_level_2', 'soil_temperature_level_3',
            'soil_temperature_level_4', 'soil_type', 'surface_pressure',
            'total_precipitation', 'volumetric_soil_water_layer_1', 'volumetric_soil_water_layer_2',
            'volumetric_soil_water_layer_3', 'volumetric_soil_water_layer_4',
        ],