Symbol Description Uses
ks surface level index 461


File Line number Procedure Code
./pkg/exf/exf_filter_rl.F 36 EXF_FILTER_RL 
      INTEGER i,j, ks
./pkg/exf/exf_filter_rl.F 43 EXF_FILTER_RL 
        ks = 1
./pkg/exf/exf_filter_rl.F 44 EXF_FILTER_RL 
        IF ( usingPCoords ) ks = Nr
./pkg/exf/exf_filter_rl.F 64 EXF_FILTER_RL 
               IF ( maskC(i,j,ks,bi,bj) .EQ. 0. ) THEN
./pkg/exf/exf_filter_rl.F 72 EXF_FILTER_RL 
               IF ( maskW(i,j,ks,bi,bj) .EQ. 0. ) THEN
./pkg/exf/exf_filter_rl.F 80 EXF_FILTER_RL 
               IF ( maskS(i,j,ks,bi,bj) .EQ. 0. ) THEN
./pkg/exf/exf_mapfields.F 190 EXF_MAPFIELDS 
     &               *( theta(i,j,ks,bi,bj) - temp_EvPrRn )
./pkg/exf/exf_mapfields.F 200 EXF_MAPFIELDS 
     &               *( theta(i,j,ks,bi,bj) - temp_EvPrRn )
./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_mapfields.F 215 EXF_MAPFIELDS 
     &              *( theta(i,j,ks,bi,bj) - runoftemp(i,j,bi,bj) )
./pkg/exf/exf_mapfields.F 255 EXF_MAPFIELDS 
     &              *exf_half*maskW(i,j,ks,bi,bj)
./pkg/exf/exf_mapfields.F 293 EXF_MAPFIELDS 
     &              *exf_half*maskS(i,j,ks,bi,bj)
./pkg/exf/exf_mapfields.F 82 EXF_MAPFIELDS 
      INTEGER i,j,ks
./pkg/exf/exf_mapfields.F 90 EXF_MAPFIELDS 
      ks = 1
./pkg/exf/exf_radiation.F 100 EXF_RADIATION 
     &          ((theta(i,j,ks,bi,bj)+cen2kel)**4)
./pkg/exf/exf_radiation.F 50 EXF_RADIATION 
      INTEGER ks, kl
./pkg/exf/exf_radiation.F 64 EXF_RADIATION 
      ks = 1
./pkg/exf/exf_radiation.F 75 EXF_RADIATION 
            Tsf = theta(i,j,ks,bi,bj) + cen2kel
./pkg/exf/exf_radiation.F 77 EXF_RADIATION 
     &             *( theta(i,j,ks,bi,bj)-theta(i,j,kl,bi,bj) )
./pkg/gchem/gchem_surfmean.F 33 GCHEM_SURFMEAN 
      INTEGER ks
./pkg/gchem/gchem_surfmean.F 39 GCHEM_SURFMEAN 
       ks = Nr
./pkg/gchem/gchem_surfmean.F 41 GCHEM_SURFMEAN 
       ks = 1
./pkg/gchem/gchem_surfmean.F 50 GCHEM_SURFMEAN 
     &        + tracer(i,j,ks,bi,bj)*rA(i,j,bi,bj)*maskInC(i,j,bi,bj)
./pkg/layers/layers_wsurf_tr.F 54 LAYERS_WSURF_TR 
      INTEGER i,j,bi,bj,ks
./pkg/layers/layers_wsurf_tr.F 64 LAYERS_WSURF_TR 
             ks = ksurfC(i,j,bi,bj)
./pkg/layers/layers_wsurf_tr.F 65 LAYERS_WSURF_TR 
             IF (ks.LE.Nr) THEN
./pkg/layers/layers_wsurf_tr.F 67 LAYERS_WSURF_TR 
               tcor = -wVelFld(i,j,ks,bi,bj) * thetaFld(i,j,ks,bi,bj)
./pkg/layers/layers_wsurf_tr.F 68 LAYERS_WSURF_TR 
               scor = -wVelFld(i,j,ks,bi,bj) * saltFld(i,j,ks,bi,bj)
./pkg/longstep/longstep_forcing_surf.F 121 LONGSTEP_FORCING_SURF 
     &          *( pTracer(i,j,ks,bi,bj,iTrc) - PTRACERS_EvPrRn(iTrc) )
./pkg/longstep/longstep_forcing_surf.F 145 LONGSTEP_FORCING_SURF 
     &            *( PTRACERS_ref(ks,iTrc) - PTRACERS_EvPrRn(iTrc) )
./pkg/longstep/longstep_forcing_surf.F 50 LONGSTEP_FORCING_SURF 
      INTEGER iTrc, ks
./pkg/longstep/longstep_forcing_surf.F 54 LONGSTEP_FORCING_SURF 
        ks = Nr
./pkg/longstep/longstep_forcing_surf.F 56 LONGSTEP_FORCING_SURF 
        ks = 1
./pkg/longstep/longstep_forcing_surf.F 88 LONGSTEP_FORCING_SURF 
     &          *( PTRACERS_EvPrRn(iTrc) - pTracer(i,j,ks,bi,bj,iTrc) )
./model/src/calc_surf_dr.F 103 CALC_SURF_DR 
     &          hfacC(i,j,ks,bi,bj), hFactmp, etaFld(i,j,bi,bj)
./model/src/calc_surf_dr.F 121 CALC_SURF_DR 
           hFac_surfC(i,j,bi,bj) = h0FacC(i,j,ks,bi,bj)
./model/src/calc_surf_dr.F 123 CALC_SURF_DR 
     &              )*recip_drF(ks)*maskC(i,j,ks,bi,bj)
./model/src/calc_surf_dr.F 131 CALC_SURF_DR 
     &         ' i,j,k,bi,bj,Thid,Iter=',i,j,ks,bi,bj,myThid,myIter
./model/src/calc_surf_dr.F 133 CALC_SURF_DR 
     &         'hFac_n-1,hFac_n,eta =', hfacC(i,j,ks,bi,bj),
./model/src/calc_surf_dr.F 147 CALC_SURF_DR 
          ks = kSurfW(i,j,bi,bj)
./model/src/calc_surf_dr.F 148 CALC_SURF_DR 
          IF (ks.LE.Nr) THEN
./model/src/calc_surf_dr.F 158 CALC_SURF_DR 
            hFac_surfW(i,j,bi,bj) = h0FacW(i,j,ks,bi,bj)
./model/src/calc_surf_dr.F 161 CALC_SURF_DR 
     &              )*recip_drF(ks)*maskW(i,j,ks,bi,bj)
./model/src/calc_surf_dr.F 168 CALC_SURF_DR 
          ks = kSurfS(i,j,bi,bj)
./model/src/calc_surf_dr.F 169 CALC_SURF_DR 
          IF (ks.LE.Nr) THEN
./model/src/calc_surf_dr.F 179 CALC_SURF_DR 
            hFac_surfS(i,j,bi,bj) = h0FacS(i,j,ks,bi,bj)
./model/src/calc_surf_dr.F 182 CALC_SURF_DR 
     &              )*recip_drF(ks)*maskS(i,j,ks,bi,bj)
./model/src/calc_surf_dr.F 52 CALC_SURF_DR 
      INTEGER ks, numbWrite, numbWrMax
./model/src/calc_surf_dr.F 84 CALC_SURF_DR 
          ks = kSurfC(i,j,bi,bj)
./model/src/calc_surf_dr.F 85 CALC_SURF_DR 
          IF (ks.LE.Nr) THEN
./model/src/calc_surf_dr.F 90 CALC_SURF_DR 
             hFactmp = h0FacC(i,j,ks,bi,bj)
./model/src/calc_surf_dr.F 91 CALC_SURF_DR 
     &         + ( rSurftmp(i,j) - Ro_surf(i,j,bi,bj) )*recip_drF(ks)
./model/src/calc_surf_dr.F 95 CALC_SURF_DR 
     &         ' i,j,k,bi,bj,Thid,Iter=',i,j,ks,bi,bj,myThid,myIter
./model/src/calc_surf_dr.F 99 CALC_SURF_DR 
     &         ' i,j,k,bi,bj,Thid,Iter=',i,j,ks,bi,bj,myThid,myIter
./model/src/calc_wsurf_tr.F 50 CALC_WSURF_TR 
      INTEGER i,j,bi,bj,ks
./model/src/calc_wsurf_tr.F 70 CALC_WSURF_TR 
             ks = ksurfC(i,j,bi,bj)
./model/src/calc_wsurf_tr.F 71 CALC_WSURF_TR 
             IF (ks.LE.Nr) THEN
./model/src/calc_wsurf_tr.F 73 CALC_WSURF_TR 
     &           + rA(i,j,bi,bj)*wVelFld(i,j,ks,bi,bj)
./model/src/calc_wsurf_tr.F 74 CALC_WSURF_TR 
     &                          *thetaFld(i,j,ks,bi,bj)
./model/src/calc_wsurf_tr.F 76 CALC_WSURF_TR 
     &           + rA(i,j,bi,bj)*wVelFld(i,j,ks,bi,bj)
./model/src/calc_wsurf_tr.F 77 CALC_WSURF_TR 
     &                          *saltFld(i,j,ks,bi,bj)
./model/src/cg3d.F 113 CG3D 
      INTEGER ks
./model/src/cg3d.F 192 CG3D 
             ks = kSurfC(i,j,bi,bj)
./model/src/cg3d.F 195 CG3D 
     &              *rA(i,j,bi,bj)*deepFac2F(ks)
./model/src/cg3d.F 374 CG3D 
             ks = kSurfC(i,j,bi,bj)
./model/src/cg3d.F 377 CG3D 
     &              *rA(i,j,bi,bj)*deepFac2F(ks)
./model/src/cg3d_ex0.F 105 CG3D_EX0 
      INTEGER ks
./model/src/cg3d_ex0.F 185 CG3D_EX0 
             ks = kSurfC(i,j,bi,bj)
./model/src/cg3d_ex0.F 188 CG3D_EX0 
     &              *rA(i,j,bi,bj)*deepFac2F(ks)
./model/src/cg3d_ex0.F 358 CG3D_EX0 
             ks = kSurfC(i,j,bi,bj)
./model/src/cg3d_ex0.F 361 CG3D_EX0 
     &              *rA(i,j,bi,bj)*deepFac2F(ks)
./pkg/monitor/mon_ke.F 297 MON_KE 
            ks = kSurfC(i,j,bi,bj)
./pkg/monitor/mon_ke.F 299 MON_KE 
     &             *maskInC(i,j,bi,bj)*deepFac2F(ks)
./pkg/monitor/mon_ke.F 300 MON_KE 
     &             *rA(i,j,bi,bj)*deepFac2F(ks)*rhoFacF(ks)
./pkg/monitor/mon_ke.F 35 MON_KE 
      INTEGER ks, kp1
./pkg/monitor/mon_surfcor.F 108 MON_SURFCOR 
             ks = ksurfC(i,j,bi,bj)
./pkg/monitor/mon_surfcor.F 109 MON_SURFCOR 
             IF (ks.LE.Nr) THEN
./pkg/monitor/mon_surfcor.F 110 MON_SURFCOR 
              tmpVal = salt(i,j,ks,bi,bj)
./pkg/monitor/mon_surfcor.F 111 MON_SURFCOR 
     &               + salt(i,j,Nr,bi,bj)*drF(Nr)*recip_drF(ks)
./pkg/monitor/mon_surfcor.F 112 MON_SURFCOR 
     &                *hFacC(i,j,Nr,bi,bj)*_recip_hFacC(i,j,ks,bi,bj)
./pkg/monitor/mon_surfcor.F 115 MON_SURFCOR 
     &                        *wVel(i,j,ks,bi,bj)*tmpVal
./pkg/monitor/mon_surfcor.F 54 MON_SURFCOR 
      INTEGER i,j,k,ks,bi,bj
./pkg/monitor/mon_surfcor.F 85 MON_SURFCOR 
            ks = ksurfC(i,j,bi,bj)
./pkg/monitor/mon_surfcor.F 86 MON_SURFCOR 
            IF (ks.LE.Nr) THEN
./pkg/monitor/mon_surfcor.F 90 MON_SURFCOR 
     &               *wVel(i,j,ks,bi,bj)*theta(i,j,ks,bi,bj)
./pkg/monitor/mon_surfcor.F 94 MON_SURFCOR 
     &                       *wVel(i,j,ks,bi,bj)*salt(i,j,ks,bi,bj)
./pkg/monitor/mon_surfcor.F 98 MON_SURFCOR 
     &                 + tmpVal*atm_cp*((rC(ks)/atm_po)**atm_kappa)
./pkg/obcs/obcs_apply_surf_dr.F 100 OBCS_APPLY_SURF_DR 
     &             + etaFld(i,j-1,bi,bj)*recip_drF(ks)
./pkg/obcs/obcs_apply_surf_dr.F 114 OBCS_APPLY_SURF_DR 
         ks = kSurfW(i,j,bi,bj)
./pkg/obcs/obcs_apply_surf_dr.F 115 OBCS_APPLY_SURF_DR 
         IF (ks.LE.Nr) THEN
./pkg/obcs/obcs_apply_surf_dr.F 117 OBCS_APPLY_SURF_DR 
           hFactmp = h0FacW(i,j,ks,bi,bj)
./pkg/obcs/obcs_apply_surf_dr.F 118 OBCS_APPLY_SURF_DR 
     &             + OBEeta(j,bi,bj)*recip_drF(ks)
./pkg/obcs/obcs_apply_surf_dr.F 120 OBCS_APPLY_SURF_DR 
           hFactmp = h0FacW(i,j,ks,bi,bj)
./pkg/obcs/obcs_apply_surf_dr.F 121 OBCS_APPLY_SURF_DR 
     &             + etaFld(i,j,bi,bj)*recip_drF(ks)
./pkg/obcs/obcs_apply_surf_dr.F 133 OBCS_APPLY_SURF_DR 
         ks = kSurfW(i,j,bi,bj)
./pkg/obcs/obcs_apply_surf_dr.F 134 OBCS_APPLY_SURF_DR 
         IF (ks.LE.Nr) THEN
./pkg/obcs/obcs_apply_surf_dr.F 136 OBCS_APPLY_SURF_DR 
           hFactmp = h0FacW(i,j,ks,bi,bj)
./pkg/obcs/obcs_apply_surf_dr.F 137 OBCS_APPLY_SURF_DR 
     &             + OBWeta(j,bi,bj)*recip_drF(ks)
./pkg/obcs/obcs_apply_surf_dr.F 139 OBCS_APPLY_SURF_DR 
           hFactmp = h0FacW(i,j,ks,bi,bj)
./pkg/obcs/obcs_apply_surf_dr.F 140 OBCS_APPLY_SURF_DR 
     &             + etaFld(i-1,j,bi,bj)*recip_drF(ks)
./pkg/obcs/obcs_apply_surf_dr.F 55 OBCS_APPLY_SURF_DR 
      INTEGER i,j,ks
./pkg/obcs/obcs_apply_surf_dr.F 74 OBCS_APPLY_SURF_DR 
         ks = kSurfS(i,j,bi,bj)
./pkg/obcs/obcs_apply_surf_dr.F 75 OBCS_APPLY_SURF_DR 
         IF (ks.LE.Nr) THEN
./pkg/obcs/obcs_apply_surf_dr.F 77 OBCS_APPLY_SURF_DR 
           hFactmp = h0FacS(i,j,ks,bi,bj)
./pkg/obcs/obcs_apply_surf_dr.F 78 OBCS_APPLY_SURF_DR 
     &             + OBNeta(i,bi,bj)*recip_drF(ks)
./pkg/obcs/obcs_apply_surf_dr.F 80 OBCS_APPLY_SURF_DR 
           hFactmp = h0FacS(i,j,ks,bi,bj)
./pkg/obcs/obcs_apply_surf_dr.F 81 OBCS_APPLY_SURF_DR 
     &             + etaFld(i,j,bi,bj)*recip_drF(ks)
./pkg/obcs/obcs_apply_surf_dr.F 93 OBCS_APPLY_SURF_DR 
         ks = kSurfS(i,j,bi,bj)
./pkg/obcs/obcs_apply_surf_dr.F 94 OBCS_APPLY_SURF_DR 
         IF (ks.LE.Nr) THEN
./pkg/obcs/obcs_apply_surf_dr.F 96 OBCS_APPLY_SURF_DR 
           hFactmp = h0FacS(i,j,ks,bi,bj)
./pkg/obcs/obcs_apply_surf_dr.F 97 OBCS_APPLY_SURF_DR 
     &             + OBSeta(i,bi,bj)*recip_drF(ks)
./pkg/obcs/obcs_apply_surf_dr.F 99 OBCS_APPLY_SURF_DR 
           hFactmp = h0FacS(i,j,ks,bi,bj)
./model/src/diags_oceanic_surf_flux.F 135 DIAGS_OCEANIC_SURF_FLUX 
     &       + PmEpR(i,j,bi,bj)*theta(i,j,ks,bi,bj)*HeatCapacity_Cp
./model/src/diags_oceanic_surf_flux.F 168 DIAGS_OCEANIC_SURF_FLUX 
     &       + PmEpR(i,j,bi,bj)*salt(i,j,ks,bi,bj)
./model/src/diags_oceanic_surf_flux.F 50 DIAGS_OCEANIC_SURF_FLUX 
      INTEGER ks
./model/src/diags_oceanic_surf_flux.F 56 DIAGS_OCEANIC_SURF_FLUX 
        ks=Nr
./model/src/diags_oceanic_surf_flux.F 58 DIAGS_OCEANIC_SURF_FLUX 
        ks=1
./pkg/ptracers/ptracers_calc_wsurf_tr.F 47 PTRACERS_CALC_WSURF_TR 
      INTEGER i,j,bi,bj,ks,itr
./pkg/ptracers/ptracers_calc_wsurf_tr.F 66 PTRACERS_CALC_WSURF_TR 
           ks = ksurfC(i,j,bi,bj)
./pkg/ptracers/ptracers_calc_wsurf_tr.F 67 PTRACERS_CALC_WSURF_TR 
           IF (ks.LE.Nr) THEN
./pkg/ptracers/ptracers_calc_wsurf_tr.F 69 PTRACERS_CALC_WSURF_TR 
     &       + rA(i,j,bi,bj)*wVelFld(i,j,ks,bi,bj)
./pkg/ptracers/ptracers_calc_wsurf_tr.F 70 PTRACERS_CALC_WSURF_TR 
     &                      *ptracer(i,j,ks,bi,bj,iTr)
./pkg/ptracers/ptracers_forcing_surf.F 134 PTRACERS_FORCING_SURF 
     &          *( PTRACERS_EvPrRn(iTrc) - pTracer(i,j,ks,bi,bj,iTrc) )
./pkg/ptracers/ptracers_forcing_surf.F 165 PTRACERS_FORCING_SURF 
     &          *( pTracer(i,j,ks,bi,bj,iTrc) - PTRACERS_EvPrRn(iTrc) )
./pkg/ptracers/ptracers_forcing_surf.F 187 PTRACERS_FORCING_SURF 
     &            *( PTRACERS_ref(ks,iTrc) - PTRACERS_EvPrRn(iTrc) )
./pkg/ptracers/ptracers_forcing_surf.F 53 PTRACERS_FORCING_SURF 
      INTEGER iTrc, ks
./pkg/ptracers/ptracers_forcing_surf.F 59 PTRACERS_FORCING_SURF 
        ks = Nr
./pkg/ptracers/ptracers_forcing_surf.F 61 PTRACERS_FORCING_SURF 
        ks = 1
./pkg/ptracers/ptracers_forcing_surf.F 82 PTRACERS_FORCING_SURF 
        cutoff = 0.1 _d 0 *drF(ks)/PTRACERS_dTLev(ks)
./pkg/ptracers/ptracers_forcing_surf.F 88 PTRACERS_FORCING_SURF 
            tmpVar = MAX( salt(i,j,ks,bi,bj), epsil )
./model/src/external_forcing_surf.F 262 EXTERNAL_FORCING_SURF 
     &          *( temp_EvPrRn - theta(i,j,ks,bi,bj) )
./model/src/external_forcing_surf.F 273 EXTERNAL_FORCING_SURF 
     &          *( salt_EvPrRn - salt(i,j,ks,bi,bj) )
./model/src/external_forcing_surf.F 297 EXTERNAL_FORCING_SURF 
     &           *( theta(i,j,ks,bi,bj) - temp_EvPrRn )
./model/src/external_forcing_surf.F 308 EXTERNAL_FORCING_SURF 
     &           *( salt(i,j,ks,bi,bj) - salt_EvPrRn )
./model/src/external_forcing_surf.F 323 EXTERNAL_FORCING_SURF 
     &           *( tRef(ks) - temp_EvPrRn )
./model/src/external_forcing_surf.F 62 EXTERNAL_FORCING_SURF 
      INTEGER ks
./model/src/external_forcing_surf.F 73 EXTERNAL_FORCING_SURF 
       ks        = Nr
./model/src/external_forcing_surf.F 75 EXTERNAL_FORCING_SURF 
       ks        = 1
./model/src/forcing_surf_relax.F 101 FORCING_SURF_RELAX 
     &         *(theta(i,j,ks,bi,bj)-SST(i,j,bi,bj))
./model/src/forcing_surf_relax.F 102 FORCING_SURF_RELAX 
     &         *drF(ks)*_hFacC(i,j,ks,bi,bj)
./model/src/forcing_surf_relax.F 106 FORCING_SURF_RELAX 
     &         *(salt(i,j,ks,bi,bj)-SSS(i,j,bi,bj))
./model/src/forcing_surf_relax.F 107 FORCING_SURF_RELAX 
     &         *drF(ks)*_hFacC(i,j,ks,bi,bj)
./model/src/forcing_surf_relax.F 136 FORCING_SURF_RELAX 
     &                    *dBHybSigF(ks)*recip_drF(ks)
./model/src/forcing_surf_relax.F 137 FORCING_SURF_RELAX 
     &                    *recip_hFacC(i,j,ks,bi,bj)
./model/src/forcing_surf_relax.F 141 FORCING_SURF_RELAX 
     &                    *dBHybSigF(ks)*recip_drF(ks)
./model/src/forcing_surf_relax.F 142 FORCING_SURF_RELAX 
     &                    *recip_hFacC(i,j,ks,bi,bj)
./model/src/forcing_surf_relax.F 150 FORCING_SURF_RELAX 
           IF (ks.EQ.kSurfC(i,j,bi,bj)) THEN
./model/src/forcing_surf_relax.F 152 FORCING_SURF_RELAX 
     &             *_recip_hFacC(i,j,ks,bi,bj)*hFac_surfC(i,j,bi,bj)
./model/src/forcing_surf_relax.F 154 FORCING_SURF_RELAX 
     &             *_recip_hFacC(i,j,ks,bi,bj)*hFac_surfC(i,j,bi,bj)
./model/src/forcing_surf_relax.F 56 FORCING_SURF_RELAX 
      INTEGER ks
./model/src/forcing_surf_relax.F 67 FORCING_SURF_RELAX 
       ks        = Nr
./model/src/forcing_surf_relax.F 69 FORCING_SURF_RELAX 
       ks        = 1
./model/src/forcing_surf_relax.F 85 FORCING_SURF_RELAX 
     &         *(theta(i,j,ks,bi,bj)-SST(i,j,bi,bj))
./model/src/forcing_surf_relax.F 86 FORCING_SURF_RELAX 
     &         *drF(ks)*_hFacC(i,j,ks,bi,bj)
./model/src/forcing_surf_relax.F 90 FORCING_SURF_RELAX 
     &         *(salt(i,j,ks,bi,bj)-SSS(i,j,bi,bj))
./model/src/forcing_surf_relax.F 91 FORCING_SURF_RELAX 
     &         *drF(ks)*_hFacC(i,j,ks,bi,bj)
./pkg/salt_plume/salt_plume_forcing_surf.F 48 SALT_PLUME_FORCING_SURF 
      INTEGER i, j, ks
./pkg/salt_plume/salt_plume_forcing_surf.F 50 SALT_PLUME_FORCING_SURF 
         ks = Nr
./pkg/salt_plume/salt_plume_forcing_surf.F 52 SALT_PLUME_FORCING_SURF 
         ks = 1
./pkg/seaice/seaice_advdiff.F 109 SEAICE_ADVDIFF 
      ks = 1
./pkg/seaice/seaice_advdiff.F 128 SEAICE_ADVDIFF 
          xA(i,j,bi,bj) = _dyG(i,j,bi,bj)*_maskW(i,j,ks,bi,bj)
./pkg/seaice/seaice_advdiff.F 129 SEAICE_ADVDIFF 
          yA(i,j,bi,bj) = _dxG(i,j,bi,bj)*_maskS(i,j,ks,bi,bj)
./pkg/seaice/seaice_advdiff.F 74 SEAICE_ADVDIFF 
      INTEGER ks
./pkg/seaice/seaice_get_dynforcing.F 131 SEAICE_GET_DYNFORCING 
     &                            *_maskW(i,j,ks,bi,bj)
./pkg/seaice/seaice_get_dynforcing.F 133 SEAICE_GET_DYNFORCING 
     &                            *_maskS(i,j,ks,bi,bj)
./pkg/seaice/seaice_get_dynforcing.F 148 SEAICE_GET_DYNFORCING 
     &          + 0.5 _d 0 * (uVel(i,j,ks,bi,bj)+uVel(i+1,j,ks,bi,bj))
./pkg/seaice/seaice_get_dynforcing.F 151 SEAICE_GET_DYNFORCING 
     &          + 0.5 _d 0 * (vVel(i,j,ks,bi,bj)+vVel(i,j+1,ks,bi,bj))
./pkg/seaice/seaice_get_dynforcing.F 192 SEAICE_GET_DYNFORCING 
     &          +0.5 _d 0*(uVel(i,j,ks,bi,bj)+uVel(i+1,j,ks,bi,bj))
./pkg/seaice/seaice_get_dynforcing.F 196 SEAICE_GET_DYNFORCING 
     &          +0.5 _d 0*(vVel(i,j,ks,bi,bj)+vVel(i,j+1,ks,bi,bj))
./pkg/seaice/seaice_get_dynforcing.F 200 SEAICE_GET_DYNFORCING 
     &          +0.5 _d 0*(uVel(i-1,j,ks,bi,bj)+uVel(i,j,ks,bi,bj))
./pkg/seaice/seaice_get_dynforcing.F 204 SEAICE_GET_DYNFORCING 
     &          +0.5 _d 0*(vVel(i-1,j,ks,bi,bj)+vVel(i-1,j+1,ks,bi,bj))
./pkg/seaice/seaice_get_dynforcing.F 206 SEAICE_GET_DYNFORCING 
     &         )*_maskW(i,j,ks,bi,bj)
./pkg/seaice/seaice_get_dynforcing.F 211 SEAICE_GET_DYNFORCING 
     &          +0.5 _d 0*(uVel(i,j,ks,bi,bj)+uVel(i+1,j,ks,bi,bj))
./pkg/seaice/seaice_get_dynforcing.F 215 SEAICE_GET_DYNFORCING 
     &          +0.5 _d 0*(vVel(i,j,ks,bi,bj)+vVel(i,j+1,ks,bi,bj))
./pkg/seaice/seaice_get_dynforcing.F 219 SEAICE_GET_DYNFORCING 
     &          +0.5 _d 0*(uVel(i,j-1,ks,bi,bj)+uVel(i+1,j-1,ks,bi,bj))
./pkg/seaice/seaice_get_dynforcing.F 223 SEAICE_GET_DYNFORCING 
     &          +0.5 _d 0*(vVel(i,j-1,ks,bi,bj)+vVel(i,j,ks,bi,bj))
./pkg/seaice/seaice_get_dynforcing.F 225 SEAICE_GET_DYNFORCING 
     &         )*_maskS(i,j,ks,bi,bj)
./pkg/seaice/seaice_get_dynforcing.F 239 SEAICE_GET_DYNFORCING 
     &         )*_maskW(i,j,ks,bi,bj)
./pkg/seaice/seaice_get_dynforcing.F 248 SEAICE_GET_DYNFORCING 
     &         )*_maskS(i,j,ks,bi,bj)
./pkg/seaice/seaice_get_dynforcing.F 268 SEAICE_GET_DYNFORCING 
     &                     *_maskW(i,j,ks,bi,bj)
./pkg/seaice/seaice_get_dynforcing.F 270 SEAICE_GET_DYNFORCING 
     &                     *_maskS(i,j,ks,bi,bj)
./pkg/seaice/seaice_get_dynforcing.F 71 SEAICE_GET_DYNFORCING 
      INTEGER ks
./pkg/seaice/seaice_get_dynforcing.F 91 SEAICE_GET_DYNFORCING 
      ks = 1
./model/src/ini_cg2d.F 207 INI_CG2D 
           ks = ksurfC(i,j,bi,bj)
./model/src/ini_cg2d.F 211 INI_CG2D 
     &      +freeSurfFac*myNorm*recip_Bo(i,j,bi,bj)*deepFac2F(ks)
./model/src/ini_cg2d.F 45 INI_CG2D 
      INTEGER i, j, k, ks
./model/src/ini_cg3d.F 214 INI_CG3D 
           ks = 1
./model/src/ini_cg3d.F 217 INI_CG3D 
             aC3d(i,j,ks,bi,bj) = aC3d(i,j,ks,bi,bj)
./model/src/ini_cg3d.F 219 INI_CG3D 
     &          *rA(i,j,bi,bj)*deepFac2F(ks)/deltaTMom/deltaTFreeSurf
./model/src/ini_cg3d.F 227 INI_CG3D 
           ks = kSurfC(i,j,bi,bj)
./model/src/ini_cg3d.F 228 INI_CG3D 
           IF ( ks.LE.Nr ) THEN
./model/src/ini_cg3d.F 229 INI_CG3D 
             aC3d(i,j,ks,bi,bj) = aC3d(i,j,ks,bi,bj)
./model/src/ini_cg3d.F 231 INI_CG3D 
     &          *rA(i,j,bi,bj)*deepFac2F(ks)/deltaTMom/deltaTFreeSurf
./model/src/ini_cg3d.F 47 INI_CG3D 
      INTEGER i, j, k, ks
./pkg/seaice/seaice_tracer_phys.F 213 SEAICE_TRACER_PHYS 
     &      *SPsalFRAC*(salt(I,j,ks,bi,bj)-SItrFromOcean(i,j))
./pkg/seaice/seaice_tracer_phys.F 42 SEAICE_TRACER_PHYS 
      INTEGER iTr, jTh, I, J, bi, bj, ks
./pkg/seaice/seaice_tracer_phys.F 58 SEAICE_TRACER_PHYS 
      ks=1
./pkg/seaice/seaice_tracer_phys.F 78 SEAICE_TRACER_PHYS 
         SItrFromOcean(i,j)=SItrFromOceanFrac(iTr)*salt(I,j,ks,bi,bj)
./pkg/seaice/seaice_tracer_phys.F 79 SEAICE_TRACER_PHYS 
         SItrFromFlood(i,j)=SItrFromFloodFrac(iTr)*salt(I,j,ks,bi,bj)
./model/src/ini_nh_fields.F 45 INI_NH_FIELDS 
      INTEGER ks
./model/src/ini_nh_fields.F 89 INI_NH_FIELDS 
             ks = kSurfC(i,j,bi,bj)
./model/src/ini_nh_fields.F 90 INI_NH_FIELDS 
             IF ( ks.LE.Nr ) THEN
./model/src/ini_nh_fields.F 91 INI_NH_FIELDS 
              dPhiNH(i,j,bi,bj) = phi_nh(i,j,ks,bi,bj)
./pkg/thsice/thsice_get_ocean.F 109 THSICE_GET_OCEAN 
     &              ( uVel(i,j,ks,bi,bj) * uVel(i,j,ks,bi,bj)
./pkg/thsice/thsice_get_ocean.F 110 THSICE_GET_OCEAN 
     &              + uVel(i+1,j,ks,bi,bj)*uVel(i+1,j,ks,bi,bj)
./pkg/thsice/thsice_get_ocean.F 111 THSICE_GET_OCEAN 
     &              + vVel(i,j+1,ks,bi,bj)*vVel(i,j+1,ks,bi,bj)
./pkg/thsice/thsice_get_ocean.F 112 THSICE_GET_OCEAN 
     &              + vVel(i,j,ks,bi,bj) * vVel(i,j,ks,bi,bj)
./pkg/thsice/thsice_get_ocean.F 120 THSICE_GET_OCEAN 
           uRel1 = uVel( i, j,ks,bi,bj)-uIce( i, j,bi,bj)
./pkg/thsice/thsice_get_ocean.F 121 THSICE_GET_OCEAN 
           uRel2 = uVel(i+1,j,ks,bi,bj)-uIce(i+1,j,bi,bj)
./pkg/thsice/thsice_get_ocean.F 122 THSICE_GET_OCEAN 
           vRel1 = vVel(i, j, ks,bi,bj)-vIce(i, j, bi,bj)
./pkg/thsice/thsice_get_ocean.F 123 THSICE_GET_OCEAN 
           vRel2 = vVel(i,j+1,ks,bi,bj)-vIce(i,j+1,bi,bj)
./pkg/thsice/thsice_get_ocean.F 56 THSICE_GET_OCEAN 
      INTEGER ks
./pkg/thsice/thsice_get_ocean.F 62 THSICE_GET_OCEAN 
        ks = 1
./pkg/thsice/thsice_get_ocean.F 68 THSICE_GET_OCEAN 
             hOceMxL(i,j,bi,bj) = drF(ks)*h0FacC(i,j,ks,bi,bj)
./pkg/thsice/thsice_get_ocean.F 76 THSICE_GET_OCEAN 
             hOceMxL(i,j,bi,bj) = drF(ks)*hFac_surfC(i,j,bi,bj)
./pkg/thsice/thsice_get_ocean.F 78 THSICE_GET_OCEAN 
             hOceMxL(i,j,bi,bj) = drF(ks)*hFacC(i,j,ks,bi,bj)
./pkg/thsice/thsice_get_ocean.F 89 THSICE_GET_OCEAN 
             hOceMxL(i,j,bi,bj) = drF(ks)*hFacC(i,j,ks,bi,bj)
./pkg/thsice/thsice_get_ocean.F 96 THSICE_GET_OCEAN 
          tOceMxL(i,j,bi,bj) = theta(i,j,ks,bi,bj)
./pkg/thsice/thsice_get_ocean.F 97 THSICE_GET_OCEAN 
          sOceMxL(i,j,bi,bj) = salt (i,j,ks,bi,bj)
./model/src/ini_nlfs_vars.F 140 INI_NLFS_VARS 
           ks = kSurfC(i,j,bi,bj)
./model/src/ini_nlfs_vars.F 141 INI_NLFS_VARS 
           IF (ks.LE.Nr) THEN
./model/src/ini_nlfs_vars.F 142 INI_NLFS_VARS 
             Rmin_tmp = rF(ks+1)
./model/src/ini_nlfs_vars.F 143 INI_NLFS_VARS 
             IF ( ks.EQ.kSurfW(i,j,bi,bj))
./model/src/ini_nlfs_vars.F 145 INI_NLFS_VARS 
             IF ( ks.EQ.kSurfW(i+1,j,bi,bj))
./model/src/ini_nlfs_vars.F 147 INI_NLFS_VARS 
             IF ( ks.EQ.kSurfS(i,j,bi,bj))
./model/src/ini_nlfs_vars.F 149 INI_NLFS_VARS 
             IF ( ks.EQ.kSurfS(i,j+1,bi,bj))
./model/src/ini_nlfs_vars.F 153 INI_NLFS_VARS 
     &        MAX( MAX(rF(ks+1),R_low(i,j,bi,bj)) + hFacInf*drF(ks),
./model/src/ini_nlfs_vars.F 154 INI_NLFS_VARS 
     &                                Rmin_tmp + hFacInfMOM*drF(ks)
./model/src/ini_nlfs_vars.F 39 INI_NLFS_VARS 
      INTEGER k, ks
./model/src/ini_p_ground.F 107 INI_P_GROUND 
           ks = 1
./model/src/ini_p_ground.F 109 INI_P_GROUND 
             IF (Hfld(i,j,bi,bj).GE.hRef(k)) ks = k
./model/src/ini_p_ground.F 111 INI_P_GROUND 
           Po_surf = rHalf(ks) + (rHalf(ks+1)-rHalf(ks))*
./model/src/ini_p_ground.F 112 INI_P_GROUND 
     &       (Hfld(i,j,bi,bj)-hRef(ks))/(hRef(ks+1)-hRef(ks))
./model/src/ini_p_ground.F 273 INI_P_GROUND 
          ks = kSurfC(i,j,bi,bj)
./model/src/ini_p_ground.F 274 INI_P_GROUND 
          IF (ks.LE.Nr) THEN
./model/src/ini_p_ground.F 278 INI_P_GROUND 
               IF ( Pfld(i,j,bi,bj).LT.rF(k) ) ks = k
./model/src/ini_p_ground.F 281 INI_P_GROUND 
           IF ( Pfld(i,j,bi,bj).GE.rC(ks) ) THEN
./model/src/ini_p_ground.F 282 INI_P_GROUND 
            phiLoc = hRef(2*ks)
./model/src/ini_p_ground.F 283 INI_P_GROUND 
     &       + (hRef(2*ks-1)-hRef(2*ks))
./model/src/ini_p_ground.F 284 INI_P_GROUND 
     &        *(Pfld(i,j,bi,bj)-rC(ks))/(rHalf(2*ks-1)-rHalf(2*ks))
./model/src/ini_p_ground.F 286 INI_P_GROUND 
            phiLoc = hRef(2*ks)
./model/src/ini_p_ground.F 287 INI_P_GROUND 
     &       + (hRef(2*ks+1)-hRef(2*ks))
./model/src/ini_p_ground.F 288 INI_P_GROUND 
     &        *(Pfld(i,j,bi,bj)-rC(ks))/(rHalf(2*ks+1)-rHalf(2*ks))
./model/src/ini_p_ground.F 62 INI_P_GROUND 
      INTEGER bi,bj,i,j,k, ks
./model/src/integr_continuity.F 161 INTEGR_CONTINUITY 
            ks = kSurfC(I,J,bi,bj)
./model/src/integr_continuity.F 164 INTEGR_CONTINUITY 
     &                                         *recip_deepFac2F(ks)
./model/src/integr_continuity.F 177 INTEGR_CONTINUITY 
            ks = kSurfC(i,j,bi,bj)
./model/src/integr_continuity.F 179 INTEGR_CONTINUITY 
     &                                         *recip_deepFac2F(ks)
./model/src/integr_continuity.F 241 INTEGR_CONTINUITY 
           ks = kSurfC(i,j,bi,bj)
./model/src/integr_continuity.F 243 INTEGR_CONTINUITY 
     &                     *deepFac2F(ks)*rhoFacF(ks)
./model/src/integr_continuity.F 58 INTEGR_CONTINUITY 
      INTEGER ks
./verification/fizhi-gridalt-hs/code/do_fizhi.F 153 DO_FIZHI 
      _RL kF,sigma_b,ks,ka,deg2rad,pi,atm_po,atm_kappa,termp,kv,kT
./verification/fizhi-gridalt-hs/code/do_fizhi.F 164 DO_FIZHI 
       ks=1. _d 0/(4. _d 0 *86400. _d 0)
./verification/fizhi-gridalt-hs/code/do_fizhi.F 190 DO_FIZHI 
         kT=ka+(ks-ka)
./verification/hs94.128x64x5/code/external_forcing.F 185 EXTERNAL_FORCING_T 
      _RL thetaLim,kT,ka,ks,sigma_b,term1,term2,thetaEq,termP
./verification/hs94.128x64x5/code/external_forcing.F 189 EXTERNAL_FORCING_T 
      ks=1. _d 0/(4. _d 0 *86400. _d 0)
./verification/hs94.128x64x5/code/external_forcing.F 201 EXTERNAL_FORCING_T 
         kT=ka+(ks-ka)
./verification/hs94.1x64x5/code/apply_forcing.F 248 APPLY_FORCING_T 
      _RL thetaLim, kT, ka, ks, sigma_b, term1, term2, thetaEq
./verification/hs94.1x64x5/code/apply_forcing.F 253 APPLY_FORCING_T 
      ks = 1. _d 0/(4. _d 0 *86400. _d 0)
./verification/hs94.1x64x5/code/apply_forcing.F 286 APPLY_FORCING_T 
         kT = ka+(ks-ka)
./verification/hs94.1x64x5/code_ad/apply_forcing.F 248 APPLY_FORCING_T 
      _RL thetaLim, kT, ka, ks, sigma_b, term1, term2, thetaEq
./verification/hs94.1x64x5/code_ad/apply_forcing.F 253 APPLY_FORCING_T 
      ks = 1. _d 0/(4. _d 0 *86400. _d 0)
./verification/hs94.1x64x5/code_ad/apply_forcing.F 286 APPLY_FORCING_T 
         kT = ka+(ks-ka)
./verification/hs94.1x64x5/code_oad/apply_forcing.F 248 APPLY_FORCING_T 
      _RL thetaLim, kT, ka, ks, sigma_b, term1, term2, thetaEq
./verification/hs94.1x64x5/code_oad/apply_forcing.F 253 APPLY_FORCING_T 
      ks = 1. _d 0/(4. _d 0 *86400. _d 0)
./verification/hs94.1x64x5/code_oad/apply_forcing.F 286 APPLY_FORCING_T 
         kT = ka+(ks-ka)
./verification/hs94.cs-32x32x5/code/apply_forcing.F 248 APPLY_FORCING_T 
      _RL thetaLim, kT, ka, ks, sigma_b, term1, term2, thetaEq
./verification/hs94.cs-32x32x5/code/apply_forcing.F 253 APPLY_FORCING_T 
      ks = 1. _d 0/(4. _d 0 *86400. _d 0)
./verification/hs94.cs-32x32x5/code/apply_forcing.F 286 APPLY_FORCING_T 
         kT = ka+(ks-ka)
./model/src/post_cg3d.F 51 POST_CG3D 
      INTEGER ks
./model/src/post_cg3d.F 68 POST_CG3D 
            ks = 1
./model/src/post_cg3d.F 71 POST_CG3D 
             dPhiNH(i,j,bi,bj) = ( phi_nh(i,j,ks,bi,bj)
./model/src/post_cg3d.F 75 POST_CG3D 
     &                     *( wVel(i,j,ks,bi,bj) - dPhiNH(i,j,bi,bj) )
./model/src/post_cg3d.F 91 POST_CG3D 
            ks = kSurfC(i,j,bi,bj)
./model/src/post_cg3d.F 92 POST_CG3D 
            IF ( ks.LE.Nr ) THEN
./model/src/post_cg3d.F 93 POST_CG3D 
             dPhiNH(i,j,bi,bj) = phi_nh(i,j,ks,bi,bj)
./model/src/pre_cg3d.F 100 PRE_CG3D 
              cg3d_b(i,j,ks,bi,bj) = cg3d_b(i,j,ks,bi,bj)
./model/src/pre_cg3d.F 125 PRE_CG3D 
              ks = 1
./model/src/pre_cg3d.F 128 PRE_CG3D 
               cg3d_b(i,j,ks,bi,bj) = cg3d_b(i,j,ks,bi,bj)
./model/src/pre_cg3d.F 130 PRE_CG3D 
     &                               + locGamma*wVel(i,j,ks,bi,bj) )
./model/src/pre_cg3d.F 132 PRE_CG3D 
     &                              *_rA(i,j,bi,bj)*deepFac2F(ks)
./model/src/pre_cg3d.F 142 PRE_CG3D 
              ks = kSurfC(i,j,bi,bj)
./model/src/pre_cg3d.F 143 PRE_CG3D 
              IF ( ks.LE.Nr ) THEN
./model/src/pre_cg3d.F 144 PRE_CG3D 
               cg3d_b(i,j,ks,bi,bj) = cg3d_b(i,j,ks,bi,bj)
./model/src/pre_cg3d.F 146 PRE_CG3D 
     &                              *_rA(i,j,bi,bj)*deepFac2F(ks)
./model/src/pre_cg3d.F 157 PRE_CG3D 
              ks = kSurfC(i,j,bi,bj)
./model/src/pre_cg3d.F 158 PRE_CG3D 
              IF ( ks.LE.Nr ) THEN
./model/src/pre_cg3d.F 159 PRE_CG3D 
               cg3d_b(i,j,ks,bi,bj) = cg3d_b(i,j,ks,bi,bj)
./model/src/pre_cg3d.F 161 PRE_CG3D 
     &                  *_rA(i,j,bi,bj)*deepFac2F(ks)/deltaTMom
./model/src/pre_cg3d.F 58 PRE_CG3D 
      INTEGER ks, kp1
./model/src/pre_cg3d.F 96 PRE_CG3D 
           ks = 1
./model/src/pre_cg3d.F 97 PRE_CG3D 
           IF ( usingPCoords ) ks = Nr
./model/src/set_ref_state.F 302 SET_REF_STATE 
          ks = (k+1)/2
./model/src/set_ref_state.F 305 SET_REF_STATE 
          phiRef(k+1) = phiRef(k)+ddPI*tLoc(ks)
./model/src/set_ref_state.F 41 SET_REF_STATE 
      INTEGER k, ks, stdUnit
./model/src/solve_for_pressure.F 193 SOLVE_FOR_PRESSURE 
           ks = kSurfC(i,j,bi,bj)
./model/src/solve_for_pressure.F 194 SOLVE_FOR_PRESSURE 
           IF ( ks.LE.Nr ) THEN
./model/src/solve_for_pressure.F 196 SOLVE_FOR_PRESSURE 
     &       -freeSurfFac*_rA(i,j,bi,bj)*deepFac2F(ks)
./model/src/solve_for_pressure.F 199 SOLVE_FOR_PRESSURE 
     &           +phi_nh(i,j,ks,bi,bj)*recip_Bo(i,j,bi,bj) )
./model/src/solve_for_pressure.F 200 SOLVE_FOR_PRESSURE 
            cg3d_b(i,j,ks,bi,bj) = cg3d_b(i,j,ks,bi,bj)
./model/src/solve_for_pressure.F 201 SOLVE_FOR_PRESSURE 
     &       -freeSurfFac*_rA(i,j,bi,bj)*deepFac2F(ks)
./model/src/solve_for_pressure.F 204 SOLVE_FOR_PRESSURE 
     &           +phi_nh(i,j,ks,bi,bj)*recip_Bo(i,j,bi,bj) )
./model/src/solve_for_pressure.F 215 SOLVE_FOR_PRESSURE 
           ks = kSurfC(i,j,bi,bj)
./model/src/solve_for_pressure.F 217 SOLVE_FOR_PRESSURE 
     &       -freeSurfFac*_rA(i,j,bi,bj)*deepFac2F(ks)
./model/src/solve_for_pressure.F 225 SOLVE_FOR_PRESSURE 
           ks = kSurfC(i,j,bi,bj)
./model/src/solve_for_pressure.F 227 SOLVE_FOR_PRESSURE 
     &       -freeSurfFac*_rA(i,j,bi,bj)*deepFac2F(ks)
./model/src/solve_for_pressure.F 59 SOLVE_FOR_PRESSURE 
      INTEGER ks
./verification/tutorial_global_oce_latlon/code/ptracers_forcing_surf.F 137 PTRACERS_FORCING_SURF 
     &          *( PTRACERS_EvPrRn(iTrc) - pTracer(i,j,ks,bi,bj,iTrc) )
./verification/tutorial_global_oce_latlon/code/ptracers_forcing_surf.F 168 PTRACERS_FORCING_SURF 
     &          *( pTracer(i,j,ks,bi,bj,iTrc) - PTRACERS_EvPrRn(iTrc) )
./verification/tutorial_global_oce_latlon/code/ptracers_forcing_surf.F 190 PTRACERS_FORCING_SURF 
     &            *( PTRACERS_ref(ks,iTrc) - PTRACERS_EvPrRn(iTrc) )
./verification/tutorial_global_oce_latlon/code/ptracers_forcing_surf.F 53 PTRACERS_FORCING_SURF 
      INTEGER iTrc, ks
./verification/tutorial_global_oce_latlon/code/ptracers_forcing_surf.F 59 PTRACERS_FORCING_SURF 
        ks = Nr
./verification/tutorial_global_oce_latlon/code/ptracers_forcing_surf.F 61 PTRACERS_FORCING_SURF 
        ks = 1
./verification/tutorial_global_oce_latlon/code/ptracers_forcing_surf.F 72 PTRACERS_FORCING_SURF 
     &                  * ( 0. _d 0 - pTracer(i,j,ks,bi,bj,iTrc) )
./verification/tutorial_global_oce_latlon/code/ptracers_forcing_surf.F 73 PTRACERS_FORCING_SURF 
     &                  * drF(ks) * _hFacC(i,j,ks,bi,bj)
./verification/tutorial_global_oce_latlon/code/ptracers_forcing_surf.F 85 PTRACERS_FORCING_SURF 
        cutoff = 0.1 _d 0 *drF(ks)/PTRACERS_dTLev(ks)
./verification/tutorial_global_oce_latlon/code/ptracers_forcing_surf.F 91 PTRACERS_FORCING_SURF 
            tmpVar = MAX( salt(i,j,ks,bi,bj), epsil )
./verification/tutorial_global_oce_optim/code_ad/external_forcing_surf.F 265 EXTERNAL_FORCING_SURF 
     &          *( temp_EvPrRn - theta(i,j,ks,bi,bj) )
./verification/tutorial_global_oce_optim/code_ad/external_forcing_surf.F 276 EXTERNAL_FORCING_SURF 
     &          *( salt_EvPrRn - salt(i,j,ks,bi,bj) )
./verification/tutorial_global_oce_optim/code_ad/external_forcing_surf.F 300 EXTERNAL_FORCING_SURF 
     &           *( theta(i,j,ks,bi,bj) - temp_EvPrRn )
./verification/tutorial_global_oce_optim/code_ad/external_forcing_surf.F 311 EXTERNAL_FORCING_SURF 
     &           *( salt(i,j,ks,bi,bj) - salt_EvPrRn )
./verification/tutorial_global_oce_optim/code_ad/external_forcing_surf.F 326 EXTERNAL_FORCING_SURF 
     &           *( tRef(ks) - temp_EvPrRn )
./verification/tutorial_global_oce_optim/code_ad/external_forcing_surf.F 65 EXTERNAL_FORCING_SURF 
      INTEGER ks
./verification/tutorial_global_oce_optim/code_ad/external_forcing_surf.F 76 EXTERNAL_FORCING_SURF 
       ks        = Nr
./verification/tutorial_global_oce_optim/code_ad/external_forcing_surf.F 78 EXTERNAL_FORCING_SURF 
       ks        = 1
./verification/tutorial_global_oce_optim/code_oad/external_forcing_surf.F 265 EXTERNAL_FORCING_SURF 
     &          *( temp_EvPrRn - theta(i,j,ks,bi,bj) )
./verification/tutorial_global_oce_optim/code_oad/external_forcing_surf.F 276 EXTERNAL_FORCING_SURF 
     &          *( salt_EvPrRn - salt(i,j,ks,bi,bj) )
./verification/tutorial_global_oce_optim/code_oad/external_forcing_surf.F 300 EXTERNAL_FORCING_SURF 
     &           *( theta(i,j,ks,bi,bj) - temp_EvPrRn )
./verification/tutorial_global_oce_optim/code_oad/external_forcing_surf.F 311 EXTERNAL_FORCING_SURF 
     &           *( salt(i,j,ks,bi,bj) - salt_EvPrRn )
./verification/tutorial_global_oce_optim/code_oad/external_forcing_surf.F 326 EXTERNAL_FORCING_SURF 
     &           *( tRef(ks) - temp_EvPrRn )
./verification/tutorial_global_oce_optim/code_oad/external_forcing_surf.F 65 EXTERNAL_FORCING_SURF 
      INTEGER ks
./verification/tutorial_global_oce_optim/code_oad/external_forcing_surf.F 76 EXTERNAL_FORCING_SURF 
       ks        = Nr
./verification/tutorial_global_oce_optim/code_oad/external_forcing_surf.F 78 EXTERNAL_FORCING_SURF 
       ks        = 1
./verification/tutorial_held_suarez_cs/code/apply_forcing.F 248 APPLY_FORCING_T 
      _RL thetaLim, kT, ka, ks, sigma_b, term1, term2, thetaEq
./verification/tutorial_held_suarez_cs/code/apply_forcing.F 253 APPLY_FORCING_T 
      ks = 1. _d 0/(4. _d 0 *86400. _d 0)
./verification/tutorial_held_suarez_cs/code/apply_forcing.F 286 APPLY_FORCING_T 
         kT = ka+(ks-ka)
./verification/tutorial_tracer_adjsens/code_ad/ptracers_forcing_surf.F 131 PTRACERS_FORCING_SURF 
     &          *( PTRACERS_EvPrRn(iTrc) - pTracer(i,j,ks,bi,bj,iTrc) )
./verification/tutorial_tracer_adjsens/code_ad/ptracers_forcing_surf.F 161 PTRACERS_FORCING_SURF 
     &          *( pTracer(i,j,ks,bi,bj,iTrc) - PTRACERS_EvPrRn(iTrc) )
./verification/tutorial_tracer_adjsens/code_ad/ptracers_forcing_surf.F 182 PTRACERS_FORCING_SURF 
     &            *( PTRACERS_ref(ks,iTrc) - PTRACERS_EvPrRn(iTrc) )
./verification/tutorial_tracer_adjsens/code_ad/ptracers_forcing_surf.F 52 PTRACERS_FORCING_SURF 
      INTEGER iTrc, ks
./verification/tutorial_tracer_adjsens/code_ad/ptracers_forcing_surf.F 58 PTRACERS_FORCING_SURF 
        ks = Nr
./verification/tutorial_tracer_adjsens/code_ad/ptracers_forcing_surf.F 60 PTRACERS_FORCING_SURF 
        ks = 1
./verification/tutorial_tracer_adjsens/code_ad/ptracers_forcing_surf.F 80 PTRACERS_FORCING_SURF 
        cutoff = 0.1 _d 0 *drF(ks)/PTRACERS_dTLev(ks)
./verification/tutorial_tracer_adjsens/code_ad/ptracers_forcing_surf.F 86 PTRACERS_FORCING_SURF 
            tmpVar = MAX( salt(i,j,ks,bi,bj), epsil )
./verification/tutorial_tracer_adjsens/code_oad/ptracers_forcing_surf.F 131 PTRACERS_FORCING_SURF 
     &          *( PTRACERS_EvPrRn(iTrc) - pTracer(i,j,ks,bi,bj,iTrc) )
./verification/tutorial_tracer_adjsens/code_oad/ptracers_forcing_surf.F 161 PTRACERS_FORCING_SURF 
     &          *( pTracer(i,j,ks,bi,bj,iTrc) - PTRACERS_EvPrRn(iTrc) )
./verification/tutorial_tracer_adjsens/code_oad/ptracers_forcing_surf.F 182 PTRACERS_FORCING_SURF 
     &            *( PTRACERS_ref(ks,iTrc) - PTRACERS_EvPrRn(iTrc) )
./verification/tutorial_tracer_adjsens/code_oad/ptracers_forcing_surf.F 52 PTRACERS_FORCING_SURF 
      INTEGER iTrc, ks
./verification/tutorial_tracer_adjsens/code_oad/ptracers_forcing_surf.F 58 PTRACERS_FORCING_SURF 
        ks = Nr
./verification/tutorial_tracer_adjsens/code_oad/ptracers_forcing_surf.F 60 PTRACERS_FORCING_SURF 
        ks = 1
./verification/tutorial_tracer_adjsens/code_oad/ptracers_forcing_surf.F 80 PTRACERS_FORCING_SURF 
        cutoff = 0.1 _d 0 *drF(ks)/PTRACERS_dTLev(ks)
./verification/tutorial_tracer_adjsens/code_oad/ptracers_forcing_surf.F 86 PTRACERS_FORCING_SURF 
            tmpVar = MAX( salt(i,j,ks,bi,bj), epsil )
./model/src/update_cg2d.F 137 UPDATE_CG2D 
           ks = kSurfC(i,j,bi,bj)
./model/src/update_cg2d.F 141 UPDATE_CG2D 
     &      +freeSurfFac*cg2dNorm*recip_Bo(i,j,bi,bj)*deepFac2F(ks)
./model/src/update_cg2d.F 53 UPDATE_CG2D 
      INTEGER i, j, k, ks
./model/src/update_surf_dr.F 106 UPDATE_SURF_DR 
          ks = kSurfW(i,j,bi,bj)
./model/src/update_surf_dr.F 107 UPDATE_SURF_DR 
          IF (ks.LE.Nr) THEN
./model/src/update_surf_dr.F 108 UPDATE_SURF_DR 
           hFacW(i,j,ks,bi,bj) = hFac_surfNm1W(i,j,bi,bj)
./model/src/update_surf_dr.F 109 UPDATE_SURF_DR 
           recip_hFacW(i,j,ks,bi,bj)= 1. _d 0 / hFac_surfNm1W(i,j,bi,bj)
./model/src/update_surf_dr.F 115 UPDATE_SURF_DR 
          ks = kSurfS(i,j,bi,bj)
./model/src/update_surf_dr.F 116 UPDATE_SURF_DR 
          IF (ks.LE.Nr) THEN
./model/src/update_surf_dr.F 117 UPDATE_SURF_DR 
           hFacS(i,j,ks,bi,bj) = hFac_surfNm1S(i,j,bi,bj)
./model/src/update_surf_dr.F 118 UPDATE_SURF_DR 
           recip_hFacS(i,j,ks,bi,bj)= 1. _d 0 / hFac_surfNm1S(i,j,bi,bj)
./model/src/update_surf_dr.F 44 UPDATE_SURF_DR 
      INTEGER ks
./model/src/update_surf_dr.F 57 UPDATE_SURF_DR 
          ks = kSurfC(i,j,bi,bj)
./model/src/update_surf_dr.F 58 UPDATE_SURF_DR 
          IF (ks.LE.Nr) THEN
./model/src/update_surf_dr.F 59 UPDATE_SURF_DR 
           hFacC(i,j,ks,bi,bj) = hFac_surfC(i,j,bi,bj)
./model/src/update_surf_dr.F 60 UPDATE_SURF_DR 
           recip_hFacC(i,j,ks,bi,bj)= 1. _d 0 / hFac_surfC(i,j,bi,bj)
./model/src/update_surf_dr.F 70 UPDATE_SURF_DR 
          ks = kSurfW(i,j,bi,bj)
./model/src/update_surf_dr.F 71 UPDATE_SURF_DR 
          IF (ks.LE.Nr) THEN
./model/src/update_surf_dr.F 72 UPDATE_SURF_DR 
           hFacW(i,j,ks,bi,bj) = hFac_surfW(i,j,bi,bj)
./model/src/update_surf_dr.F 73 UPDATE_SURF_DR 
           recip_hFacW(i,j,ks,bi,bj)= 1. _d 0 / hFac_surfW(i,j,bi,bj)
./model/src/update_surf_dr.F 79 UPDATE_SURF_DR 
          ks = kSurfS(i,j,bi,bj)
./model/src/update_surf_dr.F 80 UPDATE_SURF_DR 
          IF (ks.LE.Nr) THEN
./model/src/update_surf_dr.F 81 UPDATE_SURF_DR 
           hFacS(i,j,ks,bi,bj) = hFac_surfS(i,j,bi,bj)
./model/src/update_surf_dr.F 82 UPDATE_SURF_DR 
           recip_hFacS(i,j,ks,bi,bj)= 1. _d 0 / hFac_surfS(i,j,bi,bj)
./model/src/update_surf_dr.F 93 UPDATE_SURF_DR 
          ks = kSurfC(i,j,bi,bj)
./model/src/update_surf_dr.F 94 UPDATE_SURF_DR 
          IF (ks.LE.Nr) THEN
./model/src/update_surf_dr.F 95 UPDATE_SURF_DR 
           hFacC(i,j,ks,bi,bj) = hFac_surfNm1C(i,j,bi,bj)
./model/src/update_surf_dr.F 96 UPDATE_SURF_DR 
           recip_hFacC(i,j,ks,bi,bj)= 1. _d 0 / hFac_surfNm1C(i,j,bi,bj)
./pkg/atm_compon_interf/atm_store_aim_wndstr.F 104 ATM_STORE_AIM_WNDSTR 
     &      )* 0.5 _d 0 * vVel(i,j,ks,bi,bj)
./pkg/atm_compon_interf/atm_store_aim_wndstr.F 117 ATM_STORE_AIM_WNDSTR 
          IF ( hFacW(i,j,ks,bi,bj) .NE. 0. ) THEN
./pkg/atm_compon_interf/atm_store_aim_wndstr.F 120 ATM_STORE_AIM_WNDSTR 
     &       * 0.5 _d 0 * uVel(i,j,ks,bi,bj)
./pkg/atm_compon_interf/atm_store_aim_wndstr.F 129 ATM_STORE_AIM_WNDSTR 
          IF ( hFacS(i,j,ks,bi,bj) .NE. 0. ) THEN
./pkg/atm_compon_interf/atm_store_aim_wndstr.F 132 ATM_STORE_AIM_WNDSTR 
     &       * 0.5 _d 0 * vVel(i,j,ks,bi,bj)
./pkg/atm_compon_interf/atm_store_aim_wndstr.F 66 ATM_STORE_AIM_WNDSTR 
      INTEGER i, j, ks
./pkg/atm_compon_interf/atm_store_aim_wndstr.F 73 ATM_STORE_AIM_WNDSTR 
       ks = 1
./pkg/atm_compon_interf/atm_store_aim_wndstr.F 83 ATM_STORE_AIM_WNDSTR 
          IF ( hFacW(i,j,ks,bi,bj) .NE. 0. ) THEN
./pkg/atm_compon_interf/atm_store_aim_wndstr.F 89 ATM_STORE_AIM_WNDSTR 
     &      )* 0.5 _d 0 * uVel(i,j,ks,bi,bj)
./pkg/atm_compon_interf/atm_store_aim_wndstr.F 98 ATM_STORE_AIM_WNDSTR 
          IF ( hFacS(i,j,ks,bi,bj) .NE. 0. ) THEN
./pkg/bulk_force/bulkf_forcing.F 111 BULKF_FORCING 
     I            cloud(i,j,bi,bj),theta(i,j,ks,bi,bj),
./pkg/bulk_force/bulkf_forcing.F 127 BULKF_FORCING 
     I            Qair(i,j,bi,bj), theta(i,j,ks,bi,bj),
./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 256 BULKF_FORCING 
     &              + ( theta(i,j,ks,bi,bj) - temp_EvPrRn )
./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 64 BULKF_FORCING 
      INTEGER ks, iceornot
./pkg/bulk_force/bulkf_forcing.F 73 BULKF_FORCING 
      ks = 1
./pkg/bulk_force/bulkf_forcing.F 83 BULKF_FORCING 
           IF ( maskC(i,j,ks,bi,bj).NE.0. _d 0 ) THEN
./pkg/bulk_force/bulkf_forcing.F 88 BULKF_FORCING 
     I            theta(i,j,ks,bi,bj), flwdwn(i,j,bi,bj),
./pkg/diagnostics/diagnostics_calc_phivel.F 128 DIAGNOSTICS_CALC_PHIVEL 
     &                          *qtmp1(i,j,ks,bi,bj)*maskInW(i,j,bi,bj)
./pkg/diagnostics/diagnostics_calc_phivel.F 130 DIAGNOSTICS_CALC_PHIVEL 
     &                          *qtmp2(i,j,ks,bi,bj)*maskInS(i,j,bi,bj)
./pkg/diagnostics/diagnostics_calc_phivel.F 322 DIAGNOSTICS_CALC_PHIVEL 
              qtmp1(i,j,ks,bi,bj) =  x2d(i,j,bi,bj)
./pkg/diagnostics/diagnostics_calc_phivel.F 323 DIAGNOSTICS_CALC_PHIVEL 
              qtmp2(i,j,ks,bi,bj) =  psiVel(i,j,bi,bj)
./pkg/diagnostics/diagnostics_calc_phivel.F 69 DIAGNOSTICS_CALC_PHIVEL 
      INTEGER ks
./pkg/diagnostics/diagnostics_calc_phivel.F 93 DIAGNOSTICS_CALC_PHIVEL 
      DO ks = 1,nlevels(listId)
./pkg/diagnostics/diagnostics_calc_phivel.F 94 DIAGNOSTICS_CALC_PHIVEL 
        k = NINT(levs(ks,listId))