/[lmdze]/trunk/dyn3d/Vlsplt/vly.f
ViewVC logotype

Diff of /trunk/dyn3d/Vlsplt/vly.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 31 by guez, Thu Apr 1 14:59:19 2010 UTC revision 32 by guez, Tue Apr 6 17:52:58 2010 UTC
# Line 74  Line 74 
74        ENDIF        ENDIF
75    
76  !  !
 !PRINT*,'CALCUL EN LATITUDE'  
   
77        DO l = 1, llm        DO l = 1, llm
78  !  !
79  !   --------------------------------  !   --------------------------------
# Line 155  Line 153 
153           dyq(ip1jm+ij,l)=0.           dyq(ip1jm+ij,l)=0.
154        ENDDO        ENDDO
155    
 !CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC  
 !  En memoire de dIFferents tests sur la  
 !  limitation des pentes aux poles.  
 !CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC  
 !     PRINT*,dyq(1)  
 !     PRINT*,dyqv(iip1+1)  
 !     apn=abs(dyq(1)/dyqv(iip1+1))  
 !     PRINT*,dyq(ip1jm+1)  
 !     PRINT*,dyqv(ip1jm-iip1+1)  
 !     aps=abs(dyq(ip1jm+1)/dyqv(ip1jm-iip1+1))  
 !     DO ij=2,iim  
 !        apn=amax1(abs(dyq(ij)/dyqv(ij)),apn)  
 !        aps=amax1(abs(dyq(ip1jm+ij)/dyqv(ip1jm-iip1+ij)),aps)  
 !     ENDDO  
 !     apn=min(pente_max/apn,1.)  
 !     aps=min(pente_max/aps,1.)  
 !  
 !  
 !   cas ou on a un extremum au pole  
 !  
 !     IF(dyqv(ismin(iim,dyqv,1))*dyqv(ismax(iim,dyqv,1)).le.0.)  
 !    &   apn=0.  
 !     IF(dyqv(ismax(iim,dyqv(ip1jm-iip1+1),1)+ip1jm-iip1+1)*  
 !    &   dyqv(ismin(iim,dyqv(ip1jm-iip1+1),1)+ip1jm-iip1+1).le.0.)  
 !    &   aps=0.  
 !  
 !   limitation des pentes aux poles  
 !     DO ij=1,iip1  
 !        dyq(ij)=apn*dyq(ij)  
 !        dyq(ip1jm+ij)=aps*dyq(ip1jm+ij)  
 !     ENDDO  
 !  
 !   test  
 !      DO ij=1,iip1  
 !         dyq(iip1+ij)=0.  
 !         dyq(ip1jm+ij-iip1)=0.  
 !      ENDDO  
 !      DO ij=1,ip1jmp1  
 !         dyq(ij)=dyq(ij)*cos(rlatu((ij-1)/iip1+1))  
 !      ENDDO  
 !  
 ! changement 10 07 96  
 !     IF(dyqv(ismin(iim,dyqv,1))*dyqv(ismax(iim,dyqv,1)).le.0.)  
 !    &   THEN  
 !        DO ij=1,iip1  
 !           dyqmax(ij)=0.  
 !        ENDDO  
 !     ELSE  
 !        DO ij=1,iip1  
 !           dyqmax(ij)=pente_max*abs(dyqv(ij))  
 !        ENDDO  
 !     ENDIF  
 !  
 !     IF(dyqv(ismax(iim,dyqv(ip1jm-iip1+1),1)+ip1jm-iip1+1)*  
 !    & dyqv(ismin(iim,dyqv(ip1jm-iip1+1),1)+ip1jm-iip1+1).le.0.)  
 !    &THEN  
 !        DO ij=ip1jm+1,ip1jmp1  
 !           dyqmax(ij)=0.  
 !        ENDDO  
 !     ELSE  
 !        DO ij=ip1jm+1,ip1jmp1  
 !           dyqmax(ij)=pente_max*abs(dyqv(ij-iip1))  
 !        ENDDO  
 !     ENDIF  
 !   fin changement 10 07 96  
 !CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC  
   
156  !   calcul des pentes limitees  !   calcul des pentes limitees
157    
158        DO ij=iip2,ip1jm        DO ij=iip2,ip1jm
Legend:
Removed from v.31  
changed lines
  Added in v.32
  ViewVC Help
Powered by ViewVC 1.1.21