anchor.diff:  ! subroutine dcmip_ics(testcase,dry)
   !   ! DCMIP 2012 test case initial conditions (test cases 4 and 5)
   !   ! (idealized baroclinic instability or trop cyclone test case initial conditions)
  -!   use dcmip_initial_conditions_test_4, only: test4_baroclinic_wave 
  -!   use dcmip_initial_conditions_test_5, only: test5_tropical_cyclone 
  +!   use dcmip_initial_conditions_test_4, only: test4_baroclinic_wave
  +!   use dcmip_initial_conditions_test_5, only: test5_tropical_cyclone
   !   integer, intent(in) :: testcase
   !   logical, intent(in) :: dry
   !   real(8) lon,lat,z,X,p,u,v,w,t,phi_s,ps,rho,q,q1,q2
  
positive.code:     enddo
           ! compute sigma coordinate quantities (bottom to top).
      do k=1,nlevs+1
          si(nlevs+2-k)= ak(k)/101300.0+bk(k) ! si are now sigmas
      enddo
      do k=1,nlevs
          sl(k) = 0.5*(si(k)+si(k+1))
  
negative.code: !          call grdtospec(tracerg(:,:,k,nt),tracerspec(:,k,nt))
  !       enddo
  !   enddo
  !   call grdtospec(phis/grav, topospec)
  !
  ! end subroutine dcmip_ics
  
  