Symbol Description Uses
ks surface level index 191



File Line number Procedure Code
./pkg/ptracers/ptracers_forcing_surf.F 48 PTRACERS_FORCING_SURF
      INTEGER iTrc, ks
./pkg/ptracers/ptracers_forcing_surf.F 52 PTRACERS_FORCING_SURF
        ks = Nr
./pkg/ptracers/ptracers_forcing_surf.F 54 PTRACERS_FORCING_SURF
        ks = 1
./pkg/ptracers/ptracers_forcing_surf.F 84 PTRACERS_FORCING_SURF
     &          *( PTRACERS_EvPrRn(iTrc) - pTracer(i,j,ks,bi,bj,iTrc) )
./model/src/ini_cg3d.F 234 INI_CG3D
           ks = ksurfC(I,J,bi,bj)
./model/src/ini_cg3d.F 235 INI_CG3D
           IF ( ks.LE.Nr ) THEN
./model/src/ini_cg3d.F 236 INI_CG3D
             aC3d(I,J,ks,bi,bj) = aC3d(I,J,ks,bi,bj)
./model/src/ini_cg3d.F 50 INI_CG3D
      INTEGER I, J, K, ks
./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 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),
./model/src/ini_surf_dr.F 39 INI_SURF_DR
      INTEGER ks
./model/src/ini_surf_dr.F 65 INI_SURF_DR
           ks = ksurfC(i,j,bi,bj)
./model/src/ini_surf_dr.F 66 INI_SURF_DR
           IF (ks.LE.Nr) THEN
./model/src/ini_surf_dr.F 67 INI_SURF_DR
             Rmin_tmp = rF(ks+1)
./model/src/ini_surf_dr.F 68 INI_SURF_DR
             IF ( ks.EQ.ksurfW(i,j,bi,bj))
./model/src/ini_surf_dr.F 70 INI_SURF_DR
             IF ( ks.EQ.ksurfW(i+1,j,bi,bj))
./model/src/ini_surf_dr.F 72 INI_SURF_DR
             IF ( ks.EQ.ksurfS(i,j,bi,bj))
./model/src/ini_surf_dr.F 74 INI_SURF_DR
             IF ( ks.EQ.ksurfS(i,j+1,bi,bj))
./model/src/ini_surf_dr.F 78 INI_SURF_DR
     &        MAX( MAX(rF(ks+1),R_low(i,j,bi,bj)) + hFacInf*drF(ks),
./model/src/ini_surf_dr.F 79 INI_SURF_DR
     &                                Rmin_tmp + hFacInfMOM*drF(ks)
./verification/fizhi-gridalt-hs/code/do_fizhi.F 152 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 163 DO_FIZHI
       ks=1. _d 0/(4. _d 0 *86400. _d 0)
./verification/fizhi-gridalt-hs/code/do_fizhi.F 189 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/global_ocean.90x40x15/code/ptracers_forcing_surf.F 48 PTRACERS_FORCING_SURF
      INTEGER iTrc, ks
./verification/global_ocean.90x40x15/code/ptracers_forcing_surf.F 52 PTRACERS_FORCING_SURF
        ks = Nr
./verification/global_ocean.90x40x15/code/ptracers_forcing_surf.F 54 PTRACERS_FORCING_SURF
        ks = 1
./verification/global_ocean.90x40x15/code/ptracers_forcing_surf.F 84 PTRACERS_FORCING_SURF
     &          *( PTRACERS_EvPrRn(iTrc) - pTracer(i,j,ks,bi,bj,iTrc) )
./verification/tutorial_tracer_adjsens/code_ad/ptracers_forcing_surf.F 48 PTRACERS_FORCING_SURF
      INTEGER iTrc, ks
./verification/tutorial_tracer_adjsens/code_ad/ptracers_forcing_surf.F 52 PTRACERS_FORCING_SURF
        ks = Nr
./verification/tutorial_tracer_adjsens/code_ad/ptracers_forcing_surf.F 54 PTRACERS_FORCING_SURF
        ks = 1
./verification/tutorial_tracer_adjsens/code_ad/ptracers_forcing_surf.F 84 PTRACERS_FORCING_SURF
     &          *( PTRACERS_EvPrRn(iTrc) - pTracer(i,j,ks,bi,bj,iTrc) )
./verification/hs94.1x64x5/code/external_forcing.F 185 EXTERNAL_FORCING_T
      _RL thetaLim,kT,ka,ks,sigma_b,term1,term2,thetaEq,termP
./verification/hs94.1x64x5/code/external_forcing.F 189 EXTERNAL_FORCING_T
      ks=1. _d 0/(4. _d 0 *86400. _d 0)
./verification/hs94.1x64x5/code/external_forcing.F 201 EXTERNAL_FORCING_T
         kT=ka+(ks-ka)
./verification/hs94.1x64x5/code_ad/external_forcing.F 185 EXTERNAL_FORCING_T
      _RL thetaLim,kT,ka,ks,sigma_b,term1,term2,thetaEq,termP
./verification/hs94.1x64x5/code_ad/external_forcing.F 189 EXTERNAL_FORCING_T
      ks=1. _d 0/(4. _d 0 *86400. _d 0)
./verification/hs94.1x64x5/code_ad/external_forcing.F 201 EXTERNAL_FORCING_T
         kT=ka+(ks-ka)
./verification/hs94.cs-32x32x5/code/external_forcing.F 185 EXTERNAL_FORCING_T
      _RL thetaLim,kT,ka,ks,sigma_b,term1,term2,thetaEq,termP
./verification/hs94.cs-32x32x5/code/external_forcing.F 189 EXTERNAL_FORCING_T
      ks=1. _d 0/(4. _d 0 *86400. _d 0)
./verification/hs94.cs-32x32x5/code/external_forcing.F 201 EXTERNAL_FORCING_T
         kT=ka+(ks-ka)
./model/src/update_surf_dr.F 43 UPDATE_SURF_DR
      INTEGER ks
./model/src/update_surf_dr.F 62 UPDATE_SURF_DR
          ks = ksurfC(i,j,bi,bj)
./model/src/update_surf_dr.F 63 UPDATE_SURF_DR
          IF (ks.LE.Nr) THEN
./model/src/update_surf_dr.F 64 UPDATE_SURF_DR
           hFacC(i,j,ks,bi,bj) = hFac_surfC(i,j,bi,bj)
./model/src/update_surf_dr.F 65 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 75 UPDATE_SURF_DR
          ks = ksurfW(i,j,bi,bj)
./model/src/update_surf_dr.F 76 UPDATE_SURF_DR
          IF (ks.LE.Nr) THEN
./model/src/update_surf_dr.F 77 UPDATE_SURF_DR
           hFacW(i,j,ks,bi,bj) = hFac_surfW(i,j,bi,bj)
./model/src/update_surf_dr.F 78 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 84 UPDATE_SURF_DR
          ks = ksurfS(i,j,bi,bj)
./model/src/update_surf_dr.F 85 UPDATE_SURF_DR
          IF (ks.LE.Nr) THEN
./model/src/update_surf_dr.F 86 UPDATE_SURF_DR
           hFacS(i,j,ks,bi,bj) = hFac_surfS(i,j,bi,bj)
./model/src/update_surf_dr.F 87 UPDATE_SURF_DR
           recip_hFacS(i,j,ks,bi,bj)= 1. _d 0 / hFac_surfS(i,j,bi,bj)
./verification/tutorial_held_suarez_cs/code/external_forcing.F 185 EXTERNAL_FORCING_T
      _RL thetaLim,kT,ka,ks,sigma_b,term1,term2,thetaEq,termP
./verification/tutorial_held_suarez_cs/code/external_forcing.F 189 EXTERNAL_FORCING_T
      ks=1. _d 0/(4. _d 0 *86400. _d 0)
./verification/tutorial_held_suarez_cs/code/external_forcing.F 201 EXTERNAL_FORCING_T
         kT=ka+(ks-ka)
./verification/tutorial_global_oce_latlon/code/ptracers_forcing_surf.F 48 PTRACERS_FORCING_SURF
      INTEGER iTrc, ks
./verification/tutorial_global_oce_latlon/code/ptracers_forcing_surf.F 52 PTRACERS_FORCING_SURF
        ks = Nr
./verification/tutorial_global_oce_latlon/code/ptracers_forcing_surf.F 54 PTRACERS_FORCING_SURF
        ks = 1
./verification/tutorial_global_oce_latlon/code/ptracers_forcing_surf.F 84 PTRACERS_FORCING_SURF
     &          *( PTRACERS_EvPrRn(iTrc) - pTracer(i,j,ks,bi,bj,iTrc) )
./model/src/external_forcing_surf.F 120 EXTERNAL_FORCING_SURF
           IF (ks.EQ.ksurfC(i,j,bi,bj)) THEN
./model/src/external_forcing_surf.F 122 EXTERNAL_FORCING_SURF
     &             *_recip_hFacC(i,j,ks,bi,bj)*hFac_surfC(i,j,bi,bj)
./model/src/external_forcing_surf.F 124 EXTERNAL_FORCING_SURF
     &             *_recip_hFacC(i,j,ks,bi,bj)*hFac_surfC(i,j,bi,bj)
./model/src/external_forcing_surf.F 227 EXTERNAL_FORCING_SURF
     &         *( temp_EvPrRn - theta(i,j,ks,bi,bj) )
./model/src/external_forcing_surf.F 238 EXTERNAL_FORCING_SURF
     &         *( salt_EvPrRn - salt(i,j,ks,bi,bj) )
./model/src/external_forcing_surf.F 258 EXTERNAL_FORCING_SURF
     &      + EmPmR(i,j,bi,bj)*salt(i,j,ks,bi,bj)
./model/src/external_forcing_surf.F 55 EXTERNAL_FORCING_SURF
      INTEGER ks
./model/src/external_forcing_surf.F 62 EXTERNAL_FORCING_SURF
       ks        = Nr 
./model/src/external_forcing_surf.F 64 EXTERNAL_FORCING_SURF
       ks        = 1
./model/src/external_forcing_surf.F 79 EXTERNAL_FORCING_SURF
     &         *(theta(i,j,ks,bi,bj)-SST(i,j,bi,bj))
./model/src/external_forcing_surf.F 80 EXTERNAL_FORCING_SURF
     &         *drF(ks)*_hFacC(i,j,ks,bi,bj)
./model/src/external_forcing_surf.F 84 EXTERNAL_FORCING_SURF
     &         *(salt(i,j,ks,bi,bj)-SSS(i,j,bi,bj))
./model/src/external_forcing_surf.F 85 EXTERNAL_FORCING_SURF
     &         *drF(ks)*_hFacC(i,j,ks,bi,bj)
./model/src/external_forcing_surf.F 90 EXTERNAL_FORCING_SURF
     &         *(theta(i,j,ks,bi,bj)-SST(i,j,bi,bj))
./model/src/external_forcing_surf.F 91 EXTERNAL_FORCING_SURF
     &         *drF(ks)*_hFacC(i,j,ks,bi,bj)
./model/src/external_forcing_surf.F 95 EXTERNAL_FORCING_SURF
     &         *(salt(i,j,ks,bi,bj)-SSS(i,j,bi,bj))
./model/src/external_forcing_surf.F 96 EXTERNAL_FORCING_SURF
     &         *drF(ks)*_hFacC(i,j,ks,bi,bj)
./model/src/solve_for_pressure.F 192 SOLVE_FOR_PRESSURE
           ks = ksurfC(i,j,bi,bj)
./model/src/solve_for_pressure.F 193 SOLVE_FOR_PRESSURE
           IF ( ks.LE.Nr ) THEN
./model/src/solve_for_pressure.F 197 SOLVE_FOR_PRESSURE
            cg3d_b(i,j,ks,bi,bj) = cg3d_b(i,j,ks,bi,bj)
./model/src/solve_for_pressure.F 206 SOLVE_FOR_PRESSURE
           ks = ksurfC(i,j,bi,bj)
./model/src/solve_for_pressure.F 207 SOLVE_FOR_PRESSURE
           IF ( ks.LE.Nr ) THEN
./model/src/solve_for_pressure.F 211 SOLVE_FOR_PRESSURE
     &           +phi_nh(i,j,ks,bi,bj)*horiVertRatio/gravity )
./model/src/solve_for_pressure.F 212 SOLVE_FOR_PRESSURE
            cg3d_b(i,j,ks,bi,bj) = cg3d_b(i,j,ks,bi,bj)
./model/src/solve_for_pressure.F 215 SOLVE_FOR_PRESSURE
     &           +phi_nh(i,j,ks,bi,bj)*horiVertRatio/gravity )
./model/src/solve_for_pressure.F 389 SOLVE_FOR_PRESSURE
              ks = ksurfC(i,j,bi,bj)
./model/src/solve_for_pressure.F 390 SOLVE_FOR_PRESSURE
              IF ( ks.LE.Nr ) THEN
./model/src/solve_for_pressure.F 391 SOLVE_FOR_PRESSURE
               cg3d_b(i,j,ks,bi,bj) = cg3d_b(i,j,ks,bi,bj)
./model/src/solve_for_pressure.F 517 SOLVE_FOR_PRESSURE
            ks = ksurfC(i,j,bi,bj)
./model/src/solve_for_pressure.F 518 SOLVE_FOR_PRESSURE
            IF ( ks.LE.Nr ) THEN
./model/src/solve_for_pressure.F 520 SOLVE_FOR_PRESSURE
     &                       *(cg2d_x(i,j,bi,bj) + phi_nh(i,j,ks,bi,bj))
./model/src/solve_for_pressure.F 523 SOLVE_FOR_PRESSURE
     &                            - phi_nh(i,j,ks,bi,bj)
./model/src/solve_for_pressure.F 67 SOLVE_FOR_PRESSURE
      INTEGER ks, kp1
./model/src/ini_p_ground.F 272 INI_P_GROUND
          ks = ksurfC(i,j,bi,bj)
./model/src/ini_p_ground.F 273 INI_P_GROUND
          IF (ks.LE.Nr) THEN
./model/src/ini_p_ground.F 274 INI_P_GROUND
           IF ( Pfld(i,j,bi,bj).GE.rC(ks) ) THEN
./model/src/ini_p_ground.F 275 INI_P_GROUND
            phiLoc = hRef(2*ks)
./model/src/ini_p_ground.F 276 INI_P_GROUND
     &       + (hRef(2*ks-1)-hRef(2*ks))
./model/src/ini_p_ground.F 277 INI_P_GROUND
     &        *(Pfld(i,j,bi,bj)-rC(ks))/(rHalf(2*ks-1)-rHalf(2*ks))
./model/src/ini_p_ground.F 279 INI_P_GROUND
            phiLoc = hRef(2*ks)
./model/src/ini_p_ground.F 280 INI_P_GROUND
     &       + (hRef(2*ks+1)-hRef(2*ks))
./model/src/ini_p_ground.F 281 INI_P_GROUND
     &        *(Pfld(i,j,bi,bj)-rC(ks))/(rHalf(2*ks+1)-rHalf(2*ks))
./model/src/calc_surf_dr.F 114 CALC_SURF_DR
     &         ( rSurftmp(i,j) - MAX(rF(ks+1), R_low(i,j,bi,bj))
./model/src/calc_surf_dr.F 115 CALC_SURF_DR
     &         )*recip_drF(ks)*maskC(i,j,ks,bi,bj)
./model/src/calc_surf_dr.F 123 CALC_SURF_DR
     &         ' i,j,k,bi,bj,Thid,Iter=',i,j,ks,bi,bj,myThid,myIter
./model/src/calc_surf_dr.F 125 CALC_SURF_DR
     &         'hFac_n-1,hFac_n,eta =', hfacC(i,j,ks,bi,bj),
./model/src/calc_surf_dr.F 139 CALC_SURF_DR
          ks = ksurfW(i,j,bi,bj)
./model/src/calc_surf_dr.F 140 CALC_SURF_DR
          IF (ks.LE.Nr) THEN
./model/src/calc_surf_dr.F 141 CALC_SURF_DR
            hhm = rF(ks)
./model/src/calc_surf_dr.F 142 CALC_SURF_DR
            IF(ks.EQ.ksurfC(i-1,j,bi,bj)) hhm = rSurftmp(i-1,j)
./model/src/calc_surf_dr.F 143 CALC_SURF_DR
            hhp = rF(ks)
./model/src/calc_surf_dr.F 144 CALC_SURF_DR
            IF(ks.EQ.ksurfC(i,j,bi,bj))   hhp = rSurftmp(i,j)  
./model/src/calc_surf_dr.F 147 CALC_SURF_DR
     &          - MAX(rF(ks+1),R_low(i-1,j,bi,bj),R_low(i,j,bi,bj)) 
./model/src/calc_surf_dr.F 148 CALC_SURF_DR
     &         )*recip_drF(ks)*maskW(i,j,ks,bi,bj)
./model/src/calc_surf_dr.F 155 CALC_SURF_DR
          ks = ksurfS(i,j,bi,bj)
./model/src/calc_surf_dr.F 156 CALC_SURF_DR
          IF (ks.LE.Nr) THEN
./model/src/calc_surf_dr.F 157 CALC_SURF_DR
            hhm = rF(ks)
./model/src/calc_surf_dr.F 158 CALC_SURF_DR
            IF(ks.EQ.ksurfC(i,j-1,bi,bj)) hhm = rSurftmp(i,j-1)
./model/src/calc_surf_dr.F 159 CALC_SURF_DR
            hhp = rF(ks)
./model/src/calc_surf_dr.F 160 CALC_SURF_DR
            IF(ks.EQ.ksurfC(i,j,bi,bj))   hhp = rSurftmp(i,j)
./model/src/calc_surf_dr.F 163 CALC_SURF_DR
     &          - MAX(rF(ks+1),R_low(i,j-1,bi,bj),R_low(i,j,bi,bj)) 
./model/src/calc_surf_dr.F 164 CALC_SURF_DR
     &         )*recip_drF(ks)*maskS(i,j,ks,bi,bj)
./model/src/calc_surf_dr.F 51 CALC_SURF_DR
      INTEGER ks, numbWrite, numbWrMax
./model/src/calc_surf_dr.F 74 CALC_SURF_DR
          ks = ksurfC(i,j,bi,bj)
./model/src/calc_surf_dr.F 75 CALC_SURF_DR
          IF (ks.LE.Nr) THEN
./model/src/calc_surf_dr.F 80 CALC_SURF_DR
             hFactmp = ( rSurftmp(i,j)-MAX(rF(ks+1),R_low(i,j,bi,bj))
./model/src/calc_surf_dr.F 81 CALC_SURF_DR
     &                 )*recip_drF(ks)
./model/src/calc_surf_dr.F 85 CALC_SURF_DR
     &         ' i,j,k,bi,bj,Thid,Iter=',i,j,ks,bi,bj,myThid,myIter
./model/src/calc_surf_dr.F 89 CALC_SURF_DR
     &         ' i,j,k,bi,bj,Thid,Iter=',i,j,ks,bi,bj,myThid,myIter
./model/src/calc_surf_dr.F 93 CALC_SURF_DR
     &          hfacC(i,j,ks,bi,bj), hFactmp, etaFld(i,j,bi,bj)
./model/src/ini_phiref.F 165 INI_PHIREF
          ks = (k+1)/2
./model/src/ini_phiref.F 168 INI_PHIREF
          phiRef(k+1) = phiRef(k)+ddPI*tRef(ks)
./model/src/ini_phiref.F 40 INI_PHIREF
      INTEGER k, ks, stdUnit
./model/src/diags_oceanic_surf_flux.F 137 DIAGS_OCEANIC_SURF_FLUX
     &                        *salt(i,j,ks,bi,bj)
./model/src/diags_oceanic_surf_flux.F 48 DIAGS_OCEANIC_SURF_FLUX
      INTEGER ks
./model/src/diags_oceanic_surf_flux.F 54 DIAGS_OCEANIC_SURF_FLUX
        ks=Nr
./model/src/diags_oceanic_surf_flux.F 56 DIAGS_OCEANIC_SURF_FLUX
        ks=1
./model/src/diags_oceanic_surf_flux.F 88 DIAGS_OCEANIC_SURF_FLUX
     &                         *theta(i,j,ks,bi,bj)*HeatCapacity_Cp
./pkg/monitor/mon_surfcor.F 102 MON_SURFCOR
             ks = ksurfC(i,j,bi,bj)
./pkg/monitor/mon_surfcor.F 103 MON_SURFCOR
             IF (ks.LE.Nr) THEN 
./pkg/monitor/mon_surfcor.F 104 MON_SURFCOR
              tmpVal = salt(i,j,ks,bi,bj)
./pkg/monitor/mon_surfcor.F 105 MON_SURFCOR
     &               + salt(i,j,Nr,bi,bj)*drF(Nr)*recip_drF(ks)
./pkg/monitor/mon_surfcor.F 106 MON_SURFCOR
     &                *hFacC(i,j,Nr,bi,bj)*_recip_hFacC(i,j,ks,bi,bj)
./pkg/monitor/mon_surfcor.F 108 MON_SURFCOR
     &         + rA(i,j,bi,bj)*wVel(i,j,ks,bi,bj)*tmpVal
./pkg/monitor/mon_surfcor.F 53 MON_SURFCOR
      INTEGER i,j,k,ks,bi,bj 
./pkg/monitor/mon_surfcor.F 81 MON_SURFCOR
            ks = ksurfC(i,j,bi,bj)
./pkg/monitor/mon_surfcor.F 82 MON_SURFCOR
            IF (ks.LE.Nr) THEN
./pkg/monitor/mon_surfcor.F 85 MON_SURFCOR
     &          rA(i,j,bi,bj)*wVel(i,j,ks,bi,bj)*theta(i,j,ks,bi,bj)
./pkg/monitor/mon_surfcor.F 88 MON_SURFCOR
     &        + rA(i,j,bi,bj)*wVel(i,j,ks,bi,bj)* salt(i,j,ks,bi,bj)
./pkg/monitor/mon_surfcor.F 92 MON_SURFCOR
     &                 + tmpVal*atm_cp*((rC(ks)/atm_po)**atm_kappa) 
./pkg/obcs/obcs_apply_surf_dr.F 35 OBCS_APPLY_SURF_DR
      INTEGER i,j,ks
./pkg/obcs/obcs_apply_surf_dr.F 47 OBCS_APPLY_SURF_DR
        ks = ksurfS(i,j,bi,bj)
./pkg/obcs/obcs_apply_surf_dr.F 48 OBCS_APPLY_SURF_DR
        IF (ks.LE.Nr) THEN
./pkg/obcs/obcs_apply_surf_dr.F 49 OBCS_APPLY_SURF_DR
         hFactmp = OBNhFac0(i,bi,bj) + OBNeta(i,bi,bj)*recip_drF(ks)
./pkg/obcs/obcs_apply_surf_dr.F 56 OBCS_APPLY_SURF_DR
        ks = ksurfS(i,j+1,bi,bj)
./pkg/obcs/obcs_apply_surf_dr.F 57 OBCS_APPLY_SURF_DR
        IF (ks.LE.Nr) THEN
./pkg/obcs/obcs_apply_surf_dr.F 58 OBCS_APPLY_SURF_DR
         hFactmp = OBShFac0(i,bi,bj) + OBSeta(i,bi,bj)*recip_drF(ks)
./pkg/obcs/obcs_apply_surf_dr.F 69 OBCS_APPLY_SURF_DR
        ks = ksurfW(i,j,bi,bj)
./pkg/obcs/obcs_apply_surf_dr.F 70 OBCS_APPLY_SURF_DR
        IF (ks.LE.Nr) THEN
./pkg/obcs/obcs_apply_surf_dr.F 71 OBCS_APPLY_SURF_DR
         hFactmp = OBEhFac0(j,bi,bj) + OBEeta(j,bi,bj)*recip_drF(ks)
./pkg/obcs/obcs_apply_surf_dr.F 78 OBCS_APPLY_SURF_DR
        ks = ksurfW(i+1,j,bi,bj)
./pkg/obcs/obcs_apply_surf_dr.F 79 OBCS_APPLY_SURF_DR
        IF (ks.LE.Nr) THEN
./pkg/obcs/obcs_apply_surf_dr.F 80 OBCS_APPLY_SURF_DR
         hFactmp = OBWhFac0(j,bi,bj) + OBWeta(j,bi,bj)*recip_drF(ks)