anchor.diff:  !   else
   !      print *,'dcmip must be 40 something or 50 something'
   !      stop
  -!   endif 
  +!   endif
   !   zcoords = 0 ! p, not z, is specified.
   !   if (dry) then
   !     moist = 0
  
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:        nn = nn + 2
      enddo
   end subroutine copyspecout
  
  
   subroutine twodtooned(data2,data1)
     real(r_kind), intent(in) :: data2(nlons,nlats)
  