#include <misc.h>#include <params.h>subroutine inti ()!----------------------------------------------------------------------- ! ! Purpose: ! Set constants and call initialization procedures for time independent! physics routines! ! Method: ! <Describe the algorithm(s) used in the routine.> ! <Also include any applicable external references.> ! ! Author: J. Rosinski! !-----------------------------------------------------------------------   use precision,          only: r8   use pmgrid,             only: plev, plevp          ! Needed for hypm passed to vd_inti   use chemistry,          only: chem_initialize   use gw_drag,            only: gw_inti   use vertical_diffusion, only: vd_inti   use moistconvection,    only: mfinti   use cldwat,             only: inimc   use zm_conv,            only: zm_convi   use shr_const_mod,      only: shr_const_zvir, shr_const_cpwv, shr_const_rwv   use physconst,          only: rair, cpair, cpwv, gravit, stebol, epsilo, tmelt, &                                 latvap, latice, rh2o, zvir, cpvir, rhoh2o, pstd,  &                                 karman, rhodair   implicit none#include <comctl.h>#include <comhyb.h>!!-----------------------------------------------------------------------!! Initialize physconst variables! In adiabatic case, set zvir and cpvir explicitly to zero instead of ! computing as (rh2o/rair - 1.) and (cpwv/cpair - 1.) respectively, in order ! to guarantee an identical zero.!   if (adiabatic .or. ideal_phys) then      rh2o  = rair      zvir  = 0.      cpwv  = cpair      cpvir = 0.   else      rh2o  = shr_const_rwv      zvir  = shr_const_zvir      cpwv  = shr_const_cpwv      cpvir = cpwv/cpair - 1.   end if!! Call time independent initialization routines for parameterizations.!   if (trace_gas) call chem_initialize   call gw_inti (cpair   ,cpwv    ,gravit  ,rair    ,hypi    )   call vd_inti (cpair   ,cpwv    ,gravit  ,rair    ,zvir   , &                 hypm    ,karman    )   call tsinti  (tmelt   ,latvap  ,rair    ,stebol  ,latice  )   call radini  (gravit  ,cpair   ,epsilo  ,stebol  ,pstd*10.0 )   call esinti  (epsilo  ,latvap  ,latice  ,rh2o    ,cpair  , &                 tmelt   )   call mfinti  (rair    ,cpair   ,gravit  ,latvap  ,rhoh2o  )   call zm_convi( tmelt, epsilo, latvap, cpair )   call cldinti ()!! initialization routine for prognostic cloud water!   call inimc (tmelt, rhodair/1000.0, gravit, rh2o )   returnend subroutine inti