anchor.diff:         write(ADM_LOG_FID,*) '*** test case   : ', trim(test_case)
          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, &
  +       call jbw_moist_init( ADM_gall, ADM_kall, ADM_lall, test_case, nicamcore, prs_rebuild, &
                               chemtracer, DIAG_var(:,:,:,:) )
   
       case ('Supercell')
   
          write(ADM_LOG_FID,*) '*** test case   : ', trim(test_case)
          write(ADM_LOG_FID,*) '*** nicamcore   = ', nicamcore
  -       call sc_init( ADM_gall, ADM_kall, ADM_lall, test_case, nicamcore, DIAG_var(:,:,:,:) )
  +       call sc_init( ADM_gall, ADM_kall, ADM_lall, test_case, nicamcore, prs_rebuild, DIAG_var(:,:,:,:) )
   
       case ('Tropical-Cyclone')
   
          write(ADM_LOG_FID,*) '*** nicamcore   = ', nicamcore
  -       call tc_init( ADM_gall, ADM_kall, ADM_lall, nicamcore, DIAG_var(:,:,:,:) )
  +       call tc_init( ADM_gall, ADM_kall, ADM_lall, nicamcore, prs_rebuild, DIAG_var(:,:,:,:) )
   
       case ('Traceradvection')
   
  
positive.code:     integer,          intent(in)  :: lall
      character(len=*), intent(in)  :: test_case
      logical,          intent(in)  :: nicamcore
      real(RP),         intent(out) :: DIAG_var(ijdim,kdim,lall,6+TRC_VMAX)
  
      real(DP) :: lat, lon               ! latitude, longitude on Icosahedral grid
  
negative.code:     integer,          intent(in)  :: kdim
      integer,          intent(in)  :: lall
      logical,          intent(in)  :: nicamcore
      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
  