A numerical cross-shore model to simulate two features of wave-mud interaction on fine-grained shore profiles, i.e. fluidization of mud layer and wave height attenuation, is presented. The analytical model of Foda et al. (1993) is used to determine the fluidized depth of mud layer. The wave height transformation is computed from the energy flux conservation law combining the effects of mud bed, shoaling and wave breaking. The rheological constitutive equations of visco-elastic (VE) material are used for numerical simulation of fluid mud behavior. A preliminary application of the model shows its capability to predict wave height transformation along soft muddy profiles.