source: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/adaptdt.F @ 5100

Last change on this file since 5100 was 5099, checked in by abarral, 4 months ago

Replace most uses of CPP_DUST by the corresponding logical defined in lmdz_cppkeys_wrapper.F90
Convert several files from .F to .f90 to allow Dust to compile w/o rrtm/ecrad
Create lmdz_yoerad.f90
(lint) Remove "!" on otherwise empty line
  • 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 5099 2024-07-22 19:29:09Z abarral $
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.