MODULE constants
      REAL,PARAMETER  :: cp=1004 
      REAL,PARAMETER  :: kappa=287./1004.
     END MODULE