/[lmdze]/trunk/phylmd/suphec.f
ViewVC logotype

Diff of /trunk/phylmd/suphec.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 169 by guez, Mon Sep 14 17:13:16 2015 UTC revision 171 by guez, Tue Sep 29 19:48:59 2015 UTC
# Line 11  module suphec_m Line 11  module suphec_m
11    real, parameter:: NAVO = 6.0221367E23 ! Avogadro number, in mol-1    real, parameter:: NAVO = 6.0221367E23 ! Avogadro number, in mol-1
12    
13    ! A1.1 Astronomical constants    ! A1.1 Astronomical constants
14    REAL RSIYEA, RSIDAY, ROMEGA  
15      REAL ROMEGA
16    real, parameter:: RDAY = 86400.    real, parameter:: RDAY = 86400.
17    real, parameter:: REA = 149597870000.  
18    real, parameter:: REPSM = 0.409093    REAL, parameter:: RSIYEA = 365.25 * RDAY * 2. * PI / 6.283076
19      ! sideral year, in s
20    
21      REAL, parameter:: RSIDAY = RDAY / (1. + RDAY / RSIYEA) ! sideral day, in s
22    
23    ! A1.2 Geoide    ! A1.2 Geoide
24    real, parameter:: RG = 9.80665 ! acceleration of gravity, in m s-2    real, parameter:: RG = 9.80665 ! acceleration of gravity, in m s-2
# Line 77  contains Line 81  contains
81    
82      ! 2. DEFINE ASTRONOMICAL CONSTANTS      ! 2. DEFINE ASTRONOMICAL CONSTANTS
83    
     RSIYEA = 365.25*RDAY*2.*PI/6.283076  
     RSIDAY = RDAY/(1.+RDAY/RSIYEA)  
84      ROMEGA = 2.*PI/RSIDAY      ROMEGA = 2.*PI/RSIDAY
85    
86      print *, 'Astronomical constants '      print *, 'Astronomical constants '
     print '('' day = '', E13.7, '' s'')', RDAY  
     print '('' half g. axis = '', E13.7, '' m'')', REA  
     print '('' mean anomaly = '', E13.7, '' -'')', REPSM  
     print '('' sideral year = '', E13.7, '' s'')', RSIYEA  
     print '('' sideral day = '', E13.7, '' s'')', RSIDAY  
87      print '('' omega = '', E13.7, '' s-1'')', ROMEGA      print '('' omega = '', E13.7, '' s-1'')', ROMEGA
88    
89      ! 3. DEFINE GEOIDE.      ! 3. DEFINE GEOIDE.
Legend:
Removed from v.169  
changed lines
  Added in v.171
  ViewVC Help
Powered by ViewVC 1.1.21