Symbol Description Uses
maskC cell Center land mask 777



File Line number Procedure Code
./pkg/ecco/ecco_cost_weights.F 1640 ECCO_COST_WEIGHTS
                 wetan(i,j,bi,bj)=frame(i,j)*maskC(i,j,1,bi,bj)/
./pkg/ecco/ecco_cost_weights.F 2031 ECCO_COST_WEIGHTS
      CALL ACTIVE_WRITE_GEN_RS( 'maskCtrlC', maskC, 'XY', Nr,
./pkg/ecco/ecco_cost_weights.F 392 ECCO_COST_WEIGHTS
                 wsaltLev(i,j,k,bi,bj)=frame(i,j)*maskC(i,j,k,bi,bj)/
./pkg/ecco/ecco_cost_weights.F 407 ECCO_COST_WEIGHTS
     $             wsalt(k,bi,bj)*frame(i,j)*maskC(i,j,k,bi,bj)
./pkg/ecco/ecco_cost_weights.F 435 ECCO_COST_WEIGHTS
                 wthetaLev(i,j,k,bi,bj)=frame(i,j)*maskC(i,j,k,bi,bj)/
./pkg/ecco/ecco_cost_weights.F 450 ECCO_COST_WEIGHTS
     $             wtheta(k,bi,bj)*frame(i,j)*maskC(i,j,k,bi,bj)
./pkg/ecco/ecco_cost_weights.F 482 ECCO_COST_WEIGHTS
                 wsalt2(i,j,k,bi,bj)=frame(i,j)*maskC(i,j,k,bi,bj)/
./pkg/ecco/ecco_cost_weights.F 497 ECCO_COST_WEIGHTS
     $             wsalt(k,bi,bj)*frame(i,j)*maskC(i,j,k,bi,bj)
./pkg/ecco/ecco_cost_weights.F 529 ECCO_COST_WEIGHTS
                 wtheta2(i,j,k,bi,bj)= frame(i,j)*maskC(i,j,k,bi,bj)/
./pkg/ecco/ecco_cost_weights.F 547 ECCO_COST_WEIGHTS
     $              wtheta(k,bi,bj)*frame(i,j)*maskC(i,j,k,bi,bj)
./pkg/ecco/ecco_cost_weights.F 581 ECCO_COST_WEIGHTS
     &               wsst(i,j,bi,bj)= frame(i,j)*maskC(i,j,k,bi,bj)/
./pkg/ecco/ecco_cost_weights.F 591 ECCO_COST_WEIGHTS
     &               wtheta(k,bi,bj)*frame(i,j)*maskC(i,j,k,bi,bj)
./pkg/ecco/ecco_cost_weights.F 601 ECCO_COST_WEIGHTS
     &               wsss(i,j,bi,bj)= frame(i,j)*maskC(i,j,k,bi,bj)/
./pkg/ecco/ecco_cost_weights.F 612 ECCO_COST_WEIGHTS
     &               wsalt(k,bi,bj)*frame(i,j)*maskC(i,j,k,bi,bj)
./pkg/ecco/ecco_cost_weights.F 688 ECCO_COST_WEIGHTS
               if (maskC(i,j,k,bi,bj) .eq. 0.) then
./pkg/ecco/ecco_cost_weights.F 725 ECCO_COST_WEIGHTS
             if (maskC(i,j,k,bi,bj) .eq. 0.) then
./pkg/ecco/ecco_cost_weights.F 814 ECCO_COST_WEIGHTS
        wsshv4tmp(i,j,bi,bj)= frame(i,j)*maskC(i,j,k,bi,bj)/
./pkg/ecco/ecco_cost_weights.F 842 ECCO_COST_WEIGHTS
     &               wbp(i,j,bi,bj)= frame(i,j)*maskC(i,j,k,bi,bj)/
./pkg/ecco/ecco_cost_weights.F 863 ECCO_COST_WEIGHTS
     &               wies(i,j,bi,bj)= frame(i,j)*maskC(i,j,k,bi,bj)/
./pkg/ecco/ecco_phys.F 286 ECCO_PHYS
          tmpmsk=maskC(i,j,1,bi,bj)*gencost_mskCsurf(i,j,bi,bj,kgen)
./model/inc/GRID.h 424 NO PROCEDURE
     &  maskC, maskW, maskS,
./model/inc/GRID.h 476 NO PROCEDURE
      _RS maskC          (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
./model/inc/HFACC_MACROS.h 36 NO PROCEDURE
# define _hFacC(i,j,k,bi,bj) hFacC(i,j,k,bi,bj)*maskC(i,j,k,bi,bj)
./pkg/exf/exf_bulkformulae.F 294 EXF_BULKFORMULAE
     &              *  maskC(i,j,ksrfp1,bi,bj)
./pkg/exf/exf_check_range.F 123 EXF_CHECK_RANGE
     &         maskC(i,j,1,bi,bj) .NE. 0. ) THEN
./pkg/exf/exf_check_range.F 134 EXF_CHECK_RANGE
     &         maskC(i,j,1,bi,bj) .NE. 0. ) THEN
./pkg/exf/exf_check_range.F 144 EXF_CHECK_RANGE
     &         maskC(i,j,1,bi,bj) .NE. 0. ) THEN
./pkg/exf/exf_check_range.F 154 EXF_CHECK_RANGE
     &         maskC(i,j,1,bi,bj) .NE. 0. ) THEN
./pkg/exf/exf_check_range.F 167 EXF_CHECK_RANGE
     &         maskC(i,j,1,bi,bj) .NE. 0. ) THEN
./pkg/exf/exf_check_range.F 179 EXF_CHECK_RANGE
     &         maskC(i,j,1,bi,bj) .NE. 0. ) THEN
./pkg/exf/exf_check_range.F 197 EXF_CHECK_RANGE
     &         maskC(i,j,1,bi,bj) .NE. 0. ) THEN
./pkg/exf/exf_check_range.F 55 EXF_CHECK_RANGE
     &         maskC(i,j,1,bi,bj) .NE. 0. ) THEN
./pkg/exf/exf_check_range.F 64 EXF_CHECK_RANGE
     &         maskC(i,j,1,bi,bj) .NE. 0. ) THEN
./pkg/exf/exf_filter_rl.F 64 EXF_FILTER_RL
               IF ( maskC(i,j,ks,bi,bj) .EQ. 0. ) THEN
./pkg/exf/exf_getclim.F 192 EXF_GETCLIM
              IF ( maskC(i,j,1,bi,bj) .NE. 0. .AND.
./pkg/exf/exf_getforcing.F 231 EXF_GETFORCING
              hflux(i,j,bi,bj) = hflux(i,j,bi,bj)*maskC(i,j,k,bi,bj)
./pkg/exf/exf_getforcing.F 232 EXF_GETFORCING
              sflux(i,j,bi,bj) = sflux(i,j,bi,bj)*maskC(i,j,k,bi,bj)
./pkg/exf/exf_init_varia.F 477 EXF_INIT_VARIA
           IF ( maskC(i,j,1,bi,bj) .NE. 0. .AND.
./pkg/exf/exf_mapfields.F 202 EXF_MAPFIELDS
              Qnet(i,j,bi,bj) = Qnet(i,j,bi,bj)*maskC(i,j,ks,bi,bj)
./pkg/exf/exf_radiation.F 78 EXF_RADIATION
     &             *  maskC(i,j,kl,bi,bj)
./pkg/fizhi/fizhi_tendency_apply.F 226 FIZHI_TENDENCY_APPLY_T
     &    + ( maskC(i,j,kLev,bi,bj)*gthphy(i,j,kLev,bi,bj)
./pkg/fizhi/fizhi_tendency_apply.F 244 FIZHI_TENDENCY_APPLY_T
     &      ( maskC(i,j,kLev,bi,bj)*gthphy(i,j,kLev,bi,bj)
./pkg/fizhi/fizhi_tendency_apply.F 325 FIZHI_TENDENCY_APPLY_S
     &    + maskC(i,j,kLev,bi,bj)*gsphy(i,j,kLev,bi,bj)
./pkg/fizhi/fizhi_tendency_apply.F 332 FIZHI_TENDENCY_APPLY_S
        tmpdiag(i,j) = ( maskC(i,j,kLev,bi,bj) * gsphy(i,j,kLev,bi,bj) )
./pkg/fizhi/fizhi_wrapper.F 331 FIZHI_WRAPPER
       call ATOC(myThid,guphy,gvphy,maskC,idim1,idim2,jdim1,jdim2,Nr,
./pkg/fizhi/step_fizhi_corr.F 183 STEP_FIZHI_CORR
       call ATOC(myThid,udyntemp,vdyntemp,maskC,im1,im2,jm1,jm2,Nr,
./pkg/flt/flt_runga2.F 223 FLT_RUNGA2
     &              maskC(ic,jc,1,bi,bj).EQ.1. ) THEN
./pkg/flt/flt_runga4.F 330 FLT_RUNGA4
     &              maskC(ic,jc,1,bi,bj).EQ.1. ) THEN
./pkg/frazil/frazil_calc_rhs.F 69 FRAZIL_CALC_RHS
           IF ( maskC(i,j,k-1,bi,bj) .NE. 0. _d 0 .AND.
./pkg/frazil/frazil_calc_rhs.F 70 FRAZIL_CALC_RHS
     &          maskC(i,j,k,  bi,bj) .NE. 0. _d 0 ) THEN
./pkg/generic_advdiff/gad_advection.F 303 GAD_ADVECTION
     &                 + ( oneRS - maskC(i,j,k,bi,bj) )
./pkg/generic_advdiff/gad_advection.F 856 GAD_ADVECTION
     &                      *maskC(i,j,k-1,bi,bj)
./pkg/generic_advdiff/gad_advection.F 885 GAD_ADVECTION
     &                 *maskC(i,j,k-1,bi,bj)
./pkg/generic_advdiff/gad_advection.F 941 GAD_ADVECTION
     &                 *maskC(i,j,k-1,bi,bj)
./pkg/generic_advdiff/gad_biharm_r.F 72 GAD_BIHARM_R
     &                *tmpFac*maskC(i,j,kl,bi,bj)*maskC(i,j,km,bi,bj)
./pkg/generic_advdiff/gad_c2_adv_r.F 63 GAD_C2_ADV_R
         wT(i,j) = maskC(i,j,km1,bi,bj)*
./model/inc/RECIP_HFACC_MACROS.h 36 NO PROCEDURE
# define _recip_hFacC(i,j,k,bi,bj) recip_hFacC(i,j,k,bi,bj)*maskC(i,j,k,bi,bj)
./pkg/generic_advdiff/gad_c4_adv_r.F 76 GAD_C4_ADV_R
         maskBound = maskPM*maskC(i,j,km2,bi,bj)*maskC(i,j,kp1,bi,bj)
./pkg/generic_advdiff/gad_c4_adv_r.F 77 GAD_C4_ADV_R
         Rjp = (tracer(i,j,kp1)-tracer(i,j, k ))*maskC(i,j,kp1,bi,bj)
./pkg/generic_advdiff/gad_c4_adv_r.F 79 GAD_C4_ADV_R
         Rjm = (tracer(i,j,km1)-tracer(i,j,km2))*maskC(i,j,km1,bi,bj)
./pkg/generic_advdiff/gad_c4_adv_r.F 83 GAD_C4_ADV_R
         wT(i,j) = maskC(i,j,km1,bi,bj)*(
./pkg/generic_advdiff/gad_dst2u1_adv_r.F 90 GAD_DST2U1_ADV_R
         wT(i,j) = maskC(i,j,km1,bi,bj)*(
./pkg/generic_advdiff/gad_dst3_adv_r.F 80 GAD_DST3_ADV_R
     &         *maskC(i,j,kp1,bi,bj)
./pkg/generic_advdiff/gad_dst3_adv_r.F 82 GAD_DST3_ADV_R
     &         *maskC(i,j,k,bi,bj)*maskC(i,j,km1,bi,bj)
./pkg/generic_advdiff/gad_dst3_adv_r.F 84 GAD_DST3_ADV_R
     &         *maskC(i,j,km1,bi,bj)
./pkg/generic_advdiff/gad_dst3fl_adv_r.F 74 GAD_DST3FL_ADV_R
     &         *maskC(i,j,kp1,bi,bj)
./pkg/generic_advdiff/gad_dst3fl_adv_r.F 76 GAD_DST3FL_ADV_R
     &         *maskC(i,j,k,bi,bj)*maskC(i,j,km1,bi,bj)
./pkg/generic_advdiff/gad_dst3fl_adv_r.F 78 GAD_DST3FL_ADV_R
     &         *maskC(i,j,km1,bi,bj)
./pkg/generic_advdiff/gad_dst3fl_impl_r.F 112 GAD_DST3FL_IMPL_R
             Rjh = (tFld(i,j,k-1)-tFld(i,j,km2))*maskC(i,j,km2,bi,bj)
./pkg/generic_advdiff/gad_dst3fl_impl_r.F 114 GAD_DST3FL_IMPL_R
             Rjh = (tFld(i,j,kp1)-tFld(i,j,k)  )*maskC(i,j,kp1,bi,bj)
./pkg/generic_advdiff/gad_dst3fl_impl_r.F 146 GAD_DST3FL_IMPL_R
           mskM   = maskC(i,j,km2,bi,bj)*maskM2
./pkg/generic_advdiff/gad_dst3fl_impl_r.F 147 GAD_DST3FL_IMPL_R
           mskP   = maskC(i,j,kp1,bi,bj)*maskP1
./pkg/generic_advdiff/gad_fluxlimit_adv_r.F 107 GAD_FLUXLIMIT_ADV_R
         wT(i,j) = maskC(i,j,kM1,bi,bj)*(
./pkg/generic_advdiff/gad_fluxlimit_adv_r.F 88 GAD_FLUXLIMIT_ADV_R
     &        *maskC(i,j,kp1,bi,bj)
./pkg/generic_advdiff/gad_fluxlimit_adv_r.F 91 GAD_FLUXLIMIT_ADV_R
     &        *maskC(i,j,km2,bi,bj)
./pkg/generic_advdiff/gad_fluxlimit_impl_r.F 90 GAD_FLUXLIMIT_IMPL_R
         Rjp=(tFld(i,j,kp1)-tFld(i,j,k)  )*maskC(i,j,kp1,bi,bj)
./pkg/generic_advdiff/gad_fluxlimit_impl_r.F 92 GAD_FLUXLIMIT_IMPL_R
         Rjm=(tFld(i,j,km1)-tFld(i,j,km2))*maskC(i,j,km2,bi,bj)
./pkg/generic_advdiff/gad_implicit_r.F 128 GAD_IMPLICIT_R
           b5d(i,j,k) = -deltaTLev(k)*maskC(i,j,k-1,bi,bj)
./pkg/generic_advdiff/gad_implicit_r.F 140 GAD_IMPLICIT_R
           d5d(i,j,k) = -deltaTLev(k)*maskC(i,j,k+1,bi,bj)
./pkg/generic_advdiff/gad_implicit_r.F 200 GAD_IMPLICIT_R
     &                               *maskC(i,j,k-1,bi,bj)
./pkg/generic_advdiff/gad_implicit_r.F 318 GAD_IMPLICIT_R
     &            * maskC(i,j,k,bi,bj)
./pkg/generic_advdiff/gad_implicit_r.F 319 GAD_IMPLICIT_R
     &            * maskC(i,j,k-1,bi,bj)
./pkg/generic_advdiff/gad_os7mp_adv_r.F 74 GAD_OS7MP_ADV_R
         MskIpp  = maskC(i,j,kp2,bi,bj) * float(kp2-kp1)
./pkg/generic_advdiff/gad_os7mp_adv_r.F 75 GAD_OS7MP_ADV_R
         MskIp   = maskC(i,j,kp1,bi,bj) * float(kp1-k)
./pkg/generic_advdiff/gad_os7mp_adv_r.F 76 GAD_OS7MP_ADV_R
         MskI    = maskC(i,j,k,bi,bj)   * float(k-km1)
./pkg/generic_advdiff/gad_os7mp_adv_r.F 77 GAD_OS7MP_ADV_R
         MskIm   = maskC(i,j,km1,bi,bj) * float(km1-km2)
./pkg/generic_advdiff/gad_os7mp_adv_r.F 78 GAD_OS7MP_ADV_R
         MskImm  = maskC(i,j,km2,bi,bj) * float(km2-km3)
./pkg/generic_advdiff/gad_os7mp_adv_r.F 79 GAD_OS7MP_ADV_R
         MskImmm = maskC(i,j,km3,bi,bj) * float(km3-km4)
./pkg/generic_advdiff/gad_os7mp_adv_r.F 89 GAD_OS7MP_ADV_R
         MskIpp  = maskC(i,j,km2,bi,bj) * float(km2-km3)
./pkg/generic_advdiff/gad_os7mp_adv_r.F 90 GAD_OS7MP_ADV_R
         MskIp   = maskC(i,j,km1,bi,bj) * float(km1-km2)
./pkg/generic_advdiff/gad_os7mp_adv_r.F 91 GAD_OS7MP_ADV_R
         MskI    = maskC(i,j,k,bi,bj)   * float(k-km1)
./pkg/generic_advdiff/gad_os7mp_adv_r.F 92 GAD_OS7MP_ADV_R
         MskIm   = maskC(i,j,kp1,bi,bj) * float(kp1-k)
./pkg/generic_advdiff/gad_os7mp_adv_r.F 93 GAD_OS7MP_ADV_R
         MskImm  = maskC(i,j,kp2,bi,bj) * float(kp2-kp1)
./pkg/generic_advdiff/gad_os7mp_adv_r.F 94 GAD_OS7MP_ADV_R
         MskImmm = maskC(i,j,kp3,bi,bj) * float(kp3-kp2)
./pkg/generic_advdiff/gad_ppm_adv_r.F 113 GAD_PPM_ADV_R
     &            maskC(ix,iy,ir,bi,bj)
./pkg/generic_advdiff/gad_ppm_adv_x.F 108 GAD_PPM_ADV_X
     &          maskC(ix,iy,kk,bi,bj)
./pkg/generic_advdiff/gad_ppm_adv_y.F 108 GAD_PPM_ADV_Y
     &          maskC(ix,iy,kk,bi,bj)
./pkg/generic_advdiff/gad_pqm_adv_r.F 108 GAD_PQM_ADV_R
     &             maskC(ix,iy,ir,bi,bj)
./pkg/generic_advdiff/gad_pqm_adv_x.F 104 GAD_PQM_ADV_X
     &          maskC(ix,iy,kk,bi,bj)
./pkg/generic_advdiff/gad_pqm_adv_y.F 104 GAD_PQM_ADV_Y
     &          maskC(ix,iy,kk,bi,bj)
./pkg/generic_advdiff/gad_som_advect.F 270 GAD_SOM_ADVECT
     &                + (1. _d 0 - maskC(i,j,k,bi,bj))
./pkg/generic_advdiff/gad_som_advect.F 533 GAD_SOM_ADVECT
     &                 *maskC(i,j,k-1,bi,bj)
./pkg/generic_advdiff/gad_som_advect.F 545 GAD_SOM_ADVECT
           maskUp(i,j) = maskC(i,j,km1,bi,bj)*maskC(i,j,k,bi,bj)
./pkg/generic_advdiff/gad_u3_adv_r.F 71 GAD_U3_ADV_R
         Rjp = (tracer(i,j,kp1)-tracer(i,j, k ))*maskC(i,j,kp1,bi,bj)
./pkg/generic_advdiff/gad_u3_adv_r.F 73 GAD_U3_ADV_R
         Rjm = (tracer(i,j,km1)-tracer(i,j,km2))*maskC(i,j,km2,bi,bj)
./pkg/generic_advdiff/gad_u3_adv_r.F 76 GAD_U3_ADV_R
         wT(i,j) = maskC(i,j,km1,bi,bj)*(
./pkg/generic_advdiff/gad_u3c4_impl_r.F 104 GAD_U3C4_IMPL_R
         mskM   = maskC(i,j,km2,bi,bj)*maskM2
./pkg/generic_advdiff/gad_u3c4_impl_r.F 105 GAD_U3C4_IMPL_R
         mskP   = maskC(i,j,kp1,bi,bj)*maskP1
./pkg/generic_advdiff/gad_u3c4_impl_r.F 136 GAD_U3C4_IMPL_R
           mskM   = maskC(i,j,km2,bi,bj)*maskM2
./pkg/generic_advdiff/gad_u3c4_impl_r.F 137 GAD_U3C4_IMPL_R
           mskP   = maskC(i,j,kp1,bi,bj)*maskP1
./pkg/ggl90/ggl90_calc.F 426 GGL90_CALC
     &                            * maskC(i,j,k,bi,bj)
./pkg/ggl90/ggl90_calc.F 429 GGL90_CALC
         KappaM(i,j) = MAX(KappaM(i,j),viscArNr(k)) * maskC(i,j,k,bi,bj)
./pkg/ggl90/ggl90_calc.F 500 GGL90_CALC
         KappaE(i,j,k) = GGL90alpha * KappaM(i,j) * maskC(i,j,k,bi,bj)
./pkg/ggl90/ggl90_calc.F 566 GGL90_CALC
     &        *recip_drC(k)*maskC(i,j,k,bi,bj)
./pkg/ggl90/ggl90_calc.F 586 GGL90_CALC
     &        *recip_drC(k)*maskC(i,j,k-1,bi,bj)
./pkg/ggl90/ggl90_calc.F 622 GGL90_CALC
     &        * maskC(i,j,k,bi,bj)*maskC(i,j,km1,bi,bj)
./pkg/ggl90/ggl90_calc.F 640 GGL90_CALC
        GGL90TKE(i,j,kSurf,bi,bj) = maskC(i,j,kSurf,bi,bj)
./pkg/ggl90/ggl90_calc.F 672 GGL90_CALC
         GGL90TKE(i,j,k,bi,bj) = maskC(i,j,k,bi,bj)
./pkg/ggl90/ggl90_calc.F 697 GGL90_CALC
     &    +p8 *( (  maskC(i-1,j  ,k,bi,bj)*mskCor(i-1,j  ,bi,bj)
./pkg/ggl90/ggl90_calc.F 698 GGL90_CALC
     &           +  maskC(i+1,j  ,k,bi,bj)*mskCor(i+1,j  ,bi,bj) )
./pkg/ggl90/ggl90_calc.F 699 GGL90_CALC
     &         + (  maskC(i  ,j-1,k,bi,bj)*mskCor(i  ,j-1,bi,bj)
./pkg/ggl90/ggl90_calc.F 700 GGL90_CALC
     &           +  maskC(i  ,j+1,k,bi,bj)*mskCor(i  ,j+1,bi,bj) ) )
./pkg/ggl90/ggl90_calc.F 701 GGL90_CALC
     &    +p16*( (  maskC(i+1,j+1,k,bi,bj)* mskCor(i+1,j+1,bi,bj)
./pkg/ggl90/ggl90_calc.F 702 GGL90_CALC
     &           +  maskC(i-1,j-1,k,bi,bj)*mskCor(i-1,j-1,bi,bj) )
./pkg/ggl90/ggl90_calc.F 703 GGL90_CALC
     &         + (  maskC(i+1,j-1,k,bi,bj)*mskCor(i+1,j-1,bi,bj)
./pkg/ggl90/ggl90_calc.F 704 GGL90_CALC
     &           +  maskC(i-1,j+1,k,bi,bj)*mskCor(i-1,j+1,bi,bj) ) )
./pkg/ggl90/ggl90_calc.F 705 GGL90_CALC
     &               )*maskC(i,j,k,bi,bj)*mskCor(i,j,bi,bj)
./pkg/ggl90/ggl90_calc.F 728 GGL90_CALC
     &    +p8 *( (  maskC(i-1,j-1,k,bi,bj)*mskCor(i-1,j-1,bi,bj)
./pkg/ggl90/ggl90_calc.F 729 GGL90_CALC
     &           +  maskC(i  ,j-1,k,bi,bj)*mskCor(i  ,j-1,bi,bj) )
./pkg/ggl90/ggl90_calc.F 730 GGL90_CALC
     &         + (  maskC(i-1,j+1,k,bi,bj)*mskCor(i-1,j+1,bi,bj)
./pkg/ggl90/ggl90_calc.F 731 GGL90_CALC
     &           +  maskC(i  ,j+1,k,bi,bj)*mskCor(i  ,j+1,bi,bj) ) )
./pkg/ggl90/ggl90_calc.F 732 GGL90_CALC
     &               )*maskC(i-1,j,k,bi,bj)*mskCor(i-1,j,bi,bj)
./pkg/ggl90/ggl90_calc.F 733 GGL90_CALC
     &                *maskC(i  ,j,k,bi,bj)*mskCor(i  ,j,bi,bj)
./pkg/ggl90/ggl90_calc.F 758 GGL90_CALC
     &    +p8 *( (  maskC(i-1,j-1,k,bi,bj)*mskCor(i-1,j-1,bi,bj)
./pkg/ggl90/ggl90_calc.F 759 GGL90_CALC
     &           +  maskC(i-1,j  ,k,bi,bj)*mskCor(i-1,j  ,bi,bj) )
./pkg/ggl90/ggl90_calc.F 760 GGL90_CALC
     &         + (  maskC(i+1,j-1,k,bi,bj)*mskCor(i+1,j-1,bi,bj)
./pkg/ggl90/ggl90_calc.F 761 GGL90_CALC
     &           +  maskC(i+1,j  ,k,bi,bj)*mskCor(i+1,j  ,bi,bj) ) )
./pkg/ggl90/ggl90_calc.F 762 GGL90_CALC
     &               )*maskC(i,j-1,k,bi,bj)*mskCor(i,j-1,bi,bj)
./pkg/ggl90/ggl90_calc.F 763 GGL90_CALC
     &                *maskC(i,j  ,k,bi,bj)*mskCor(i,j  ,bi,bj)
./pkg/ggl90/ggl90_idemix.F 119 GGL90_IDEMIX
     &                 *maskC(i,j,k,bi,bj)*maskC(i,j,k-1,bi,bj)
./pkg/ggl90/ggl90_idemix.F 244 GGL90_IDEMIX
     &        IDEMIX_V0(i-1,j,k,bi,bj)*maskC(i-1,j,k,bi,bj)
./pkg/ggl90/ggl90_idemix.F 245 GGL90_IDEMIX
     &       +IDEMIX_V0(i  ,j,k,bi,bj)*maskC(i  ,j,k,bi,bj))
./pkg/ggl90/ggl90_idemix.F 261 GGL90_IDEMIX
     &        IDEMIX_V0(i,j  ,k,bi,bj)*maskC(i,j  ,k,bi,bj)
./pkg/ggl90/ggl90_idemix.F 262 GGL90_IDEMIX
     &       +IDEMIX_V0(i,j-1,k,bi,bj)*maskC(i,j-1,k,bi,bj) )
./pkg/ggl90/ggl90_idemix.F 281 GGL90_IDEMIX
     &         *maskC(i,j,k,bi,bj) ! paranoia setting
./pkg/ggl90/ggl90_idemix.F 333 GGL90_IDEMIX
     &        *maskC(i,j,k,bi,bj)!*maskC(i,j,km1,bi,bj)
./pkg/ggl90/ggl90_idemix.F 346 GGL90_IDEMIX
     &        *maskC(i,j,k-1,bi,bj)
./pkg/ggl90/ggl90_idemix.F 371 GGL90_IDEMIX
     &         *maskC(i,j,k,bi,bj)*maskC(i,j,km1,bi,bj)
./pkg/ggl90/ggl90_idemix.F 375 GGL90_IDEMIX
     &        *maskC(i,j,km1,bi,bj)
./pkg/ggl90/ggl90_idemix.F 379 GGL90_IDEMIX
     &         *maskC(i,j,k,bi,bj)
./pkg/ggl90/ggl90_idemix.F 391 GGL90_IDEMIX
     &          *maskC(i,j,k,bi,bj)*maskC(i,j,km1,bi,bj)
./pkg/ggl90/ggl90_idemix.F 395 GGL90_IDEMIX
     &          *maskC(i,j,k,bi,bj)
./pkg/ggl90/ggl90_idemix.F 399 GGL90_IDEMIX
     &          *maskC(i,j,k,bi,bj)*maskC(i,j,km1,bi,bj)
./pkg/ggl90/ggl90_idemix.F 403 GGL90_IDEMIX
     &          *maskC(i,j,km1,bi,bj)
./pkg/ggl90/ggl90_idemix.F 415 GGL90_IDEMIX
     &        *maskC(i,j,k,bi,bj)
./pkg/ggl90/ggl90_idemix.F 420 GGL90_IDEMIX
     &     *maskC(i,j,k,bi,bj)
./pkg/ggl90/ggl90_idemix.F 445 GGL90_IDEMIX
     &     /max(1. _d -12,Nsquare(i,j,k))*maskC(i,j,k,bi,bj)
./pkg/ggl90/ggl90_init_varia.F 151 GGL90_INIT_VARIA
     &            GGL90TKEmin)*maskC(I,J,K,bi,bj)
./pkg/ggl90/ggl90_init_varia.F 54 GGL90_INIT_VARIA
            GGL90TKE(I,J,K,bi,bj)=GGL90eps*maskC(I,J,K,bi,bj)
./pkg/ggl90/ggl90_init_varia.F 57 GGL90_INIT_VARIA
            GGL90TKE(I,J,K,bi,bj)=GGL90TKEmin*maskC(I,J,K,bi,bj)
./pkg/gmredi/gmredi_calc_tensor.F 1044 GMREDI_CALC_TENSOR
     &       ( maskC(i,j-1,k,bi,bj)*
./pkg/gmredi/gmredi_calc_tensor.F 1046 GMREDI_CALC_TENSOR
     &        +maskC(i, j ,k,bi,bj)*
./pkg/gmredi/gmredi_calc_tensor.F 1050 GMREDI_CALC_TENSOR
     &       ( maskC(i,j-1,kp1,bi,bj)*
./pkg/gmredi/gmredi_calc_tensor.F 1052 GMREDI_CALC_TENSOR
     &        +maskC(i, j ,kp1,bi,bj)*
./pkg/gmredi/gmredi_calc_tensor.F 316 GMREDI_CALC_TENSOR
     &                      )*maskC(i,j,k,bi,bj)
./pkg/gmredi/gmredi_calc_tensor.F 319 GMREDI_CALC_TENSOR
     &                      )*maskC(i,j,k,bi,bj)
./pkg/gmredi/gmredi_calc_tensor.F 345 GMREDI_CALC_TENSOR
          IF ( maskC(i,j,k,bi,bj).NE.0. ) THEN
./pkg/gmredi/gmredi_calc_tensor.F 369 GMREDI_CALC_TENSOR
     &                 + maskC(i-1,j,k,bi,bj)
./pkg/gmredi/gmredi_calc_tensor.F 370 GMREDI_CALC_TENSOR
     &                 + maskC(i+1,j,k,bi,bj)
./pkg/gmredi/gmredi_calc_tensor.F 371 GMREDI_CALC_TENSOR
     &                 + maskC(i,j-1,k,bi,bj)
./pkg/gmredi/gmredi_calc_tensor.F 372 GMREDI_CALC_TENSOR
     &                 + maskC(i,j+1,k,bi,bj)
./pkg/gmredi/gmredi_calc_tensor.F 435 GMREDI_CALC_TENSOR
         SlopeX(i,j)=SlopeX(i,j)*maskC(i,j,k,bi,bj)
./pkg/gmredi/gmredi_calc_tensor.F 436 GMREDI_CALC_TENSOR
         SlopeY(i,j)=SlopeY(i,j)*maskC(i,j,k,bi,bj)
./pkg/gmredi/gmredi_calc_tensor.F 437 GMREDI_CALC_TENSOR
         SlopeSqr(i,j)=SlopeSqr(i,j)*maskC(i,j,k,bi,bj)
./pkg/gmredi/gmredi_calc_tensor.F 826 GMREDI_CALC_TENSOR
     &       ( maskC(i-1,j,k,bi,bj)*
./pkg/gmredi/gmredi_calc_tensor.F 828 GMREDI_CALC_TENSOR
     &        +maskC( i ,j,k,bi,bj)*
./pkg/gmredi/gmredi_calc_tensor.F 832 GMREDI_CALC_TENSOR
     &       ( maskC(i-1,j,kp1,bi,bj)*
./pkg/gmredi/gmredi_calc_tensor.F 834 GMREDI_CALC_TENSOR
     &        +maskC( i ,j,kp1,bi,bj)*
./pkg/gmredi/gmredi_diagnostics_impl.F 62 GMREDI_DIAGNOSTICS_IMPL
     &            * maskC(i,j,k,bi,bj)*maskC(i,j,k-1,bi,bj)
./pkg/gmredi/gmredi_k3d.F 383 GMREDI_K3D
         tmpU(i,j) = tmpU(i,j) * maskC(i,j,k,bi,bj)
./pkg/gmredi/gmredi_k3d.F 384 GMREDI_K3D
         tmpV(i,j) = tmpV(i,j) * maskC(i,j,k,bi,bj)
./pkg/gmredi/gmredi_k3d.F 424 GMREDI_K3D
         N2(i,j,k) = MAX(N2(i,j,k),GM_K3D_minN2)*maskC(i,j,k,bi,bj)
./pkg/gmredi/gmredi_k3d.F 457 GMREDI_K3D
     I       kLow_C, maskC(:,:,:,bi,bj),
./pkg/gmredi/gmredi_k3d.F 622 GMREDI_K3D
        IF(maskC(i,j,1,bi,bj).EQ.0.0) surfkz(i,j)=0.0
./pkg/gmredi/gmredi_residual_flow.F 87 GMREDI_RESIDUAL_FLOW
     &              )*maskC(i,j,k,bi,bj)
./pkg/gmredi/gmredi_xtransport.F 158 GMREDI_XTRANSPORT
     &       ( maskC(i-1,j,k,bi,bj)*
./pkg/gmredi/gmredi_xtransport.F 160 GMREDI_XTRANSPORT
     &        +maskC( i ,j,k,bi,bj)*
./pkg/gmredi/gmredi_xtransport.F 164 GMREDI_XTRANSPORT
     &       ( maskC(i-1,j,kp1,bi,bj)*
./pkg/gmredi/gmredi_xtransport.F 166 GMREDI_XTRANSPORT
     &        +maskC( i ,j,kp1,bi,bj)*
./pkg/gmredi/gmredi_ytransport.F 157 GMREDI_YTRANSPORT
     &       ( maskC(i,j-1,k,bi,bj)*
./pkg/gmredi/gmredi_ytransport.F 159 GMREDI_YTRANSPORT
     &        +maskC(i, j ,k,bi,bj)*
./pkg/gmredi/gmredi_ytransport.F 163 GMREDI_YTRANSPORT
     &       ( maskC(i,j-1,kp1,bi,bj)*
./pkg/gmredi/gmredi_ytransport.F 165 GMREDI_YTRANSPORT
     &        +maskC(i, j ,kp1,bi,bj)*
./pkg/grdchk/grdchk_get_position.F 145 GRDCHK_GET_POSITION
                     if ( maskC(i,j,k,bi,bj) .gt. 0.) then
./pkg/grdchk/grdchk_get_position.F 148 GRDCHK_GET_POSITION
                     wetlocal = maskC(i,j,k,bi,bj)
./pkg/grdchk/grdchk_loc.F 193 GRDCHK_LOC
                     if ( maskC(i,j,k,bi,bj) .gt. 0.) then
./pkg/icefront/icefront_thermodynamics.F 146 ICEFRONT_THERMODYNAMICS
            FreshWaterFlux(I,J) = maskC(I,J,K,bi,bj) *
./pkg/icefront/icefront_thermodynamics.F 150 ICEFRONT_THERMODYNAMICS
            HeatFlux(I,J) = maskC(I,J,K,bi,bj) * HeatCapacity_Cp *
./pkg/kl10/kl10_calc.F 105 KL10_CALC
               IF ( (rhoS(K).LT.rhoS(K-1)).AND.(maskC(I,J,K,bi
./model/src/apply_forcing.F 592 APPLY_FORCING_T
     &                   )*maskC(i,j,km,bi,bj)*0.25 _d 0
./model/src/apply_forcing.F 605 APPLY_FORCING_T
     &                   )*maskC(i,j,kp,bi,bj)*0.25 _d 0
./model/src/apply_forcing.F 690 APPLY_FORCING_T
     &   -Qsw(i,j,bi,bj)*(swfracb(1)*maskC(i,j,k,bi,bj)
./model/src/apply_forcing.F 691 APPLY_FORCING_T
     &                   -swfracb(2)*maskC(i,j,kp1, bi,bj))
./pkg/kpp/kpp_calc.F 336 KPP_CALC
               dbloc(i,j,k) = dbloc(i,j,k) * maskC(i,j,k,bi,bj)
./pkg/kpp/kpp_calc.F 337 KPP_CALC
     &              * maskC(i,j,km1,bi,bj)
./pkg/kpp/kpp_calc.F 338 KPP_CALC
               ghat(i,j,k)  = ghat(i,j,k)  * maskC(i,j,k,bi,bj)
./pkg/kpp/kpp_calc.F 339 KPP_CALC
     &              * maskC(i,j,km1,bi,bj)
./pkg/kpp/kpp_calc.F 340 KPP_CALC
               Ritop(i,j,k) = Ritop(i,j,k) * maskC(i,j,k,bi,bj)
./pkg/kpp/kpp_calc.F 341 KPP_CALC
     &              * maskC(i,j,km1,bi,bj)
./pkg/kpp/kpp_calc.F 549 KPP_CALC
     I     , maskC(1-OLx,1-OLy,1,bi,bj)
./pkg/kpp/kpp_calc.F 574 KPP_CALC
         KPPviscAz(i,j,k,bi,bj) = vddiff(i,j,k-1,1) * maskC(i,j,k,bi,bj)
./pkg/kpp/kpp_calc.F 575 KPP_CALC
     &        * maskC(i,j,km1,bi,bj)
./pkg/kpp/kpp_calc.F 576 KPP_CALC
         KPPdiffKzS(i,j,k,bi,bj)= vddiff(i,j,k-1,2) * maskC(i,j,k,bi,bj)
./pkg/kpp/kpp_calc.F 577 KPP_CALC
     &        * maskC(i,j,km1,bi,bj)
./pkg/kpp/kpp_calc.F 578 KPP_CALC
         KPPdiffKzT(i,j,k,bi,bj)= vddiff(i,j,k-1,3) * maskC(i,j,k,bi,bj)
./pkg/kpp/kpp_calc.F 579 KPP_CALC
     &        * maskC(i,j,km1,bi,bj)
./pkg/kpp/kpp_calc.F 580 KPP_CALC
         KPPghat(i,j,k,bi,bj)   = ghat(i,j,k)       * maskC(i,j,k,bi,bj)
./pkg/kpp/kpp_calc.F 581 KPP_CALC
     &        * maskC(i,j,km1,bi,bj)
./pkg/kpp/kpp_calc.F 587 KPP_CALC
        KPPhbl(i,j,bi,bj) = hbl(i,j) * maskC(i,j,k,bi,bj)
./pkg/kpp/kpp_forcing_surf.F 307 KPP_FORCING_SURF
     &        maskC(I,J,k,bi,bj) .GT. 0. .AND.
./model/src/calc_adv_flow.F 129 CALC_ADV_FLOW
           maskUp(i,j) = maskC(i,j,k-1,bi,bj)*maskC(i,j,k,bi,bj)
./model/src/calc_adv_flow.F 96 CALC_ADV_FLOW
     &                  * maskC(i,j,k,bi,bj)*maskC(i,j,k+1,bi,bj)
./pkg/kpp/kpp_routines.F 1254 SMOOTH_HORIZ
     &           p25   *   maskC(i  ,j  ,k,bi,bj)   +
./pkg/kpp/kpp_routines.F 1255 SMOOTH_HORIZ
     &           p125  * ( maskC(im1,j  ,k,bi,bj)   +
./pkg/kpp/kpp_routines.F 1256 SMOOTH_HORIZ
     &                     maskC(ip1,j  ,k,bi,bj)   +
./pkg/kpp/kpp_routines.F 1257 SMOOTH_HORIZ
     &                     maskC(i  ,jm1,k,bi,bj)   +
./pkg/kpp/kpp_routines.F 1258 SMOOTH_HORIZ
     &                     maskC(i  ,jp1,k,bi,bj) ) +
./pkg/kpp/kpp_routines.F 1259 SMOOTH_HORIZ
     &           p0625 * ( maskC(im1,jm1,k,bi,bj)   +
./pkg/kpp/kpp_routines.F 1260 SMOOTH_HORIZ
     &                     maskC(im1,jp1,k,bi,bj)   +
./pkg/kpp/kpp_routines.F 1261 SMOOTH_HORIZ
     &                     maskC(ip1,jm1,k,bi,bj)   +
./pkg/kpp/kpp_routines.F 1262 SMOOTH_HORIZ
     &                     maskC(ip1,jp1,k,bi,bj) )
./pkg/kpp/kpp_routines.F 1265 SMOOTH_HORIZ
     &              p25  * fld(i  ,j  )*maskC(i  ,j  ,k,bi,bj) +
./pkg/kpp/kpp_routines.F 1266 SMOOTH_HORIZ
     &              p125 *(fld(im1,j  )*maskC(im1,j  ,k,bi,bj) +
./pkg/kpp/kpp_routines.F 1267 SMOOTH_HORIZ
     &                     fld(ip1,j  )*maskC(ip1,j  ,k,bi,bj) +
./pkg/kpp/kpp_routines.F 1268 SMOOTH_HORIZ
     &                     fld(i  ,jm1)*maskC(i  ,jm1,k,bi,bj) +
./pkg/kpp/kpp_routines.F 1269 SMOOTH_HORIZ
     &                     fld(i  ,jp1)*maskC(i  ,jp1,k,bi,bj))+
./pkg/kpp/kpp_routines.F 1270 SMOOTH_HORIZ
     &              p0625*(fld(im1,jm1)*maskC(im1,jm1,k,bi,bj) +
./pkg/kpp/kpp_routines.F 1271 SMOOTH_HORIZ
     &                     fld(im1,jp1)*maskC(im1,jp1,k,bi,bj) +
./pkg/kpp/kpp_routines.F 1272 SMOOTH_HORIZ
     &                     fld(ip1,jm1)*maskC(ip1,jm1,k,bi,bj) +
./pkg/kpp/kpp_routines.F 1273 SMOOTH_HORIZ
     &                     fld(ip1,jp1)*maskC(ip1,jp1,k,bi,bj)))
./pkg/layers/layers_fluxcalc.F 529 LAYERS_DIAPYCNAL
         IF (maskC(i,j,kp1,bi,bj).EQ.zeroRS) kp1=k
./pkg/layers/layers_thermodynamics.F 153 LAYERS_CALC_RHS
     &            -Qsw(i,j,bi,bj)*(swfracb(1)*maskC(i,j,k,bi,bj)
./pkg/layers/layers_thermodynamics.F 154 LAYERS_CALC_RHS
     &                      -swfracb(2)*maskC(i,j,kdown,bi,bj))
./model/src/calc_gw.F 194 CALC_GW
           IF ( maskC(i,j,k-1,bi,bj).EQ.0. .OR.
./model/src/calc_gw.F 195 CALC_GW
     &          maskC(i,j, k ,bi,bj).EQ.0. ) THEN
./model/src/calc_gw.F 244 CALC_GW
           IF ( maskC(i,j,k,bi,bj).EQ.0. ) THEN
./model/src/calc_oce_mxlayer.F 202 CALC_OCE_MXLAYER
     &                       tmpFac * ( maskC(i-1,j,1,bi,bj) + 
./model/src/calc_oce_mxlayer.F 203 CALC_OCE_MXLAYER
     &                                  maskC(i+1,j,1,bi,bj) +
./model/src/calc_oce_mxlayer.F 204 CALC_OCE_MXLAYER
     &                                  maskC(i,j-1,1,bi,bj) + 
./model/src/calc_oce_mxlayer.F 205 CALC_OCE_MXLAYER
     &                                  maskC(i,j+1,1,bi,bj) )
./model/src/calc_oce_mxlayer.F 206 CALC_OCE_MXLAYER
     &                  ) * maskC(i,j,1,bi,bj)
./model/src/calc_phi_hyd.F 194 CALC_PHI_HYD
           alphaRho(i,j) = alphaRho(i,j)*maskC(i,j,k,bi,bj)
./model/src/calc_phi_hyd.F 349 CALC_PHI_HYD
            alphaRho(i,j)=maskC(i,j,k,bi,bj)*
./model/src/calc_phi_hyd.F 449 CALC_PHI_HYD
     &                      - thetaRef )*maskC(i,j,k,bi,bj)
./model/src/calc_surf_dr.F 123 CALC_SURF_DR
     &              )*recip_drF(ks)*maskC(i,j,ks,bi,bj)
./pkg/mom_common/mom_u_implicit_r.F 150 MOM_U_IMPLICIT_R
     &                                 *maskC( i ,j,k-1,bi,bj)
./pkg/mom_common/mom_u_implicit_r.F 152 MOM_U_IMPLICIT_R
     &                                 *maskC(i-1,j,k-1,bi,bj)
./pkg/mom_common/mom_v_implicit_r.F 150 MOM_V_IMPLICIT_R
     &                                 *maskC(i, j ,k-1,bi,bj)
./pkg/mom_common/mom_v_implicit_r.F 152 MOM_V_IMPLICIT_R
     &                                 *maskC(i,j-1,k-1,bi,bj)
./pkg/mom_fluxform/mom_u_adv_wu.F 100 MOM_U_ADV_WU
     &          (maskC( i ,j,k,bi,bj)-maskC( i ,j,k-1,bi,bj))
./pkg/mom_fluxform/mom_u_adv_wu.F 102 MOM_U_ADV_WU
     &          (maskC(i-1,j,k,bi,bj)-maskC(i-1,j,k-1,bi,bj))
./pkg/mom_fluxform/mom_v_adv_wv.F 101 MOM_V_ADV_WV
     &          (maskC(i, j ,k,bi,bj)-maskC(i, j ,k-1,bi,bj))
./pkg/mom_fluxform/mom_v_adv_wv.F 103 MOM_V_ADV_WV
     &          (maskC(i,j-1,k,bi,bj)-maskC(i,j-1,k-1,bi,bj))
./model/src/cg3d.F 144 CG3D
     &                          * maskC(i,j,k,bi,bj)
./pkg/mom_vecinv/mom_vi_u_vertshear.F 69 MOM_VI_U_VERTSHEAR
     &    wFld(I,J,K,bi,bj)*rA(i,j,bi,bj)*maskC(I,J,Km1,bi,bj)
./pkg/mom_vecinv/mom_vi_u_vertshear.F 70 MOM_VI_U_VERTSHEAR
     &   +wFld(I-1,J,K,bi,bj)*rA(i-1,j,bi,bj)*maskC(I-1,J,Km1,bi,bj)
./pkg/mom_vecinv/mom_vi_u_vertshear.F 83 MOM_VI_U_VERTSHEAR
     &    wFld(I,J,K,bi,bj)*maskC(I,J,Km1,bi,bj)
./pkg/mom_vecinv/mom_vi_u_vertshear.F 84 MOM_VI_U_VERTSHEAR
     &   +wFld(I-1,J,K,bi,bj)*maskC(I-1,J,Km1,bi,bj)
./pkg/mom_vecinv/mom_vi_v_vertshear.F 69 MOM_VI_V_VERTSHEAR
     &    wFld(I,J,K,bi,bj)*rA(i,j,bi,bj)*maskC(i,j,Km1,bi,bj)
./pkg/mom_vecinv/mom_vi_v_vertshear.F 70 MOM_VI_V_VERTSHEAR
     &   +wFld(I,J-1,K,bi,bj)*rA(i,j-1,bi,bj)*maskC(i,j-1,Km1,bi,bj)
./pkg/mom_vecinv/mom_vi_v_vertshear.F 83 MOM_VI_V_VERTSHEAR
     &    wFld(I,J,K,bi,bj)*maskC(i,j,Km1,bi,bj)
./pkg/mom_vecinv/mom_vi_v_vertshear.F 84 MOM_VI_V_VERTSHEAR
     &   +wFld(I,J-1,K,bi,bj)*maskC(i,j-1,Km1,bi,bj)
./model/src/cg3d_ex0.F 137 CG3D_EX0
     &                          * maskC(i,j,k,bi,bj)
./pkg/monitor/mon_ke.F 117 MON_KE
     &        )*maskC(i,j,k,bi,bj)*maskInC(i,j,bi,bj)
./pkg/monitor/mon_surfcor.F 134 MON_SURFCOR
     &           *maskC(i,j,k-1,bi,bj)*maskC(i,j,k,bi,bj)
./pkg/monitor/monitor.F 116 MONITOR
     &            maskC, maskInC, rA , thickFacF, dummyRL, myThid )
./pkg/monitor/monitor_ad.F 134 ADMONITOR
     &             maskC, maskInC, rA , drC, dummyRL, myThid )
./pkg/monitor/monitor_ad.F 164 ADMONITOR
     &             maskC, maskInC, rA , drC, dummyRL, myThid )
./pkg/monitor/monitor_g.F 121 G_MONITOR
     &           maskC, maskInC, rA , drC, dummyRL, myThid )
./pkg/my82/my82_calc.F 177 MY82_CALC
     &        * maskC(I,J,K,bi,bj)
./pkg/my82/my82_calc.F 179 MY82_CALC
     &        * maskC(I,J,K,bi,bj)
./pkg/mypackage/mypackage_tendency_apply.F 195 MYPACKAGE_TENDENCY_APPLY_T
     &     + maskC(i,j,k,bi,bj) * myPa_TendScal1(i,j,k,bi,bj)
./pkg/mypackage/mypackage_tendency_apply.F 252 MYPACKAGE_TENDENCY_APPLY_S
     &     + maskC(i,j,k,bi,bj) * myPa_TendScal2(i,j,k,bi,bj)
./pkg/obcs/obcs_apply_w.F 111 OBCS_APPLY_W
            obc_mask = maskC(Iobc,j,km1,bi,bj)*maskC(Iobc,j,k,bi,bj)
./pkg/obcs/obcs_apply_w.F 115 OBCS_APPLY_W
            obc_mask = maskC(Iobc,j,km1,bi,bj)*maskC(Iobc,j,k,bi,bj)
./pkg/obcs/obcs_apply_w.F 130 OBCS_APPLY_W
            obc_mask = maskC(Iobc,j,km1,bi,bj)*maskC(Iobc,j,k,bi,bj)
./pkg/obcs/obcs_apply_w.F 134 OBCS_APPLY_W
            obc_mask = maskC(Iobc,j,km1,bi,bj)*maskC(Iobc,j,k,bi,bj)
./pkg/obcs/obcs_apply_w.F 72 OBCS_APPLY_W
            obc_mask = maskC(i,Jobc,km1,bi,bj)*maskC(i,Jobc,k,bi,bj)
./pkg/obcs/obcs_apply_w.F 76 OBCS_APPLY_W
            obc_mask = maskC(i,Jobc,km1,bi,bj)*maskC(i,Jobc,k,bi,bj)
./pkg/obcs/obcs_apply_w.F 91 OBCS_APPLY_W
            obc_mask = maskC(i,Jobc,km1,bi,bj)*maskC(i,Jobc,k,bi,bj)
./pkg/obcs/obcs_apply_w.F 95 OBCS_APPLY_W
            obc_mask = maskC(i,Jobc,km1,bi,bj)*maskC(i,Jobc,k,bi,bj)
./pkg/obcs/obcs_calc_stevens.F 1055 OBCS_STEVENS_SAVE_TRACERS
     &           *maskC(i,Jobc+1,k,bi,bj)
./pkg/obcs/obcs_calc_stevens.F 1057 OBCS_STEVENS_SAVE_TRACERS
     &           *maskC(i,Jobc+1,k,bi,bj)
./pkg/obcs/obcs_calc_stevens.F 1071 OBCS_STEVENS_SAVE_TRACERS
     &           *maskC(i,Jobc+1,k,bi,bj)
./pkg/obcs/obcs_calc_stevens.F 1073 OBCS_STEVENS_SAVE_TRACERS
     &           *maskC(i,Jobc+1,k,bi,bj)
./pkg/obcs/obcs_calc_stevens.F 1087 OBCS_STEVENS_SAVE_TRACERS
     &           *maskC(Iobc-1,j,k,bi,bj)
./pkg/obcs/obcs_calc_stevens.F 1089 OBCS_STEVENS_SAVE_TRACERS
     &           *maskC(Iobc-1,j,k,bi,bj)
./pkg/obcs/obcs_calc_stevens.F 1103 OBCS_STEVENS_SAVE_TRACERS
     &           *maskC(Iobc+1,j,k,bi,bj)
./pkg/obcs/obcs_calc_stevens.F 1105 OBCS_STEVENS_SAVE_TRACERS
     &           *maskC(Iobc+1,j,k,bi,bj)
./pkg/obcs/obcs_seaice_sponge.F 125 OBCS_SEAICE_SPONGE_A
     &           - maskC(i,j,1,bi,bj) * lambda_obcs
./pkg/obcs/obcs_seaice_sponge.F 152 OBCS_SEAICE_SPONGE_A
     &           - maskC(i,j,1,bi,bj) * lambda_obcs
./pkg/obcs/obcs_seaice_sponge.F 228 OBCS_SEAICE_SPONGE_H
     &           - maskC(i,j,1,bi,bj) * lambda_obcs
./pkg/obcs/obcs_seaice_sponge.F 255 OBCS_SEAICE_SPONGE_H
     &           - maskC(i,j,1,bi,bj) * lambda_obcs
./pkg/obcs/obcs_seaice_sponge.F 282 OBCS_SEAICE_SPONGE_H
     &           - maskC(i,j,1,bi,bj) * lambda_obcs
./pkg/obcs/obcs_seaice_sponge.F 309 OBCS_SEAICE_SPONGE_H
     &           - maskC(i,j,1,bi,bj) * lambda_obcs
./pkg/obcs/obcs_seaice_sponge.F 385 OBCS_SEAICE_SPONGE_SL
     &           - maskC(i,j,1,bi,bj) * lambda_obcs
./pkg/obcs/obcs_seaice_sponge.F 412 OBCS_SEAICE_SPONGE_SL
     &           - maskC(i,j,1,bi,bj) * lambda_obcs
./pkg/obcs/obcs_seaice_sponge.F 439 OBCS_SEAICE_SPONGE_SL
     &           - maskC(i,j,1,bi,bj) * lambda_obcs
./pkg/obcs/obcs_seaice_sponge.F 466 OBCS_SEAICE_SPONGE_SL
     &           - maskC(i,j,1,bi,bj) * lambda_obcs
./pkg/obcs/obcs_seaice_sponge.F 542 OBCS_SEAICE_SPONGE_SN
     &           - maskC(i,j,1,bi,bj) * lambda_obcs
./pkg/obcs/obcs_seaice_sponge.F 569 OBCS_SEAICE_SPONGE_SN
     &           - maskC(i,j,1,bi,bj) * lambda_obcs
./pkg/obcs/obcs_seaice_sponge.F 596 OBCS_SEAICE_SPONGE_SN
     &           - maskC(i,j,1,bi,bj) * lambda_obcs
./pkg/obcs/obcs_seaice_sponge.F 623 OBCS_SEAICE_SPONGE_SN
     &           - maskC(i,j,1,bi,bj) * lambda_obcs
./pkg/obcs/obcs_seaice_sponge.F 71 OBCS_SEAICE_SPONGE_A
     &           - maskC(i,j,1,bi,bj) * lambda_obcs
./pkg/obcs/obcs_seaice_sponge.F 98 OBCS_SEAICE_SPONGE_A
     &           - maskC(i,j,1,bi,bj) * lambda_obcs
./pkg/obcs/obcs_sponge.F 496 OBCS_SPONGE_T
     &            -  maskC(i,j,k,bi,bj) * lambda_obcs_t
./pkg/obcs/obcs_sponge.F 536 OBCS_SPONGE_T
     &            - maskC(i,j,k,bi,bj) * lambda_obcs_t
./pkg/obcs/obcs_sponge.F 576 OBCS_SPONGE_T
     &            - maskC(i,j,k,bi,bj) * lambda_obcs_t
./pkg/obcs/obcs_sponge.F 616 OBCS_SPONGE_T
     &            - maskC(i,j,k,bi,bj) * lambda_obcs_t
./pkg/obcs/obcs_sponge.F 728 OBCS_SPONGE_S
     &            - maskC(i,j,k,bi,bj) * lambda_obcs_s
./pkg/obcs/obcs_sponge.F 768 OBCS_SPONGE_S
     &           - maskC(i,j,k,bi,bj) * lambda_obcs_s
./pkg/obcs/obcs_sponge.F 808 OBCS_SPONGE_S
     &           - maskC(i,j,k,bi,bj) * lambda_obcs_s
./pkg/obcs/obcs_sponge.F 848 OBCS_SPONGE_S
     &           - maskC(i,j,k,bi,bj) * lambda_obcs_s
./model/src/do_atmospheric_phys.F 95 DO_ATMOSPHERIC_PHYS
     &         - thetaRef )*maskC(i,j,k,bi,bj)
./pkg/opps/opps_calc.F 485 OPPS_CALC
     &                  )*maskC(i,j,kRef,bi,bj)
./pkg/opps/opps_calc.F 496 OPPS_CALC
     &                  )*maskC(i,j,kRef,bi,bj)
./pkg/opps/opps_calc.F 501 OPPS_CALC
        pLoc = pRef4EOS(kRef)*maskC(i,j,kRef,bi,bj)
./pkg/opps/opps_calc.F 503 OPPS_CALC
        pLoc = rhoConst*phiRef(2*kRef)*maskC(i,j,kRef,bi,bj)
./pkg/opps/opps_calc.F 507 OPPS_CALC
       pLoc = rC(kRef)* maskC(i,j,kRef,bi,bj)
./pkg/pp81/pp81_calc.F 115 PP81_CALC
     &        * maskC(I,J,K,bi,bj)
./pkg/pp81/pp81_calc.F 117 PP81_CALC
     &        * maskC(I,J,K,bi,bj)
./model/src/do_oceanic_phys.F 816 DO_OCEANIC_PHYS
     &         - thetaRef )*maskC(i,j,k,bi,bj)
./pkg/profiles/profiles_interp.F 107 PROFILES_INTERP
        tab_coeffs3(q)=maskC(i_cur(q),j_cur(q),k,bi,bj)
./pkg/ptracers/ptracers_init_varia.F 117 PTRACERS_INIT_VARIA
            IF (maskC(i,j,k,bi,bj).EQ.0.)
./pkg/ptracers/ptracers_reset.F 98 PTRACERS_RESET
             IF (maskC(i,j,k,bi,bj).EQ.0.)
./pkg/rbcs/rbcs_init_fixed.F 244 RBCS_INIT_FIXED
     &                                    * maskC(i,j,k,bi,bj)
./model/src/external_forcing.F 464 EXTERNAL_FORCING_T
     &                   )*maskC(i,j,km,bi,bj)*0.25 _d 0
./model/src/external_forcing.F 477 EXTERNAL_FORCING_T
     &                   )*maskC(i,j,kp,bi,bj)*0.25 _d 0
./model/src/external_forcing.F 547 EXTERNAL_FORCING_T
     &   -Qsw(i,j,bi,bj)*(swfracb(1)*maskC(i,j,kLev,bi,bj)
./model/src/external_forcing.F 548 EXTERNAL_FORCING_T
     &                   -swfracb(2)*maskC(i,j,kp1, bi,bj))
./pkg/salt_plume/salt_plume_tendency_apply_s.F 106 SALT_PLUME_TENDENCY_APPLY_S
     &                     *maskC(i,j,k,bi,bj)
./pkg/salt_plume/salt_plume_tendency_apply_s.F 144 SALT_PLUME_TENDENCY_APPLY_S
         plumefrac(I,J) = (plumekb(2)-plumekb(1))*maskC(i,j,k,bi,bj)
./pkg/sbo/sbo_calc.F 230 SBO_CALC
            darea = rA(i,j,bi,bj)*maskC(i,j,1,bi,bj)
./pkg/sbo/sbo_calc.F 320 SBO_CALC
           IF ( maskC(i,j,1,bi,bj) .NE. 0. ) THEN
./pkg/sbo/sbo_calc.F 323 SBO_CALC
              darea = rA(i,j,bi,bj)*maskC(i,j,1,bi,bj)
./pkg/sbo/sbo_calc.F 348 SBO_CALC
     &             * maskC(i,j,k,bi,bj)*hFacC(i,j,k,bi,bj)
./pkg/seaice/seaice_bottomdrag_coeffs.F 153 SEAICE_BOTTOMDRAG_COEFFS
     &           * maskC(I,J,kSrf,bi,bj)
./pkg/seaice/seaice_budget_ocean.F 148 SEAICE_BUDGET_OCEAN
        sflux(i,j,bi,bj) = sflux(i,j,bi,bj)*maskC(i,j,1,bi,bj)
./model/src/grad_sigma.F 102 GRAD_SIGMA
         sigmaR(i,j,k)= maskC(i,j,k,bi,bj)*maskC(i,j,k-1,bi,bj)
./pkg/seaice/seaice_calc_strainrates.F 157 SEAICE_CALC_STRAINRATES
     &         *maskC(i  ,j  ,k,bi,bj)*maskC(i-1,j  ,k,bi,bj)
./pkg/seaice/seaice_calc_strainrates.F 158 SEAICE_CALC_STRAINRATES
     &         *maskC(i  ,j-1,k,bi,bj)*maskC(i-1,j-1,k,bi,bj)
./pkg/seaice/seaice_calc_viscosities.F 199 SEAICE_CALC_VISCOSITIES
          sumNorm  = maskC(I,J,  k,bi,bj)+maskC(I-1,J,  k,bi,bj)
./pkg/seaice/seaice_calc_viscosities.F 200 SEAICE_CALC_VISCOSITIES
     &             + maskC(I,J-1,k,bi,bj)+maskC(I-1,J-1,k,bi,bj)
./pkg/seaice/seaice_calc_viscosities.F 216 SEAICE_CALC_VISCOSITIES
           maskZ = maskC(I,J,  k,bi,bj)*maskC(I-1,J,  k,bi,bj)
./pkg/seaice/seaice_calc_viscosities.F 217 SEAICE_CALC_VISCOSITIES
     &          *  maskC(I,J-1,k,bi,bj)*maskC(I-1,J-1,k,bi,bj)
./pkg/seaice/seaice_cost_driver.F 103 SEAICE_COST_DRIVER
     &     maskC, wsmrarea,
./pkg/seaice/seaice_cost_driver.F 114 SEAICE_COST_DRIVER
     &     maskC, wsmrarea,
./pkg/seaice/seaice_cost_driver.F 92 SEAICE_COST_DRIVER
     &     maskC, wsmrarea,
./model/src/impldiff.F 341 IMPLDIFF
     &            * maskC(i,j,k,bi,bj)
./model/src/impldiff.F 342 IMPLDIFF
     &            * maskC(i,j,k-1,bi,bj)
./pkg/seaice/seaice_evp.F 438 SEAICE_EVP
     &            * recip_rA(I,J,bi,bj) ) * maskC(I,J,kSrf,bi,bj)
./pkg/seaice/seaice_evp.F 448 SEAICE_EVP
            sumNorm = maskC(I,J,  kSrf,bi,bj)+maskC(I-1,J,  kSrf,bi,bj)
./pkg/seaice/seaice_evp.F 449 SEAICE_EVP
     &              + maskC(I,J-1,kSrf,bi,bj)+maskC(I-1,J-1,kSrf,bi,bj)
./pkg/seaice/seaice_fake.F 62 SEAICE_FAKE
            EmPmR(i,j,bi,bj)  = maskC(i,j,1,bi,bj)*(
./pkg/seaice/seaice_growth.F 1061 SEAICE_GROWTH
     &         * SEAICE_deltaTtherm * maskC(i,j,kSurface,bi,bj)
./pkg/seaice/seaice_growth.F 2213 SEAICE_GROWTH
     &           ) * maskC(I,J,kSurface,bi,bj)
./pkg/seaice/seaice_growth.F 2262 SEAICE_GROWTH
           tmpscal3=rhoConstFresh*maskC(I,J,kSurface,bi,bj)*(
./pkg/seaice/seaice_growth.F 2315 SEAICE_GROWTH
     &            maskC(I,J,kSurface,bi,bj)*convertHI2Q*(
./pkg/seaice/seaice_growth.F 2321 SEAICE_GROWTH
           tmpscal1 = rhoConstFresh*maskC(I,J,kSurface,bi,bj)
./pkg/seaice/seaice_growth.F 2325 SEAICE_GROWTH
           tmpscal2=rhoConstFresh*maskC(I,J,kSurface,bi,bj)*
./pkg/seaice/seaice_growth.F 2374 SEAICE_GROWTH
          EmPmR(I,J,bi,bj)  = maskC(I,J,kSurface,bi,bj)*(
./pkg/seaice/seaice_growth.F 2386 SEAICE_GROWTH
     &           * maskC(I,J,kSurface,bi,bj)
./pkg/seaice/seaice_growth.F 2389 SEAICE_GROWTH
           SIatmFW(I,J,bi,bj) = maskC(I,J,kSurface,bi,bj)*(
./pkg/seaice/seaice_growth.F 2405 SEAICE_GROWTH
            frWtrAtm(I,J,bi,bj) = maskC(I,J,kSurface,bi,bj)*(
./pkg/seaice/seaice_growth.F 2423 SEAICE_GROWTH
          frWtrAtm(I,J,bi,bj) = maskC(I,J,kSurface,bi,bj)*(
./pkg/seaice/seaice_growth.F 2544 SEAICE_GROWTH
            DIAGarrayA(I,J) = maskC(I,J,kSurface,bi,bj)
./pkg/seaice/seaice_growth.F 2555 SEAICE_GROWTH
           DIAGarrayB(I,J) = maskC(I,J,kSurface,bi,bj) *
./pkg/seaice/seaice_growth.F 2565 SEAICE_GROWTH
           DIAGarrayA(I,J) = maskC(I,J,kSurface,bi,bj)
./pkg/seaice/seaice_growth.F 2569 SEAICE_GROWTH
           DIAGarrayC(I,J) = maskC(I,J,kSurface,bi,bj)
./pkg/seaice/seaice_growth.F 2578 SEAICE_GROWTH
           DIAGarrayB(I,J) = -maskC(I,J,kSurface,bi,bj)*rhoConstFresh
./pkg/seaice/seaice_lsr.F 1537 SEAICE_LSR_RHSU
     &       *maskC(I  ,J  ,kSrf,bi,bj)*maskC(I-1,J  ,kSrf,bi,bj)
./pkg/seaice/seaice_lsr.F 1538 SEAICE_LSR_RHSU
     &       *maskC(I  ,J-1,kSrf,bi,bj)*maskC(I-1,J-1,kSrf,bi,bj)
./pkg/seaice/seaice_lsr.F 1560 SEAICE_LSR_RHSU
     &        *maskC(I  ,J  ,kSrf,bi,bj)*maskC(I-1,J  ,kSrf,bi,bj)
./pkg/seaice/seaice_lsr.F 1561 SEAICE_LSR_RHSU
     &        *maskC(I  ,J-1,kSrf,bi,bj)*maskC(I-1,J-1,kSrf,bi,bj)
./pkg/seaice/seaice_lsr.F 1669 SEAICE_LSR_RHSV
     &       *maskC(I  ,J  ,kSrf,bi,bj)*maskC(I-1,J  ,kSrf,bi,bj)
./pkg/seaice/seaice_lsr.F 1670 SEAICE_LSR_RHSV
     &       *maskC(I  ,J-1,kSrf,bi,bj)*maskC(I-1,J-1,kSrf,bi,bj)
./pkg/seaice/seaice_lsr.F 1692 SEAICE_LSR_RHSV
     &        *maskC(I  ,J  ,kSrf,bi,bj)*maskC(I-1,J  ,kSrf,bi,bj)
./pkg/seaice/seaice_lsr.F 1693 SEAICE_LSR_RHSV
     &        *maskC(I  ,J-1,kSrf,bi,bj)*maskC(I-1,J-1,kSrf,bi,bj)
./model/src/ini_cg3d.F 148 INI_CG3D
           faceArea = _rA(i,j,bi,bj)*maskC(i,j, k ,bi,bj)
./model/src/ini_cg3d.F 149 INI_CG3D
     &                              *maskC(i,j,k-1,bi,bj)
./pkg/seaice/seaice_oceandrag_coeffs.F 114 SEAICE_OCEANDRAG_COEFFS
          CwatC(I,J,bi,bj) = CwatC(I,J,bi,bj) * maskC(I,J,kSrf,bi,bj)
./pkg/shelfice/shelfice_thermodynamics.F 407 SHELFICE_THERMODYNAMICS
            shelfIceHeatFlux(I,J,bi,bj) = maskC(I,J,K,bi,bj)
./model/src/ini_eos.F 467 EOS_CHECK
        mskSave = maskC(i,j,k,bi,bj)
./model/src/ini_eos.F 469 EOS_CHECK
        maskC(i,j,k,bi,bj) = 1.
./model/src/ini_eos.F 548 EOS_CHECK
        maskC(i,j,k,bi,bj) = mskSave
./pkg/smooth/smooth_filtervar2d.F 179 SMOOTH_FILTERVAR2D
           if (maskC(i,j,1,bi,bj).NE.0) then
./pkg/smooth/smooth_filtervar2d.F 216 SMOOTH_FILTERVAR2D
           if (maskC(i,j,1,bi,bj).NE.0) then
./pkg/smooth/smooth_filtervar3d.F 219 SMOOTH_FILTERVAR3D
           if (maskC(i,j,k,bi,bj).NE.0) then
./pkg/smooth/smooth_filtervar3d.F 260 SMOOTH_FILTERVAR3D
           if (maskC(i,j,k,bi,bj).NE.0) then
./pkg/smooth/smooth_rhs.F 104 SMOOTH_RHS
     &        maskC(i,j,k-1,bi,bj)*maskC(i,j,k,bi,bj)
./pkg/smooth/smooth_rhs.F 125 SMOOTH_RHS
     &       ( maskC(i-1,j,k,bi,bj)*
./pkg/smooth/smooth_rhs.F 127 SMOOTH_RHS
     &       +maskC( i ,j,k,bi,bj)*
./pkg/smooth/smooth_rhs.F 131 SMOOTH_RHS
     &       ( maskC(i-1,j,MIN(k+1,Nr),bi,bj)*
./pkg/smooth/smooth_rhs.F 133 SMOOTH_RHS
     &       +maskC( i ,j,MIN(k+1,Nr),bi,bj)*
./pkg/smooth/smooth_rhs.F 193 SMOOTH_RHS
     &       ( maskC(i,j-1,k,bi,bj)*
./pkg/smooth/smooth_rhs.F 195 SMOOTH_RHS
     &         +maskC(i, j ,k,bi,bj)*
./pkg/smooth/smooth_rhs.F 199 SMOOTH_RHS
     &       ( maskC(i,j-1,MIN(k+1,Nr),bi,bj)*
./pkg/smooth/smooth_rhs.F 201 SMOOTH_RHS
     &         +maskC(i, j ,MIN(k+1,Nr),bi,bj)*
./model/src/ini_masks_etc.F 482 INI_MASKS_ETC
            maskC(i,j,k,bi,bj) = 1.
./model/src/ini_masks_etc.F 485 INI_MASKS_ETC
            maskC(i,j,k,bi,bj) = 0.
./pkg/thsice/thsice_test_addfluid.h 20 NO PROCEDURE
     &                            *maskC(i,j,1,bi,bj)
./model/src/ini_pressure.F 125 INI_PRESSURE
     &             * maskC(i,j,k,bi,bj)
./model/src/ini_pressure.F 126 INI_PRESSURE
              sumTile = sumTile + maskC(i,j,k,bi,bj)
./model/src/ini_salt.F 101 INI_SALT
            IF ( maskC(i,j,k,bi,bj).NE.0.
./model/src/ini_salt.F 94 INI_SALT
            IF (maskC(i,j,k,bi,bj).EQ.0.) salt(i,j,k,bi,bj) = 0.
./model/src/ini_theta.F 102 INI_THETA
            IF (  maskC(i,j,k,bi,bj).NE.0.
./model/src/ini_theta.F 95 INI_THETA
            IF (maskC(i,j,k,bi,bj).EQ.0.) theta(i,j,k,bi,bj) = 0.
./verification/advect_xy/code/ini_salt.F 100 INI_SALT
            IF ( maskC(i,j,k,bi,bj).NE.0.
./verification/advect_xy/code/ini_salt.F 93 INI_SALT
            IF (maskC(i,j,k,bi,bj).EQ.0.) salt(i,j,k,bi,bj) = 0.
./verification/advect_xy/code/ini_theta.F 101 INI_THETA
            IF (  maskC(i,j,k,bi,bj).NE.0.
./verification/advect_xy/code/ini_theta.F 94 INI_THETA
            IF (maskC(i,j,k,bi,bj).EQ.0.) theta(i,j,k,bi,bj) = 0.
./model/src/integr_continuity.F 122 INTEGR_CONTINUITY
     &       +maskC(i,j,k,bi,bj)*( uTrans(i+1,j)-uTrans(i,j)
./model/src/integr_continuity.F 271 INTEGR_CONTINUITY
     &               +mass2rUnit*PmEpR(i,j,bi,bj)*maskC(i,j,k,bi,bj)
./model/src/integrate_for_w.F 106 INTEGRATE_FOR_W
     &         *maskC(i,j,k,bi,bj)*maskC(i,j,k-1,bi,bj)
./model/src/integrate_for_w.F 116 INTEGRATE_FOR_W
     &        )*maskC(i,j,k,bi,bj)*maskC(i,j,k-1,bi,bj)
./model/src/integrate_for_w.F 132 INTEGRATE_FOR_W
     &        )*maskC(i,j,k,bi,bj)
./model/src/integrate_for_w.F 143 INTEGRATE_FOR_W
     &        )*maskC(i,j,k,bi,bj)
./model/src/integrate_for_w.F 158 INTEGRATE_FOR_W
     &        )*maskC(i,j,k,bi,bj)
./model/src/integrate_for_w.F 168 INTEGRATE_FOR_W
     &        )*maskC(i,j,k,bi,bj)
./model/src/integrate_for_w.F 182 INTEGRATE_FOR_W
     &         *maskC(i,j,k,bi,bj)
./model/src/integrate_for_w.F 192 INTEGRATE_FOR_W
     &        )*maskC(i,j,k,bi,bj)
./verification/fizhi-gridalt-hs/code/fizhi_tendency_apply.F 138 FIZHI_TENDENCY_APPLY_T
        gT_arr(i,j) = maskC(i,j,kLev,bi,bj)
./verification/fizhi-gridalt-hs/code/fizhi_tendency_apply.F 184 FIZHI_TENDENCY_APPLY_S
        gS_arr(i,j) = maskC(i,j,kLev,bi,bj)
./verification/global_ocean.cs32x15/code_ad/cost_test.F 82 COST_TEST
         if (maskC(i,j,k,bi,bj).NE.0.) then
./verification/halfpipe_streamice/code_ad/ctrl_map_gentim2d.F 136 CTRL_MAP_GENTIM2D
     I      maskC,
./verification/halfpipe_streamice/code_ad/ctrl_map_ini_gentim2d.F 145 CTRL_MAP_INI_GENTIM2D
           call SMOOTH_CORREL2D(xx_gen,maskC,smoothOpNb,mythid)
./verification/halfpipe_streamice/code_ad/ctrl_map_ini_gentim2d.F 151 CTRL_MAP_INI_GENTIM2D
               if ((maskC(i,j,1,bi,bj).NE.0.).AND.
./verification/halfpipe_streamice/code_oad/ctrl_map_gentim2d.F 136 CTRL_MAP_GENTIM2D
     I      maskC,
./verification/hs94.128x64x5/code/external_forcing.F 207 EXTERNAL_FORCING_T
     &            *maskC(i,j,kLev,bi,bj)
./verification/hs94.1x64x5/code/apply_forcing.F 291 APPLY_FORCING_T
     &                *maskC(i,j,k,bi,bj)
./verification/hs94.1x64x5/code/ini_theta.F 101 INI_THETA
           IF (maskC(I,J,K,bi,bj).EQ.0.) theta(I,J,K,bi,bj) = 0.
./verification/hs94.1x64x5/code/ini_theta.F 107 INI_THETA
            IF (  maskC(I,J,K,bi,bj).NE.0.
./verification/hs94.1x64x5/code_ad/apply_forcing.F 291 APPLY_FORCING_T
     &                *maskC(i,j,k,bi,bj)
./verification/hs94.1x64x5/code_ad/ctrl_map_gentim2d.F 131 CTRL_MAP_GENTIM2D
     &            maskC(i,j,1,bi,bj)*rA(i,j,bi,bj)
./verification/hs94.1x64x5/code_ad/ctrl_map_gentim2d.F 145 CTRL_MAP_GENTIM2D
     &            LOCsumGlob/globalArea*maskC(i,j,1,bi,bj)
./verification/hs94.1x64x5/code_ad/ctrl_map_gentim2d.F 98 CTRL_MAP_GENTIM2D
     I      maskC,
./verification/hs94.1x64x5/code_ad/ini_theta.F 101 INI_THETA
           IF (maskC(I,J,K,bi,bj).EQ.0.) theta(I,J,K,bi,bj) = 0.
./verification/hs94.1x64x5/code_ad/ini_theta.F 107 INI_THETA
            IF (  maskC(I,J,K,bi,bj).NE.0.
./verification/hs94.1x64x5/code_oad/apply_forcing.F 291 APPLY_FORCING_T
     &                *maskC(i,j,k,bi,bj)
./verification/hs94.1x64x5/code_oad/ctrl_map_gentim2d.F 131 CTRL_MAP_GENTIM2D
     &            maskC(i,j,1,bi,bj)*rA(i,j,bi,bj)
./verification/hs94.1x64x5/code_oad/ctrl_map_gentim2d.F 145 CTRL_MAP_GENTIM2D
     &            LOCsumGlob/globalArea*maskC(i,j,1,bi,bj)
./verification/hs94.1x64x5/code_oad/ctrl_map_gentim2d.F 98 CTRL_MAP_GENTIM2D
     I      maskC,
./verification/hs94.1x64x5/code_oad/ini_theta.F 101 INI_THETA
           IF (maskC(I,J,K,bi,bj).EQ.0.) theta(I,J,K,bi,bj) = 0.
./verification/hs94.1x64x5/code_oad/ini_theta.F 107 INI_THETA
            IF (  maskC(I,J,K,bi,bj).NE.0.
./verification/hs94.cs-32x32x5/code/apply_forcing.F 291 APPLY_FORCING_T
     &                *maskC(i,j,k,bi,bj)
./verification/hs94.cs-32x32x5/code/ini_theta.F 101 INI_THETA
           IF (maskC(I,J,K,bi,bj).EQ.0.) theta(I,J,K,bi,bj) = 0.
./verification/hs94.cs-32x32x5/code/ini_theta.F 107 INI_THETA
            IF (  maskC(I,J,K,bi,bj).NE.0.
./model/src/post_cg3d.F 122 POST_CG3D
     &                            )*maskC(i,j,k,bi,bj)
./model/src/pre_cg3d.F 221 PRE_CG3D
     &       +( wVel(i,j, k ,bi,bj)*wFacKm*maskC(i,j,k-1,bi,bj)
./verification/isomip/code_ad/cost_test.F 82 COST_TEST
         if (maskC(i,j,k,bi,bj).NE.0.) then
./verification/isomip/code_oad/cost_test.F 82 COST_TEST
         if (maskC(i,j,k,bi,bj).NE.0.) then
./model/src/rotate_uv2en.F 103 ROTATE_UV2EN_RL
            tmpU(i,j) = tmpU(i,j) * maskC(i,j,kk,bi,bj)
./model/src/rotate_uv2en.F 104 ROTATE_UV2EN_RL
            tmpV(i,j) = tmpV(i,j) * maskC(i,j,kk,bi,bj)
./model/src/rotate_uv2en.F 115 ROTATE_UV2EN_RL
              tmpU(i,j) = tmpU(i,j) * maskC(i,j,kk,bi,bj)
./model/src/rotate_uv2en.F 116 ROTATE_UV2EN_RL
              tmpV(i,j) = tmpV(i,j) * maskC(i,j,kk,bi,bj)
./model/src/rotate_uv2en.F 176 ROTATE_UV2EN_RL
              uFldX(i,j,k,bi,bj)=uFldX(i,j,k,bi,bj)*maskC(i,j,kk,bi,bj)
./model/src/rotate_uv2en.F 177 ROTATE_UV2EN_RL
              vFldY(i,j,k,bi,bj)=vFldY(i,j,k,bi,bj)*maskC(i,j,kk,bi,bj)
./model/src/rotate_uv2en.F 284 ROTATE_UV2EN_RS
            tmpU(i,j) = tmpU(i,j) * maskC(i,j,kk,bi,bj)
./model/src/rotate_uv2en.F 285 ROTATE_UV2EN_RS
            tmpV(i,j) = tmpV(i,j) * maskC(i,j,kk,bi,bj)
./model/src/rotate_uv2en.F 296 ROTATE_UV2EN_RS
              tmpU(i,j) = tmpU(i,j) * maskC(i,j,kk,bi,bj)
./model/src/rotate_uv2en.F 297 ROTATE_UV2EN_RS
              tmpV(i,j) = tmpV(i,j) * maskC(i,j,kk,bi,bj)
./model/src/rotate_uv2en.F 357 ROTATE_UV2EN_RS
              uFldX(i,j,k,bi,bj)=uFldX(i,j,k,bi,bj)*maskC(i,j,kk,bi,bj)
./model/src/rotate_uv2en.F 358 ROTATE_UV2EN_RS
              vFldY(i,j,k,bi,bj)=vFldY(i,j,k,bi,bj)*maskC(i,j,kk,bi,bj)
./verification/rotating_tank/code/apply_forcing.F 595 APPLY_FORCING_T
     &                   )*maskC(i,j,km,bi,bj)*0.25 _d 0
./verification/rotating_tank/code/apply_forcing.F 608 APPLY_FORCING_T
     &                   )*maskC(i,j,kp,bi,bj)*0.25 _d 0
./verification/rotating_tank/code/apply_forcing.F 678 APPLY_FORCING_T
     &   -Qsw(i,j,bi,bj)*(swfracb(1)*maskC(i,j,k,bi,bj)
./verification/rotating_tank/code/apply_forcing.F 679 APPLY_FORCING_T
     &                   -swfracb(2)*maskC(i,j,kp1, bi,bj))
./verification/tidal_basin_2d/code/external_forcing.F 276 EXTERNAL_FORCING_T
     &   -Qsw(i,j,bi,bj)*(swfracb(1)*maskC(i,j,klev,bi,bj)
./verification/tidal_basin_2d/code/external_forcing.F 277 EXTERNAL_FORCING_T
     &                   -swfracb(2)*maskC(i,j,kp1, bi,bj))
./verification/tutorial_global_oce_biogeo/code_ad/cost_tracer.F 55 COST_TRACER
            locfc = locfc + maskC(i,j,k,bi,bj)*
./verification/tutorial_global_oce_biogeo/code_oad/cost_tracer.F 55 COST_TRACER
            locfc = locfc + maskC(i,j,k,bi,bj)*
./verification/tutorial_global_oce_latlon/code/ptracers_apply_forcing.F 92 PTRACERS_APPLY_FORCING
     &                  + 1. _d 0 * maskC(i,j,k,bi,bj)
./verification/tutorial_global_oce_optim/code_ad/cost_hflux.F 45 COST_HFLUX
             tmpC = tmpC + maskC(i,j,1,bi,bj)
./verification/tutorial_global_oce_optim/code_ad/cost_hflux.F 59 COST_HFLUX
            locfc = locfc + tmpC*maskC(i,j,1,bi,bj)*
./verification/tutorial_global_oce_optim/code_ad/cost_temp.F 50 COST_TEMP
             tmp = tmp + maskC(i,j,k,bi,bj)
./verification/tutorial_global_oce_optim/code_ad/cost_temp.F 66 COST_TEMP
              locfc = locfc + tmp*maskC(i,j,k,bi,bj)*
./verification/tutorial_global_oce_optim/code_oad/cost_hflux.F 45 COST_HFLUX
             tmpC = tmpC + maskC(i,j,1,bi,bj)
./verification/tutorial_global_oce_optim/code_oad/cost_hflux.F 59 COST_HFLUX
            locfc = locfc + tmpC*maskC(i,j,1,bi,bj)*
./verification/tutorial_global_oce_optim/code_oad/cost_temp.F 50 COST_TEMP
             tmp = tmp + maskC(i,j,k,bi,bj)
./verification/tutorial_global_oce_optim/code_oad/cost_temp.F 66 COST_TEMP
              locfc = locfc + tmp*maskC(i,j,k,bi,bj)*
./verification/tutorial_held_suarez_cs/code/apply_forcing.F 291 APPLY_FORCING_T
     &                *maskC(i,j,k,bi,bj)
./model/src/timestep_wvel.F 92 TIMESTEP_WVEL
     &             *maskC(i,j,k,bi,bj)*maskC(i,j,km1,bi,bj)
./model/src/tracers_iigw_correction.F 70 TRACERS_IIGW_CORRECTION
     &                  - gW(i,j,kp1,bi,bj) )*maskC(i,j,k,bi,bj)
./model/src/update_r_star.F 102 UPDATE_R_STAR
            recip_hFacC(i,j,k,bi,bj) = maskC(i,j,k,bi,bj)
./model/src/update_r_star.F 103 UPDATE_R_STAR
     &        / ( _hFacC(i,j,k,bi,bj) + (1.-maskC(i,j,k,bi,bj)) )
./model/src/update_r_star.F 109 UPDATE_R_STAR
           IF (maskC(i,j,k,bi,bj).NE.0.)
./model/src/update_r_star.F 67 UPDATE_R_STAR
            recip_hFacC(i,j,k,bi,bj) = maskC(i,j,k,bi,bj)
./model/src/update_r_star.F 68 UPDATE_R_STAR
     &        / ( _hFacC(i,j,k,bi,bj) + (1.-maskC(i,j,k,bi,bj)) )
./model/src/update_r_star.F 74 UPDATE_R_STAR
           IF (maskC(i,j,k,bi,bj).NE.0.)
./model/src/update_sigma.F 117 UPDATE_SIGMA
            hFacC(i,j,k,bi,bj) = maskC(i,j,k,bi,bj)
./model/src/update_sigma.F 121 UPDATE_SIGMA
            recip_hFacC(i,j,k,bi,bj) = maskC(i,j,k,bi,bj)*drF(k)
./pkg/aim_v23/aim_dyn2aim.F 168 AIM_DYN2AIM
         IF ( maskC(i,j,k,bi,bj).EQ.oneRS ) THEN
./pkg/aim_v23/aim_tendency_apply.F 133 AIM_TENDENCY_APPLY_U
     &      -( (1.-maskC(i-1,j,k-1,bi,bj))*aim_drag(i-1,j,bi,bj)
./pkg/aim_v23/aim_tendency_apply.F 134 AIM_TENDENCY_APPLY_U
     &        +(1.-maskC( i ,j,k-1,bi,bj))*aim_drag( i ,j,bi,bj)
./pkg/aim_v23/aim_tendency_apply.F 321 AIM_TENDENCY_APPLY_V
     &     -( (1.-maskC(i,j-1,k-1,bi,bj))*aim_drag(i,j-1,bi,bj)
./pkg/aim_v23/aim_tendency_apply.F 322 AIM_TENDENCY_APPLY_V
     &       +(1.-maskC(i, j ,k-1,bi,bj))*aim_drag(i, j ,bi,bj)
./pkg/aim_v23/aim_tendency_apply.F 441 AIM_TENDENCY_APPLY_T
        gT_arr(i,j) = maskC(i,j,k,bi,bj)
./pkg/aim_v23/aim_tendency_apply.F 503 AIM_TENDENCY_APPLY_S
        gS_arr(i,j) = maskC(i,j,k,bi,bj)
./pkg/atm2d/calc_1dto2d.F 63 CALC_1DTO2D
          IF (maskC(i,j,1,1,1).EQ.1.) THEN
./pkg/atm2d/calc_zonal_means.F 73 CALC_ZONAL_MEANS
         IF (maskC(i,j,1,1,1).EQ.1.) THEN
./pkg/atm2d/fixed_flux_add.F 127 FIXED_FLUX_ADD
         IF (maskC(i,j,1,1,1).EQ.1.) THEN
./pkg/atm2d/init_atm2d.F 158 INIT_ATM2D
          IF (maskC(i,j,1,1,1).EQ.1.) THEN
./pkg/atm2d/init_atm2d.F 249 INIT_ATM2D
            if ( (maskC(i,j,1,1,1).EQ.1.) .AND.
./pkg/atm2d/init_atm2d.F 250 INIT_ATM2D
     &          ( (maskC(i-1,j,1,1,1).EQ.0.).OR.
./pkg/atm2d/init_atm2d.F 251 INIT_ATM2D
     &            (maskC(i+1,j,1,1,1).EQ.0.).OR.
./pkg/atm2d/init_atm2d.F 252 INIT_ATM2D
     &            (maskC(i,j-1,1,1,1).EQ.0.).OR.
./pkg/atm2d/init_atm2d.F 253 INIT_ATM2D
     &            (maskC(i,j+1,1,1,1).EQ.0.).OR.
./pkg/atm2d/init_atm2d.F 254 INIT_ATM2D
     &            (maskC(i+1,j+1,1,1,1).EQ.0.).OR.
./pkg/atm2d/init_atm2d.F 255 INIT_ATM2D
     &            (maskC(i-1,j-1,1,1,1).EQ.0.).OR.
./pkg/atm2d/init_atm2d.F 256 INIT_ATM2D
     &            (maskC(i+1,j-1,1,1,1).EQ.0.).OR.
./pkg/atm2d/init_atm2d.F 257 INIT_ATM2D
     &            (maskC(i-1,j+1,1,1,1).EQ.0.) ) ) THEN
./pkg/atm2d/init_atm2d.F 272 INIT_ATM2D
            totrun_b(ib)=totrun_b(ib)+runoffVal(i,j)*maskC(i,j,1,1,1)
./pkg/atm2d/init_atm2d.F 278 INIT_ATM2D
            runoffVal(i,j)=runoffVal(i,j)*maskC(i,j,1,1,1)/totrun_b(ib)
./pkg/atm2d/norm_ocn_fluxes.F 42 NORM_OCN_FLUXES
          IF (maskC(i,j,1,1,1) .EQ. 1.) THEN
./pkg/atm2d/relax_add.F 100 RELAX_ADD
          IF (maskC(i,j,1,1,1) .EQ. 1.) THEN
./pkg/atm2d/sum_ocn_fluxes.F 41 SUM_OCN_FLUXES
          IF (maskC(i,j,1,1,1) .EQ. 1.) THEN
./pkg/atm2d/sum_yr_end_diags.F 60 SUM_YR_END_DIAGS
          IF (maskC(i,j,1,1,1) .EQ. 1.) THEN
./pkg/atm_phys/atm_phys_tendency_apply.F 235 ATM_PHYS_TENDENCY_APPLY_T
     &     + maskC(i,j,k,bi,bj) * atmPhys_dT(i,j,k,bi,bj)
./pkg/atm_phys/atm_phys_tendency_apply.F 290 ATM_PHYS_TENDENCY_APPLY_S
     &     + maskC(i,j,k,bi,bj) * atmPhys_dQ(i,j,k,bi,bj)
./pkg/bbl/bbl_tendency_apply.F 109 BBL_TENDENCY_APPLY_S
     &       + maskC(i,j,k,bi,bj) * bbl_TendSalt(i,j,bi,bj)
./pkg/bbl/bbl_tendency_apply.F 56 BBL_TENDENCY_APPLY_T
     &       + maskC(i,j,k,bi,bj) * bbl_TendTheta(i,j,bi,bj)
./pkg/bling/bling_airseaflux.F 118 BLING_AIRSEAFLUX
     &                          * maskC(i,j,1,bi,bj)
./pkg/bling/bling_airseaflux.F 120 BLING_AIRSEAFLUX
     &                          * maskC(i,j,1,bi,bj)
./pkg/bling/bling_airseaflux.F 123 BLING_AIRSEAFLUX
             surfsi(i,j)   = SILICA(i,j,bi,bj) * maskC(i,j,1,bi,bj)
./pkg/bling/bling_airseaflux.F 162 BLING_AIRSEAFLUX
          IF ( maskC(i,j,klev,bi,bj).NE.0. _d 0 ) THEN
./pkg/bling/bling_airseaflux.F 186 BLING_AIRSEAFLUX
          IF ( maskC(i,j,1,bi,bj).NE.0. _d 0 ) THEN
./pkg/bling/bling_airseaflux.F 249 BLING_AIRSEAFLUX
            IF (maskC(i,j,1,bi,bj).NE.0.) THEN
./pkg/bling/bling_carbonate_init.F 170 BLING_CARBONATE_INIT
     &                          * maskC(i,j,k,bi,bj)
./pkg/bling/bling_carbonate_init.F 172 BLING_CARBONATE_INIT
     &                          * maskC(i,j,k,bi,bj)
./pkg/bling/bling_carbonate_init.F 175 BLING_CARBONATE_INIT
             sitmp(i,j) =  Silica(i,j,bi,bj) * maskC(i,j,k,bi,bj)
./pkg/bling/bling_carbonate_init.F 177 BLING_CARBONATE_INIT
             sitmp(i,j) =  0.03 * maskC(i,j,k,bi,bj)
./pkg/bling/bling_carbonate_init.F 180 BLING_CARBONATE_INIT
     &                          * maskC(i,j,k,bi,bj)
./pkg/bling/bling_carbonate_init.F 196 BLING_CARBONATE_INIT
            IF ( maskC(i,j,k,bi,bj) .NE. 0. _d 0) THEN
./pkg/bling/bling_ini_forcing.F 79 BLING_INI_FORCING
             WIND(i,j,bi,bj)   = 5. _d 0*maskC(i,j,1,bi,bj)
./pkg/bling/bling_ini_forcing.F 80 BLING_INI_FORCING
             AtmosP(i,j,bi,bj) = 1. _d 0*maskC(i,j,1,bi,bj)
./pkg/bling/bling_ini_forcing.F 81 BLING_INI_FORCING
             SILICA(i,j,bi,bj) = 7.6838 _d -3*maskC(i,j,1,bi,bj)
./pkg/bling/bling_ini_forcing.F 84 BLING_INI_FORCING
             InputFe(i,j,bi,bj)= 1. _d -11*maskC(i,j,1,bi,bj)
./pkg/bling/bling_init_varia.F 142 BLING_INIT_VARIA
            IF (maskC(i,j,k,bi,bj).EQ.0.) THEN
./pkg/bling/bling_light.F 155 BLING_LIGHT
     &                 -parfrac*Qsw(i,j,bi,bj)*maskC(i,j,1,bi,bj))
./pkg/bling/bling_min_val.F 35 BLING_MIN_VAL
           PTR(i,j,k) = max( min_val*maskC(i,j,k,bi,bj), PTR(i,j,k) )
./pkg/bling/bling_production.F 428 BLING_PROD
     &           * maskC(i,j,k,bi,bj)
./pkg/bling/bling_production.F 433 BLING_PROD
     &           * maskC(i,j,k,bi,bj)
./pkg/bling/bling_production.F 437 BLING_PROD
          Pc_m      = max(Pc_m     ,maskC(i,j,k,bi,bj)*1. _d -15)
./pkg/bling/bling_production.F 438 BLING_PROD
          Pc_m_diaz = max(Pc_m_diaz,maskC(i,j,k,bi,bj)*1. _d -15)
./pkg/bling/bling_production.F 713 BLING_PROD
          DON_remin(i,j,k) = MAX(maskC(i,j,k,bi,bj)*gamma_DON 
./pkg/bling/bling_production.F 715 BLING_PROD
          DOP_remin(i,j,k) = MAX(maskC(i,j,k,bi,bj)*gamma_DOP
./pkg/bling/bling_production.F 718 BLING_PROD
          DON_remin(i,j,k) = maskC(i,j,k,bi,bj)*gamma_DON
./pkg/bling/bling_production.F 720 BLING_PROD
          DOP_remin(i,j,k) = maskC(i,j,k,bi,bj)*gamma_DOP
./pkg/bulk_force/bulkf_forcing.F 176 BULKF_FORCING
            IF ( maskC(i,j,ks,bi,bj).NE.0. _d 0 ) THEN
./pkg/bulk_force/bulkf_forcing.F 259 BULKF_FORCING
              Qnet(i,j,bi,bj) = Qnet(i,j,bi,bj)*maskC(i,j,ks,bi,bj)
./pkg/bulk_force/bulkf_forcing.F 83 BULKF_FORCING
           IF ( maskC(i,j,ks,bi,bj).NE.0. _d 0 ) THEN
./pkg/cfc/cfc11_surfforcing.F 70 CFC11_SURFFORCING
           IF ( maskC(i,j,k,bi,bj).NE.zeroRS ) THEN
./pkg/cfc/cfc12_surfforcing.F 70 CFC12_SURFFORCING
           IF ( maskC(i,j,k,bi,bj).NE.zeroRS ) THEN
./pkg/cfc/cfc_fields_load.F 192 CFC_FIELDS_LOAD
             locWind(i,j) = 5. _d 0*maskC(i,j,1,bi,bj)
./pkg/cfc/cfc_fields_load.F 216 CFC_FIELDS_LOAD
             ATMOSP(i,j,bi,bj) = maskC(i,j,1,bi,bj)
./pkg/cheapaml/cheapaml.F 372 CHEAPAML
     &                          )*maskC(i,j,1,bi,bj)
./pkg/cheapaml/cheapaml.F 425 CHEAPAML
     &            )*xef*maskC(i,j,1,bi,bj)
./pkg/cheapaml/cheapaml.F 429 CHEAPAML
     &            )*xef*maskC(i,j,1,bi,bj)
./pkg/cheapaml/cheapaml.F 465 CHEAPAML
     &                        )*maskC(i,j,1,bi,bj)
./pkg/cheapaml/cheapaml.F 504 CHEAPAML
     &                        /rhoa*maskC(i,j,1,bi,bj)
./pkg/cheapaml/cheapaml.F 507 CHEAPAML
     &                         )*maskC(i,j,1,bi,bj)
./pkg/cheapaml/cheapaml_init_fixed.F 131 CHEAPAML_INIT_FIXED
               IF( maskC(i,j,1,bi,bj).EQ.0. _d 0) THEN
./pkg/cheapaml/cheapaml_init_fixed.F 213 CHEAPAML_INIT_FIXED
     &                               (oneRL - maskC(i,j,1,bi,bj)) )
./pkg/cheapaml/cheapaml_lanl_flux.F 98 CHEAPAML_LANL_FLUX
             if (maskC(i,j,1,bi,bj).ne.0. _d 0) then
./pkg/cost/cost_accumulate_mean.F 63 COST_ACCUMULATE_MEAN
     &                 *maskW(i,j,k,bi,bj)*maskC(i,j,k,bi,bj)
./pkg/cost/cost_accumulate_mean.F 69 COST_ACCUMULATE_MEAN
     &                 *maskS(i,j,k,bi,bj)*maskC(i,j,k,bi,bj)
./pkg/cost/cost_atlantic_heat.F 107 COST_ATLANTIC_HEAT
     &                 *maskS(i,j,k,bi,bj)*maskC(i,j,k,bi,bj)
./pkg/cost/cost_atlantic_heat.F 110 COST_ATLANTIC_HEAT
     &            maskS(i,j,k,bi,bj)*maskC(i,j,k,bi,bj)
./pkg/cost/cost_atlantic_heat.F 156 COST_ATLANTIC_HEAT
     &                      *maskC(i,j,k,bi,bj)*maskC(i-1,j,k,bi,bj)
./pkg/cost/cost_atlantic_heat.F 161 COST_ATLANTIC_HEAT
     &                 *maskC(i,j,k,bi,bj)*maskC(i-1,j,k,bi,bj)
./pkg/cost/cost_atlantic_heat.F 163 COST_ATLANTIC_HEAT
     &                   *maskC(i,j,k,bi,bj)*maskC(i-1,j,k,bi,bj)
./pkg/cost/cost_atlantic_heat.F 165 COST_ATLANTIC_HEAT
     &                   *maskC(i,j,k,bi,bj)*maskC(i-1,j,k,bi,bj)
./pkg/cost/cost_depth.F 44 COST_DEPTH
          IF ( maskC(I,J,1,BI,BJ) .NE. 0. ) THEN
./pkg/ctrl/ctrl_cost_driver.F 115 CTRL_COST_DRIVER
     &     maskC, myThid )
./pkg/ctrl/ctrl_cost_driver.F 138 CTRL_COST_DRIVER
     &     maskC, myThid )
./pkg/ctrl/ctrl_cost_driver.F 88 CTRL_COST_DRIVER
     &     maskC, myThid )
./pkg/ctrl/ctrl_init.F 953 CTRL_INIT
          whflux    (i,j,bi,bj) = maskC(i,j,1,bi,bj)
./pkg/ctrl/ctrl_init.F 954 CTRL_INIT
          wsflux    (i,j,bi,bj) = maskC(i,j,1,bi,bj)
./pkg/ctrl/ctrl_init.F 957 CTRL_INIT
          watemp    (i,j,bi,bj) = maskC(i,j,1,bi,bj)
./pkg/ctrl/ctrl_init.F 958 CTRL_INIT
          waqh      (i,j,bi,bj) = maskC(i,j,1,bi,bj)
./pkg/ctrl/ctrl_init.F 959 CTRL_INIT
          wprecip   (i,j,bi,bj) = maskC(i,j,1,bi,bj)
./pkg/ctrl/ctrl_init.F 960 CTRL_INIT
          wswflux   (i,j,bi,bj) = maskC(i,j,1,bi,bj)
./pkg/ctrl/ctrl_init.F 961 CTRL_INIT
          wswdown   (i,j,bi,bj) = maskC(i,j,1,bi,bj)
./pkg/ctrl/ctrl_init.F 962 CTRL_INIT
          wuwind    (i,j,bi,bj) = maskC(i,j,1,bi,bj)
./pkg/ctrl/ctrl_init.F 963 CTRL_INIT
          wvwind    (i,j,bi,bj) = maskC(i,j,1,bi,bj)
./pkg/ctrl/ctrl_init.F 964 CTRL_INIT
          wlwflux   (i,j,bi,bj) = maskC(i,j,1,bi,bj)
./pkg/ctrl/ctrl_init.F 965 CTRL_INIT
          wlwdown   (i,j,bi,bj) = maskC(i,j,1,bi,bj)
./pkg/ctrl/ctrl_init.F 966 CTRL_INIT
          wevap     (i,j,bi,bj) = maskC(i,j,1,bi,bj)
./pkg/ctrl/ctrl_init.F 967 CTRL_INIT
          wsnowprecip(i,j,bi,bj) = maskC(i,j,1,bi,bj)
./pkg/ctrl/ctrl_init.F 968 CTRL_INIT
          wapressure(i,j,bi,bj) = maskC(i,j,1,bi,bj)
./pkg/ctrl/ctrl_init.F 969 CTRL_INIT
          wrunoff   (i,j,bi,bj) = maskC(i,j,1,bi,bj)
./pkg/ctrl/ctrl_init.F 970 CTRL_INIT
          wsst      (i,j,bi,bj) = maskC(i,j,1,bi,bj)
./pkg/ctrl/ctrl_init.F 971 CTRL_INIT
          wsss      (i,j,bi,bj) = maskC(i,j,1,bi,bj)
./pkg/ctrl/ctrl_init_wet.F 117 CTRL_INIT_WET
                if (maskC(i,j,k,bi,bj) .ne. 0.) then
./pkg/ctrl/ctrl_init_wet.F 603 CTRL_INIT_WET
      CALL ACTIVE_WRITE_GEN_RS( 'maskCtrlC', maskC, 'XY', Nr,
./pkg/ctrl/ctrl_map_gentim2d.F 124 CTRL_MAP_GENTIM2D
     &            maskC(i,j,1,bi,bj)*rA(i,j,bi,bj)
./pkg/ctrl/ctrl_map_gentim2d.F 138 CTRL_MAP_GENTIM2D
     &            LOCsumGlob/globalArea*maskC(i,j,1,bi,bj)
./pkg/ctrl/ctrl_map_gentim2d.F 91 CTRL_MAP_GENTIM2D
     I      maskC,
./pkg/ctrl/ctrl_map_ini_ecco.F 275 CTRL_MAP_INI_ECCO
     & xx_kapgm_dummy, boundsVec, kapGM, maskC, paramSmooth,
./pkg/ctrl/ctrl_map_ini_ecco.F 288 CTRL_MAP_INI_ECCO
     & xx_kapredi_dummy, boundsVec, kapRedi, maskC, paramSmooth,
./pkg/ctrl/ctrl_map_ini_genarr.F 271 CTRL_MAP_GENARR2D
       IF (dowc01) call SMOOTH_CORREL2D(xx_gen,maskC,numsmo,mythid)
./pkg/ctrl/ctrl_map_ini_genarr.F 272 CTRL_MAP_GENARR2D
       IF (dosmooth) call SMOOTH2D(xx_gen,maskC,numsmo,mythid)
./pkg/ctrl/ctrl_map_ini_genarr.F 282 CTRL_MAP_GENARR2D
                if ( (maskC(i,j,1,bi,bj).NE.0.).AND.
./pkg/ctrl/ctrl_map_ini_genarr.F 298 CTRL_MAP_GENARR2D
      CALL CTRL_BOUND_2D(fld,maskC,xx_genarr2d_bounds(1,iarr),myThid)
./pkg/ctrl/ctrl_map_ini_genarr.F 442 CTRL_MAP_GENARR3D
        call ECCO_CPRSRL(maskC,nr,localmask,nr,myThid)
./pkg/ctrl/ctrl_map_ini_genarr.F 456 CTRL_MAP_GENARR3D
                if ( (maskC(i,j,k,bi,bj).NE.0.).AND.
./pkg/ctrl/ctrl_map_ini_genarr.F 478 CTRL_MAP_GENARR3D
      CALL CTRL_BOUND_3D(fld,maskC,xx_genarr3d_bounds(1,iarr),myThid)
./pkg/ctrl/ctrl_map_ini_gentim2d.F 316 CTRL_MAP_INI_GENTIM2D
          IF (dowc01) call SMOOTH_CORREL2D(xx_gen,maskC,numsmo,mythid)
./pkg/ctrl/ctrl_map_ini_gentim2d.F 317 CTRL_MAP_INI_GENTIM2D
          IF (dosmooth) call SMOOTH2D(xx_gen,maskC,numsmo,mythid)
./pkg/ctrl/ctrl_map_ini_gentim2d.F 326 CTRL_MAP_INI_GENTIM2D
             if ((maskC(i,j,1,bi,bj).NE.0.).AND.
./pkg/ctrl/ctrl_map_ini_gentim2d.F 340 CTRL_MAP_INI_GENTIM2D
         CALL CTRL_BOUND_2D(xx_gen,maskC,
./pkg/diagnostics/diagnostics_interp_vert.F 175 DIAGNOSTICS_INTERP_VERT
                IF (maskC(i,j,kLev,bi,bj).NE.0.) THEN
./pkg/diagnostics/diagnostics_mnc_out.F 370 DIAGNOSTICS_MNC_OUT
                IF ( maskC(i,j,klev,bi,bj) .EQ. 0. )
./pkg/dic/alk_surfforcing.F 54 ALK_SURFFORCING
            IF (maskC(i,j,1,bi,bj).NE.0. _d 0) THEN
./pkg/dic/bio_export.F 149 BIO_EXPORT
     &        lit(i,j)/(lit(i,j)+lit0)*maskC(i,j,k,bi,bj)*
./pkg/dic/bio_export.F 86 BIO_EXPORT
        lit(i,j)=-parfrac*Qsw(i,j,bi,bj)*maskC(i,j,1,bi,bj)
./pkg/dic/car_flux.F 111 CAR_FLUX
     &          *maskC(i,j,kop1,bi,bj)
./pkg/dic/dic_atmos.F 108 DIC_ATMOS
     &                         *maskC(i,j,1,bi,bj)*dTtracerLev(1)
./pkg/dic/dic_biotic_forcing.F 276 DIC_BIOTIC_FORCING
           RDOP(i,j,k)= MAX(maskC(i,j,k,bi,bj)*KDOPRemin*PTR_DOP(i,j,k)
./pkg/dic/dic_biotic_forcing.F 279 DIC_BIOTIC_FORCING
           RDOP(i,j,k)= maskC(i,j,k,bi,bj)*KDOPRemin*PTR_DOP(i,j,k)
./pkg/dic/dic_cost.F 52 DIC_COST
     &                       *maskC(i,j,1,bi,bj)*dTtracerLev(1)
./pkg/dic/dic_ini_forcing.F 79 DIC_INI_FORCING
             WIND(i,j,bi,bj)   = 5. _d 0*maskC(i,j,1,bi,bj)
./pkg/dic/dic_ini_forcing.F 80 DIC_INI_FORCING
             AtmosP(i,j,bi,bj) = 1. _d 0*maskC(i,j,1,bi,bj)
./pkg/dic/dic_ini_forcing.F 81 DIC_INI_FORCING
             SILICA(i,j,bi,bj) = 7.6838 _d -3*maskC(i,j,1,bi,bj)
./pkg/dic/dic_ini_forcing.F 85 DIC_INI_FORCING
             PAR(i,j,bi,bj)    = 100. _d 0*maskC(i,j,1,bi,bj)
./pkg/dic/dic_ini_forcing.F 89 DIC_INI_FORCING
             CHL(i,j,bi,bj)    = 1. _d -2*maskC(i,j,1,bi,bj)
./pkg/dic/dic_ini_forcing.F 92 DIC_INI_FORCING
             InputFe(i,j,bi,bj)= 1. _d -11*maskC(i,j,1,bi,bj)
./pkg/dic/dic_surfforcing.F 101 DIC_SURFFORCING
     &                          * maskC(i,j,kLev,bi,bj)
./pkg/dic/dic_surfforcing.F 105 DIC_SURFFORCING
     &                          * maskC(i,j,kLev,bi,bj)
./pkg/dic/dic_surfforcing.F 106 DIC_SURFFORCING
             surfphos(i,j)  = 5.1225 _d -4 * maskC(i,j,kLev,bi,bj)
./pkg/dic/dic_surfforcing.F 109 DIC_SURFFORCING
             surfsi(i,j)   = SILICA(i,j,bi,bj) * maskC(i,j,kLev,bi,bj)
./pkg/dic/dic_surfforcing.F 158 DIC_SURFFORCING
          IF ( maskC(i,j,kLev,bi,bj).NE.0. _d 0 ) THEN
./pkg/dic/dic_surfforcing.F 181 DIC_SURFFORCING
          IF ( maskC(i,j,kLev,bi,bj).NE.0. _d 0 ) THEN
./pkg/dic/dic_surfforcing.F 229 DIC_SURFFORCING
            IF (maskC(i,j,kLev,bi,bj).NE.0. _d 0) THEN
./pkg/dic/dic_surfforcing.F 93 DIC_SURFFORCING
     &                          * maskC(i,j,kLev,bi,bj)
./pkg/dic/dic_surfforcing.F 96 DIC_SURFFORCING
     &                          * maskC(i,j,kLev,bi,bj)
./pkg/dic/dic_surfforcing.F 99 DIC_SURFFORCING
     &                          * maskC(i,j,kLev,bi,bj)
./pkg/dic/dic_surfforcing_init.F 176 DIC_SURFFORCING_INIT
     &                          * maskC(i,j,kLev,bi,bj)
./pkg/dic/dic_surfforcing_init.F 179 DIC_SURFFORCING_INIT
     &                          * maskC(i,j,kLev,bi,bj)
./pkg/dic/dic_surfforcing_init.F 182 DIC_SURFFORCING_INIT
     &                          * maskC(i,j,kLev,bi,bj)
./pkg/dic/dic_surfforcing_init.F 184 DIC_SURFFORCING_INIT
     &                          * maskC(i,j,kLev,bi,bj)
./pkg/dic/dic_surfforcing_init.F 188 DIC_SURFFORCING_INIT
     &                          * maskC(i,j,kLev,bi,bj)
./pkg/dic/dic_surfforcing_init.F 189 DIC_SURFFORCING_INIT
             surfphos(i,j)= 5.1225 _d -4 * maskC(i,j,kLev,bi,bj)
./pkg/dic/dic_surfforcing_init.F 192 DIC_SURFFORCING_INIT
             surfsi(i,j)   =  Silica(i,j,bi,bj) * maskC(i,j,kLev,bi,bj)
./pkg/dic/dic_surfforcing_init.F 204 DIC_SURFFORCING_INIT
     &                    * maskC(i,j,kLev,bi,bj)
./pkg/dic/dic_surfforcing_init.F 224 DIC_SURFFORCING_INIT
            IF ( maskC(i,j,kLev,bi,bj) .NE. 0. _d 0) THEN
./pkg/dic/fe_chem.F 73 FE_CHEM
          IF ( maskC(i,j,k,bi,bj).GT.0. ) THEN
./pkg/dic/o2_surfforcing.F 72 O2_SURFFORCING
            IF (maskC(i,j,k,bi,bj).NE.0.) THEN
./pkg/dic/phos_flux.F 117 PHOS_FLUX
     &          *maskC(i,j,kop1,bi,bj)
./pkg/ebm/ebm_area_t.F 79 EBM_AREA_T
     &                       + maskC(i,j,k,bi,bj)*theta(i,j,k,bi,bj)
./pkg/ebm/ebm_area_t.F 81 EBM_AREA_T
     &                       + maskC(i,j,k,bi,bj)
./pkg/ebm/ebm_area_t.F 85 EBM_AREA_T
     &                       + maskC(i,j,k,bi,bj)*theta(i,j,k,bi,bj)
./pkg/ebm/ebm_area_t.F 87 EBM_AREA_T
     &                       + maskC(i,j,k,bi,bj)
./pkg/ebm/ebm_area_t.F 91 EBM_AREA_T
     &                       + maskC(i,j,k,bi,bj)*theta(i,j,k,bi,bj)
./pkg/ebm/ebm_area_t.F 93 EBM_AREA_T
     &                       + maskC(i,j,k,bi,bj)
./pkg/ebm/ebm_area_t.F 97 EBM_AREA_T
     &                       + maskC(i,j,k,bi,bj)*theta(i,j,k,bi,bj)
./pkg/ebm/ebm_area_t.F 99 EBM_AREA_T
     &                       + maskC(i,j,k,bi,bj)
./pkg/ebm/ebm_zonalmean.F 65 EBM_ZONALMEAN
     &                         + maskC(i,j,k,bi,bj)
./pkg/ebm/ebm_zonalmean.F 67 EBM_ZONALMEAN
     &                         + maskC(i,j,k,bi,bj)*theta(i,j,k,bi,bj)
./pkg/ecco/cost_forcing.F 112 COST_FORCING
     &     maskC, mythid )
./pkg/ecco/cost_forcing.F 132 COST_FORCING
     &     maskC, mythid )
./pkg/ecco/cost_forcing.F 150 COST_FORCING
     &     maskC, mythid )
./pkg/ecco/cost_forcing.F 192 COST_FORCING
     &     maskC, mythid )
./pkg/ecco/cost_forcing.F 234 COST_FORCING
     &     maskC, mythid )
./pkg/ecco/cost_forcing.F 255 COST_FORCING
     &     maskC, mythid )
./pkg/ecco/cost_forcing.F 275 COST_FORCING
     &     maskC, mythid )
./pkg/ecco/cost_forcing.F 295 COST_FORCING
     &     maskC, mythid )
./pkg/ecco/cost_forcing.F 315 COST_FORCING
     &     maskC, mythid )
./pkg/ecco/cost_forcing.F 335 COST_FORCING
     &     maskC, mythid )
./pkg/ecco/cost_forcing.F 355 COST_FORCING
     &     maskC, mythid )
./pkg/ecco/cost_forcing.F 374 COST_FORCING
     &     maskC, mythid )
./pkg/ecco/cost_forcing.F 393 COST_FORCING
     &     maskC, mythid )
./pkg/ecco/cost_forcing.F 412 COST_FORCING
     &     maskC, mythid )
./pkg/ecco/cost_forcing.F 94 COST_FORCING
     &     maskC, mythid )
./pkg/ecco/cost_gencost_customize.F 125 COST_GENCOST_CUSTOMIZE
     &      m_eta(i,j,bi,bj)*maskC(i,j,1,bi,bj)
./pkg/ecco/cost_gencost_customize.F 134 COST_GENCOST_CUSTOMIZE
     &      THETA(i,j,1,bi,bj)*maskC(i,j,1,bi,bj)
./pkg/ecco/cost_gencost_customize.F 137 COST_GENCOST_CUSTOMIZE
     &      SALT(i,j,1,bi,bj)*maskC(i,j,1,bi,bj)
./pkg/ecco/cost_gencost_customize.F 146 COST_GENCOST_CUSTOMIZE
     &      )*maskC(i,j,1,bi,bj)
./pkg/ecco/cost_gencost_customize.F 151 COST_GENCOST_CUSTOMIZE
     &      geothermalFlux(i,j,bi,bj)*maskC(i,j,1,bi,bj)
./pkg/ecco/cost_gencost_customize.F 156 COST_GENCOST_CUSTOMIZE
     &      zontau(i,j,bi,bj)*maskC(i,j,1,bi,bj)
./pkg/ecco/cost_gencost_customize.F 159 COST_GENCOST_CUSTOMIZE
     &      mertau(i,j,bi,bj)*maskC(i,j,1,bi,bj)
./pkg/ecco/cost_gencost_customize.F 162 COST_GENCOST_CUSTOMIZE
     &      zonwind(i,j,bi,bj)*maskC(i,j,1,bi,bj)
./pkg/ecco/cost_gencost_customize.F 165 COST_GENCOST_CUSTOMIZE
     &      merwind(i,j,bi,bj)*maskC(i,j,1,bi,bj)
./pkg/ecco/cost_gencost_customize.F 169 COST_GENCOST_CUSTOMIZE
     &      atemp(i,j,bi,bj)*maskC(i,j,1,bi,bj)
./pkg/ecco/cost_gencost_customize.F 172 COST_GENCOST_CUSTOMIZE
     &      aqh(i,j,bi,bj)*maskC(i,j,1,bi,bj)
./pkg/ecco/cost_gencost_customize.F 175 COST_GENCOST_CUSTOMIZE
     &      precip(i,j,bi,bj)*maskC(i,j,1,bi,bj)
./pkg/ecco/cost_gencost_customize.F 180 COST_GENCOST_CUSTOMIZE
     &      swdown(i,j,bi,bj)*maskC(i,j,1,bi,bj)
./pkg/ecco/cost_gencost_customize.F 183 COST_GENCOST_CUSTOMIZE
     &      lwdown(i,j,bi,bj)*maskC(i,j,1,bi,bj)
./pkg/ecco/cost_gencost_customize.F 187 COST_GENCOST_CUSTOMIZE
     &      wspeed(i,j,bi,bj)*maskC(i,j,1,bi,bj)
./pkg/ecco/cost_gencost_customize.F 193 COST_GENCOST_CUSTOMIZE
     &      bottomDragFld(i,j,bi,bj)*maskC(i,j,1,bi,bj)
./pkg/ecco/cost_gencost_customize.F 200 COST_GENCOST_CUSTOMIZE
     &      area(i,j,bi,bj)*maskC(i,j,1,bi,bj) 
./pkg/ecco/cost_gencost_customize.F 203 COST_GENCOST_CUSTOMIZE
     &      theta(i,j,1,bi,bj)*maskC(i,j,1,bi,bj)
./pkg/ecco/cost_gencost_customize.F 207 COST_GENCOST_CUSTOMIZE
     &      heff(i,j,bi,bj)*maskC(i,j,1,bi,bj)
./pkg/ecco/cost_gencost_customize.F 210 COST_GENCOST_CUSTOMIZE
     &      hsnow(i,j,bi,bj)*maskC(i,j,1,bi,bj)
./pkg/ecco/cost_gencost_customize.F 217 COST_GENCOST_CUSTOMIZE
     &       theta(i,j,k2,bi,bj)*maskC(i,j,k2,bi,bj)
./pkg/ecco/cost_gencost_customize.F 223 COST_GENCOST_CUSTOMIZE
     &       salt(i,j,k2,bi,bj)*maskC(i,j,k2,bi,bj)
./pkg/ecco/cost_gencost_customize.F 230 COST_GENCOST_CUSTOMIZE
     &       pTracer(i,j,k2,bi,bj,itr)*maskC(i,j,k2,bi,bj)
./pkg/ecco/cost_gencost_customize.F 237 COST_GENCOST_CUSTOMIZE
     &       m_UE(i,j,k2,bi,bj)*maskC(i,j,k2,bi,bj)
./pkg/ecco/cost_gencost_customize.F 243 COST_GENCOST_CUSTOMIZE
     &       m_VN(i,j,k2,bi,bj)*maskC(i,j,k2,bi,bj)
./pkg/ecco/cost_gencost_customize.F 265 COST_GENCOST_CUSTOMIZE
     &       diffkr(i,j,k2,bi,bj)*maskC(i,j,k2,bi,bj)
./pkg/ecco/cost_gencost_customize.F 274 COST_GENCOST_CUSTOMIZE
     &       kapgm(i,j,k2,bi,bj)*maskC(i,j,k2,bi,bj)
./pkg/ecco/cost_gencost_customize.F 282 COST_GENCOST_CUSTOMIZE
     &       kapredi(i,j,k2,bi,bj)*maskC(i,j,k2,bi,bj)
./pkg/ecco/cost_gencost_seaicev4.F 190 COST_GENCOST_SEAICEV4
      call ECCO_CPRSRL(maskC,Nr,localmask,Nr,myThid)
./pkg/ecco/cost_gencost_sshv4.F 351 COST_GENCOST_SSHV4
              if ( (maskC(i,j,1,bi,bj) .eq. 0.).OR.
./pkg/ecco/cost_generic.F 315 COST_GENLOOP
        call ECCO_CPRSRL(maskC,nr,localmask,nr,myThid)
./pkg/ecco/cost_hyd.F 287 COST_HYD
     &     num_temp0, objf_temp0, maskC, myThid )
./pkg/ecco/cost_hyd.F 298 COST_HYD
     &     num_salt0, objf_salt0, maskC, myThid )
./pkg/ecco/cost_hyd.F 371 COST_HYD
     &     maskC, myThid )
./pkg/ecco/cost_hyd.F 411 COST_HYD
     &     maskC, wsmrarea,
./pkg/ecco/cost_internal_params.F 100 COST_INTERNAL_PARAMS
     &     num_diffkr, objf_diffkr, maskC, myThid )
./pkg/ecco/cost_internal_params.F 117 COST_INTERNAL_PARAMS
     &     maskC, myThid )
./pkg/ecco/cost_internal_params.F 65 COST_INTERNAL_PARAMS
     &     num_zero_mean, objf_eddytau, maskC, myThid )
./pkg/ecco/cost_internal_params.F 68 COST_INTERNAL_PARAMS
     &     num_zero_mean, objf_eddytau, maskC, myThid )
./pkg/ecco/cost_internal_params.F 80 COST_INTERNAL_PARAMS
     &     num_kapgm, objf_kapgm, maskC, myThid )
./pkg/ecco/cost_internal_params.F 90 COST_INTERNAL_PARAMS
     &     num_kapredi, objf_kapredi, maskC, myThid )
./pkg/ecco/cost_ssh.F 150 COST_SSH
              mdtmask(i,j,bi,bj)=maskC(i,j,1,bi,bj)