libf/dyn3d/adaptdt.f

!
! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/adaptdt.F,v 1.1.1.1 2004/05/19 12:53:05 lmdzadmin Exp $
!
      subroutine adaptdt(nadv,dtbon,n,pbaru,
     c                   masse)

      use dimens_m
      use paramet_m
      use comconst
      use comvert
      use conf_gcm_m
      use logic
      use comgeom
      use temps
      use ener
            use comdissip
      IMPLICIT NONE


c----------------------------------------------------------
c     Arguments
c----------------------------------------------------------
      INTEGER n,nadv
      REAL dtbon 
      REAL pbaru(iip1,jjp1,llm)
      REAL masse(iip1,jjp1,llm)
c----------------------------------------------------------    
c     Local
c----------------------------------------------------------
      INTEGER i,j,l
      REAL CFLmax,aaa,bbb
      
        CFLmax=0.
        do l=1,llm
         do j=2,jjm
          do i=1,iim
             aaa=pbaru(i,j,l)*dtvr/masse(i,j,l)
             CFLmax=max(CFLmax,aaa)
             bbb=-pbaru(i,j,l)*dtvr/masse(i+1,j,l)
             CFLmax=max(CFLmax,bbb)
          enddo
         enddo
        enddo              
        n=int(CFLmax)+1
c pour reproduire cas VL du code qui appele x,y,z,y,x
c        if (nadv.eq.30) n=n/2   ! Pour Prather
        dtbon=dtvr/n
        
       return
       end


1	!
2	! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/adaptdt.F,v 1.1.1.1 2004/05/19 12:53:05 lmdzadmin Exp $
3	!
4	subroutine adaptdt(nadv,dtbon,n,pbaru,
5	c masse)
6
7	use dimens_m
8	use paramet_m
9	use comconst
10	use comvert
11	use conf_gcm_m
12	use logic
13	use comgeom
14	use temps
15	use ener
16	use comdissip
17	IMPLICIT NONE
18
19
20	c----------------------------------------------------------
21	c Arguments
22	c----------------------------------------------------------
23	INTEGER n,nadv
24	REAL dtbon
25	REAL pbaru(iip1,jjp1,llm)
26	REAL masse(iip1,jjp1,llm)
27	c----------------------------------------------------------
28	c Local
29	c----------------------------------------------------------
30	INTEGER i,j,l
31	REAL CFLmax,aaa,bbb
32
33	CFLmax=0.
34	do l=1,llm
35	do j=2,jjm
36	do i=1,iim
37	aaa=pbaru(i,j,l)*dtvr/masse(i,j,l)
38	CFLmax=max(CFLmax,aaa)
39	bbb=-pbaru(i,j,l)*dtvr/masse(i+1,j,l)
40	CFLmax=max(CFLmax,bbb)
41	enddo
42	enddo
43	enddo
44	n=int(CFLmax)+1
45	c pour reproduire cas VL du code qui appele x,y,z,y,x
46	c if (nadv.eq.30) n=n/2 ! Pour Prather
47	dtbon=dtvr/n
48
49	return
50	end
51
52
53
54
55
56
57