在cesm中寻找科氏参数，或者纬度的技术思路看fortran代码(4)

最新推荐文章于 2023-03-06 20:18:32 发布

mayubins

最新推荐文章于 2023-03-06 20:18:32 发布

阅读量1.1k

点赞数 1

分类专栏：地球系统模式文章标签： fortran linux

本文链接：https://blog.csdn.net/mayubins/article/details/123799368

版权

地球系统模式专栏收录该内容

28 篇文章 13 订阅

订阅专栏

因为计算边界层高度的公式要用到科氏参数，我没想到这居然是最难找的一个变量。

我们先看看能不能简单计算

$\LARGE f=2\Omega sin\varphi$

其中\varphi是纬度

\Omega是地球自转速度。

\Omega是常量

也就是说实际上我只需要找到纬度传进来即可

再次使用grep大法！

我们看到实际上有很多latvals

我们搜索一下latvals

in radians 表示弧度制

in degrees 表示度数制

我们用弧度制好一点

我们看到这个函数用了latvals

我们看到这个文件里应用了这个函数

我们看到这前面没有定义从哪来的latvals

再往前找

我在pbl_util 里写了关于计算f的代码

结果编译出了问题

应该是因为不能在这写运算语言

数组表达式的形状不符合。 [F]

改了之后还有这个问题，应该是数组定义没定义好

/public/home/chengxl/cesm/components/cam/src/physics/cam/vertical_diffusion.F90(988): error #6404: This name does not have a type, and must have an explicit type.   [LATVALS]

/public/home/chengxl/cesm/components/cam/src/physics/cam/vertical_diffusion.F90(988): error #6514: Substring or array slice notation requires CHARACTER type or array.   [LATVALS]

/public/home/chengxl/cesm/components/cam/src/physics/cam/vertical_diffusion.F90(988): error #6633: The type of the actual argument differs from the type of the dummy argument.   [LATVALS]

/public/home/chengxl/cesm/components/cam/src/physics/cam/vertical_diffusion.F90(1029): error #6514: Substring or array slice notation requires CHARACTER type or array.   [LATVALS]

/public/home/chengxl/cesm/components/cam/src/physics/cam/vertical_diffusion.F90(1029): error #6633: The type of the actual argument differs from the type of the dummy argument.   [LATVALS]

接下来这些是调用的问题了，分别是：需要定义，切片的话需要数组，实参和虚参类型不同。

也就是说，vertical——diffusion里面没有latvals这个变量。

我得从别的文件中找到这个变量。

我试一下在vertical_diffusion 里use const_init,这样不知道行不行呀

这样包了很多错，看起来不行的

我加一个only latvals,这样感觉不行，我需要的是输出那个气柱的纬度，但是这样似乎不是那样的