CALL SMOOTH_DIFF2D(smoothTmpFld,maskc,nbt_in,mythid)

           if (maskc(i,j,1,bi,bj).NE.0) then

      call SMOOTH_DIFF2D(smoothTmpFld,maskc,nbt_in,mythid)

     & xx_theta_dummy, boundsVec, theta, maskc, paramSmooth,

     & xx_salt_dummy, boundsVec, salt, maskc, paramSmooth, 

     & xx_etan_dummy, boundsVec, etan, maskc, paramSmooth, 

     & lambdaThetaClimRelax, maskc, paramSmooth, mythid )

     & lambdaSaltClimRelax, maskc, paramSmooth, mythid )

     & xx_bottomdrag_dummy, boundsVec, bottomdragfld, maskc, 

     & xx_diffkr_dummy, boundsVec, diffkr, maskc, paramSmooth, 

     & xx_kapgm_dummy, boundsVec, kapgm, maskc, paramSmooth, 

     & xx_kapredi_dummy, boundsVec, kapredi, maskc, paramSmooth, 

     &              ( fac*maskc(i,j,1,bi,bj)*bpbar(i,j,bi,bj) 

     &                      ( fac*maskc(i,j,1,bi,bj)*bpbar(i,j,bi,bj) 

     &                      *bpmask(i,j,bi,bj)*maskc(i,j,1,bi,bj)

      call SMOOTH_BASIC2D(mean_slaobs,maskc,150000. _d 0,1000,mythid)

               if ( maskc(i,j,1,bi,bj) .NE. 0. )

      call SMOOTH_BASIC2D(mean_psMssh_all,maskc,

     &              ( maskc(i,j,1,bi,bj) .NE. 0. ) ) then

      call SMOOTH_BASIC2D(anom_psMslaobs,maskc,300000. _d 0,2000,mythid)

      call SMOOTH_BASIC2D(anom_slaobs,maskc,300000. _d 0,2000,mythid)

     &    + junk*junk*wsshv4(i,j,5,bi,bj)*maskc(i,j,1,bi,bj)

     &         (maskc(i,j,1,bi,bj) .ne. 0.) )

              wvwind(i,j,bi,bj) = wvwind(i,j,bi,bj)*maskc(i,j,k,bi,bj)

              wtauu(i,j,bi,bj) = wtauu(i,j,bi,bj)*maskc(i,j,k,bi,bj)

              wtauu2(i,j,bi,bj) = wtau0*maskc(i,j,k,bi,bj)*frame(i,j)

              wuwind(i,j,bi,bj) = wuwind(i,j,bi,bj)*maskc(i,j,k,bi,bj)

              wtauv(i,j,bi,bj) = wtauv(i,j,bi,bj)*maskc(i,j,k,bi,bj)

              wtauv2(i,j,bi,bj) = wtau0*maskc(i,j,k,bi,bj)*frame(i,j)

     &     maskc, climsst, xx_sst0, xx_sst1, xx_sst_dummy,

     &     maskc, atemp, xx_atemp0, xx_atemp1, xx_atemp_dummy,

     &     maskc, aqh, xx_aqh0, xx_aqh1, xx_aqh_dummy,

     &     maskc, precip, xx_precip0, xx_precip1, xx_precip_dummy,

     &     maskc, swflux, xx_swflux0, xx_swflux1, xx_swflux_dummy,

     &     maskc, swdown, xx_swdown0, xx_swdown1, xx_swdown_dummy,

     &     maskc, lwflux, xx_lwflux0, xx_lwflux1, xx_lwflux_dummy,

     &     maskc, lwdown, xx_lwdown0, xx_lwdown1, xx_lwdown_dummy,

     &     maskc, evap, xx_evap0, xx_evap1, xx_evap_dummy,

     &     maskc, snowprecip, xx_snowprecip0, xx_snowprecip1, 

     &     maskc, apressure, xx_apressure0, xx_apressure1, 

     &     maskc, uwind, xx_uwind0, xx_uwind1, xx_uwind_dummy,

     &     maskc, vwind, xx_vwind0, xx_vwind1, xx_vwind_dummy,

     &     maskc, tmpUE, xx_uwind0, xx_uwind1, xx_uwind_dummy,

     &     maskc, tmpVN, xx_vwind0, xx_vwind1, xx_vwind_dummy,

     &     maskc, tmpUE, xx_tauu0, xx_tauu1, xx_tauu_dummy,

     &     maskc, tmpVN, xx_tauv0, xx_tauv1, xx_tauv_dummy,

     &     maskc, hflux, xx_hflux0, xx_hflux1, xx_hflux_dummy,

     &     maskc, sflux, xx_sflux0, xx_sflux1, xx_sflux_dummy,