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

Diff of /trunk/Sources/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 66 by guez, Thu Sep 20 13:00:41 2012 UTC
# Line 13  Line 13 
13        use dimens_m        use dimens_m
14        use paramet_m        use paramet_m
15        use comconst        use comconst
16        use comvert        use disvert_m
17        use logic        use conf_gcm_m
18        use comgeom        use comgeom
19          USE nr_util, ONLY : pi
20        IMPLICIT NONE        IMPLICIT NONE
21  !  !
22  !  !
# Line 74  Line 75 
75        ENDIF        ENDIF
76    
77  !  !
 !PRINT*,'CALCUL EN LATITUDE'  
   
78        DO l = 1, llm        DO l = 1, llm
79  !  !
80  !   --------------------------------  !   --------------------------------
# Line 155  Line 154 
154           dyq(ip1jm+ij,l)=0.           dyq(ip1jm+ij,l)=0.
155        ENDDO        ENDDO
156    
 !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  
   
157  !   calcul des pentes limitees  !   calcul des pentes limitees
158    
159        DO ij=iip2,ip1jm        DO ij=iip2,ip1jm

Legend:
Removed from v.31  
changed lines
  Added in v.66

  ViewVC Help
Powered by ViewVC 1.1.21