16 |
use dimens_m |
use dimens_m |
17 |
use paramet_m |
use paramet_m |
18 |
use comconst |
use comconst |
19 |
use comvert |
use disvert_m |
20 |
use logic |
use conf_gcm_m |
21 |
use comgeom |
use comgeom |
22 |
|
USE nr_util, ONLY : pi |
23 |
IMPLICIT NONE |
IMPLICIT NONE |
24 |
c |
c |
25 |
c |
c |
115 |
|
|
116 |
c calcul des pentes limites aux poles |
c calcul des pentes limites aux poles |
117 |
|
|
|
c print*,dyqv(iip1+1) |
|
|
c apn=abs(dyq(1)/dyqv(iip1+1)) |
|
|
c print*,dyq(ip1jm+1) |
|
|
c print*,dyqv(ip1jm-iip1+1) |
|
|
c aps=abs(dyq(ip1jm+1)/dyqv(ip1jm-iip1+1)) |
|
|
c do ij=2,iim |
|
|
c apn=amax1(abs(dyq(ij)/dyqv(ij)),apn) |
|
|
c aps=amax1(abs(dyq(ip1jm+ij)/dyqv(ip1jm-iip1+ij)),aps) |
|
|
c enddo |
|
|
c apn=min(pente_max/apn,1.) |
|
|
c aps=min(pente_max/aps,1.) |
|
|
|
|
|
|
|
118 |
c cas ou on a un extremum au pole |
c cas ou on a un extremum au pole |
119 |
|
|
120 |
c if(dyqv(ismin(iim,dyqv,1))*dyqv(ismax(iim,dyqv,1)).le.0.) |
c if(dyqv(ismin(iim,dyqv,1))*dyqv(ismax(iim,dyqv,1)).le.0.) |