LMDZ
adaptdt.F
Go to the documentation of this file.
1 !
2 ! $Id: adaptdt.F 1952 2014-01-28 13:05:47Z lguez $
3 !
4  subroutine adaptdt(nadv,dtbon,n,pbaru,
5  c masse)
6 
7  USE control_mod
8  IMPLICIT NONE
9 
10 #include "dimensions.h"
11 c#include "paramr2.h"
12 #include "paramet.h"
13 #include "comconst.h"
14 #include "comdissip.h"
15 #include "comvert.h"
16 #include "comgeom2.h"
17 #include "logic.h"
18 #include "temps.h"
19 #include "ener.h"
20 #include "description.h"
21 
22 c----------------------------------------------------------
23 c Arguments
24 c----------------------------------------------------------
25  INTEGER n,nadv
26  REAL dtbon
27  REAL pbaru(iip1,jjp1,llm)
28  REAL masse(iip1,jjp1,llm)
29 c----------------------------------------------------------
30 c Local
31 c----------------------------------------------------------
32  INTEGER i,j,l
33  REAL CFLmax,aaa,bbb
34 
35  cflmax=0.
36  do l=1,llm
37  do j=2,jjm
38  do i=1,iim
39  aaa=pbaru(i,j,l)*dtvr/masse(i,j,l)
40  cflmax=max(cflmax,aaa)
41  bbb=-pbaru(i,j,l)*dtvr/masse(i+1,j,l)
42  cflmax=max(cflmax,bbb)
43  enddo
44  enddo
45  enddo
46  n=int(cflmax)+1
47 c pour reproduire cas VL du code qui appele x,y,z,y,x
48 c if (nadv.eq.30) n=n/2 ! Pour Prather
49  dtbon=dtvr/n
50 
51  return
52  end
53 
54 
55 
56 
57 
58 
59 
subroutine adaptdt(nadv, dtbon, n, pbaru, masse)
Definition: adaptdt.F:6
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vq_toga w_prof vq_prof llm day day1 day day1 *dt_dice endif!time annee_ref dt_dice swup_dice vg_dice omega_dice tg_prof vg_profd w_profd omega_profd!do llm!print llm l llm
!$Header jjp1
Definition: paramet.h:14
!$Id mode_top_bound COMMON comconstr dtvr
Definition: comconst.h:7
c c zjulian c cym CALL iim cym klev iim
Definition: ini_bilKP_ave.h:24