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