trunk/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(dtbon, n, pbaru, masse)

  USE dimens_m
  USE paramet_m
  USE comconst
  USE disvert_m
  USE conf_gcm_m
  USE conf_gcm_m
  USE comgeom
  USE temps
  IMPLICIT NONE


  ! ----------------------------------------------------------
  ! Arguments
  ! ----------------------------------------------------------
  INTEGER n
  REAL dtbon
  REAL, INTENT (IN) :: pbaru(iip1, jjp1, llm)
  REAL masse(iip1, jjp1, llm)
  ! ----------------------------------------------------------
  ! Local
  ! ----------------------------------------------------------
  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)
      END DO
    END DO
  END DO
  n = int(cflmax) + 1
  dtbon = dtvr/n

  RETURN
END SUBROUTINE adaptdt


1	guez	3
2	guez	81	! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/adaptdt.F,v 1.1.1.1 2004/05/19
3			! 12:53:05 lmdzadmin Exp $
4	guez	3
5	guez	104	SUBROUTINE adaptdt(dtbon, n, pbaru, masse)
6	guez	3
7	guez	81	USE dimens_m
8			USE paramet_m
9			USE comconst
10			USE disvert_m
11			USE conf_gcm_m
12			USE conf_gcm_m
13			USE comgeom
14			USE temps
15			IMPLICIT NONE
16	guez	3
17
18	guez	81	! ----------------------------------------------------------
19			! Arguments
20			! ----------------------------------------------------------
21	guez	104	INTEGER n
22	guez	81	REAL dtbon
23			REAL, INTENT (IN) :: pbaru(iip1, jjp1, llm)
24			REAL masse(iip1, jjp1, llm)
25			! ----------------------------------------------------------
26			! Local
27			! ----------------------------------------------------------
28			INTEGER i, j, l
29			REAL cflmax, aaa, bbb
30	guez	3
31	guez	81	cflmax = 0.
32			DO l = 1, llm
33			DO j = 2, jjm
34			DO i = 1, iim
35			aaa = pbaru(i, j, l)*dtvr/masse(i, j, l)
36			cflmax = max(cflmax, aaa)
37			bbb = -pbaru(i, j, l)*dtvr/masse(i+1, j, l)
38			cflmax = max(cflmax, bbb)
39			END DO
40			END DO
41			END DO
42			n = int(cflmax) + 1
43			dtbon = dtvr/n
44	guez	3
45	guez	81	RETURN
46			END SUBROUTINE adaptdt
47	guez	3
48
49
50	guez	81
51
52
53