! ! $Header$ ! SUBROUTINE suphel !L. Fita, LMD. September 2013 IMPLICIT NONE !C #include "YOMCST.h" #include "YOETHF.h" #include "../share/module_model_constants.F.h" !IM cf. JLD LOGICAL firstcall SAVE firstcall !$OMP THREADPRIVATE(firstcall) DATA firstcall /.TRUE./ IF (firstcall) THEN PRINT*, 'suphel initialise les constantes du GCM' firstcall = .FALSE. ELSE PRINT*, 'suphel DEJA APPELE ' RETURN ENDIF !C ----------------------------------------------------------------- !C !C* 1. DEFINE FUNDAMENTAL CONSTANTS. !C ----------------------------- !C WRITE(UNIT=6,FMT='(''0*** Constants of the ICM ***'')') RPI=2.*ASIN(1.) RCLUM=299792458. RHPLA=6.6260755E-34 RKBOL=1.380658E-23 RNAVO=6.0221367E+23 WRITE(UNIT=6,FMT='('' *** Fundamental constants ***'')') WRITE(UNIT=6,FMT='('' PI = '',E13.7,'' -'')')RPI WRITE(UNIT=6,FMT='('' c = '',E13.7,''m s-1'')') & & RCLUM WRITE(UNIT=6,FMT='('' h = '',E13.7,''J s'')') & & RHPLA WRITE(UNIT=6,FMT='('' K = '',E13.7,''J K-1'')') & & RKBOL WRITE(UNIT=6,FMT='('' N = '',E13.7,''mol-1'')') & & RNAVO !C !C ---------------------------------------------------------------- !C !C* 2. DEFINE ASTRONOMICAL CONSTANTS. !C ------------------------------ !C RDAY=86400. REA=149597870000. REPSM=0.409093 !C RSIYEA=365.25*RDAY*2.*RPI/6.283076 RSIDAY=RDAY/(1.+RDAY/RSIYEA) ROMEGA=2.*RPI/RSIDAY !c !c exp1 R_ecc = 0.05 !c exp1 R_peri = 102.04 !c exp1 R_incl = 22.5 !c exp1 print*, 'Parametres orbitaux modifies' !c ref R_ecc = 0.016724 !c ref R_peri = 102.04 !c ref R_incl = 23.5 !c !IM 161002 : pour avoir les ctes AMIP II !IM 161002 R_ecc = 0.016724 !IM 161002 R_peri = 102.04 !IM 161002 R_incl = 23.5 !IM on mets R_ecc, R_peri, R_incl dans conf_phys.F90 !c R_ecc = 0.016715 !c R_peri = 102.7 !c R_incl = 23.441 !c WRITE(UNIT=6,FMT='('' *** Astronomical constants ***'')') WRITE(UNIT=6,FMT='('' day = '',E13.7,'' s'')')RDAY WRITE(UNIT=6,FMT='('' half g. axis = '',E13.7,'' m'')')REA WRITE(UNIT=6,FMT='('' mean anomaly = '',E13.7,'' -'')')REPSM WRITE(UNIT=6,FMT='('' sideral year = '',E13.7,'' s'')')RSIYEA WRITE(UNIT=6,FMT='('' sideral day = '',E13.7,'' s'')')RSIDAY WRITE(UNIT=6,FMT='('' omega = '',E13.7,'' s-1'')') & & ROMEGA !c write(unit=6,fmt='('' excentricite = '',e13.7,''-'')')R_ecc !c write(unit=6,fmt='('' equinoxe = '',e13.7,''-'')')R_peri !c write(unit=6,fmt='('' inclinaison = '',e13.7,''-'')')R_incl !C !C ------------------------------------------------------------------ !C !C* 3. DEFINE GEOIDE. !C -------------- !C !L. Fita, LMD. September 2013 ! RG=9.80665 ! RA=6371229. RG=g RA=1.d0/reradius R1SA=SNGL(1.D0/DBLE(RA)) WRITE(UNIT=6,FMT='('' *** Geoide ***'')') WRITE(UNIT=6,FMT='('' Gravity = '',E13.7,'' m s-2'')') & & RG WRITE(UNIT=6,FMT='('' Earth radius = '',E13.7,'' m'')')RA WRITE(UNIT=6,FMT='('' Inverse E.R. = '',E13.7,'' m'')')R1SA !C !C ----------------------------------------------------------------- !C !C* 4. DEFINE RADIATION CONSTANTS. !C --------------------------- !C !c z.x.li RSIGMA=2. * RPI**5 * RKBOL**4 /(15.* RCLUM**2 * RHPLA**3) !L. Fita, LMD. September 2013 ! rsigma = 2.*rpi**5 * (rkbol/rhpla)**3 * rkbol/rclum/rclum/15. rsigma = STBOLT !IM init. dans conf_phys.F90 RI0=1365. WRITE(UNIT=6,FMT='('' *** Radiation ***'')') WRITE(UNIT=6,FMT='('' Stefan-Bol. = '',E13.7,'' W m-2 K-4'' & & )') RSIGMA !IM init. dans conf_phys.F90 WRITE(UNIT=6,FMT='('' Solar const. = '',E13.7,'' W m-2'')') !IM init. dans conf_phys.F90 S RI0 !C !C ----------------------------------------------------------------- !C !C* 5. DEFINE THERMODYNAMIC CONSTANTS, GAS PHASE. !C ------------------------------------------ !C R=RNAVO*RKBOL !L. Fita, LMD. September 2013 ! RMD=28.9644 RMD = mwdry RMO3=47.9942 RMV=18.0153 !L. Fita, LMD. September 2013 ! RD=1000.*R/RMD ! RV=1000.*R/RMV ! RCPD=3.5*RD RD = r_d RV = r_v RCPD = cp RCVD=RCPD-RD RCPV=4. *RV RCVV=RCPV-RV RKAPPA=RD/RCPD RETV=RV/RD-1. WRITE(UNIT=6,FMT='('' *** Thermodynamic, gas ***'')') WRITE(UNIT=6,FMT='('' Perfect gas = '',e13.7)') R WRITE(UNIT=6,FMT='('' Dry air mass = '',e13.7)') RMD WRITE(UNIT=6,FMT='('' Ozone mass = '',e13.7)') RMO3 WRITE(UNIT=6,FMT='('' Vapour mass = '',e13.7)') RMV WRITE(UNIT=6,FMT='('' Dry air cst. = '',e13.7)') RD WRITE(UNIT=6,FMT='('' Vapour cst. = '',e13.7)') RV WRITE(UNIT=6,FMT='('' Cpd = '',e13.7)') RCPD WRITE(UNIT=6,FMT='('' Cvd = '',e13.7)') RCVD WRITE(UNIT=6,FMT='('' Cpv = '',e13.7)') RCPV WRITE(UNIT=6,FMT='('' Cvv = '',e13.7)') RCVV WRITE(UNIT=6,FMT='('' Rd/Cpd = '',e13.7)') RKAPPA WRITE(UNIT=6,FMT='('' Rv/Rd-1 = '',e13.7)') RETV !C !C ---------------------------------------------------------------- !C !C* 6. DEFINE THERMODYNAMIC CONSTANTS, LIQUID PHASE. !C --------------------------------------------- !C RCW=RCPV WRITE(UNIT=6,FMT='('' *** Thermodynamic, liquid ***'')') WRITE(UNIT=6,FMT='('' Cw = '',E13.7)') RCW !C !C ---------------------------------------------------------------- !C !C* 7. DEFINE THERMODYNAMIC CONSTANTS, SOLID PHASE. !C -------------------------------------------- !C RCS=RCPV WRITE(UNIT=6,FMT='('' *** thermodynamic, solid ***'')') WRITE(UNIT=6,FMT='('' Cs = '',E13.7)') RCS !C !C ---------------------------------------------------------------- !C !C* 8. DEFINE THERMODYNAMIC CONSTANTS, TRANSITION OF PHASE. !C ---------------------------------------------------- !C !L. Fita, LMD. September 2013 ! RTT=273.16 RTT = SVPT0 RLVTT=2.5008E+6 RLSTT=2.8345E+6 RLMLT=RLSTT-RLVTT RATM=100000. WRITE(UNIT=6,FMT='('' *** Thermodynamic, trans. ***'')') WRITE(UNIT=6,FMT='('' Fusion point = '',E13.7)') RTT WRITE(UNIT=6,FMT='('' RLvTt = '',E13.7)') RLVTT WRITE(UNIT=6,FMT='('' RLsTt = '',E13.7)') RLSTT WRITE(UNIT=6,FMT='('' RLMlt = '',E13.7)') RLMLT WRITE(UNIT=6,FMT='('' Normal press. = '',E13.7)') RATM WRITE(UNIT=6,FMT='('' Latent heat : '')') !C !C ---------------------------------------------------------------- !C !C* 9. SATURATED VAPOUR PRESSURE. !C -------------------------- !C !L. Fita, LMD. September 2013 ! RESTT=611.14 RESTT = SVP1*1000. RGAMW=(RCW-RCPV)/RV RBETW=RLVTT/RV+RGAMW*RTT RALPW=LOG(RESTT)+RBETW/RTT+RGAMW*LOG(RTT) RGAMS=(RCS-RCPV)/RV RBETS=RLSTT/RV+RGAMS*RTT RALPS=LOG(RESTT)+RBETS/RTT+RGAMS*LOG(RTT) RGAMD=RGAMS-RGAMW RBETD=RBETS-RBETW RALPD=RALPS-RALPW !C !C ------------------------------------------------------------------ !c !c calculer les constantes pour les fonctions thermodynamiques !c RVTMP2=RCPV/RCPD-1. RHOH2O=RATM/100. R2ES=RESTT*RD/RV R3LES=17.269 R3IES=21.875 R4LES=35.86 R4IES=7.66 R5LES=R3LES*(RTT-R4LES) R5IES=R3IES*(RTT-R4IES) !C RETURN END