anchor.diff:  
       !--- calculate other parameters
       CNST_PI    = 4.D0 * atan( 1.D0 )
  +    CNST_D2R   = CNST_PI / 180.D0
   
       CNST_CV    = CNST_CP - CNST_RAIR
       CNST_GAMMA = CNST_CP / CNST_CV
  
positive.code:        ADM_proc_stop
      implicit none
  
      !--- Parameters to be controled
      real(8) :: earth_radius      !--- Earth radius
      real(8) :: earth_angvel      !--- Anguler velocity of the earth
      real(8) :: earth_gravity     !--- Gravitational accelaration
      real(8) :: gas_cnst          !--- Gas constant of dry air
      real(8) :: gas_cnst_vap      !--- Gas constant of water vapour
      real(8) :: specific_heat_pre !--- Specific heat of air( const pre )
      real(8) :: specific_heat_pre_vap !--- Specific heat of water vapour ( const pre )
      real(8) :: latent_heat_vap   !--- latent heat of vaporization LH0 ( 0 deg )
      real(8) :: latent_heat_sub   !--- latent heat of sublimation LHS0 ( 0 deg )
  
      namelist / CNSTPARAM / &
           earth_radius,     & !--- earth radius
           earth_angvel,     & !--- anguler velocity of the earth
           earth_gravity,    & !--- gravitational accelaration
           gas_cnst,         & !--- gas constant of dry air ( Rd )
           gas_cnst_vap,     & !--- gas constant of water vapour
           specific_heat_pre,& !--- specific heat of air ( Cp )
           specific_heat_pre_vap,&  !--- specific heat of water vapour
           latent_heat_vap,  &
           latent_heat_sub
  
      integer :: ierr
      !---------------------------------------------------------------------------
  
negative.code:     endif
      write(ADM_LOG_FID,CNSTPARAM)
  
      CNST_ERADIUS = earth_radius / small_planet_factor
      CNST_EOHM    = earth_angvel * small_planet_factor
      CNST_EGRAV   = earth_gravity
      CNST_RAIR    = gas_cnst
      CNST_RVAP    = gas_cnst_vap
  