   !# 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**: <http://brams.cptec.inpe.br/documentation/>
   !#
   !# **Author**: Rodrigues, L.F. **&#9993;**<mailto:luiz.rodrigues@inpe.br>
   !#
   !# **Date**: 2018Sep
   !# @endnote
   !#
   !# @changes
   !#
   !# +
   !# @endchanges
   !# @bug
   !# No active bugs reported now
   !# @endbug
   !#
   !# @todo
   !#  &#9744; <br/>
   !# @endtodo
   !#
   !# @warning
   !# Now is under CC-GPL License, please see
   !# &copy; <https://creativecommons.org/licenses/GPL/2.0/legalcode.pt>
   !# @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

   !# 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_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

   !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
