76 |
REAL eps, x1, xo1, f, df, xdm, y1, yo1, ydm |
REAL eps, x1, xo1, f, df, xdm, y1, yo1, ydm |
77 |
REAL coslatm, coslatp, radclatm, radclatp |
REAL coslatm, coslatp, radclatm, radclatp |
78 |
REAL, dimension(iip1, jjp1):: cuij1, cuij2, cuij3, cuij4 ! in m |
REAL, dimension(iip1, jjp1):: cuij1, cuij2, cuij3, cuij4 ! in m |
79 |
REAL, dimension(iip1, jjp1):: cvij1, cvij2, cvij3, cvij4 |
REAL, dimension(iip1, jjp1):: cvij1, cvij2, cvij3, cvij4 ! in m |
80 |
REAL rlatu1(jjm), yprimu1(jjm), rlatu2(jjm), yprimu2(jjm) |
REAL rlatu1(jjm), yprimu1(jjm), rlatu2(jjm), yprimu2(jjm) |
81 |
real yprimv(jjm), yprimu(jjp1) |
real yprimv(jjm), yprimu(jjp1) |
82 |
REAL gamdi_gdiv, gamdi_grot, gamdi_h |
REAL gamdi_gdiv, gamdi_grot, gamdi_h |
332 |
+ cvij3(i, j)) |
+ cvij3(i, j)) |
333 |
cuv(i, j) = 0.5 * (cuij2(i, j) + cuij3(i, j) + cuij1(i, j + 1) & |
cuv(i, j) = 0.5 * (cuij2(i, j) + cuij3(i, j) + cuij1(i, j + 1) & |
334 |
+ cuij4(i, j + 1)) |
+ cuij4(i, j + 1)) |
335 |
unscv2_2d(i, j) = 1. / (cv_2d(i, j) * cv_2d(i, j)) |
unscv2_2d(i, j) = 1. / cv_2d(i, j)**2 |
336 |
END DO |
END DO |
337 |
DO i = 1, iim |
DO i = 1, iim |
338 |
cuvsurcv_2d(i, j) = airev_2d(i, j) * unscv2_2d(i, j) |
cuvsurcv_2d(i, j) = airev_2d(i, j) * unscv2_2d(i, j) |
356 |
DO i = 1, iim |
DO i = 1, iim |
357 |
cu_2d(i, j) = 0.5 * (cuij1(i, j) + cuij4(i + 1, j) + cuij2(i, j) & |
cu_2d(i, j) = 0.5 * (cuij1(i, j) + cuij4(i + 1, j) + cuij2(i, j) & |
358 |
+ cuij3(i + 1, j)) |
+ cuij3(i + 1, j)) |
359 |
unscu2_2d(i, j) = 1. / (cu_2d(i, j) * cu_2d(i, j)) |
unscu2_2d(i, j) = 1. / cu_2d(i, j)**2 |
360 |
cvusurcu_2d(i, j) = aireu_2d(i, j) * unscu2_2d(i, j) |
cvusurcu_2d(i, j) = aireu_2d(i, j) * unscu2_2d(i, j) |
361 |
cusurcvu_2d(i, j) = 1. / cvusurcu_2d(i, j) |
cusurcvu_2d(i, j) = 1. / cvusurcu_2d(i, j) |
362 |
cvuscugam1_2d(i, j) = cvusurcu_2d(i, j)**(-gamdi_gdiv) |
cvuscugam1_2d(i, j) = cvusurcu_2d(i, j)**(-gamdi_gdiv) |