e t, q, u, v, tr are defined
             shape must conform to that of xi
             units = meters

	 dtau_datmos  -- real, dimension(:,:)
	     shape must conform to first two dimensions of input fields
	     Set = 0.0 if the surface flux is treated explicitly. If
	     dependence of flux on lowest level of atmosphere is
	     treated implicitly, set = derivative of surface stress
	     with respect to wind component in lowest model layer.  
             Units = kg/(m**2 s)

        dt_t, dt_q -- real,