Two methods for computational reconstruction of fluctuating flows in the ocean based on field measurements are suggested. Here full-3D CFD treatment with the spatial scale of O(100m) - O(1000m) is taken into account. In the first method, the fluctuating flow based on measurement data is generated at the inflow boundary and introduced into the computational domain. In the second method, small wavenumber components of fluctuation in the ocean are forced in the computational domain, and large wavenumber components are evolved freely in the computation. The results of each method are discussed with comparing with the measurements, and their availability is suggested.