1 |
module tau2alpha_m |
module tau2alpha_m |
2 |
|
|
3 |
|
USE paramet_m, ONLY : iip1, jjp1 |
4 |
|
USE dimens_m, ONLY : jjm |
5 |
|
|
6 |
IMPLICIT NONE |
IMPLICIT NONE |
7 |
|
|
8 |
REAL lat_min_guide, lat_max_guide |
private iip1, jjp1, jjm |
9 |
|
|
10 |
|
REAL dxdys(iip1, jjp1), dxdyu(iip1, jjp1), dxdyv(iip1, jjm) |
11 |
|
|
12 |
contains |
contains |
13 |
|
|
14 |
SUBROUTINE tau2alpha(type, pim, pjm, factt, taumin, taumax, alpha) |
SUBROUTINE tau2alpha(type, pim, pjm, factt, taumin, taumax, alpha) |
15 |
|
|
|
USE dimens_m, ONLY : iim, jjm |
|
|
USE paramet_m, ONLY : iip1, jjp1 |
|
|
USE nr_util, ONLY : pi |
|
16 |
USE comgeom, ONLY : cu_2d, cv_2d, rlatu, rlatv |
USE comgeom, ONLY : cu_2d, cv_2d, rlatu, rlatv |
17 |
|
use conf_guide_m, only: lat_min_guide, lat_max_guide |
18 |
|
use dump2d_m, only: dump2d |
19 |
|
USE dimens_m, ONLY : iim |
20 |
|
USE nr_util, ONLY : pi |
21 |
USE serre, ONLY : clat, clon, grossismx, grossismy |
USE serre, ONLY : clat, clon, grossismx, grossismy |
22 |
|
|
23 |
! arguments : |
! arguments : |
24 |
INTEGER type |
INTEGER type |
25 |
INTEGER pim, pjm |
INTEGER pim, pjm |
26 |
REAL factt, taumin, taumax |
REAL, intent(in):: factt, taumin, taumax |
27 |
REAL dxdy_, alpha(pim, pjm) |
REAL dxdy_, alpha(pim, pjm) |
28 |
REAL dxdy_min, dxdy_max |
REAL dxdy_min, dxdy_max |
29 |
|
|
37 |
DATA first/ .TRUE./ |
DATA first/ .TRUE./ |
38 |
|
|
39 |
REAL zdx(iip1, jjp1), zdy(iip1, jjp1) |
REAL zdx(iip1, jjp1), zdy(iip1, jjp1) |
|
|
|
40 |
REAL zlat |
REAL zlat |
41 |
REAL dxdys(iip1, jjp1), dxdyu(iip1, jjp1), dxdyv(iip1, jjm) |
|
42 |
COMMON /comdxdy/dxdys, dxdyu, dxdyv |
!------------------------------------------------------------ |
43 |
|
|
44 |
IF (first) THEN |
IF (first) THEN |
45 |
DO j = 2, jjm |
DO j = 2, jjm |
72 |
END DO |
END DO |
73 |
DO j = 1, jjm |
DO j = 1, jjm |
74 |
DO i = 1, iip1 |
DO i = 1, iip1 |
75 |
dxdyv(i, j) = 0.5*(dxdys(i, j)+dxdys(i+1, j)) |
dxdyv(i, j) = 0.5*(dxdys(i, j)+dxdys(i, j + 1)) |
76 |
END DO |
END DO |
77 |
END DO |
END DO |
78 |
|
|