!# Include for all codes with physics and othes constants !# !# @note !# ![](http://brams.cptec.inpe.br/wp-content/uploads/2015/11/logo-brams.navigation.png "") !# !# **Brief**: Physics & other constants !# !# **Documentation**: !# !# **Author**: Rodrigues, L.F. **✉** !# !# **Date**: 2018Sep !# @endnote !# !# @changes !# !# + !# @endchanges !# @bug !# No active bugs reported now !# @endbug !# !# @todo !# ☐
!# @endtodo !# !# @warning !# Now is under CC-GPL License, please see !# © !# @endwarning !# !#--- ---------------------------------------------------------------------------------------- ! !integer kinds integer, parameter :: kind_ib = selected_int_kind(13) !# 8 byte integer integer, parameter :: kind_im = selected_int_kind(6) !# 4 byte integer integer, parameter :: kind_in = kind(1) !# native integer !real kinds integer, parameter :: kind_rb = selected_real_kind(12) !# 8 byte real integer, parameter :: kind_rm = selected_real_kind(6) !# 4 byte real integer, parameter :: kind_rn = kind(1.0) !# native real character(len=*), parameter :: modelVersion='- Rev. 5.4 -' !# Last version of model !# Phys & others constants real(kind=kind_rb), parameter :: c_pi = 3.1415926535897932384626433 real(kind=kind_rb), parameter :: c_rgas = 287. real(kind=kind_rb), parameter :: c_cp = 1004. real(kind=kind_rb), parameter :: c_cv = 717. real(kind=kind_rb), parameter :: c_rm = 461. real(kind=kind_rb), parameter :: c_p00 = 1.e5 real(kind=kind_rb), parameter :: c_t00 = 273.16 !Absolute temperature real(kind=kind_rb), parameter :: c_pi180 = c_pi / 180. real(kind=kind_rb), parameter :: c_i_pi180 = 1./c_pi180 real(kind=kind_rb), parameter :: c_pi4 = c_pi * 4. real(kind=kind_rb), parameter :: c_spcon = 111120. real(kind=kind_rb), parameter :: c_erad = 6367000. real(kind=kind_rb), parameter :: c_vonk = 0.40 real(kind=kind_rb), parameter :: c_tkmin = 5.e-4 !# Minimum TKE [J/kg] real(kind=kind_rb), parameter :: c_grav = 9.80665 !# Gravity acceleration [m/s] real(kind=kind_rb), parameter :: c_alvl = 2.50e6 real(kind=kind_rb), parameter :: c_alvi = 2.834e6 real(kind=kind_rb), parameter :: c_alli = 0.334e6 real(kind=kind_rb), parameter :: c_alvl2 = 6.25e12 real(kind=kind_rb), parameter :: c_alvi2 = 8.032e12 real(kind=kind_rb), parameter :: c_solar = 1.3533e3 real(kind=kind_rb), parameter :: c_stefan = 5.6696e-8 real(kind=kind_rb), parameter :: c_cww = 4218. real(kind=kind_rb), parameter :: c_c0 = 752.55 * 4.18684e4 real(kind=kind_rb), parameter :: c_viscos = .15e-4 real(kind=kind_rb), parameter :: c_rowt = 1.e3 real(kind=kind_rb), parameter :: c_dlat = 111120. real(kind=kind_rb), parameter :: c_omega = 7.292e-5 real(kind=kind_rb), parameter :: c_rocp = c_rgas / c_cp real(kind=kind_rb), parameter :: c_p00i = 1. / c_p00 real(kind=kind_rb), parameter :: c_cpor = c_cp / c_rgas real(kind=kind_rb), parameter :: c_rocv = c_rgas / c_cv real(kind=kind_rb), parameter :: c_cpi = 1. / c_cp real(kind=kind_rb), parameter :: c_cpi4 = 4. * c_cpi real(kind=kind_rb), parameter :: c_cp253i = c_cpi / 253. real(kind=kind_rb), parameter :: c_allii = 1. / c_alli real(kind=kind_rb), parameter :: c_aklv = c_alvl / c_cp real(kind=kind_rb), parameter :: c_akiv = c_alvi / c_cp real(kind=kind_rb), parameter :: c_gama = c_cp / c_cv real(kind=kind_rb), parameter :: c_gg = .5 * c_grav real(kind=kind_rb), parameter :: c_ep = c_rgas / c_rm real(kind=kind_rb), parameter :: c_p00k = 26.870941 !# = p00 ** rocp real(kind=kind_rb), parameter :: c_p00ki = 1. / c_p00k real(kind=kind_rb), parameter :: c_onethird = 1./3. !# 1/3 real(kind=kind_rb), parameter :: c_half = 1./2. !# 1/2 ! Lower bounds for turbulence-related variables ! real(kind=kind_rb), parameter :: c_sigwmin = 1.e-4 !# Minimum sigma-w [m/s] real(kind=kind_rb), parameter :: c_abslmomin = 1.e-4 !# Minimum abs value of Obukhov length [m] real(kind=kind_rb), parameter :: c_ltscalemax = 1.e5 !# Maximum Lagrangian timescale [s] real(kind=kind_rb), parameter :: c_abswltlmin = 1.e-4 !# Minimum abs value of Theta* [K m/s] real(kind=kind_rb), parameter :: c_lturbmin = 1.e-3 !# Minimum abs value of turb. lenght [m] real,parameter :: c_scday=86400.0 !# Seconds by day !Tiny numbers real(kind=kind_rb), parameter :: tinyReal = 1.17549435E-38 !# Tiny real number !real(kind=kind_rb), parameter :: tinyDouble = 2.2250738585072014E-308 !# Tiny double precision number !Errors type for dump functions integer, parameter :: c_noError=0 !# No Error, just write integer, parameter :: c_notice=1 !# Notice Error integer, parameter :: c_warning=2 !# warning Error integer, parameter :: c_fatal=3 !# Fatal error logical, parameter :: c_yes=.true. !# Yes is the .true. value logical, parameter :: c_no=.false. !# No is the false value character(len=*), parameter :: c_empty='' !# empty string !Colors strings for print and write commands character(len=5), parameter :: c_darkGrey =achar(27)//'[90m' character(len=5), parameter :: c_peach =achar(27)//'[91m' character(len=5), parameter :: c_lightGreen =achar(27)//'[92m' character(len=5), parameter :: c_lightYellow=achar(27)//'[93m' character(len=5), parameter :: c_lightBlue =achar(27)//'[94m' character(len=5), parameter :: c_pink =achar(27)//'[95m' character(len=5), parameter :: c_lightAqua =achar(27)//'[96m' character(len=5), parameter :: c_pearlWhite =achar(27)//'[97m' character(len=5), parameter :: c_black =achar(27)//'[30m' character(len=5), parameter :: c_red =achar(27)//'[31m' character(len=5), parameter :: c_green =achar(27)//'[32m' character(len=5), parameter :: c_yellow =achar(27)//'[33m' character(len=5), parameter :: c_blue =achar(27)//'[34m' character(len=5), parameter :: c_purple =achar(27)//'[35m' character(len=5), parameter :: c_aqua =achar(27)//'[36m' character(len=*), parameter :: c_blink =achar(27)//'[31;5;95;38;5;214m' character(len=4), parameter :: c_noColor =achar(27)//'[0m' integer, parameter :: c_tty=6 !# Default TTY (terminal) output !#ifdef splitlog ! integer, parameter :: logUnit=69 !#else integer, parameter :: logUnit=6 !#endif !# Number of file to log integer :: iErrNumber !# integer to use with dumps real(kind=kind_rb), parameter :: c_adjust=0.01 !# to convert from Pascal to mbar character(len=3), parameter, dimension(12) :: month_name=(/'jan','feb','mar' & , 'apr','may','jun' & , 'jul','aug','sep' & , 'oct','nov','dec'/)