!MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier !MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence !MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt !MNH_LIC for details. version 1. !----------------------------------------------------------------- !--------------- special set of characters for RCS information !----------------------------------------------------------------- ! $Source: /srv/cvsroot/MNH-VX-Y-Z/src/MNH/ini_cst.f90,v $ $Revision: 1.1.10.1.2.1.2.2.2.3 $ ! MASDEV4_7 init 2006/05/18 13:07:25 !----------------------------------------------------------------- ! ################### MODULE MODI_INI_CST ! ################### ! INTERFACE ! SUBROUTINE INI_CST END SUBROUTINE INI_CST ! END INTERFACE ! END MODULE MODI_INI_CST ! ! ! ! ################## SUBROUTINE INI_CST ! ################## ! !!**** *INI_CST * - routine to initialize the module MODD_CST !! !! PURPOSE !! ------- ! The purpose of this routine is to initialize the physical constants ! stored in module MODD_CST. ! ! !!** METHOD !! ------ !! The physical constants are set to their numerical values !! !! !! EXTERNAL !! -------- !! FMLOOK : to retrieve logical unit number associated to a file !! !! IMPLICIT ARGUMENTS !! ------------------ !! Module MODD_CST : contains physical constants !! !! REFERENCE !! --------- !! Book2 of the documentation (module MODD_CST, routine INI_CST) !! !! !! AUTHOR !! ------ !! V. Ducrocq * Meteo France * !! !! MODIFICATIONS !! ------------- !! Original 18/05/94 !! J. Stein 02/01/95 add the volumic mass of liquid water !! J.-P. Pinty 13/12/95 add the water vapor pressure over solid ice !! J. Stein 29/06/97 add XTH00 !! V. Masson 05/10/98 add XRHOLI !! C. Mari 31/10/00 add NDAYSEC !! V. Masson 01/03/03 add XCONDI !! J. Escobar 28/03/2014 for pb with emissivity/aerosol reset XMNH_TINY=1.0e-80 in real8 case !! !------------------------------------------------------------------------------- ! !* 0. DECLARATIONS ! ------------ ! USE MODD_CST ! IMPLICIT NONE ! !------------------------------------------------------------------------------- ! !* 1. FUNDAMENTAL CONSTANTS ! --------------------- ! XPI = 2.*ASIN(1.) XKARMAN = 0.38!0.4 XLIGHTSPEED = 299792458. XPLANCK = 6.6260755E-34 XBOLTZ = 1.380658E-23 XAVOGADRO = 6.0221367E+23 ! !------------------------------------------------------------------------------- ! !* 2. ASTRONOMICAL CONSTANTS ! ---------------------- ! XDAY = 86400. XSIYEA = 365.25*XDAY*2.*XPI/ 6.283076 XSIDAY = XDAY/(1.+XDAY/XSIYEA) XOMEGA = 2.*XPI/XSIDAY XOMEGA = 0. !===>FA NDAYSEC = 24*3600 ! Number of seconds in a day ! !-------------------------------------------------------------------------------! ! ! !* 3. TERRESTRIAL GEOIDE CONSTANTS ! ---------------------------- ! XRADIUS = 6371229. XG = 9.80665!*1e-6 !XG =0.0058 !===>FA ! !------------------------------------------------------------------------------- ! !* 4. REFERENCE PRESSURE ! ------------------- ! !XP00 = 0.952327E+05 !XTH00 = 304.070 XP00 = 1.E5 XTH00 = 300. !------------------------------------------------------------------------------- ! !* 5. RADIATION CONSTANTS ! ------------------- ! !JUAN OVERFLOW XSTEFAN = 2.* XPI**5 * XBOLTZ**4 / (15.* XLIGHTSPEED**2 * XPLANCK**3) XSTEFAN = ( 2.* XPI**5 / 15. ) * ( (XBOLTZ / XPLANCK) * XBOLTZ ) * (XBOLTZ/(XLIGHTSPEED*XPLANCK))**2 XI0 = 1370. ! !------------------------------------------------------------------------------- ! !* 6. THERMODYNAMIC CONSTANTS ! ----------------------- ! XMD = 28.9644E-3 XMV = 18.0153E-3 XRD = XAVOGADRO * XBOLTZ / XMD !XRD = 287. !===>FA XRV = XAVOGADRO * XBOLTZ / XMV XCPD = 7.* XRD /2. !XCPD = 1004. !===>FA XCPV = 4.* XRV !XCPV = 717. !===>FA XRHOLW = 1000. XRHOLI = 900. XCONDI = 2.22 XCL = 4.218E+3 XCI = 2.106E+3 XTT = 273.16 XLVTT = 2.5008E+6 XLSTT = 2.8345E+6 XLMTT = XLSTT - XLVTT XESTT = 611.14 XGAMW = (XCL - XCPV) / XRV XBETAW = (XLVTT/XRV) + (XGAMW * XTT) XALPW = LOG(XESTT) + (XBETAW /XTT) + (XGAMW *LOG(XTT)) XGAMI = (XCI - XCPV) / XRV XBETAI = (XLSTT/XRV) + (XGAMI * XTT) XALPI = LOG(XESTT) + (XBETAI /XTT) + (XGAMI *LOG(XTT)) ! ! Some machine precision value depending of real4/8 use ! XMNH_EPSILON = EPSILON (XMNH_EPSILON ) XMNH_HUGE = HUGE (XMNH_HUGE ) #ifdef MNH_MPI_DOUBLE_PRECISION XMNH_TINY = 1.0e-80 XEPS_DT = 1.0e-5 XRES_FLAT_CART = 1.0e-12 XRES_OTHER = 1.0e-9 XRES_PREP = 1.0e-8 #else XMNH_TINY = TINY (XMNH_TINY ) XEPS_DT = 1.0e-4 XRES_FLAT_CART = 1.0e-12 XRES_OTHER = 1.0e-7 XRES_PREP = 1.0e-4 #endif XMNH_TINY_12 = SQRT (XMNH_TINY ) ! !------------------------------------------------------------------------------- ! END SUBROUTINE INI_CST