--- trunk/libf/dyn3d/limy.f 2008/02/27 13:16:39 3 +++ trunk/libf/dyn3d/limy.f 2010/04/06 17:52:58 32 @@ -52,7 +52,6 @@ REAL SSUM integer ismax,ismin EXTERNAL SSUM, convflu,ismin,ismax - EXTERNAL filtreg data first/.true./ @@ -115,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.)