| 114 |
|
|
| 115 |
c calcul des pentes limites aux poles |
c calcul des pentes limites aux poles |
| 116 |
|
|
|
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.) |
|
|
|
|
|
|
|
| 117 |
c cas ou on a un extremum au pole |
c cas ou on a un extremum au pole |
| 118 |
|
|
| 119 |
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.) |
Parent Directory
| 
Revision Log
| 
Patch