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

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

Cleanup in the dynamics: turn temps.h into module temps_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 2601 2016-07-24 09:51:55Z emillour $
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 "logic.h"
15      include "ener.h"
16      include "description.h"
17
18c----------------------------------------------------------
19c     Arguments
20c----------------------------------------------------------
21      INTEGER n,nadv
22      REAL dtbon
23      REAL pbaru(iip1,jjp1,llm)
24      REAL masse(iip1,jjp1,llm)
25c----------------------------------------------------------   
26c     Local
27c----------------------------------------------------------
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          enddo
40         enddo
41        enddo             
42        n=int(CFLmax)+1
43c pour reproduire cas VL du code qui appele x,y,z,y,x
44c        if (nadv.eq.30) n=n/2   ! Pour Prather
45        dtbon=dtvr/n
46       
47       return
48       end
49
50
51
52
53
54
55
Note: See TracBrowser for help on using the repository browser.