anchor.diff:      logical,          intent(in)  :: prs_rebuild
       real(RP),         intent(out) :: DIAG_var(ijdim,kdim,lall,6+TRC_VMAX)
   
  +    real(DP), parameter :: Mvap = 0.608d0 ! Ratio of molar mass of dry air/water (imported from tropical_cyclone_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
  
positive.code:                                       rho(k),    &  ! [OUT] density (kg m^-3)
                                        q(k)       )  ! [OUT] water vapor mixing ratio (kg/kg)
         enddo
  
         call diag_pressure( kdim, z, rho, tmp, q, ps, prs, prs_rebuild, prs_dry )
         call conv_vxvyvz ( kdim, lat, lon, wix, wiy, vx_local, vy_local, vz_local )
  
negative.code:                                      q(k)        )  ! [OUT] water vapor mixing ratio (kg/kg)
         enddo
  
         ! Re-Evaluation of temperature from virtual temperature
         tmp(:) = tmp(:) * ( (1.d0+Mvap*q(:)) / (1.d0+Mvap2*q(:)) )
  
         call diag_pressure( kdim, z, rho, tmp, q, ps, prs, prs_rebuild, prs_dry )
         call conv_vxvyvz  ( kdim, lat, lon, wix, wiy, vx_local, vy_local, vz_local )
  
         do k = 1, kdim
            DIAG_var(n,k,l,1     ) = real(prs(k),kind=RP)
  