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
      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	IMPLICIT NONE
17
18
19	c----------------------------------------------------------
20	c Arguments
21	c----------------------------------------------------------
22	INTEGER n,nadv
23	REAL dtbon
24	REAL pbaru(iip1,jjp1,llm)
25	REAL masse(iip1,jjp1,llm)
26	c----------------------------------------------------------
27	c Local
28	c----------------------------------------------------------
29	INTEGER i,j,l
30	REAL CFLmax,aaa,bbb
31
32	CFLmax=0.
33	do l=1,llm
34	do j=2,jjm
35	do i=1,iim
36	aaa=pbaru(i,j,l)*dtvr/masse(i,j,l)
37	CFLmax=max(CFLmax,aaa)
38	bbb=-pbaru(i,j,l)*dtvr/masse(i+1,j,l)
39	CFLmax=max(CFLmax,bbb)
40	enddo
41	enddo
42	enddo
43	n=int(CFLmax)+1
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
46	dtbon=dtvr/n
47
48	return
49	end
50
51
52
53
54
55
56