source: LMDZ6/branches/blowing_snow/libf/dyn3d_common/adaptdt.F @ 5420

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

Some code tidying: turn ener.h into ener_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 2622 2016-09-04 06:12:02Z fhourdin $
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 "description.h"
15
16c----------------------------------------------------------
17c     Arguments
18c----------------------------------------------------------
19      INTEGER n,nadv
20      REAL dtbon
21      REAL pbaru(iip1,jjp1,llm)
22      REAL masse(iip1,jjp1,llm)
23c----------------------------------------------------------   
24c     Local
25c----------------------------------------------------------
26      INTEGER i,j,l
27      REAL CFLmax,aaa,bbb
28     
29        CFLmax=0.
30        do l=1,llm
31         do j=2,jjm
32          do i=1,iim
33             aaa=pbaru(i,j,l)*dtvr/masse(i,j,l)
34             CFLmax=max(CFLmax,aaa)
35             bbb=-pbaru(i,j,l)*dtvr/masse(i+1,j,l)
36             CFLmax=max(CFLmax,bbb)
37          enddo
38         enddo
39        enddo             
40        n=int(CFLmax)+1
41c pour reproduire cas VL du code qui appele x,y,z,y,x
42c        if (nadv.eq.30) n=n/2   ! Pour Prather
43        dtbon=dtvr/n
44       
45       return
46       end
47
48
49
50
51
52
53
Note: See TracBrowser for help on using the repository browser.