/[lmdze]/trunk/dyn3d/comconst.f90
ViewVC logotype

Annotation of /trunk/dyn3d/comconst.f90

Parent Directory Parent Directory | Revision Log Revision Log


Revision 328 - (hide annotations)
Thu Jun 13 14:40:06 2019 UTC (5 years ago) by guez
File size: 1193 byte(s)
Change all `.f` suffixes to `.f90`. (The opposite was done in revision
82.)  Because of change of philosopy in GNUmakefile: we already had a
rewritten rule for `.f`, so it does not make the makefile longer to
replace it by a rule for `.f90`. And it spares us options of
makedepf90 and of the compiler. Also we prepare the way for a simpler
`CMakeLists.txt`.

1 guez 3 module comconst
2    
3 guez 297 use nr_util, only: twopi
4 guez 3
5     implicit none
6    
7     real, parameter:: daysec = 86400. ! number of seconds per day
8 guez 79
9 guez 297 REAL, protected:: dtvr ! time step for dynamics, in s
10     REAL, protected:: dtphys ! time step for physics, in s
11 guez 79
12 guez 3 real, parameter:: rad = 6371229. ! radius of the Earth (in m)
13 guez 79
14     real, parameter:: cpp = 1004.70885
15     ! specific heat capacity at constant pressure of dry air, in J K-1 kg-1
16    
17     real, parameter:: kappa = 0.2857143 ! r / cpp
18    
19     real, parameter:: r = cpp * kappa
20     ! specific ideal gas constant for dry air, in J K-1 kg-1
21    
22 guez 3 real, parameter:: g = 9.8 ! acceleration of gravity (in m s-2)
23    
24 guez 297 real, parameter:: omeg = twopi / daysec
25 guez 39 ! angular speed of rotation of the Earth (in rad s-1)
26 guez 3
27 guez 297 private twopi
28 guez 3
29 guez 79 contains
30    
31     SUBROUTINE iniconst
32    
33     ! From dyn3d/iniconst.F,v 1.1.1.1 2004/05/19 12:53:05
34     ! P. Le Van
35    
36     USE conf_gcm_m, ONLY: day_step, iphysiq
37    
38     IMPLICIT NONE
39    
40     !-----------------------------------------------------------------------
41    
42     dtvr = daysec / real(day_step)
43     dtphys = iphysiq * dtvr
44    
45     print *, 'dtvr = ', dtvr
46     print *, 'dtphys = ', dtphys
47     PRINT *, 'R = ', r
48    
49     END SUBROUTINE iniconst
50    
51 guez 3 end module comconst

  ViewVC Help
Powered by ViewVC 1.1.21