我不希望发布但不确定这个问题是否或不属于这里https://codereview.stackexchange.com/questions/185368/numerical-differentiation-on-sphere-with-python。这是一个直接的例子的问题我问,大气科学fortran代码评审吗?似乎是无关的气象学SE编写的,但它实际上是非常相关的气象当你考虑它在风的空间导数在等压面或潜在温度纬向的方向。显然我重述了这个网站。
应对危险风暴的要求
这个问题可以这样
这是一个算法的空间导数为全球直线网格等压纬向和经向风
#格丽扫描顺序向东走然后向北#空间衍生品内部点#杜(区域)/ dx = u(东)- u(西方)# dv(经向)/ dy = v(北)- v(南)环南到北环西向东dudx = u (i + 1) - u(张)/ 2 * dlambda dvdy - v (j + 1) v (j - 1) / 2 * dphi结束循环结束循环#节处理波兰人(南北)类似