下面的页面概述了一些已与FVCOM一起使用的有用数据集。其中包括有关如何使用FVCOM matlab工具箱获取和处理这些数据的说明。
NCEP再分析II
网站:http : //www.esrl.noaa.gov/psd/data/gridded/data.ncep.reanalysis2.html
MATLAB FVCOM工具箱功能:get_NCEP_forcing.m
。有关示例用法get_NCEP_forcing.m
,请参见NCEP配置部分。
符合Office统一模型
网站:http : //badc.nerc.ac.uk/browse/badc/ukmo-um。需要用户名/密码(每个项目每个人一个)。
文档:http : //badc.nerc.ac.uk/view/badc.nerc.ac.uk__ATOM__dataent_UM
2011年前的模型输出对具有唯一ID(存储代码-参见下表)的每个变量使用单独的文件。对于2011年以后的版本,每个模型运行的输出均存储为单个文件。这样,隐藏代码在很大程度上被变量名称所取代。一些文档在这里:
http://badc.nerc.ac.uk/view/badc.nerc.ac.uk__ATOM__DE_5b70af60-2097-11e3-a39b-00163e251233(NAE输出)
http://badc.nerc.ac.uk/view/badc.nerc.ac.uk__ATOM__DE_f693d60c-208d-11e3-ad34-00163e251233(全局输出)
我(皮埃尔)我目前加载使用的下列变量HEATING_ON
或HEATING_CALCULATED
为2011年以后的数据:
变量名 | 长名 | 外型尺寸 | 单位 |
longitude | – | latitude | 度数 |
latitude | – | latitude | 度数 |
x-wind | x风分量(相对于网格) | t_1 ,ht ,latitude ,longitude | 毫秒-1 |
y-wind | y风分量(相对于电网) | t_1 ,ht ,latitude ,longitude | 毫秒-1 |
solar | 净短波辐射通量 | t_3 ,surface ,latitude ,longitude | 宽m -2 |
longwave | 净长波辐射通量 | t_3 ,surface ,latitude ,longitude | 宽m -2 |
sh | 显热通量 | t_1 ,surface ,latitude ,longitude | Wm -2 |
lh | 潜热通量 | t_1 ,surface ,latitude ,longitude | Wm -2 |
rh | 相对湿度 | t_1 ,ht ,latitude ,longitude | % |
precip | 总降水率 | t_1 ,surface ,latitude ,longitude | 千克米-2 s -1 |
p_4 要么 p_3 | 压力 | t_1 ,msl ,latitude ,longitude | 帕 |
temp | 温度T | t_1 ,surface ,latitude ,longitude | 摄氏温度 |
降水需要转换为m / s(除以淡水密度1000kg m -3)。与NCEP数据一样,我正在计算潜热的蒸发量,并从降水量中减去它。然后将蒸发量设置为零,并使用沉淀蒸发量作为沉淀变量。
变量清单
缩略语 | 长名 |
美联社 | 分析,压力水平 |
上午 | 分析,模型级别 |
FP | 预测,压力水平 |
调频 | 预测,模型级别 |
美联社
藏 | STASHmaster说明 | 单位 | 有效期自 | 有效 | CF标准名称 | 用过的? |
23 | TIMESTEP后打雪量KG / M2 | 千克米2 | 503 | snowfall_amount | ||
23 | 陆上雪量TSTP KG / M2 | 千克米2 | 505 | snowfall_amount | ||
24 | 时间戳后的表面温度 | ķ | surface_temperature | |||
25 | 时间戳后的边界层深度 | 米 | atmosphere_boundary_layer_thickness | |||
26 | 时间戳后的粗糙度 | 米 | surface_roughness_length | |||
30 | LAND MASK(逻辑:LAND = TRUE) | 1个 | 407 | land_binary_mask | ||
30 | LAND MASK(无光晕)(LAND = TRUE) | 1个 | 500 | land_binary_mask | ||
31 | 跨海后海冰的分形 | 1个 | 503 | sea_ice_area_fraction | ||
31 | 时间戳后的海冰分数 | 1个 | 502 | sea_ice_area_fraction | ||
33 | 地形(/ STRAT LOWER BC) | 米 | surface_altitude | |||
1201 | 净落面SW流量:仅SW TS | 宽m-2 | surface_net_downward_shortwave_flux | 是啊 | ||
1235 | 向下总表面通量 | 宽m-2 | surface_downwelling_shortwave_flux_in_air | 是啊 | ||
2201 | 净落面LW RAD通量 | 宽m-2 | surface_net_downward_longwave_flux | 是啊 | ||
2207 | 向下LW RAD通量:表面 | 宽m-2 | 407 | surface_downwelling_longwave_flux_in_air | 是啊 | |
3217 | 带M2的表面热通量 | 宽m-2 | 501 | surface_upward_sensible_heat_flux | 是啊 | |
3217 | 带M2的表面和B层热流 | 宽m-2 | 500 | |||
3223 | SURF&BL TOTL水分通量KG / M2 / S | 千克m-2 s-1 | 500 | upward_water_vapor_flux_in_air | ||
3223 | 表面总水分通量KG / M2 / S | 千克m-2 s-1 | 501 | surface_upward_water_flux | ||
3224 | WIND MIX EN'GY FL到SEA:SEA MN W / M2 | 宽m-2 | 503 | wind_mixing_energy_flux_into_sea_water | ||
3224 | W / M2混合风入海 | 宽m-2 | 502 | wind_mixing_energy_flux_into_sea_water | ||
3225 | 10米风量U-COMP | 1秒 | 500 | eastward_wind | 是啊 | |
3225 | 10米风量U-COMP | 1秒 | 500 | x_wind | 是啊 | |
3225 | 10米风U-COMP B网格 | 1秒 | 501 | eastward_wind | 是啊 | |
3225 | 10米风U-COMP B网格 | 1秒 | 501 | x_wind | 是啊 | |
3226 | 10米风速V-COMP | 1秒 | 500 | northward_wind | 是啊 | |
3226 | 10米风速V-COMP | 1秒 | 500 | y_wind | 是啊 | |
3226 | 10米风速V-COMP B网格 | 1秒 | 501 | northward_wind | 是啊 | |
3226 | 10米风速V-COMP B网格 | 1秒 | 501 | y_wind | 是啊 | |
3234 | W / M2表面潜热通量 | 宽m-2 | 407 | surface_upward_latent_heat_flux | 是啊 | |
3236 | 温度为150万 | ķ | air_temperature | 是啊 | ||
3237 | 150万的特定湿度 | 1个 | 407 | specific_humidity | ||
3245 | 相对湿度150万 | % | relative_humidity | 是啊 | ||
3248 | 雾度为150万 | 407 | ||||
3250 | 露点为1.5M(K) | ķ | 407 | dew_point_temperature | ||
3281 | TS S / M之后用于SO4的抗性S | 401 | 407 | |||
3281 | VIS AT 1.5M(包括沉淀)M | 606 | ||||
3281 | 可见度为1.5MM | 503 | 601 | |||
3316 | SNOW-ADJSTD标题K上的表面温度 | 405 | 407 | |||
3316 | 标题K上的表面温度 | 503 | ||||
3321 | 无冰砖上的冠层水KG / M2 | 405 | 407 | |||
3321 | 瓷砖雨水KG / M2 | 503 | ||||
3463 | 风 | 606 | ||||
3476 | 合并边界层类型 | 606 | ||||
4201 | 大雨量KG / M2 / TS | 千克米2 | stratiform_rainfall_amount | |||
4202 | 大雪量KG / M2 / TS | 千克米2 | stratiform_snowfall_amount | |||
4203 | 大型降雨率KG / M2 / S | 千克m-2 s-1 | stratiform_rainfall_rate | |||
4204 | 大空速KG / M2 / S | 千克m-2 s-1 | stratiform_snowfall_rate | |||
5201 | 对流雨量KG / M2 / TS | 千克米2 | convective_rainfall_amount | |||
5202 | 对流雪量KG / M2 / TS | 千克米2 | convective_snowfall_amount | |||
5205 | 对流降雨率KG / M2 / S | 千克m-2 s-1 | 407 | convective_rainfall_rate | ||
5206 | 对流降雪率KG / M2 / S | 千克m-2 s-1 | 407 | convective_snowfall_rate | ||
5207 | 对流云基的压力 | 帕 | 407 | air_pressure_at_convective_cloud_base | ||
5208 | 对流云顶的压力 | 帕 | 407 | air_pressure_at_convective_cloud_top | ||
5210 | 国际民航组织对流云基础 | 407 | ||||
5211 | 国际民航组织对流云顶 | 407 | ||||
5216 | 总沉淀率KG / M2 / S | 千克m-2 s-1 | 407 | precipitation_flux | 是啊 | |
5226 | 总降水量KG / M2 / TS | |||||
5262 | 2D对流云量 | 503 | ||||
5270 | 浅层对流指示器 | 505 | ||||
8208 | 土壤水分含量 | |||||
8209 | 冠层含水量 | 千克米2 | 407 | canopy_water_amount | ||
8223 | 一层中的土壤水分含量 | 千克米2 | moisture_content_of_soil_layer | |||
8225 | 深层土壤温度 水文学后 | ķ | soil_temperature | |||
8236 | 瓷砖上的雪量KG / M2 | 503 | ||||
9203 | 低云量 | 407 | ||||
9204 | 中型云量 | 407 | ||||
9205 | 高云量 | 407 | ||||
9210 | GT.2.5 OCTA KFT的云基础 | 407 | ||||
9212 | GT.4.5 OCTA KFT的云基础 | 407 | ||||
9216 | 总云量-随机重叠 | |||||
9217 | 总云量最大/随机重叠 | 407 | ||||
9218 | 1000英尺ASL以下的云层分数 | 407 | ||||
9221 | 湿球冻结高度(M) | 407 | ||||
15201 | 压力水平风的U分量 | 1秒 | 500 | eastward_wind | ||
15201 | U风压等级B网格 | 1秒 | 501 | eastward_wind | ||
15202 | 压力水平风的V分量 | 1秒 | 500 | northward_wind | ||
15202 | V风压等级B网格 | 1秒 | 501 | northward_wind | ||
15212 | 50米风量U分量 | 1秒 | 500 | eastward_wind | ||
15212 | 50米风量U分量 | 1秒 | 500 | x_wind | ||
15212 | 50米风量U分量B网格 | 1秒 | 501 | eastward_wind | ||
15212 | 50米风量U分量B网格 | 1秒 | 501 | x_wind | ||
15213 | 50米风速V分量 | 1秒 | 500 | northward_wind | ||
15213 | 50米风速V分量 | 1秒 | 500 | y_wind | ||
15213 | 50米风速V分量B网格 | 1秒 | 501 | northward_wind | ||
15213 | 50米风速V分量B网格 | 1秒 | 501 | y_wind | ||
15242 | W级(风向)在压力水平上 | 1秒 | 500 | upward_air_velocity | ||
16202 | P LEV / P GRID上的地势高度 | 米 | 606 | geopotential_height | ||
16202 | 地势高度:压力水平 | 米 | 601 | geopotential_height | ||
16203 | P LEV / P GRID上的温度 | ķ | 606 | air_temperature | ||
16203 | 压力水平的温度 | ķ | 601 | air_temperature | ||
16204 | 压力水平相对湿度 | % | 503 | relative_humidity | ||
16204 | P LVS上的相对湿度WRT冰 | 505 | 601 | |||
16204 | 左/右网格上的RH WRT冰 | 606 | ||||
16205 | 湿泡势温度 | 407 | ||||
16205 | 湿泡势温度K | 502 | ||||
16222 | 海平面压力 | 帕 | air_pressure_at_sea_level | 是啊 | ||
20001 | 河流蓄水M2 | 505 | 505 | |||
20001 | 厚度:1000-500 MB | 601 | ||||
20002 | GRIDBOX出水量KG / S | 505 | 505 | |||
20002 | 厚度:1000-850 MB | 601 | ||||
20003 | GRIDBOX流入KG / S | 505 | 505 | |||
20003 | 风速:10米 | 1秒 | 601 | wind_speed | ||
20004 | 河流流量KG / M2 / S | 505 | 505 | 是啊 | ||
20004 | 风速:压力等级 | 1秒 | 601 | wind_speed | ||
20012 | 对流云深度 | 601 | ||||
20014 | 降水符号 | 601 | ||||
20015 | 当前的天气代码 | 601 | ||||
20020 | 最大风量U分量 | 601 | ||||
20021 | 最大风量V分量 | 601 | ||||
20022 | 最大风压 | 601 | ||||
20023 | 最大风速ICAO高度 | 601 | ||||
20024 | 对流线压力 | 帕 | 601 | tropopause_air_pressure | ||
20025 | 对流线温度 | ķ | 601 | tropopause_air_temperature | ||
20026 | 对虾高度 | 米 | 601 | tropopause_altitude | ||
20027 | TROOPAUSE ICAO身高 | 601 | ||||
20028 | 下雪概率 | 601 | ||||
20033 | 冻结地势高度 | 601 | ||||
20034 | 冻结水平压力 | 帕 | 601 | air_pressure_at_freezing_level | ||
20035 | 冻结ICAO高度 | 601 | ||||
20036 | -20C等温地势高度 | 601 | ||||
20037 | -20C等温压力 | 601 | ||||
20038 | -20C ISOTHERM国际民航组织高度 | 601 |
上午
藏 | STASHmaster说明 | 单位 | 有效期自 | 有效 | CF标准名称 | 用过的? |
2 | 时间戳后的风量 | 毫秒-1 | eastward_wind | 是啊 | ||
2 | 时间戳后的风量 | 毫秒-1 | x_wind | 是啊 | ||
3 | 时间戳后的风速V分量 | 毫秒-1 | northward_wind | 是啊 | ||
3 | 时间戳后的风速V分量 | 毫秒-1 | y_wind | 是啊 | ||
10 | 计时后的特定湿度 | 1个 | specific_humidity | |||
30 | LAND MASK(逻辑:LAND = TRUE) | 1个 | 407 | land_binary_mask | ||
30 | LAND MASK(无光晕)(LAND = TRUE) | 1个 | 500 | land_binary_mask | ||
33 | 地形(/ STRAT LOWER BC) | 米 | surface_altitude | |||
90 | 总气溶胶(可见度) | |||||
150 | W COMP OF WIND CP TH LEVS:VAR DUMMY | 401 | 407 | |||
150 | 时间戳后W的风分量 | 毫秒-1 | 500 | upward_air_velocity | ||
266 | 大块云分数 | 500 | 500 | |||
266 | 每层的大块云分数 | 501 | ||||
407 | TS后RHO水平的压力 | 帕 | 500 | air_pressure | 是啊 | |
408 | TS之后的THETA压力水平 | 帕 | 500 | air_pressure | 是啊 | |
409 | 时间戳后的表面压力 | 帕 | 500 | surface_air_pressure | 是啊 | |
1004 | SW RAD增量后的温度 | ķ | air_temperature | 是啊 | ||
2004年 | LW辐射增加后的温度 | ķ | air_temperature | 是啊 | ||
3216 | 带M2的边界层热流 | 501 | ||||
3216 | 肋面层(海冰) | 400.2 | ||||
3219 | 冲浪和BL风应力的X-COMP N / M2 | 帕 | 407 | surface_downward_eastward_stress | ||
3220 | 冲浪和BL风应力的Y-COMP N / M2 | 帕 | 407 | surface_downward_northward_stress | ||
3222 | B.总水分通量KG / M2 / S | 501 | ||||
4205 | LS PRECIP之后的云水 | |||||
4206 | LS PREPIP之后的云冰含量 | |||||
4222 | 降雨率超出模型水平 | 405 | 407 | |||
5181 | 温度INCR:对流 | 钾-1 | 501 | tendency_of_air_temperature_due_to_convection | ||
5182 | 特别湿度INCR:对流 | s -1 | 501 | tendency_of_specific_humidity_due_to_convection | ||
5212 | 转换 每个模型级别上的云量 | 1个 | 407 | convective_cloud_area_fraction_in_atmosphere_layer | ||
5250 | 超级质量通量(PA / S) | 404.1 | 407 | |||
5250 | 超级质量通量(Pa / s) | 503 | ||||
6201 | 重力波应力的U分量 | 帕 | 501 | atmosphere_eastward_stress_due_to_gravity_wave_drag | ||
6201 | X重力波应力分量 | 帕 | 502 | atmosphere_eastward_stress_due_to_gravity_wave_drag | ||
6202 | 重力波应力的V分量 | 帕 | 501 | atmosphere_northward_stress_due_to_gravity_wave_drag | ||
6202 | 重力波Y分量 | 帕 | 502 | atmosphere_northward_stress_due_to_gravity_wave_drag | ||
6223 | GW饱和应力的X分量 | 502 | ||||
6224 | GW饱和应力的Y分量 | 502 | ||||
9229 | 主云后的相对湿度 | 502 | 是啊 | |||
16004 | 模型水平上的温度 | ķ | 501 | 503 | air_temperature | 是啊 |
16004 | THETA水平上的温度 | ķ | 505 | air_temperature | 是啊 | |
16201 | THETA水平的地势高度 | 米 | 503 | geopotential_height | ||
16201 | 地势高度:层边界 | 米 | 407 | geopotential_height |
将Met Office Unified Model输出文件转换为NetCDF文件
使用xconv / convsh实用程序
有一个方便的实用程序,称为Xconv/convsh
,可将Met Office UM输出文件转换为不同的格式。这是将.pp文件转换为NetCDF文件的便捷步骤。
制备
-
下载
Xconv
并convsh
从这里:http://badc.nerc.ac.uk/help/software/xconv/index.html。如果使用Windows,请同时获取Xconv
和convsh
。如果使用Linux,请获取Xconv
并创建指向convsh
(ln -s xconv1.91 convsh1.91
)的符号链接。 -
创建名为的文本文件
conv2nc.tcl
,然后从此处复制代码:http : //www.met.rdg.ac.uk/~jeff/xconv/example2.html。省略第一行#!/home/jeff/bin/convsh
。您可能已经注意到,此脚本会将多个.pp文件转换为单个NetCDF文件。请随意修改下面的批处理文件以执行相同操作,然后将其上传到Wiki。 -
Windows用户:请按照此处的说明分别处理文件:http ://badc.nerc.ac.uk/help/software/xconv/xconv_windows.html(Xconv)。这也应适用于Linux用户。
-
Windows用户:您还可以使用批处理文件处理多个.pp文件。我无法在可用的时间内找到一个优雅的解决方案,但是这是一个可以解决问题的优雅方法(只要您的文件夹路径中没有空格):
示例脚本
-
创建一个名为的文本文件
run_Xconv.bat
。 -
将以下脚本复制到其中:
@echo off
rem Batch file to convert Met Office .pp files to NetCDF format
rem Karen Amoudry, 18 June 2013, National Oceanography Centre, Liverpool
setlocal EnableDelayedExpansion
call :treeProcess
:treeProcess
for %%f in (*.pp) do (
set outname=%%~nf.nc
set "cmd_str=convsh1.91.exe < conv2nc.tcl -i %%~ff -o !outname!"
echo !cmd_str!
)
echo ^@echo off
for /D %%d in (*) do (
cd %%d
call :treeProcess
cd ..
)
exit /b
-
广场
run_Xconv.bat
,convsh1.91.exe
并conv2nc.tcl
在你要处理的文件.PP上面的父文件夹。 -
打开命令窗口,导航到包含步骤©中的文件的文件夹,然后键入以下内容:
run_Xconv.bat > out.bat
-
该脚本完成后,在命令窗口中键入以下内容:
out.bat > log.txt
-
现在,您应该将所有NetCDF文件与步骤©中的文件放在同一目录中。
-
为了获得额外的信誉,编写一个批处理文件以将NetCDF文件放回原始文件夹结构中,并将其上传到Wiki。我有空的时候无法完成这项工作!
Linux用户:原则上,您应该能够创建一个不错的Shell脚本来执行与步骤(4)相同的操作。我没有,因为我的文件在Windows中。如果您制作一个,请放在Wiki上!
使用MATLAB
FVCOM工具箱包含一个函数pp2nc.m
,该函数调用convsh进行许多文件的转换。因此,它仍然需要convsh安装,但是可以更轻松地将其包含在现有FVCOM工作流程中。该功能依赖于使用FVCOM工具箱附带的convsh脚本(该脚本最初从此处获取)。下面显示了一个简单的用例。
files = {'/path/to/file1.pp', '/path/to/file2.pp', '/path/to/file3.pp'};
convsh = '/users/modellers/pica/bin/convsh';
tcl = '/users/modellers/pica/MATLAB/fvcom-toolbox/utilities/pp2nc.tcl';
pp2nc(files, convsh, tcl);
注意:convsh二进制文件的路径将需要在pp2nc.tcl脚本中进行编辑(将第一个列表更改为指向您的convsh二进制文件)。