ens_datmos, devap_datmos,    &
                    sens, evap, tau_u, tau_v, flux_tr,               &
                    dt_u, dt_v, dt_t, dt_q, dt_tr, dissipative_heat, &
		    kbot)
 
    input: 

         delt -- real
           time step , seconds
           if using leapfrog step, set equal to 2*dt = time(i+1) - time(i-1)

         u -- real, dimension(:,:,:)   zonal wind (m/sec)
         v -- real, dimension(:,:,:)   meridional wind (m/sec)
         T -- real, dimension(:,:,:)   temperature  (K)
    