59 |
REAL, intent(out):: zpic(:, :) ! Maximum altitude |
REAL, intent(out):: zpic(:, :) ! Maximum altitude |
60 |
real, intent(out):: zval(:, :) ! Minimum altitude |
real, intent(out):: zval(:, :) ! Minimum altitude |
61 |
|
|
62 |
real, intent(out):: mask(:, :) |
real, intent(out):: mask(:, :) ! fraction of land |
63 |
|
|
64 |
! Variables local to the procedure: |
! Variables local to the procedure: |
65 |
|
|
141 |
zusn(i, jusn+2)=zusn(i+iusn/2, jusn+1) |
zusn(i, jusn+2)=zusn(i+iusn/2, jusn+1) |
142 |
zusn(i+iusn/2+iext, jusn+2)=zusn(i, jusn+1) |
zusn(i+iusn/2+iext, jusn+2)=zusn(i, jusn+1) |
143 |
ENDDO |
ENDDO |
144 |
! |
|
145 |
! COMPUTE LIMITS OF MODEL GRIDPOINT AREA |
! COMPUTE LIMITS OF MODEL GRIDPOINT AREA |
146 |
! ( REGULAR GRID) |
! ( REGULAR GRID) |
147 |
|
|
190 |
ENDDO |
ENDDO |
191 |
|
|
192 |
! SUMMATION OVER GRIDPOINT AREA |
! SUMMATION OVER GRIDPOINT AREA |
193 |
! |
|
194 |
zleny=pi/real(jusn)*rad |
zleny=pi/real(jusn)*rad |
195 |
xincr=pi/2./real(jusn) |
xincr=pi/2./real(jusn) |
196 |
DO ii = 1, iim+1 |
DO ii = 1, iim+1 |
324 |
zllmval=AMAX1(zval(ii, jj), zllmval) |
zllmval=AMAX1(zval(ii, jj), zllmval) |
325 |
ENDDO |
ENDDO |
326 |
ENDDO |
ENDDO |
327 |
print *, ' MEAN ORO:', zllmmea |
print *, 'MEAN ORO: ', zllmmea |
328 |
print *, ' ST. DEV.:', zllmstd |
print *, 'ST. DEV.: ', zllmstd |
329 |
print *, ' PENTE:', zllmsig |
print *, 'PENTE: ', zllmsig |
330 |
print *, ' ANISOTROP:', zllmgam |
print *, 'ANISOTROP: ', zllmgam |
331 |
print *, ' ANGLE:', zminthe, zllmthe |
print *, 'ANGLE: ', zminthe, zllmthe |
332 |
print *, ' pic:', zllmpic |
print *, 'pic: ', zllmpic |
333 |
print *, ' val:', zllmval |
print *, 'val: ', zllmval |
334 |
|
|
335 |
! gamma and theta a 1. and 0. at poles |
! gamma and theta a 1. and 0. at poles |
336 |
zmea(iim+1, :)=zmea(1, :) |
zmea(iim+1, :)=zmea(1, :) |