wavesp0 (i,j,bi,bj) = 0.
CALL READ_REC_XY_RL( WavePFile,wavesp0,intime0,
_EXCH_XY_RL( wavesp0 , myThid )
wavesp(i,j,bi,bj) = bWght*wavesp0(i,j,bi,bj)
& wavesp0, wavesp1,
_RL wavesp0 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)