/[lmdze]/trunk/Sources/dyn3d/adaptdt.f
ViewVC logotype

Diff of /trunk/Sources/dyn3d/adaptdt.f

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

trunk/libf/dyn3d/adaptdt.f revision 3 by guez, Wed Feb 27 13:16:39 2008 UTC trunk/Sources/dyn3d/adaptdt.f revision 134 by guez, Wed Apr 29 15:47:56 2015 UTC
# Line 1  Line 1 
1  !  
2  ! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/adaptdt.F,v 1.1.1.1 2004/05/19 12:53:05 lmdzadmin Exp $  ! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/adaptdt.F,v 1.1.1.1 2004/05/19
3  !  ! 12:53:05 lmdzadmin Exp $
4        subroutine adaptdt(nadv,dtbon,n,pbaru,  
5       c                   masse)  SUBROUTINE adaptdt(dtbon, n, pbaru, masse)
6    
7        use dimens_m    USE dimens_m
8        use paramet_m    USE paramet_m
9        use comconst    USE comconst
10        use comvert    USE disvert_m
11        use conf_gcm_m    USE conf_gcm_m
12        use logic    USE conf_gcm_m
13        use comgeom    USE comgeom
14        use temps    USE temps
15        use ener    IMPLICIT NONE
16              use comdissip  
17        IMPLICIT NONE  
18      ! ----------------------------------------------------------
19      ! Arguments
20  c----------------------------------------------------------    ! ----------------------------------------------------------
21  c     Arguments    INTEGER n
22  c----------------------------------------------------------    REAL dtbon
23        INTEGER n,nadv    REAL, INTENT (IN) :: pbaru(iip1, jjp1, llm)
24        REAL dtbon    REAL masse(iip1, jjp1, llm)
25        REAL pbaru(iip1,jjp1,llm)    ! ----------------------------------------------------------
26        REAL masse(iip1,jjp1,llm)    ! Local
27  c----------------------------------------------------------        ! ----------------------------------------------------------
28  c     Local    INTEGER i, j, l
29  c----------------------------------------------------------    REAL cflmax, aaa, bbb
30        INTEGER i,j,l  
31        REAL CFLmax,aaa,bbb    cflmax = 0.
32            DO l = 1, llm
33          CFLmax=0.      DO j = 2, jjm
34          do l=1,llm        DO i = 1, iim
35           do j=2,jjm          aaa = pbaru(i, j, l)*dtvr/masse(i, j, l)
36            do i=1,iim          cflmax = max(cflmax, aaa)
37               aaa=pbaru(i,j,l)*dtvr/masse(i,j,l)          bbb = -pbaru(i, j, l)*dtvr/masse(i+1, j, l)
38               CFLmax=max(CFLmax,aaa)          cflmax = max(cflmax, bbb)
39               bbb=-pbaru(i,j,l)*dtvr/masse(i+1,j,l)        END DO
40               CFLmax=max(CFLmax,bbb)      END DO
41            enddo    END DO
42           enddo    n = int(cflmax) + 1
43          enddo                  dtbon = dtvr/n
44          n=int(CFLmax)+1  
45  c pour reproduire cas VL du code qui appele x,y,z,y,x    RETURN
46  c        if (nadv.eq.30) n=n/2   ! Pour Prather  END SUBROUTINE adaptdt
         dtbon=dtvr/n  
           
        return  
        end  
47    
48    
49    
Legend:
Removed from v.3  
changed lines
  Added in v.134
  ViewVC Help
Powered by ViewVC 1.1.21