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