关键词不能为空

当前您在: 主页 > 英语 >

ROMSTOOLS 中文翻译

作者:高考题库网
来源:https://www.bjmy2z.cn/gaokao
2021-03-02 16:53
tags:

-

2021年3月2日发(作者:sterile)


ROMSTOOLS


文件



表的内容



介绍



安装




系统要求




获取文件



教程︰南部格拉示例




入门︰处理输入的文件




构建网格




让风和其他表面通量




初始和侧边界条件




得到的径流强迫



编译模型




配置参数︰



param.h



cppdefs.h




数值模拟和物理选项︰




编译脚本︰



jobcomp



运行模式




长时间的模拟




得到的结果




诊断程序




潮汐



生物地球化学



年度间模拟




表面的强迫数据得到






CFSR


(推荐)



NCEP2


再分析获取表面强迫数据



每月气候学数据




QuikSCAT



获取侧边界条件




运行年际运行的模型



嵌套




嵌套网格的制备




编译和运行模型



?



海洋预报︰



Forecast_tools



?



参考书目



ROMSTOOLSUser


的指南



-ROMSTOOLS-



吉尔达斯由布莱特



Penven



Patrick Marchesiello




Laurent




研究所倒



le D



é



veloppement


(税务局)




44


大道



de


敦刻尔克,


CS 90009


13572


马赛楼



cedex 02


法国



介绍



区域海洋建模系统



(ROM)


是新的一代海洋环流模式(


Shchepetkin


和威廉姆斯,


2005


年),已


专为区域海洋系统的准确的模拟。读者可参考



Shchepetkin


和威廉姆斯(


2003


年)、



Shchepetkin


和威廉姆斯(


2005


年)为模型的完整描述。


ROM


已经 申请区域模拟研究各种不


同地区的世界海洋


(例如



Marchesiello et al.



2003



;Penven et al.



2001



;


麦克格雷迪



et al.



2002



;


海德福格尔等人,


2000



;Di Lorenzo et al.



2003



;Blanke


等人,

< p>
2002


年)。



若要执行区域模拟



ROM

< p>
,建模者需要提供特定格式的多个数据文件︰水平网格,海底地形表面


强迫 、侧边界条件


......


他还需要对模型输出分析。这里描述 的工具原本是为了执行这些任务。目


标是要能够建立一个标准的区域模型配置在最短的时 间。



安装



系统要求



利用



matlab


设计了这个工具箱。它需要至少



2 Gbites


的磁盘空间。它已经过测试的几个



Matlab


版本从



Matlab6




Matlab2012


。主要测试在



Linux


工作站上,但它可以使用在任何平


台上如果



NetCDF




LoadDAP Matlab Mex


文件提供。


NetCDF Matlab Mex


文件需要读取和


写入到



NetCDF


文件,它可以在



web


位置找到︰



/



LoadDAP


Matlab Mex


文件用于从年度间和预测模拟的简化服务器下载数据。它可以在



web


位置找到︰



/download/



Matlab LoadDAP Mex


文件提供任何简化可


访问的数据读入



Matlab


方法。请注意,


LibDAP


图书馆之前必须先安装您的系统上安装



LoadDAP



详细信息可以在



web


位置找到︰





MexCDF




LoadDAP


mex


文件为



Linux


提供



(FEDORA 32


位系统︰



mexcdf




Opendap_tools/FEDORA;


系统



CENTOS




FEDORA 64bits: mexnc




Opendap_t ools/FEDORA_X64)


,但他们不在所有平


台上工 作。





ROMSTOOLS


包中包含所有其他必要



Matlab


工具箱(即海


-


气、面具、

< p>


netcdf



< /p>


m_map......



< p>
全球数据集,


例如地形


(史密斯和桑德,


1997


年)



水文< /p>



Conkright


等人,


2002


年)或地表通量(


Da Silva et al.



1994


年),也被包 括在内。



获取文件



All the necessary compressed tar files () containing the Matlab programs, several


datasets and other toolboxes and softwares needed by ROMSTOOLS are located at


in the download section.


下载所有压缩的



tar


文件。解压缩并解压所有


文件



(gunzip


和焦油


-xvf)< /p>




定义的不同的目录︰



?



Aforc_NCEP


︰用于恢复的表面迫使年度间模拟数据(基于



CFSR


再分析)的脚本。



?



Aforc_QuikSCAT


︰用于风应力回收卫星散射计数据



(QuickSCAT)


的脚本。



?



CARS2009: CSIRO


区域海洋地图集数据库。每 月、季、年气候的温度、盐度、硝酸盐、


磷酸盐和氧气



?



COADS05


︰目录表面通量全球每月气候学在决议(


Da Silva et al.



1994


年)。



?



编译︰空暂存目录



ROM


编译。



?



Diagnostic_tools


︰几个



Matlab


脚本动画和基本的统计分析。



?



ROMSTOOLS


用户指南文档︰位置。



?



Forecast_tools


︰脚本生成的业务海洋预报系统



?



掩码︰土地掩模版工具箱开发的



A.Y.


谢尔比纳。



?



mex60: Matlab NetCDF


接口为



32




64




Linux


体系结构和旧的



matlab


版本︰



6


之前。



?



mexcdf/mexnc: Matlab NetCDF


接口为



32




64




Linux


体系结构中,直到



R2008a


MatlabR14sp1 (/downloads/)


。< /p>


下一发行


版的



Matlab



R2008b



R2009a


,它是



mexcf


的更简单,要么使用本机的



NetCDF




matlab


或使用最后一次发布版本的相同



url


上后



R2008a



()


/downloads/



?



mexcdf/netcdf_toolbox: Matlab NetCDF


工具箱可用相同的



mexcdf


包中。



?



m_map: Matlab


映射工具箱



(/rich/)




?



Nesting_tools


︰预处理工具用于编写嵌套的模型。



?



netcdf_g77: NetCDF Fortran


库对于



Linux


,编译使用



77


国集团



(/packages/netcdf/)




?



netcdf_ifc: NetCDF Fortran


库对于



Linux

< br>,




ifort


编译。


英特尔



Fortran


编译器



(ifort)



可在



/s oftware/products/compilers/flin/




?



netcdf_x86_64: NetCDF Fortran


库对于



Linux


,用



ifort




64


位体系结构上编译。



?



Oforc_OGCM


︰用于初始和侧边界条件从全球



OGCMs


(汽水(纸箱等人,


2005


年)或



ECCO


(口吃



et al.



1999


年))的年度间模拟恢复脚本。

< p>


?



Opendap_tools: LoadDAP mexcdf


和几个脚本,自动通过互联网下载数据。



?



Preprocessing_tools


︰预处理



Matlab


脚本



(make_grid.m




make_forcing


等。


.. ....)




?



河流︰脚本编写时变径流迫使文件并计算径流位置。



?



RUNOFF_DAI


︰戴伯斯全球河流流量和欧式放电每月气候学。


()


/cas/catalog/surface/dai- runoff/



?



Roms_Agrif: ROM Fortran


来源。



?



运行︰工作目录。这是



ROM


输入的文件生成的位置和运行模型的位置。



?



SeaWifs


︰表面叶绿素



a


气候学



SeaWifs


观察的基础。



?



SST_pathfinder: SST


全球每月气候学比



COADS05


,计算从



AVHRR


探路者观察从



1985


年到



1997


年的精细分辨率(


9.28


公里)。(凯西和科尔尼永,


1999


年)



?



潮汐︰



Matlab


例程



ROM


潮汐模拟做准备。潮汐数据来自俄勒冈州立大学全球模型的海


洋潮汐

< br>


TPXO6




TPXO7


(埃格伯特和



Erofeeva



2002

< br>年)︰



/research/po/research /tide/




?



Topo


︰位置



2


上的全球地形数据集°分辨率(史密斯和桑德,


1997


年)。 原始数据可以


发现在︰



/cgi- bin/get_



?



TPX06


︰目录全球模型的海洋潮汐



TPXO6


(埃格伯特和



Erofeeva



2002

< br>年)。



?



TPX07


︰目录全球模型的海洋潮汐



TPXO7


(埃格伯特和



Erofeeva



2002

< br>年)。



?



Visualization_tools: Matlab


脚本



ROM


可视化图形用户界面。



?



WOA2009


︰世界海洋地图集



2009


年全球数据库。



引用列表︰


/OC5/WOA09/



?



WOAPISCES


︰生物地球化学双鱼座数据(年度、季节 气候学)全球数据集。引用︰




铁和



DOC: Aumont et Bopp



2006





寺,


O2



NO3


,从


WOA2005 PO4




DIC


及碱度来从



Goyet et al




教程︰南部格拉示例



本节介绍了编写和运行区域的



ROM


仿真的基本步骤。


这是完成之后的南部的格拉在低分辨率和


使用气候强迫在表面和边界模型的示例。



入门︰处理输入的文件



成功安装后,启动



Matlab


会话目录中︰



~/Roms_tools/


运行。



运行


start.m


脚本来为此会话设 置



Matlab


路径。



在此步骤中此安装,您必须知道关于



matlab


安装和旅游计算机环境的事情︰



?



我机



32




64


位 体系结构是什么?为此做


uname-a.



?



我有的本机



matlab


安装库是什么?



?



如果我有我的



Matlab


版本已经本机



netcdf


例程和图书馆,不必使用



netcdf


库提供的



Roms_tools


,所以从



start.m


文件中删除它


.



?



如果我将本机



m_map


例程已经跟我的



Matlab


版本了,我不需要使用



netcdf


库提供的



Roms_tools


,所以从



start.m


文件中删除它。



?



对于这些问题,它可用于在



matlab


会话中编辑你



matlab


路径与



matlab


命令


路径






Roms_tools


,提供了



matlab


一些



NetCDF


库︰



?



mex60: matlab12


(老),


32




64


位体系结构



?



mexcdf/mexnc: matlab7



matlab2008



matlab2009



32


位和



64


位体系结构





Roms_tools

< p>
,还提供一些



简化



和图书馆︰



?



Opendap_tools/FEDORA: LibDAP




LoadDAP bin


和图书馆



Fedora Linux


发行版,


32


位体系结构



?



Opendap_tools/FEDORA_X64


︰相同的



64


位体系结构。



然而,如果您的



Linux


发行版不同于



Fedora


,最好是要编译和安装由您自己的



LibDAP



< br>LoadDAP


。现在你已经准备好创建新的配置。它是重要的是尊重以下的预处 理步骤的顺序︰



make_grid



make_forcing



mak e_clim


。所有的预处理步骤,对于只有一个要编辑的文件



~/Roms_tools/Run/romstools_param.m


。此文件包含必需的参数为



ROM


输入



NetCDF



件的生成。



Romstools_param.m


的第一节定义用于生成 网格的配置参数。例如标题、工作目录或文件名称


:



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%



%



%1-


配置参数



%


使用的


make_grid.m


(和其他人


......


< p>


%



%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

< br>


%



%Rom


标题名称和目录



%



ROMS_title =


本格拉模型


;



ROMS_config = 'Benguela_LR';



%



%


网格尺寸︰



%



lonmin = 8;%


最低


[


度东经


]



lonmax = 22;%


最大东经


[



]



latmin =-38;%


最低


latitudeF [



]



latmax =-26;%


最大纬度


[



]



%



%


网格分辨率


[



]



%



dl = 1/3;



%



%


的垂直级别数


(


!应该是相同的


param.h



)



%



N = 32;



%



%


垂直网格参数


(


!应该是相同的




)



%



theta_s = 6



;



θ



试卷


= 0



;



hc = 10



;%



拉伸正面厚度流



vtransform = 1



;%s


的坐标类型


(1


︰老






; 2


︰新坐标


)



%



在岸边


[m]


(取决于分辨率,最低深度


%



%


的经验法则︰


dl = 1



= 300



dl = 1/4


,两


= 150



...)



%


这影响过滤因为它适用于研究生


(h) / h




%




= 75;



%



最大深度在岸边


[m]


(以防止代


%



太大墙沿着海岸的


%




%



hmax_coast = 500;



%



最大深度


%[m]


(切地形以防止



外推


%


低于


WOA


数据)



%



hmax = 5000;





%




ETOPO 2


或任何其他


netcdf


文件


%


地形


netcdf


文件名称



在相同的格式


%




%



TOPODIR = '



/';



topofile = [TOPODIR,'Topo/'];



%



%


边坡参 数(


r=grad(h)/h)


最大值为地形平滑



%



rtarget = 0.25;



%



选择性的过滤,以减少孤立的通行证数目


%



%


对深海的海山。



%



n_filter_deep_topo = 4;



%



末尾的单一宁滤波器的通


%


数目



%smooting


程序,以确保在没有


2DX


噪声



%


的地形。



%



n_filter_final = 2;



%



%Gshss


用户定义的海岸线(见


m_map




%Xxx_


全分辨率数据



%Xxx_


高分辨率数据



%Xxx_


中间分辨率数据



%Xxx_


低分辨率数据



%Xxx_


粗分辨率数据



%



coastfileplot = 'coastline_';



coastfilemask = 'coastline_l_';



%



客观分析解相关规模


[m] %



%(


如果


Roa = 0


︰最近的外推法


;


粗但便宜得多


)



%



居港权


%= 300e3;



Roa = 0;



%



interp_method =


线性


;%


的插值方法


:


线性或立方


'



%



makeplot = 1;%1


︰每个预处理步骤之后创建几个图形



%



变量的描述︰



?



ROMS_title =


本格拉测试模式︰一般标题。你可 以给任何你想要为您的配置的名称。



?



ROMS_config = 'Benguela_LR':


配置的名称。这用于特定配置的



NCEP


或模式数据的存


储。



?



lonmin = 8


︰西方极限的经度在网格



[-360


°、


360


°


]


。网格是长方形的纬度


/


经度。



?



lonmax = 22


︰东部极限



[-360


°、


360


°


]


。应该是优于



lonmin




?



latmin =-38


︰网格在纬度的南限



[-90


°、



90


°


].



?



latmax =-26


︰北部限制



[-90


°,


90


°


]


。应该是优于



latmin.



?



dl = 1/3


︰网格经度决议以度为单位。纬度间距理论推导,得到各向同性的网格,使用关 系




d


Φ



= dl


,因为



(


Φ


)



?



N = 32


︰垂直的级别数。警告!


N


已在该文件中还定义︰



~/Roms_tools/Run/param.h



前编译模型。



?



theta_s = 6


。︰垂直



S


坐标表面拉伸参数


.


当建筑气候学和初始



ROM


文件,我们必须定


义垂直网格。警告!不同垂直网格参数应该是相同的和



ROM


输入文件中此文件



(




~ /


Roms_tools/Run/)



这是



bug


的一个严重原因。


可以使用



Matlab


脚本测试的



theta_s




theta_b




hc




N


的影响︰



~/Roms_tools /Preprocessing_tools/test_vgrid.m.



?



theta_b = 0


。︰拉伸参数垂直的



S


坐标底部。



?



vtransform = 1



: S


坐标类型



(1


︰老






; 2


︰新坐标


)


.


它是关联到


#NEW_S_COORD


cpp



ROMS_AGRIF


源代码中。



?



hc = 10



︰垂直



S


坐标



Hc


参数。


它给大约过渡深度水平表 面能级和底部地形跟随水平之


间。它应该是次于两。



?





= 75


︰以米为单位的最小深度。模型深度削减这水平,防止,例如,模型网格的发生细


胞没 有水。这并不影响掩蔽例程。在较低的分辨率,两应该是很大(例如为



dl 150 m = 1/2


)。


否 则,既然基于地形平滑,底部斜坡可以完全侵蚀。



?



hmax_coast = 500


︰面具下的最大深度。它可以防止选定的等深线



(


这里



500 m)


去下了


面具。如果是这 样,这可能(例如)是西边界流问题的来源。



?



hmax = 5000;



?



TOPODIR = '



/';



?



topofile = [TOPODIR,'Topo/'];


地形的默认文件。我们在这里使用



etopo2


(史密


斯和桑德,


1997

年)。



?



rtarget = 0.25


︰此变量控制的最大值的措施 的西格玛层(贝克曼和海德福格尔,


1993


年)



-


参数︰为了防止横向压力梯度误差,众所周知, 在地形跟随坐标模型(哈尼,


1991


年),切


合实际的地形要求一些平滑。实证结果表明那可靠的模型得到的结果,如果不超过



0.2




?



n_filter_deep_topo = 4


︰通宁筛选器防止噪音和孤立的海山上深地区发生的数目。



?



n_filter_final = 2


︰通宁数目筛选平滑过程结束时肯定没有噪音是目前在地形。



?



coastfileplot = 'coastline_':


二进制



GSHSS


沿海文件由



m_map

< p>
用于图形脱脂棉。



封信之前


选择海岸线决议。


f


︰全分辨率,


h


︰高分辨率,


i


︰中间决议,


l


︰低分辨率



c


︰粗分辨率。



?



coastfilemask = 'coastline_l_':


用于海岸线掩蔽工具箱中的二进制文件。



?



Roa = 0


︰解相关长度尺度在客观分析米(


300


公里是一个合理的值,就业数据集)。如果



Roa = 0



最近





Matlab


外推方法而不是客观的分析。这是便宜得多,< /p>


但结果可能会在较低的


质量。



?



interp_method = '


立方︰使用后的客观分析 水平的内插方法。它可以是线性或立方。



?



makeplot = 1


︰选择此选项可生成图像预处理的每一步之后。



第二部分在



romstools_param.m


定义的通用文件和目录的名称︰



%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%



%



%2-


通用文件和目录名称



%



%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%





%



%Romstools


目录



%



ROMSTOOLS_dir = '



/';



%



%


运行目录



%



RUN_dir = [pwd



'/'];



%



%Rom


输入


netcdf


文件目录



%



ROMS_files_dir = [RUN_dir



' ROMS_FILES /'];



%



%


全球数据目录


(etopo



coads


,数据集下载从


ftp


等。


.)



%



DATADIR = ROMSTOOLS_dir;



%



%


强迫数据目录


(ncep


,低空,数据集下载与简化等。


.)



%



FORC_DATA_DIR = [RUN_dir



'


数据


/'];



%



eval (['



mkdir '



ROMS_files_dir])



%



%Rom


文件名(网格,迫使,散装,气候学,初始)



%



grdname = [ROMS_files_dir,'roms_'];



frcname = [ROMS_files_dir,'roms_'];



blkname = [ROMS_files_dir,'roms_'];



clmname = [ROMS_files_dir,'roms_'];



bryname = [ROMS_files_dir,'roms_'];



ininame = [ROMS_files_dir,'roms_'];



bioname = [ROMS_files_dir,'roms_'];


双鱼座


%


铁粉尘迫使



rivname = [ROMS_files_dir,'roms_'];



%



%z


级的中间数据文件(不用于模拟)



%



oaname = [ROMS_ files_dir,'roms_'];


气候学数据加工


%



Zbryname = [ROMS_files_dir ,'roms_bry_'];%


的边界数据处理



%



迫使年际模拟


(NCEP/GFS)


根文件名


%


泛型



%



frc_prefix = [ROMS_files_dir



'roms_frc'];%


强制文件名称

< br>


blk_prefix = [ROMS_files_dir



'roms_blk'];%


大容量文件名称



%





?



ROMSTOOLS_dir = '



/':


目录。



?



RUN_dir = [ROMSTOOLS_dir



'


运行



/']: Roms_tool s/


运行目录。这是在哪里做了所有的


工作。

< br>


?



ROMS_files_dir = [RUN_dir



' ROMS_FILES /']: R oms_tools/


运行


/ROMS_FILES/


目录。



是在哪里



ROM


输入



NetCDF


文件的存储。



?



ROMS_files_dir = [RUN_dir



' ROMS_FILES /']:



?



DATADIR = ROMSTOOLS_dir;


︰全球数据目录



(ETOPO



COADS


,数据集下载从



ftp


等。


.)



?



FORC_DATA_DIR = [RUN_dir



'


数据



/']:


强迫数据目录



(NCEP

< p>
,低空,数据集下载与简


化等。


.)



?



grdname=[ROMS_files_dir,'roms_']: ROM


名称输入



NetCDF

< p>
网格文件。这是水平


网格参数的存储位置。一般情况下,我们遵循的风格︰



XXX_




?



frcname=[ROMS_files_dir,'roms_']::




ROM


名称输入



NetCDF

< p>
迫使文件。这是表


面的强迫变量(例如风应力)的存储位置。一般情况下, 我们遵循的风格︰



XXX_




?



blkname=[ROMS_files_dir,'roms_']: ROM


名称输入



NetCDF

< p>
大容量文件。


这是用于


大容量参数


(如空气温度)


的大气变量的存储位置。


一般情况下,


我们遵循的风格︰



XXX_




?



clmname=[ROMS_files_dir,'roms_']: ROM


名称输入



NetCDF


气候学文件。这是



ROM


的预后变量(


u




v


、温度、盐、乌巴、



vbar




zeta


)为侧边界和内部微移的存储位置。


该文件可以是大的 因为存储变量以供所有



ROM


网格 内部点。它被称为



气候学文件



,因为这是


在过去的恢复光盘使用的文件解决方案对原位气候(如莱维图 斯的例子)。一般情况下,我们遵


循的风格︰



XXX_




?



bryname=[ROMS_files_dir,'roms_']: ROM


名称输入



NetCDF

< p>
边界文件。


这是气候学


文件替代。


在这种情况下,


变量只存储为横向边界。


一般情况下,


我们遵循的风格︰



XXX_




?



ininame=[ROMS_files_dir,'roms_']: ROM


名称输入



NetCDF


初始文件。这是



ROM

< p>
预后变量(


u




v


、温度、盐、乌巴、



vbar




zeta


)为初始条件的存储位置。一般情况下,我


们遵循的风 格︰



XXX_




?



bioname=[ROMS_files_dir,'roms_']: ROM


名称输入



NetCDF

< p>
铁粉尘迫使双鱼座


生物地球化学模型文件。



?



rivname = [ROMS_files_dir,'roms_'];ROM


名称输入



NetCDF


决胜选举迫使文件



?



oaname = [ROMS_files_dir,'roms_']:


不使用的



ROM

< br>的中间文件的名称。这是相当


于气候学文件,


但在垂直的



z


协调。


首先,


变量水平内插来创建一个



roms_


文件



OA



件)。然后,他们垂直插补上气候学文件



ROM


的坐标。一般情况下,我们遵循的风格︰



XXX_




?



Zbryname=[ROMS_files_dir,'roms_bry_']:


边界文件上的



z

< br>坐标的中间文件。一般情


况下,我们遵循的风格︰



XXX_bry_




?



frc_prefix = [ROMS_files_dir



'roms_frc']:




inter_annual


模拟的情况下,迫使文件名称


的第一部分。在这种情况下,为每月创 建一个单独的文件。例如,迫使文件基于



NCEP 2000




1


月为︰



roms_frc_NCEP_



?



blk_prefix = [ROMS_files_dir



'roms_blk']:




inter_annual


模拟大容量文件名称的第一部


分。在这种情况下,为每月创建一个单 独的文件。例如,大容量文件基于



NCEP 2000




1



为︰



roms_blk_NCEP_



构建网格



保存



romstools_param.m




make_grid




Matlab


会话中运行︰



>>


make_grid



制造网格


:



/Run/ROMS_FILES/ro ms_



标题︰本格拉测试模型



分辨率︰


1/3




创建网格文件


...



法学硕士


= 41




= 42



填充网格文件


...



计算指标


......



Min dx = 29.1913


公里


-


马克斯


dx = 33.3244


公里



Min dy = 29.2434


公里


-


马克斯


dy = 33.1967


公里



填充网格文件


...



添加地形


...



ROM


的决议︰


31.3


公里



地形数据分辨率︰


3.42


公里



地形分辨率减半


4




新的地形决议︰


54.6


公里



正在处理


coastline_...



你想要使用


editmask


吗?


y



[n]



在深海要删除孤立的海山上应用筛选器︰



4


通过选择性筛选器。



应用


log(h)


减少


grad (h) /h


的选择性的过滤︰



13


迭代


-rmax = 0.24381



最后一次平滑地形以防止


2DX


噪音︰



2


通过汉宁的更顺畅。



把它写下来


......



做一个阴谋


......



你应该在过程中保持法学硕士和



MMm


的值。他们将



ROM


参数文件



/Roms_tools/Run/param.h


的必要条件。在这个测试案例,


LLm0 = 23



MMm0 = 31




在网格生成过程中,这个问题



做你想要使用



editmask



y


< p>
[n]


问。



默认回答是



n


(为



no)


。如果答案是



y


(代表



yes)


,将推出


< br>editmask



bina


, 由


开发的图形界面来手动编辑掩码


(请注意,

< br>就目前而言,


editmask


不正在与



matlab7




mexnc





否则该掩码是从未经过滤的地 形数据生成的。过程可以防止孤立的地(或海面)存在点。





(1.1)


提出了海底地形与南部格拉例如



make_grid.m


获得。请注意,在这种低分辨率



(1 /3)



地形已经强烈理顺。



.






1.1


︰结果的本格拉例如



make_grid.m



让风和其他表面通量



下一步是创建包含不同地表通量的文件。您应该编辑文件



romstools_param.m


的部分是︰



%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%



%



%3-


表面迫使参数



make_forcing.m



make_bulk.m


使用的


%



%



%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%



%



%Coads


目录(用于气候学运行)



%



coads_dir = [DATADIR



' COADS05 /'];



%



%Coads


时间(对于气候学运行)



%



coads_time=(15: 30:345)



%:


每个月的中间



coads_cycle = 360;


一典型


360


年天的


%


重复



%



%coads_time=(15.2188:30.4375:350.0313); %365.25


年天案



%coad s_cycle=365.25%



QSCAT


实验



%


气候热通量。



%



探路者


SST


数据使用的


pathfinder_sst.m %



%



pathfinder_sst_name = [DATADIR



......



]; SST_pathfinder/climato_


< /p>


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%



变量的描述︰



?



coads_dir = [ROMSTOOLS_dir



' COADS05 /']:


目录在哪里在



1/2 < /p>


的表面海洋数据全球地


图册°位于决议(


Da Silva et al.



1994

< br>年)。



?



coads_time=(15:30:345)


︰在天为每 月气候学的时间。


它对应于每个月的中间。


ROM

< p>
使用


这一次在时间线性插值的强迫的变量。



?



coads_cycle = 360


︰持续时间的循环强制 变量。在这里,为了简单起见,我们的模型上运


行重复气候



360


年天。



?



pathfin der_sst_name=[ROMSTOOLS_dir,SST_pathfinder/climato _']:


目录每


月气候学海表面温度从探路者卫星观测


(凯西和科尔尼永,


1999


年)

< p>


这可以用有一个



Da Silva


等人(


1994


年)



SST


的替代。



保存



romstools_param.m




make_forcing




Matlab


会话中运行︰





>>


make_forcing




您应该获得︰



本格拉测试模型



在网格中读


......





创建强制文件


...



Taux


获取时间指数


1



获取时间指数


1


陶伊河



...



让几个情节


......




此程序可以采取相对较长的时间来处理强迫的所有变量。图



(1.2)


提出了风应力矢量和风应力准

则从表面海洋数据全球地图册索取



1/2


°在



4


个不同阶段的年决议



Da Silva et al.



1994


年)



Da Silva


等人(


1994


年)海表温度



(SST)


用于恢复期限



(dQdSST) < /p>


在热通量计算。要提高模


型的求解就可以使用


SST


气候学在更精细的分辨率



9.28


公里)


(凯西和科尔尼永,


1999< /p>


年)



要这样做,你可以在



Matlab


会话中运行



pathfinder_sst.m:




>>


pathfinder_sst




您应该获得︰



...


月指数︰


1



...


月指数︰


2



...



对于表面的强迫,而不直接订明通量,有可能要使用批量公式在模型运行过程中从大气变量生成

< br>地表通量。在这种情况下,光盘需要重新编译与定义的



BULK_FLUX cpp


键。若要生成批量迫


使文件,您需要在



Matlab


会话中运行



make_bulk:



>>


make_bulk




您应该获得︰



本格拉测试模型





在网格中读


......





创建批量迫使文件


...



让坐时间指数


1



让坐时间指数


2



...



让几个情节


......






1.2


︰风



stress[N.m


?


]


获得本格拉例如使用



make_forcing.m.



初始和侧边界条件



最后一个预处理步骤包括在生成包含



ROM


初次和外侧的开放边界条件的必要信息的文件。


此脚


本将生成两个文件︰气候学文件



(XXX_)


给出了外侧的边界条件和初始条件文件



(XXX_)




应该由用户在文件



romstools_param.m


中编辑的部分是。




应该由用户在文件



romstools_param.m


中编辑的部分是︰



%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

< br>


%



%4-


开放边界和初始条件参数



使用


make_clim.m



make_biol.m



make_br y.m %



%make_ogcm.m



make_OGCM_frcst.m



%



%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%



%



%


开放边界交换机


(


!应符合


cppdefs.h



)



%



obc = [1 1 1 1];%


开放边界


(1 =


开,


[S E N W])



%



%


的地计算参考水平



%



zref =-1000;



%



%


初始


/


边界数据选项


(1 =


进程


)



%


(用


make_clim



make_biol



make_bry




%make_ogcm.m



make_OGCM_frcst.m




%



makeini = 1;%


的初始数据



makeclim = 1;%


气候数据(用于边界和逼近层)



makebry = 1;%


侧边界数据



makenpzd = 0;%nchlpzd



N2ChlPZD2


模型的初始条件和边界数据



makebioebus = 0;BioEBUS


模型


%


初始和边界数据



makepisces = 0;


双鱼座模型


%


初始和边界数据



%



%



makeoa = 1;%oa


数据(中间文件)



makeZbry = 1;Z %


边界数据协调(中间文件)



insitu2pot = 1;


对潜在温度


%


变换原位温度



%



psource_ts = 0;


示踪剂径流浓度


%



%



初始化为气候学实验


%



(= 0: 0 h 1



1



)



%



tini = 0;



%



%


世界海洋地图集目录


(WOA2009)



%


(盐、生物武器和临时变量)



woa_dir = [DATADIR



' WOA2009 /'];



%



%Cars2009


气候学目录


(CARS2009)



%


(盐、生物武器和临时变量)



cars2009_dir = [DATADIR



' CARS2009 /'];



%



双鱼座生物地球化学的季节性气候学


%



%



woapisces_dir = [DATADIR



' WOAPISCES /'];



%



%


气候数据


dir



t



s


、生物变量)



%



climato_dir = cars2009_dir;



%


< /p>


%


表面叶绿素季节性气候学(


WOA20 01



SeaWifs




%



chla_dir = [DATADIR



' SeaWifs /'];



%



径流每月季节性气候学(傣族和伯斯)


%



global_clim_riverdir = [DATADIR



' RUNOFF_DAI /'];



global_clim_rivername=[ global_clim_riverdir,'Dai_Trenberth_runoff_global< /p>


_'];



%



%


设置时间和周期边界条件为︰



%


每月气候学



%



woa_time=(15:30 :345)



%:


每个月的中间



woa_cycle = 360;


一典型


360


年天的


%


重复



%



%woa_time=(15.2188:30.4375:350.0313);%3 65.25


年天案



%woa_cy cle=365.25%



QSCAT


实验



%


气候边界条件



%



变量的描述︰



?



obc = [1 1 1 1]:


开关打开



(1 =



)


或关闭



(0 =



)


侧边界



[


南东部北部西部


]

< br>。这用于大


规模执法的应用。


要知道,

< br>这应该是与开放边界文件



/Roms_tools/Run/cppdefs.h


中的



CPP



换机兼容。



?



zref = 1000


年︰深度



[



]


地转速度计算没有运动的水平。



?



makeini = 1


︰开关来定义如果生成初始文件



(roms_)


。应为



1




?



makeclim = 1


︰开关来定义如果生成的气候


(侧边界条件)


文件


(roms_)



应为



1




?



makebry = 1


︰开关来定义如果生成的边界文件



(roms_)


。仅用



make_bry




?



makenpzd = 0


︰切换到



NChlPZD




N2ChlPZD2


的生物地球化学模型过程初始条件和


边界数据



?



makebioebus = 0


︰切换到



BioEBUS


生物地球化学模型过程初始条件和边界数据



?



makepisces = 0


︰切换到双鱼座生物地球化学模型过程初始条件和边界数据



?



makeoa = 1


︰开关来定义如果



OA


(客观分析


; roms_


)文件生成的。这应该是



1



OA


文件是



ROM


水平网格,而



z


垂直网格水文数据的存储位置的中间文件。


S


坐标转化再完


成。此文件不使用



ROM




?



makeZbry = 1


︰开关,如果边界上



z


坐标



(roms_bry_)


的中间文件生成定义。仅用



make_bry




?



insitu2pot = 1


︰如果它提供的原位温度开关 定义。在这种情况下,在现场温度转化为潜在


的温度。



?



tini = 0


︰初始化在气候学实验中的一天



(15 = 1




15



)




?



woa_dir = [ROMSTOOLS_dir



' WOA2009 /']:


海洋地图集



2009


年世界气候(


Conkright


等 人,


2002


年)所在的目录。



?



cars2009_dir = [DATADIR



' CARS2009 /'];


汽车



2009


年气候



(


邓恩



et al



2009


)



于输出目录。



?



woapisces_dir = [DATADIR



' WOAPISCES /'];



?



climato_dir = cars2009_dir;



?



chla_dir = [ROMSTOOLS_dir



' SeaWifs /']:


表面叶绿素季节性气候学目录。



?



woa_time=(15:30:345)


︰在世界海洋地 图集每月气候学天时间。它对应于每个月的中间。


ROM


使用此变量在时间线性插值的气候变量。



?



woa_cycle = 360


︰持续时间的循环气候变量 。在这里,为了简单起见,我们的模型上运行


重复气候



360


年天。



?



qbar_time = [15:30:365]:


时间在天为傣族和伯斯径流每月气候学。



?



qbar_cycle = 360


︰持续时间的循环气候变量。



?



psource_ts = 0


︰开关来定义或不示踪剂(< /p>


t




s




biogechemical


)的浓度。它是在



cpp


关键



#PSOURCE_NCFILE_TS


的源代码相关联



保存



romstools_param.m




make_clim




Matlab


会话中运行︰




>>


make_clim




您应该获得︰


< br>制作中


:



/Run/ROMS _FILES/roms_





标题︰本格拉测试模型





在网格中读


......



创建气候学文件


...



创建文件


:



/Run/R OMS_FILES/roms_



...



让几个情节


......





此程序还可以很长的时间来运行。图



(1.3)


提出了



4


温度为本格拉示例的初始状态文件的不同

< br>部分。部分是在



X


方向(东


-


西),第一节是为南部的域部分,最后一个是域的北部地区。






1.3


︰结果的本格拉例如



make_clim.m



使用气候学 文件的替代方法是创建一个边界文件。在这种情况下,只有边界值存储。应定义



FRC_BRY




cpp


键和重新编译的



ROM


。在



Matlab


会话中运行



make_bry:



>>


make_bry




您应该获得︰



使该文件


:



/Run/ROMS_FILES/roms_





标题︰本格拉测试模型





在网格中读


......



...



得到的径流强迫



%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


%



%4-


开放边界和初始条件参数



使用


make_clim.m



make_biol.m



make_br y.m %



%make_ogcm.m



make_OGCM_frcst.m



%



%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%



... ... ...



... ... ..



...%


集时间和周期为径流条件︰



%


每月气候学



qbar_time = [15:30:365];



qbar_cycle = 360;



psource_ts = 0;



如果


psource_ts



%


定义人工示踪剂(


t



s


和最终生物地球化学示踪剂



%


浓度



temp_src0 = [11 9 9 12 20 20 24 25 21 18 13 12];



temp_src(:,:) = [temp_src0; temp_src0 + 2; temp_src0 + 2.8];



%



salt_src0 = [2 3 5 1 5 3 2 1 4 2 1 2];



salt_src(:,:) = [salt_src0; salt_src0; salt_src0];



%



no3_src0 = [0 0 0 0 0 0 0 0 0 0 0 0];



no3_src(:,:) = [no3_src0; no3_src0 + 2; no3_src0 + 2.8];



%



temp_src_time = [15:30:365];



temp_src_cycle = 360;



salt_src_time = [15:30:365];



salt_src_cycle = 360;



其他



temp_src_time = [];temp_src_cycle = [];



salt_src_time = [];salt_src_cycle = [];



结束



%


< /p>


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%



%



>>


make_runoff



您应该获得︰



创建径流迫使从傣族和伯斯的全球每月气候径流数据集



河域中︰



河的数量


= 2



域包含河


=



-


橙色


(Senqu)


在海洋


ATL


中流动



-


在海洋


ATL


中流动的多灵



第一个猜测︰



============





-


处理河


#1


︰橙色


(Senqu)



位置是我和


J = 36 approximetly = 27



lon src


中网格(


rho


点)


~16.6667



lat src


中网格(


rho


点)


~-28.2677





-


处理河


#2


︰多灵



位置是我和


approximetly J = 21 = 36



lon src


中网格(


rho


点)


~19.6667



lat src


中网格(


rho


点)


~-32.5713



你想使用河


(



[1]



No[0])


?橙色


(Senqu)





... ... ... ...



... ... ... ...



... ... ... ...





在径流文件写入



......


河职位



......


放电






psource_ncfile


节中的



文件中输入的行︰



-------- -------------------------------------------------- -------



psource_ncfile: Nsrc Isrc Jsrc Dsrc qbardir Lsrc


橡径流文件名称



ROMS_FILES/roms_



2



24 35 0-1 30 * T 5.0 0.0



35 13 1-1 30 * T 5.0 0.0



-------------------- ---------------------------------------------

< p>




编译模型



一旦所有的



netcdf


数据文件都准备好了(即



XXX_




XXX_




XXX_




XXX_


),我们可以准备编译的



ROM


。所有都是在



~/Roms_tools/Run/


目录中完成。



配置参数︰



param.h



编辑文件



~ / Roms_too ls/Run/param.h


。需要更改的行是︰




# elif


定义



BENGUELA_LR



参数



(LLm0 = 41



MMm0 = 42



N = 32)


!南部的格拉测试用例



#


其他




这些都是模型网格的大小值︰



LLm0


点在



X


方向,


MMm0


点的



y


轴方向和



N


垂直水平。通


过运行



make_grid.m


,给出了



LLm0




MMm0




N




romstools_param.m


中定义。在第



2.4


节中详细介绍了



param.h


参数



数值模拟和物理选项︰



cppdefs.h



要编辑的第二个文件是



~ / Ro ms_tools/Run/cppdefs.h



此文件定义 通过使用



CPP




C-


预处 理


器编译



ROM


时。


C


预处理程序选择需要根据已定义的



CPP


选项编译的



Fortran


代码的不同


部分。< /p>


这些选项被分为两个部分


(基本选项键和高级的选项键)




cppdefs.h


中。


Roms_Agrif



户指南中所述的所有键和他们的组织。



编译脚本︰



jobcomp



ROM


可以编译通过运行



UNIX bash


脚本



~/Roms_to ols/Run/jobcomp



Jobcomp


应该能够识


别您的系统。它已经过测试的



Linux




IBM


、太阳和康柏公司的系统。


Linux Pc


上默认编译器是



GNU 77


国集团,但它是可能取消



jobcomp


使用



g95




ifort


中的特定行的注释。后者是强制

< br>性时,使用



AGRIF




OPEN_MP

< br>。更改编译器时,您应该提供一个相应的



NetCDF


库。一旦


编译完毕,你应该在



/Roms_tools/Run


目录中获得一个新的可执行文件



(rom)


。应重新编译的



ROM


,修改每个时间



param.h




cppdefs.h


。如果你编译使用



MPI


并行,


jobcomp


以股代息


检测到它,并将编译器设置为



OpenMPI


,所以你需要把它安装



()



运行模式



编辑输入的参数文件


: ~ / Roms_tools/Ru n/



垂直网格参数



(THETA_S



THETA_B

< br>,


HC)


应与



romstools_param.m


中的相同。否则,不应更改其他的默认值。在每个



ROM


仿真开始时


给出的所有输入的 变量定义。要运行模型,在目录



~/Roms_tools/Run/


中键入


:./ rom




如果您使用



paralelle


计算,需要一些更具体的命令︰在



OpenMP


并行化,设置环境变量



OMP_NUM_THREADS



number_of_cpu_used


(


例如出口



OMP_NUM_THREADS = 4 4 cpu


并行运行


)


然后


./ rom


。在



MPI


并行的情况下使用下面的命令︰



mpirun-np


number_of_processus_used


./ rom




名单




的描述是一时半第



2.6


节的详细信息。



在屏幕上,您应检查



Cu_max


参数︰如果它是大于



1


你就违反



CFL

< br>标准。在这种情况下,你


就应该减少的时间步长。



模型运行的示例︰




./ rom




您应该获得︰



本格拉测试模型



720 ntimes



timesteps 3D


方程的总数。



3600.00 dt


步长


[



] 3D


方程



60 ndtfast



2D timesteps


内每个


3D


的步数。



1



timesteps


之间运行诊断程序的数目。





6.000E + 00 theta_s S


坐标表面控制参数。



0.000E + 00 theta_b S


坐标底部控制参数。



1.500E + 02 Tcline S


坐标面


/


底部层宽度中使用



伸展,米的垂直坐标。



网格文件︰


ROMS_FILES/roms_





警告︰无法识别的关键字︰迫使


--> DISREGARDED






大容量数据文件︰


ROMS_FILES/roms_





警告︰无法识别的关键字︰气候学


--> DISREGARDED






边界文件︰


ROMS_FILES/roms_



初始状态文件︰


ROMS_FILES/roms_


记录︰


1



重新启动文件︰


ROMS_FILES/roms_


奥动


= 720 rec/


文件


:-1



历史文件︰


ROMS_FILES/roms_


新建︰


T nwrt = 72 rec/


文件


= 0



1 ntsavg


起始步长对输出的积累



时间平均数据。



72 navg timesteps


之间的平均时间为写作数目



到平均文件的数据。



平均︰


ROMS_FILES/roms_ rec/


文件


= 0





要将保存在历史记录中的字段


: (T/F)



T



zeta


自由表面。



T


写乌巴


2D U


势头组件。



T



VBAR 2D V


势头组件。



T



U 3D U-


动量分量。



T



V 3D V-


动量分量。




T T (1)


示踪剂的索引为


1





T T (2)


示踪指标为


2






F



RHO


密度异常。



F


写欧米茄欧米茄垂直速度。



T



W


真实的垂直速度。



F



Akv


垂直粘度。



T



Akt


垂直扩散系数的温度。



F



Aks


垂直扩散系数的盐度。



T



Hbl


的印尼模型边界层的深度。



T



Hbbl Depth


的底部行星边界层。





T



shflx_rlw [W/m2]


长波热通量。



T



shflx_lat [W/m2]


潜热通量。



T



shflx_sen [W/m2]


感热通量



T



Bostr


底应力。



T



Wstress


风应力。



T U Wstress


写工时


U


风应力。



T V Wstress


写工时


V


风应力。



T



Shflx [W/m2]


表面净热通量





T



Swflx [


厘米


/



]


地表淡水通量


(E-P)



T



Shflx_rsw [W/m2]


短波表面辐射



T



Visc3d


横向扩散系数。







警告︰无法识别的关键字︰


gls_history_fields--> DISREGARDED








字段可保存在平均文件


: (T/F)



T



zeta


自由表面。



T


写乌巴


2D U


势头组件。



T



VBAR 2D V


势头组件。



T



U 3D U-


动量分量。



T



V 3D V-


动量分量。




T T (1)


示踪剂的索引为


1





T T (2)


示踪指标为


2






F



RHO


密度异常



T


写欧米茄欧米茄垂直速度。



T



W


真实的垂直速度。



F



Akv


垂直粘度



T



Akt


垂直扩散系数的温度。



F



Aks


垂直扩散系数的盐度。



T


写宏宝深度的印尼模型边界层



T



Hbbl Depth


的底部行星边界层





T



shflx_rlw [W/m2]


长波热通量。



T



shflx_lat [W/m2] Latente


热通量。



T



shflx_sen [W/m2]


感热通量。



T



Bostr


底应力。



T



Wstr


风应力。



T U Wstress


写工时


U


风应力。



T V Wstress


写工时


V


风应力。



T



Shflx [W/m2]


表面净热通量。





T



Swflx [


厘米


/



]


地表淡水通量


(E-P)



T



Shflx_rsw [W/m2]


短波表面辐射。



T



diff3d


横向扩散系数







警告︰无法识别的关键字︰


gls_averages--> DISREGARDED






1025.0000 rho0 Boussinesq < /p>


近似平均密度,


kg/m3


< p>


0.000E + 00 visc2


水平拉普拉斯混合系数


[m2/s]



势头。





警告︰无法识别的关键字︰


lateral_visc--> DISREGARDED






0.000E + 00 tnu2



1


)水平拉普拉斯混合系数


(m2/s)



示踪剂


1




0.000E + 00 tnu2



2


)水平拉普拉斯混合系 数


(m2/s)



示踪剂


2




0.000E + 00 tnu4



1


)水平双调和混合系数


[m4/s]



示踪剂


1




0.000E + 00 tnu4



2)


水平双调和混合系数


[m4/s]



示踪剂


2






警告︰无法识别的关键字︰


vertical_mixing--> DISREGARDED






3.000E-04 rdrg


线性底部阻力系数


(m/si)




0.000E + 00 rdrg2


二次底部阻力系数。



0.000E + 00 Zob


底部粗糙度为对数律


(m)




1.000E-04 Cdb_min


最低底部阻力系数。



1.000E-01 Cdb_max


最大底部阻力系数。





1.00 gamma2


滑溜参数︰免费滑


+ 1


或无滑移


-1






SPONGE_GRID


定义︰海绵层参数



一般是在


set_nudgcof


中设置的。


F


例程





1.157E-05 tauT_in Nudging


系数


[sec ^-1]



3.215E-08 tauT_out Nudging


系数


[sec ^-1]



3.858E-06 tauM_in Nudging


系数


[sec ^-1]



3.215E-08 tauM_out Nudging


系数


[sec ^-1]





激活


C


预处理选项︰





区域



BENGUELA_LR



OBC_EAST



OBC_WEST



OBC_NORTH



OBC_SOUTH



CURVGRID



球形



掩蔽



NEW_S_COORD



SOLVE3D



UV_COR



UV_ADV



TS_HADV_RSUP3



TS_HADV_C4



TS_DIF4



海绵



SPONGE_GRID



LMD_MIXING



LMD_SKPP



LMD_SKPP2005



LMD_BKPP



LMD_RIMIX



LMD_CONVEC



LMD_NONLOCAL



盐度



NONLIN_EOS



SPLIT_EOS



BULK_FLUX



BULK_FAIRALL



BULK_LW



BULK_EP



BULK_SMFLUX



FRC_BRY



Z_FRC_BRY



M2_FRC_BRY



M3_FRC_BRY



T_FRC_BRY



ANA_BSFLUX



ANA_BTFLUX



OBC_M2CHARACT


-


-


-


-


-


-


-


-



本文更新与2021-03-02 16:53,由作者提供,不代表本网站立场,转载请注明出处:https://www.bjmy2z.cn/gaokao/691017.html

ROMSTOOLS 中文翻译的相关文章