anchor.diff:      real(DP) :: rho(kdim),   q(kdim)   ! density and water vapor mixing ratio in ICO-grid field
       real(DP) :: thetav(kdim)           ! Virtual potential temperature in ICO-grid field
       real(DP) :: p = 0.0_RP             ! dummy variable
  +    real(DP) :: Mvap2 ! Ratio of molar mass of dry air/water based on NICAM CONSTANTs
  +    real(DP) :: RdovRv
   
       real(DP) :: z(kdim)
       real(RP) :: vx_local(kdim)
  
positive.code:     logical,          intent(in)  :: prs_rebuild
      real(RP),         intent(out) :: DIAG_var(ijdim,kdim,lall,6+TRC_VMAX)
  
      real(DP) :: lat, lon               ! latitude, longitude on Icosahedral grid
      real(DP) :: prs(kdim), tmp(kdim)   ! presssure and temperature in ICO-grid field
      real(DP) :: wix(kdim),   wiy(kdim) ! zonal/meridional wind components in ICO-grid field
  
negative.code:     logical,          intent(in)  :: prs_rebuild
      real(RP),         intent(out) :: DIAG_var(ijdim,kdim,lall,6+TRC_VMAX)
  
      real(DP), parameter :: Mvap = 0.61d0 ! Ratio of molar mass of dry air/water (imported from supercell_test.f90)
      real(DP) :: lat, lon               ! latitude, longitude on Icosahedral grid
      real(DP) :: prs(kdim), tmp(kdim)   ! presssure and temperature in ICO-grid field
      real(DP) :: wix(kdim),   wiy(kdim) ! zonal/meridional wind components in ICO-grid field
  