/[lmdze]/trunk/dyn3d/adaptdt.f
ViewVC logotype

Annotation of /trunk/dyn3d/adaptdt.f

Parent Directory Parent Directory | Revision Log Revision Log

Revision 57 - (hide annotations)
Mon Jan 30 12:54:02 2012 UTC (12 years, 4 months ago) by guez
Original Path: trunk/libf/dyn3d/adaptdt.f
File size: 1307 byte(s) 
Write used namelists to file "" instead of standard output.

Avoid aliasing in "inidissip" in calls to "divgrad2", "divgrad",
"gradiv2", "gradiv", "nxgraro2" and "nxgrarot". Add a degenerate
dimension to arrays so they have rank 3, like the dummy arguments in
"divgrad2", "divgrad", "gradiv2", "gradiv", "nxgraro2" and "nxgrarot".

Extract the initialization part from "bilan_dyn" and make a separate
procedure, "init_dynzon", from it.

Move variables from modules "iniprint" and "logic" to module
"conf_gcm_m".

Promote internal procedures of "fxy" to private procedures of module
"fxy_m".

Extracted documentation from "inigeom". Removed useless "save"
attributes. Removed useless intermediate variables. Extracted
processing of poles from loop on latitudes. Write coordinates to file
"longitude_latitude.txt" instead of standard output.

Do not use ozone tracer for radiative transfer.
1 guez 3 !
2     ! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/adaptdt.F,v 1.1.1.1 2004/05/19 12:53:05 lmdzadmin Exp $
3     !
4     subroutine adaptdt(nadv,dtbon,n,pbaru,
5     c masse)
6    
7     use dimens_m
8     use paramet_m
9     use comconst
10     use comvert
11     use conf_gcm_m
12 guez 57 use conf_gcm_m
13 guez 3 use comgeom
14     use temps
15     use ener
16     IMPLICIT NONE
17    
18    
19     c----------------------------------------------------------
20     c Arguments
21     c----------------------------------------------------------
22     INTEGER n,nadv
23     REAL dtbon
24 guez 31 REAL, intent(in):: pbaru(iip1,jjp1,llm)
25 guez 3 REAL masse(iip1,jjp1,llm)
26     c----------------------------------------------------------
27     c Local
28     c----------------------------------------------------------
29     INTEGER i,j,l
30     REAL CFLmax,aaa,bbb
31    
32     CFLmax=0.
33     do l=1,llm
34     do j=2,jjm
35     do i=1,iim
36     aaa=pbaru(i,j,l)*dtvr/masse(i,j,l)
37     CFLmax=max(CFLmax,aaa)
38     bbb=-pbaru(i,j,l)*dtvr/masse(i+1,j,l)
39     CFLmax=max(CFLmax,bbb)
40     enddo
41     enddo
42     enddo
43     n=int(CFLmax)+1
44     c pour reproduire cas VL du code qui appele x,y,z,y,x
45     c if (nadv.eq.30) n=n/2 ! Pour Prather
46     dtbon=dtvr/n
47    
48     return
49     end
50    
51    
52    
53    
54    
55    
56    
  ViewVC Help
Powered by ViewVC 1.1.21