--- trunk/libf/dyn3d/limy.f 2010/03/25 14:29:07 27 +++ trunk/libf/dyn3d/limy.f 2010/04/06 17:52:58 32 @@ -114,19 +114,6 @@ c calcul des pentes limites aux poles -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.) - - c cas ou on a un extremum au pole c if(dyqv(ismin(iim,dyqv,1))*dyqv(ismax(iim,dyqv,1)).le.0.)