trunk/dyn3d/adaptdt.f90


! $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 dimensions
  USE paramet_m
  USE comconst
  USE disvert_m
  USE conf_gcm_m
  USE conf_gcm_m
  USE comgeom

  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
2	! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/adaptdt.F,v 1.1.1.1 2004/05/19
3	! 12:53:05 lmdzadmin Exp $
4
5	SUBROUTINE adaptdt(dtbon, n, pbaru, masse)
6
7	USE dimensions
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
15	IMPLICIT NONE
16
17
18	! ----------------------------------------------------------
19	! Arguments
20	! ----------------------------------------------------------
21	INTEGER n
22	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
31	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
45	RETURN
46	END SUBROUTINE adaptdt
47
48
49
50
51
52
53