FVCOM数据集

下面的页面概述了一些已与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_ONHEATING_CALCULATED为2011年以后的数据:

变量名长名外型尺寸单位
longitudelatitude度数
latitudelatitude度数
x-windx风分量(相对于网格)t_1htlatitudelongitude毫秒-1
y-windy风分量(相对于电网)t_1htlatitudelongitude毫秒-1
solar净短波辐射通量t_3surfacelatitudelongitude宽m -2
longwave净长波辐射通量t_3surfacelatitudelongitude宽m -2
sh显热通量t_1surfacelatitudelongitudeWm -2
lh潜热通量t_1surfacelatitudelongitudeWm -2
rh相对湿度t_1htlatitudelongitude
precip总降水率t_1surfacelatitudelongitude千克米-2 s -1
p_4 要么 p_3压力t_1msllatitudelongitude
temp温度Tt_1surfacelatitudelongitude摄氏温度

降水需要转换为m / s(除以淡水密度1000kg m -3)。与NCEP数据一样,我正在计算潜热的蒸发量,并从降水量中减去它。然后将蒸发量设置为零,并使用沉淀蒸发量作为沉淀变量。

变量清单

缩略语长名
美联社分析,压力水平
上午分析,模型级别
FP预测,压力水平
调频预测,模型级别

美联社

STASHmaster说明单位有效期自有效CF标准名称用过的?
23TIMESTEP后打雪量KG / M2千克米2 503snowfall_amount 
23陆上雪量TSTP KG / M2千克米2505 snowfall_amount 
24时间戳后的表面温度ķ  surface_temperature 
25时间戳后的边界层深度  atmosphere_boundary_layer_thickness 
26时间戳后的粗糙度  surface_roughness_length 
30LAND MASK(逻辑:LAND = TRUE)1个 407land_binary_mask 
30LAND MASK(无光晕)(LAND = TRUE)1个500 land_binary_mask 
31跨海后海冰的分形1个503 sea_ice_area_fraction 
31时间戳后的海冰分数1个 502sea_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 407surface_downwelling_longwave_flux_in_air是啊
3217带M2的表面热通量宽m-2501 surface_upward_sensible_heat_flux是啊
3217带M2的表面和B层热流宽m-2 500  
3223SURF&BL TOTL水分通量KG / M2 / S千克m-2 s-1 500upward_water_vapor_flux_in_air 
3223表面总水分通量KG / M2 / S千克m-2 s-1501 surface_upward_water_flux 
3224WIND MIX EN'GY FL到SEA:SEA MN W / M2宽m-2503 wind_mixing_energy_flux_into_sea_water 
3224W / M2混合风入海宽m-2 502wind_mixing_energy_flux_into_sea_water 
322510米风量U-COMP1秒 500eastward_wind是啊
322510米风量U-COMP1秒 500x_wind是啊
322510米风U-COMP B网格1秒501 eastward_wind是啊
322510米风U-COMP B网格1秒501 x_wind是啊
322610米风速V-COMP1秒 500northward_wind是啊
322610米风速V-COMP1秒 500y_wind是啊
322610米风速V-COMP B网格1秒501 northward_wind是啊
322610米风速V-COMP B网格1秒501 y_wind是啊
3234W / M2表面潜热通量宽m-2 407surface_upward_latent_heat_flux是啊
3236温度为150万ķ  air_temperature是啊
3237150万的特定湿度1个 407specific_humidity 
3245相对湿度150万  relative_humidity是啊
3248雾度为150万  407 
3250露点为1.5M(K)ķ 407dew_point_temperature 
3281TS S / M之后用于SO4的抗性S 401407  
3281VIS AT 1.5M(包括沉淀)M 606   
3281可见度为1.5MM 503601  
3316SNOW-ADJSTD标题K上的表面温度 405407  
3316标题K上的表面温度 503   
3321无冰砖上的冠层水KG / M2 405407  
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 407convective_rainfall_rate 
5206对流降雪率KG / M2 / S千克m-2 s-1 407convective_snowfall_rate 
5207对流云基的压力 407air_pressure_at_convective_cloud_base 
5208对流云顶的压力 407air_pressure_at_convective_cloud_top 
5210国际民航组织对流云基础  407  
5211国际民航组织对流云顶  407  
5216总沉淀率KG / M2 / S千克m-2 s-1 407precipitation_flux是啊
5226总降水量KG / M2 / TS     
52622D对流云量 503   
5270浅层对流指示器 505   
8208土壤水分含量     
8209冠层含水量千克米2 407canopy_water_amount 
8223一层中的土壤水分含量千克米2  moisture_content_of_soil_layer 
8225深层土壤温度 水文学后ķ  soil_temperature 
8236瓷砖上的雪量KG / M2 503   
9203低云量  407  
9204中型云量  407  
9205高云量  407  
9210GT.2.5 OCTA KFT的云基础  407  
9212GT.4.5 OCTA KFT的云基础  407  
9216总云量-随机重叠     
9217总云量最大/随机重叠  407  
92181000英尺ASL以下的云层分数  407  
9221湿球冻结高度(M)  407  
15201压力水平风的U分量1秒 500eastward_wind 
15201U风压等级B网格1秒501 eastward_wind 
15202压力水平风的V分量1秒 500northward_wind 
15202V风压等级B网格1秒501 northward_wind 
1521250米风量U分量1秒 500eastward_wind 
1521250米风量U分量1秒 500x_wind 
1521250米风量U分量B网格1秒501 eastward_wind 
1521250米风量U分量B网格1秒501 x_wind 
1521350米风速V分量1秒 500northward_wind 
1521350米风速V分量1秒 500y_wind 
1521350米风速V分量B网格1秒501 northward_wind 
1521350米风速V分量B网格1秒501 y_wind 
15242W级(风向)在压力水平上1秒500 upward_air_velocity 
16202P LEV / P GRID上的地势高度606 geopotential_height 
16202地势高度:压力水平 601geopotential_height 
16203P LEV / P GRID上的温度ķ606 air_temperature 
16203压力水平的温度ķ 601air_temperature 
16204压力水平相对湿度 503relative_humidity 
16204P LVS上的相对湿度WRT冰 505601  
16204左/右网格上的RH WRT冰 606   
16205湿泡势温度  407  
16205湿泡势温度K 502   
16222海平面压力  air_pressure_at_sea_level是啊
20001河流蓄水M2 505505  
20001厚度:1000-500 MB 601   
20002GRIDBOX出水量KG / S 505505  
20002厚度:1000-850 MB 601   
20003GRIDBOX流入KG / S 505505  
20003风速:10米1秒601 wind_speed 
20004河流流量KG / M2 / S 505505 是啊
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 
20027TROOPAUSE 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 
30LAND MASK(逻辑:LAND = TRUE)1个 407land_binary_mask 
30LAND MASK(无光晕)(LAND = TRUE)1个500 land_binary_mask 
33地形(/ STRAT LOWER BC)  surface_altitude 
90总气溶胶(可见度)     
150W COMP OF WIND CP TH LEVS:VAR DUMMY 401407  
150时间戳后W的风分量毫秒-1500 upward_air_velocity 
266大块云分数 500500  
266每层的大块云分数 501   
407TS后RHO水平的压力500 air_pressure是啊
408TS之后的THETA压力水平500 air_pressure是啊
409时间戳后的表面压力500 surface_air_pressure是啊
1004SW RAD增量后的温度ķ  air_temperature是啊
2004年LW辐射增加后的温度ķ  air_temperature是啊
3216带M2的边界层热流 501   
3216肋面层(海冰)  400.2  
3219冲浪和BL风应力的X-COMP N / M2 407surface_downward_eastward_stress 
3220冲浪和BL风应力的Y-COMP N / M2 407surface_downward_northward_stress 
3222B.总水分通量KG / M2 / S 501   
4205LS PRECIP之后的云水     
4206LS PREPIP之后的云冰含量     
4222降雨率超出模型水平 405407  
5181温度INCR:对流钾-1501 tendency_of_air_temperature_due_to_convection 
5182特别湿度INCR:对流s -1501 tendency_of_specific_humidity_due_to_convection 
5212转换 每个模型级别上的云量1个 407convective_cloud_area_fraction_in_atmosphere_layer 
5250超级质量通量(PA / S) 404.1407  
5250超级质量通量(Pa / s) 503   
6201重力波应力的U分量 501atmosphere_eastward_stress_due_to_gravity_wave_drag 
6201X重力波应力分量502 atmosphere_eastward_stress_due_to_gravity_wave_drag 
6202重力波应力的V分量 501atmosphere_northward_stress_due_to_gravity_wave_drag 
6202重力波Y分量502 atmosphere_northward_stress_due_to_gravity_wave_drag 
6223GW饱和应力的X分量 502   
6224GW饱和应力的Y分量 502   
9229主云后的相对湿度 502  是啊
16004模型水平上的温度ķ501503air_temperature是啊
16004THETA水平上的温度ķ505 air_temperature是啊
16201THETA水平的地势高度503 geopotential_height 
16201地势高度:层边界 407geopotential_height 

将Met Office Unified Model输出文件转换为NetCDF文件

 

使用xconv / convsh实用程序

有一个方便的实用程序,称为Xconv/convsh,可将Met Office UM输出文件转换为不同的格式。这是将.pp文件转换为NetCDF文件的便捷步骤。

制备

  1. 下载Xconvconvsh从这里:http://badc.nerc.ac.uk/help/software/xconv/index.html。如果使用Windows,请同时获取Xconvconvsh。如果使用Linux,请获取Xconv并创建指向convshln -s xconv1.91 convsh1.91)的符号链接。
  2. 创建名为的文本文件conv2nc.tcl,然后从此处复制代码:http : //www.met.rdg.ac.uk/~jeff/xconv/example2.html。省略第一行#!/home/jeff/bin/convsh。您可能已经注意到,此脚本会将多个.pp文件转换为单个NetCDF文件。请随意修改下面的批处理文件以执行相同操作,然后将其上传到Wiki。
  3. Windows用户:请按照此处的说明分别处理文件:http ://badc.nerc.ac.uk/help/software/xconv/xconv_windows.html(Xconv)。这也应适用于Linux用户。
  4. 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.batconvsh1.91.execonv2nc.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二进制文件)。

 

  • 1
    点赞
  • 9
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
FVCOM 3.2.1是Finite-Volume Coastal Ocean Model的版本号。FVCOM是一种用于模拟近岸海洋和海岸区域的数值模型,可以对海洋流动、海洋生态系统、气候变化等进行模拟和预测。 FVCOM采用有限体积方法对海洋进行离散化,将海洋分为许多小区域,然后利用方程描述物理过程(如流体力学、热力学、海洋化学等)在这些小区域内的交互作用。FVCOM模型还包括边界条件、初始条件和外部强迫,通过计算这些参数的变化来模拟海洋的动态变化。其优点在于能够对不同空间和时间尺度上的物理和生物过程进行高精度的模拟,并能够准确预测海洋环境中的变化。 FVCOM 3.2.1作为一个特定版本的FVCOM,在之前版本的基础上进行了优化和改进。它可能会提供更高的计算效率、更精确的数值计算、更多的功能和更稳定的模拟结果。然而,具体的改进和优化内容需要参考版本更新文档以及相关的出版物和研究论文。 FVCOM 3.2.1的应用范围广泛,可以用于研究海洋生态系统、海洋环境、沿岸工程等问题。同时,它还可以用于预测海洋变化对人类活动的影响,为渔业、航运、海岸防御等提供决策支持。此外,由于它可以模拟全球范围内的海洋现象,还可以用于研究全球气候变化和海洋循环等重要科学问题。 总之,FVCOM 3.2.1是一种用于模拟和预测近岸海洋和海岸区域的数值模型,它的改进可能包括提高计算效率和精度。它在海洋科学、环境保护和工程应用等领域具有重要作用。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值