subroutine su_watercycle

      use watercommon_h
      implicit none
#include "comcstfi.h"


!==================================================================
!
!     Purpose
!     -------
!     Set up relevant constants and parameters for the water cycle
!
!     Authors
!     -------
!     Robin Wordsworth (2010)
!
!==================================================================

      epsi   = mH2O / mugaz
      RCPD   = cpp 

      !RV = 1000.*R/mH2O
      RV = 1000.*8.314/mH2O ! caution! R is R/mugaz already!

      RCPV   = 4. *RV	 ! could be more precise...
      !RCPV   = 5/2 *8.31*1000/mH2O ! BC modif here: doesn't appear to be correct

      RVTMP2 = RCPV/RCPD-1. ! not currently used...

      end subroutine su_watercycle