source: LMDZ5/trunk/libf/dyn3d_common/adaptdt.F @ 2609

Last change on this file since 2609 was 2603, checked in by Ehouarn Millour, 8 years ago

Cleanup in the dynamics: turn logic.h into module logic_mod.F90
EM
  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 1.2 KB
Line 
1!
2! $Id: adaptdt.F 2603 2016-07-25 09:31:56Z acozic $
3!
4      subroutine adaptdt(nadv,dtbon,n,pbaru,
5     c                   masse)
6
7      USE comconst_mod, ONLY: dtvr
8      IMPLICIT NONE
9
10      include "dimensions.h"
11      include "paramet.h"
12      include "comdissip.h"
13      include "comgeom2.h"
14      include "ener.h"
15      include "description.h"
16
17c----------------------------------------------------------
18c     Arguments
19c----------------------------------------------------------
20      INTEGER n,nadv
21      REAL dtbon
22      REAL pbaru(iip1,jjp1,llm)
23      REAL masse(iip1,jjp1,llm)
24c----------------------------------------------------------   
25c     Local
26c----------------------------------------------------------
27      INTEGER i,j,l
28      REAL CFLmax,aaa,bbb
29     
30        CFLmax=0.
31        do l=1,llm
32         do j=2,jjm
33          do i=1,iim
34             aaa=pbaru(i,j,l)*dtvr/masse(i,j,l)
35             CFLmax=max(CFLmax,aaa)
36             bbb=-pbaru(i,j,l)*dtvr/masse(i+1,j,l)
37             CFLmax=max(CFLmax,bbb)
38          enddo
39         enddo
40        enddo             
41        n=int(CFLmax)+1
42c pour reproduire cas VL du code qui appele x,y,z,y,x
43c        if (nadv.eq.30) n=n/2   ! Pour Prather
44        dtbon=dtvr/n
45       
46       return
47       end
48
49
50
51
52
53
54
Note: See TracBrowser for help on using the repository browser.