fortran基本函数

FORTRAN 90标准函数（一）

(2012-07-03 17:14:57)

转载▼

分类：学习

标签：

fortran

函数

教育

符号约定：

●I代表整型;R代表实型;C代表复型;CH代表字符型;S代表字符串;L代表逻辑型;A代表数组;P代表指针;T代表派生类型;AT为任意类型。

●s:P表示s类型为P类型(任意kind值)。s:P(k)表示s类型为P类型(kind值=k)。

●[…]表示可选参数。

●*表示常用函数。

注：三角函数名前有C、D的函数为复数、双精度型函数。

注：指数函数名、平方根函数名、对数函数名前有C、D的函数为复数、双精度型函数。

表4 参数查询函数

atan2函数的值域是多少？我从网上找到一个fortran函数的日志，说此值域是-π~π，但正常反正切函数的值域应该是-π/2~π/2。对atan2函数不够了解，所以不知道你的答案对不对，我个人认为不对。我是用正常的反正切函数atan(v/u)来算的：

FORTRAN：

if (u>0..and.v>0.) dir=270-atan(v/u)*180/pi

if (u<0..and.v>0.) dir=90-atan(v/u)*180/pi

if (u<0..and.v<0.) dir=90-atan(v/u)*180/pi

if (u>0..and.v<0.) dir=270-atan(v/u)*180/pi

if (u==0..and.v>0.) dir=180

if (u==0..and.v<0.) dir=0

if (u>0..and.v==0.) dir=270

if (u<0..and.v==0.) dir=90

if (u==0..and.v==0.) dir=999

其中uv等于零的五种情况要单独挑出来，不然程序会有瑕疵。atan函数换成atand函数的话直接是度数，不用*180/pi

我四个象限和轴都试了，应该没错。

最需要注意的问题，一个是函数值域，另一个是uv矢量方向和风向是反着的，并且风向角度数是从正Y轴开始顺时针算，和三角函数里度数从正X轴开始逆时针算不一样。