/[lmdze]/trunk/libf/phylmd/suphec.f90
ViewVC logotype

Diff of /trunk/libf/phylmd/suphec.f90

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

revision 3 by guez, Wed Feb 27 13:16:39 2008 UTC revision 17 by guez, Tue Aug 5 13:31:32 2008 UTC
# Line 24  contains Line 24  contains
24      IF (firstcall) THEN      IF (firstcall) THEN
25         PRINT *, 'suphec initialise les constantes du GCM'         PRINT *, 'suphec initialise les constantes du GCM'
26         firstcall = .FALSE.         firstcall = .FALSE.
27         !  
28         !*       1.    DEFINE FUNDAMENTAL CONSTANTS.         !*       1.    DEFINE FUNDAMENTAL CONSTANTS.
29         !  
30         WRITE(UNIT=6,FMT='(''0*** Constants of the ICM   ***'')')         WRITE(UNIT=6,FMT='(''0*** Constants of the ICM   ***'')')
31         RPI=2.*ASIN(1.)         RPI=2.*ASIN(1.)
32         RCLUM=299792458.         RCLUM=299792458.
# Line 43  contains Line 43  contains
43              RKBOL              RKBOL
44         WRITE(UNIT=6,FMT='(''            N = '',E13.7,''mol-1'')') &         WRITE(UNIT=6,FMT='(''            N = '',E13.7,''mol-1'')') &
45              RNAVO              RNAVO
46         !  
        !  
47         !*       2.    DEFINE ASTRONOMICAL CONSTANTS.         !*       2.    DEFINE ASTRONOMICAL CONSTANTS.
48         !  
49         RDAY=86400.         RDAY=86400.
50         REA=149597870000.         REA=149597870000.
51         REPSM=0.409093         REPSM=0.409093
52         !  
53         RSIYEA=365.25*RDAY*2.*RPI/6.283076         RSIYEA=365.25*RDAY*2.*RPI/6.283076
54         RSIDAY=RDAY/(1.+RDAY/RSIYEA)         RSIDAY=RDAY/(1.+RDAY/RSIYEA)
55         ROMEGA=2.*RPI/RSIDAY         ROMEGA=2.*RPI/RSIDAY
56         !  
57         WRITE(UNIT=6,FMT='('' *** Astronomical constants ***'')')         WRITE(UNIT=6,FMT='('' *** Astronomical constants ***'')')
58         WRITE(UNIT=6,FMT='(''          day = '',E13.7,'' s'')')RDAY         WRITE(UNIT=6,FMT='(''          day = '',E13.7,'' s'')')RDAY
59         WRITE(UNIT=6,FMT='('' half g. axis = '',E13.7,'' m'')')REA         WRITE(UNIT=6,FMT='('' half g. axis = '',E13.7,'' m'')')REA
# Line 63  contains Line 62  contains
62         WRITE(UNIT=6,FMT='(''  sideral day = '',E13.7,'' s'')')RSIDAY         WRITE(UNIT=6,FMT='(''  sideral day = '',E13.7,'' s'')')RSIDAY
63         WRITE(UNIT=6,FMT='(''        omega = '',E13.7,'' s-1'')') &         WRITE(UNIT=6,FMT='(''        omega = '',E13.7,'' s-1'')') &
64              ROMEGA              ROMEGA
65         !  
66         !*       3.    DEFINE GEOIDE.         !*       3.    DEFINE GEOIDE.
67         !  
68         RG=9.80665         RG=9.80665
69         RA=6371229.         RA=6371229.
70         R1SA=SNGL(1.D0/DBLE(RA))         R1SA=SNGL(1.D0/DBLE(RA))
# Line 74  contains Line 73  contains
73              RG              RG
74         WRITE(UNIT=6,FMT='('' Earth radius = '',E13.7,'' m'')')RA         WRITE(UNIT=6,FMT='('' Earth radius = '',E13.7,'' m'')')RA
75         WRITE(UNIT=6,FMT='('' Inverse E.R. = '',E13.7,'' m'')')R1SA         WRITE(UNIT=6,FMT='('' Inverse E.R. = '',E13.7,'' m'')')R1SA
76         !  
        !  
77         !*       4.    DEFINE RADIATION CONSTANTS.         !*       4.    DEFINE RADIATION CONSTANTS.
78         !  
        ! z.x.li      RSIGMA=2. * RPI**5 * RKBOL**4 /(15.* RCLUM**2 * RHPLA**3)  
79         rsigma = 2.*rpi**5 * (rkbol/rhpla)**3 * rkbol/rclum/rclum/15.         rsigma = 2.*rpi**5 * (rkbol/rhpla)**3 * rkbol/rclum/rclum/15.
80         !IM init. dans conf_phys.F90   RI0=1365.         !IM init. dans conf_phys.F90   RI0=1365.
81         WRITE(UNIT=6,FMT='('' ***        Radiation       ***'')')         WRITE(UNIT=6,FMT='('' ***        Radiation       ***'')')
82         WRITE(UNIT=6,FMT='('' Stefan-Bol.  = '',E13.7,'' W m-2 K-4'')')  RSIGMA         WRITE(UNIT=6,FMT='('' Stefan-Bol.  = '',E13.7,'' W m-2 K-4'')')  RSIGMA
83    
        !  
84         !*       5.    DEFINE THERMODYNAMIC CONSTANTS, GAS PHASE.         !*       5.    DEFINE THERMODYNAMIC CONSTANTS, GAS PHASE.
85         !  
86         R=RNAVO*RKBOL         R=RNAVO*RKBOL
87         RMD=28.9644         RMD=28.9644
88         RMO3=47.9942         RMO3=47.9942
# Line 112  contains Line 108  contains
108         WRITE(UNIT=6,FMT='(''         Cvv  = '',e13.7)') RCVV         WRITE(UNIT=6,FMT='(''         Cvv  = '',e13.7)') RCVV
109         WRITE(UNIT=6,FMT='(''      Rd/Cpd  = '',e13.7)') RKAPPA         WRITE(UNIT=6,FMT='(''      Rd/Cpd  = '',e13.7)') RKAPPA
110         WRITE(UNIT=6,FMT='(''     Rv/Rd-1  = '',e13.7)') RETV         WRITE(UNIT=6,FMT='(''     Rv/Rd-1  = '',e13.7)') RETV
111         !  
        !  
112         !*       6.    DEFINE THERMODYNAMIC CONSTANTS, LIQUID PHASE.         !*       6.    DEFINE THERMODYNAMIC CONSTANTS, LIQUID PHASE.
113         !  
114         RCW=RCPV         RCW=RCPV
115         WRITE(UNIT=6,FMT='('' *** Thermodynamic, liquid  ***'')')         WRITE(UNIT=6,FMT='('' *** Thermodynamic, liquid  ***'')')
116         WRITE(UNIT=6,FMT='(''         Cw   = '',E13.7)') RCW         WRITE(UNIT=6,FMT='(''         Cw   = '',E13.7)') RCW
117         !  
        !  
118         !*       7.    DEFINE THERMODYNAMIC CONSTANTS, SOLID PHASE.         !*       7.    DEFINE THERMODYNAMIC CONSTANTS, SOLID PHASE.
119         !  
120         RCS=RCPV         RCS=RCPV
121         WRITE(UNIT=6,FMT='('' *** thermodynamic, solid   ***'')')         WRITE(UNIT=6,FMT='('' *** thermodynamic, solid   ***'')')
122         WRITE(UNIT=6,FMT='(''         Cs   = '',E13.7)') RCS         WRITE(UNIT=6,FMT='(''         Cs   = '',E13.7)') RCS
123         !  
        !  
124         !*       8.    DEFINE THERMODYNAMIC CONSTANTS, TRANSITION OF PHASE.         !*       8.    DEFINE THERMODYNAMIC CONSTANTS, TRANSITION OF PHASE.
125         !  
126         RTT=273.16         RTT=273.16
127         RLVTT=2.5008E+6         RLVTT=2.5008E+6
128         RLSTT=2.8345E+6         RLSTT=2.8345E+6
# Line 142  contains Line 135  contains
135         WRITE(UNIT=6,FMT='(''        RLMlt  = '',E13.7)') RLMLT         WRITE(UNIT=6,FMT='(''        RLMlt  = '',E13.7)') RLMLT
136         WRITE(UNIT=6,FMT='('' Normal press. = '',E13.7)') RATM         WRITE(UNIT=6,FMT='('' Normal press. = '',E13.7)') RATM
137         WRITE(UNIT=6,FMT='('' Latent heat :  '')')         WRITE(UNIT=6,FMT='('' Latent heat :  '')')
138         !  
        !  
139         !*       9.    SATURATED VAPOUR PRESSURE.         !*       9.    SATURATED VAPOUR PRESSURE.
140         !  
141         RESTT=611.14         RESTT=611.14
142         RGAMW=(RCW-RCPV)/RV         RGAMW=(RCW-RCPV)/RV
143         RBETW=RLVTT/RV+RGAMW*RTT         RBETW=RLVTT/RV+RGAMW*RTT
# Line 156  contains Line 148  contains
148         RGAMD=RGAMS-RGAMW         RGAMD=RGAMS-RGAMW
149         RBETD=RBETS-RBETW         RBETD=RBETS-RBETW
150         RALPD=RALPS-RALPW         RALPD=RALPS-RALPW
151         !  
        !  
152         ! calculer les constantes pour les fonctions thermodynamiques         ! calculer les constantes pour les fonctions thermodynamiques
153         !  
154         RVTMP2=RCPV/RCPD-1.         RVTMP2=RCPV/RCPD-1.
155         RHOH2O=RATM/100.         RHOH2O=RATM/100.
156         R2ES=RESTT*RD/RV         R2ES=RESTT*RD/RV
# Line 169  contains Line 160  contains
160         R4IES=7.66         R4IES=7.66
161         R5LES=R3LES*(RTT-R4LES)         R5LES=R3LES*(RTT-R4LES)
162         R5IES=R3IES*(RTT-R4IES)         R5IES=R3IES*(RTT-R4IES)
   
163      ELSE      ELSE
164         PRINT *, 'suphec DEJA APPELE '         PRINT *, 'suphec DEJA APPELE '
165      ENDIF      ENDIF
Legend:
Removed from v.3  
changed lines
  Added in v.17
  ViewVC Help
Powered by ViewVC 1.1.21