gfortran中缺失的函数：tand, dtand

在gfortran中，包括g77到现在的gfortran 4.1都缺失了一些函数实现，包括以角度为参数的三角函数等。
那我们在编译这些函数时要如何处理呢？
gcc的邮件列表里有一个解决方案。
我本来想用这个方案的，结果fortran语法全部忘光，这个方案里面居然没有看懂。没办法，就自己实现了一个简单的方案，不过还没有进行测试。

tand(x)

变为：

tan(x*atan(1.0)/45.0)

而

dtand(x)

则变为：

tan(x*atan(1.0d0)/45.0d0)

若有错误，还请帮忙指出。