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

Diff of /trunk/dyn3d/adaptdt.f

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

trunk/libf/dyn3d/adaptdt.f revision 27 by guez, Thu Mar 25 14:29:07 2010 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        IMPLICIT NONE  
17    
18      ! ----------------------------------------------------------
19  c----------------------------------------------------------    ! Arguments
20  c     Arguments    ! ----------------------------------------------------------
21  c----------------------------------------------------------    INTEGER n
22        INTEGER n,nadv    REAL dtbon
23        REAL dtbon    REAL, INTENT (IN) :: pbaru(iip1, jjp1, llm)
24        REAL pbaru(iip1,jjp1,llm)    REAL masse(iip1, jjp1, llm)
25        REAL masse(iip1,jjp1,llm)    ! ----------------------------------------------------------
26  c----------------------------------------------------------        ! Local
27  c     Local    ! ----------------------------------------------------------
28  c----------------------------------------------------------    INTEGER i, j, l
29        INTEGER i,j,l    REAL cflmax, aaa, bbb
30        REAL CFLmax,aaa,bbb  
31            cflmax = 0.
32          CFLmax=0.    DO l = 1, llm
33          do l=1,llm      DO j = 2, jjm
34           do j=2,jjm        DO i = 1, iim
35            do i=1,iim          aaa = pbaru(i, j, l)*dtvr/masse(i, j, l)
36               aaa=pbaru(i,j,l)*dtvr/masse(i,j,l)          cflmax = max(cflmax, aaa)
37               CFLmax=max(CFLmax,aaa)          bbb = -pbaru(i, j, l)*dtvr/masse(i+1, j, l)
38               bbb=-pbaru(i,j,l)*dtvr/masse(i+1,j,l)          cflmax = max(cflmax, bbb)
39               CFLmax=max(CFLmax,bbb)        END DO
40            enddo      END DO
41           enddo    END DO
42          enddo                  n = int(cflmax) + 1
43          n=int(CFLmax)+1    dtbon = dtvr/n
44  c pour reproduire cas VL du code qui appele x,y,z,y,x  
45  c        if (nadv.eq.30) n=n/2   ! Pour Prather    RETURN
46          dtbon=dtvr/n  END SUBROUTINE adaptdt
           
        return  
        end  
47    
48    
49    
Legend:
Removed from v.27  
changed lines
  Added in v.134
  ViewVC Help
Powered by ViewVC 1.1.21