/[lmdze]/trunk/dyn3d/limy.f
ViewVC logotype

Diff of /trunk/dyn3d/limy.f

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

revision 3 by guez, Wed Feb 27 13:16:39 2008 UTC revision 66 by guez, Thu Sep 20 13:00:41 2012 UTC
# Line 16  c   ------------------------------------ Line 16  c   ------------------------------------
16        use dimens_m        use dimens_m
17        use paramet_m        use paramet_m
18        use comconst        use comconst
19        use comvert        use disvert_m
20        use logic        use conf_gcm_m
21        use comgeom        use comgeom
22          USE nr_util, ONLY : pi
23        IMPLICIT NONE        IMPLICIT NONE
24  c  c
25  c  c
# Line 52  c Line 53  c
53        REAL      SSUM        REAL      SSUM
54        integer ismax,ismin        integer ismax,ismin
55        EXTERNAL  SSUM, convflu,ismin,ismax        EXTERNAL  SSUM, convflu,ismin,ismax
       EXTERNAL filtreg  
56    
57        data first/.true./        data first/.true./
58    
# Line 115  c   calcul des pentes aux poles Line 115  c   calcul des pentes aux poles
115    
116  c   calcul des pentes limites aux poles  c   calcul des pentes limites aux poles
117    
 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.)  
   
   
118  c   cas ou on a un extremum au pole  c   cas ou on a un extremum au pole
119    
120  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.)
Legend:
Removed from v.3  
changed lines
  Added in v.66
  ViewVC Help
Powered by ViewVC 1.1.21