anchor.diff:      real(DP) :: thetav(kdim)           ! Virtual potential temperature in ICO-grid field
       real(DP) :: p = 0.0_RP             ! dummy variable
   
  -    real(DP) :: z_local
  +    real(DP) :: z(kdim)
       real(RP) :: vx_local(kdim)
       real(RP) :: vy_local(kdim)
       real(RP) :: vz_local(kdim)
       real(DP) :: ps, phis
   
       real(RP) :: cl, cl2
  +!    real(RP) :: uave, dz, f_cf
  +    real(RP) :: tv(kdim)
   
       integer, parameter :: deep    = 0 ! deep atmosphere (1 = yes or 0 = no)
       integer, parameter :: zcoords = 1 ! 1 if z is specified, 0 if p is specified
  
positive.code: 
            ! Todo : the mixing ratios are dry
            ! i.e. the ratio between the density of the species and the density of dry air.
            DIAG_var(n,:,l,6+NQW_MAX+NCHEM_STR) = cl
            DIAG_var(n,:,l,6+NQW_MAX+NCHEM_END) = cl2
         endif
  
      enddo
      enddo
  
  !    write (ADM_LOG_FID,*) " |            Vertical Coordinate used in JBW initialization              |"
  !    write (ADM_LOG_FID,*) " |------------------------------------------------------------------------|"
  !    do k = 1, kdim
  !       write (ADM_LOG_FID,'(3X,"(k=",I3,") HGT:",F8.2," [m]",2X,"PRS: ",F9.2," [Pa]",2X,"GH: ",F8.2," [m]",2X,"ETA: ",F9.5)') &
  !       k, z_local(k), prs(k), geo(k)/g, eta(k,1)
  !    enddo
  !    write (ADM_LOG_FID,*) " |------------------------------------------------------------------------|"
  
      return
    end subroutine jbw_moist_init
  
negative.code:        write(ADM_LOG_FID,*) '*** eps_geo2prs = ', eps_geo2prs
         write(ADM_LOG_FID,*) '*** nicamcore   = ', nicamcore
         write(ADM_LOG_FID,*) '*** chemtracer  = ', chemtracer
         call jbw_moist_init( ADM_gall, ADM_kall, ADM_lall, test_case, nicamcore, &
                              chemtracer, DIAG_var(:,:,:,:) )
  
      case ('Traceradvection')
  