anchor.diff:         lall,         &
          test_case,    &
   !       eps_geo2prs,  &
  -!       nicamcore,    &
  +       nicamcore,    &
          chemtracer,   &
          DIAG_var      )
       use mod_adm, only: &
  
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:     integer,          intent(in)  :: lall
      character(len=*), intent(in)  :: test_case
  !    real(RP),         intent(in)  :: eps_geo2prs
      logical,          intent(in)  :: nicamcore
      logical,          intent(in)  :: chemtracer
      real(RP),         intent(out) :: DIAG_var(ijdim,kdim,lall,6+TRC_VMAX)
  
  