如何使用fortran编写高精度的累积分布函数

yqman

本人想在WRF中的参数化方案中加入一些东西，其中需要计算正太分布的累积分布函数，已经使用误差函数和Abramowitz and Stegun提供的标准正态分布的累积分布函数近似公式，来编写累积分布函数，代码如下图，但是得出来的结果均和python中norm.cdf得出来的结果相差甚远，因为常用python，fortran不太懂，因此特来求助。