18 |
! Local: |
! Local: |
19 |
REAL alphamin, alphamax, xi |
REAL alphamin, alphamax, xi |
20 |
INTEGER i, j, n_lon, n_lat |
INTEGER i, j, n_lon, n_lat |
|
REAL zlat |
|
21 |
|
|
22 |
!------------------------------------------------------------ |
!------------------------------------------------------------ |
23 |
|
|
31 |
alphamax = factt / taumin |
alphamax = factt / taumin |
32 |
|
|
33 |
DO j = 1, n_lat |
DO j = 1, n_lat |
34 |
zlat = rlat(j) * 180. / pi |
IF (lat_min_guide <= rlat(j) .AND. rlat(j) <= lat_max_guide) THEN |
|
IF (lat_min_guide <= zlat .AND. zlat <= lat_max_guide) THEN |
|
35 |
DO i = 1, n_lon |
DO i = 1, n_lon |
36 |
xi = min(((dxdy_max - dxdy(i, j)) & |
xi = min(((dxdy_max - dxdy(i, j)) & |
37 |
/ (dxdy_max - dxdy_min))**gamma, 1.) |
/ (dxdy_max - dxdy_min))**gamma, 1.) |