Symbol Description Uses
theta potential temperature (oC, held at pressure/tracer point) 512



File Line number Procedure Code
./pkg/exf/exf_bulkformulae.F 288 EXF_BULKFORMULAE
           Tsf = theta(i,j,ksrf,bi,bj) + cen2kel
./pkg/exf/exf_bulkformulae.F 291 EXF_BULKFORMULAE
     &              *( theta(i,j,ksrf,bi,bj)-theta(i,j,ksrfp1,bi,bj) )
./pkg/exf/exf_getclim.F 163 EXF_GETCLIM
     &          theta(i,j,1,bi,bj) = climsst(i,j,bi,bj)
./pkg/exf/exf_getclim.F 180 EXF_GETCLIM
     &             theta(i,j,1,bi,bj) .EQ. 0. ) then
./pkg/exf/exf_getclim.F 181 EXF_GETCLIM
                print *, 'ph-warn-exf-clim ', i, j, theta(i,j,1,bi,bj)
./pkg/exf/exf_init.F 400 EXF_INIT
     &          theta(i,j,1,bi,bj) = climsst(i,j,bi,bj)
./pkg/exf/exf_init.F 416 EXF_INIT
     &             theta(i,j,1,bi,bj) .EQ. 0. ) then
./pkg/exf/exf_init.F 417 EXF_INIT
                print *, 'ph-warn-exf-init ', i, j, theta(i,j,1,bi,bj)
./pkg/exf/exf_mapfields.F 187 EXF_MAPFIELDS
     &               *( theta(i,j,ks,bi,bj) - temp_EvPrRn )
./pkg/exf/exf_mapfields.F 197 EXF_MAPFIELDS
     &               *( theta(i,j,ks,bi,bj) - temp_EvPrRn )
./pkg/exf/exf_mapfields.F 212 EXF_MAPFIELDS
     &              *( theta(i,j,ks,bi,bj) - runoftemp(i,j,bi,bj) )
./pkg/exf/exf_radiation.F 72 EXF_RADIATION
            Tsf = theta(i,j,ks,bi,bj) + cen2kel
./pkg/exf/exf_radiation.F 74 EXF_RADIATION
     &             *( theta(i,j,ks,bi,bj)-theta(i,j,kl,bi,bj) )
./pkg/exf/exf_radiation.F 97 EXF_RADIATION
     &          ((theta(i,j,ks,bi,bj)+cen2kel)**4)
./pkg/fizhi/fizhi_init_vars.F 182 FIZHI_INIT_VARS
        call DYN2PHYS(theta,pedyn,im1,im2,jm1,jm2,Nr,nSx,nSy,
./pkg/fizhi/fizhi_wrapper.F 125 FIZHI_WRAPPER
         thdynbef(i,j,L,bi,bj) = theta(i,j,L,bi,bj)
./pkg/fizhi/step_fizhi_corr.F 196 STEP_FIZHI_CORR
        thdyntemp(i,j,L,bi,bj)=theta(i,j,L,bi,bj)-thdyntemp(i,j,L,bi,bj)
./pkg/fizhi/step_fizhi_corr.F 429 STEP_FIZHI_CORR
     .    (theta(i,j,L,bi,bj)-thdynbef(i,j,L,bi,bj))*dtinv -
./pkg/flt/flt_traj.F 140 FLT_TRAJ
              CALL FLT_BILINEAR  (ix,jy,tt,theta, kp,0,bi,bj,myThid)
./pkg/flt/flt_up.F 138 FLT_UP
               CALL FLT_BILINEAR  (ix,jy,tt,theta, k,0,bi,bj,myThid)
./pkg/frazil/frazil_calc_rhs.F 74 FRAZIL_CALC_RHS
            tLoc = SW_TEMP(sLoc,theta(i,j,k,bi,bj),pLoc,0. _d 0)
./pkg/gmredi/gmredi_calc_tensor.F 1001 GMREDI_CALC_TENSOR
     &           (theta(i,j-1,k,bi,bj)-theta(i,j-1,kp1,bi,bj))
./pkg/gmredi/gmredi_calc_tensor.F 1003 GMREDI_CALC_TENSOR
     &           (theta(i, j ,k,bi,bj)-theta(i, j ,kp1,bi,bj))
./pkg/gmredi/gmredi_calc_tensor.F 794 GMREDI_CALC_TENSOR
     &           (theta(i-1,j,km1,bi,bj)-theta(i-1,j,k,bi,bj))
./pkg/gmredi/gmredi_calc_tensor.F 796 GMREDI_CALC_TENSOR
     &           (theta( i ,j,km1,bi,bj)-theta( i ,j,k,bi,bj))
./pkg/gmredi/gmredi_calc_tensor.F 800 GMREDI_CALC_TENSOR
     &           (theta(i-1,j,k,bi,bj)-theta(i-1,j,kp1,bi,bj))
./pkg/gmredi/gmredi_calc_tensor.F 802 GMREDI_CALC_TENSOR
     &           (theta( i ,j,k,bi,bj)-theta( i ,j,kp1,bi,bj))
./pkg/gmredi/gmredi_calc_tensor.F 995 GMREDI_CALC_TENSOR
     &           (theta(i,j-1,km1,bi,bj)-theta(i,j-1,k,bi,bj))
./pkg/gmredi/gmredi_calc_tensor.F 997 GMREDI_CALC_TENSOR
     &           (theta(i, j ,km1,bi,bj)-theta(i, j ,k,bi,bj))
./pkg/gmredi/gmredi_diagnostics_impl.F 61 GMREDI_DIAGNOSTICS_IMPL
     &            *(theta(i,j,k,bi,bj)-theta(i,j,k-1,bi,bj))
./pkg/icefront/icefront_thermodynamics.F 113 ICEFRONT_THERMODYNAMICS
            tLoc = theta(I,J,K,bi,bj)
./pkg/kl10/kl10_calc.F 82 KL10_CALC
            CALL FIND_RHO_SCALAR(theta(I,J,K,bi,bj), salt(I,J,K,bi,bj),
./pkg/kl10/kl10_calc.F 97 KL10_CALC
               CALL FIND_RHO_SCALAR(theta(I,J,K,bi,bj), salt(I,J,K,bi
./pkg/kpp/kpp_routines.F 1630 STATEKPP
     I     theta(1-OLx,1-OLy,1,bi,bj), salt(1-OLx,1-OLy,1,bi,bj),
./pkg/kpp/kpp_routines.F 1672 STATEKPP
     I        theta(1-OLx,1-OLy,k,bi,bj), salt(1-OLx,1-OLy,k,bi,bj),
./pkg/kpp/kpp_routines.F 1684 STATEKPP
     I        theta(1-OLx,1-OLy,k-1,bi,bj),salt(1-OLx,1-OLy,k-1,bi,bj),
./pkg/kpp/kpp_routines.F 1696 STATEKPP
     I        theta(1-OLx,1-OLy,1,bi,bj), salt(1-OLx,1-OLy,1,bi,bj),
./pkg/kpp/kpp_routines.F 1843 KPP_DOUBLEDIFF
         alphaDT(I,J) = ( theta(I,J,Km1,bi,bj)-theta(I,J,K,bi,bj) )
./model/src/apply_forcing.F 506 APPLY_FORCING_T
     &          *( temp_addMass - theta(i,j,k,bi,bj) )
./model/src/apply_forcing.F 576 APPLY_FORCING_T
     &                  *( theta(i,j,km,bi,bj)*salt(i,j,km,bi,bj)
./model/src/apply_forcing.F 577 APPLY_FORCING_T
     &                   + theta(i,j,kc,bi,bj)*salt(i,j,kc,bi,bj)
./model/src/apply_forcing.F 589 APPLY_FORCING_T
     &                  *( theta(i,j,kc,bi,bj)*salt(i,j,kc,bi,bj)
./model/src/apply_forcing.F 590 APPLY_FORCING_T
     &                   + theta(i,j,kp,bi,bj)*salt(i,j,kp,bi,bj)
./pkg/layers/layers_calc.F 101 LAYERS_CALC
        CALL LAYERS_DIAPYCNAL( theta,iLa,
./pkg/layers/layers_calc.F 132 LAYERS_CALC
     &                       theta(1-OLx,1-OLy,k,bi,bj),
./pkg/layers/layers_calc.F 94 LAYERS_CALC
        CALL LAYERS_FLUXCALC( uVel,vVel,theta,iLa,
./pkg/longstep/longstep_average.F 93 LONGSTEP_AVERAGE
      CALL LONGSTEP_FILL_3D(LS_thetaCount, LS_theta, theta, Nr, myThid)
./model/src/calc_oce_mxlayer.F 112 CALC_OCE_MXLAYER
     I        theta(1-OLx,1-OLy,k,bi,bj), salt(1-OLx,1-OLy,k,bi,bj),
./model/src/calc_oce_mxlayer.F 158 CALC_OCE_MXLAYER
     I        theta(1-OLx,1-OLy,k,bi,bj), salt(1-OLx,1-OLy,k,bi,bj),
./pkg/monitor/mon_solution.F 45 MON_SOLUTION
     I                Nr, theta, hFacC, maskInC, rA, drF,
./pkg/monitor/mon_surfcor.F 133 MON_SURFCOR
     &           *(theta(i,j,k,bi,bj)+theta(i,j,k-1,bi,bj))*0.5 _d 0
./pkg/monitor/mon_surfcor.F 156 MON_SURFCOR
              tmpVal  = rStarDhCDt(i,j,bi,bj)*theta(i,j,k,bi,bj)
./pkg/monitor/mon_surfcor.F 90 MON_SURFCOR
     &               *wVel(i,j,ks,bi,bj)*theta(i,j,ks,bi,bj)
./pkg/monitor/monitor.F 117 MONITOR
         CALL MON_WRITESTATS_RL( Nr, theta,'_theta',
./pkg/monitor/monitor.F 130 MONITOR
         CALL MON_WRITESTATS_RL( 1, theta(1-OLx,1-OLy,k,1,1), '_sst',
./pkg/my82/my82_ri_number.F 91 MY82_RI_NUMBER
     I     theta(1-OLx,1-OLy,Km1,bi,bj), salt(1-OLx,1-OLy,Km1,bi,bj),
./pkg/my82/my82_ri_number.F 96 MY82_RI_NUMBER
     I     theta(1-OLx,1-OLy,K,bi,bj), salt(1-OLx,1-OLy,K,bi,bj),
./pkg/obcs/obcs_calc.F 11 OBCS_CALC
     &                      uVel, vVel, wVel, theta, salt,
./pkg/obcs/obcs_calc.F 146 OBCS_CALC
     &          uVel, vVel, wVel, theta, salt,
./pkg/obcs/obcs_calc.F 196 OBCS_CALC
     &          uVel, vVel, wVel, theta, salt,
./pkg/obcs/obcs_calc.F 230 OBCS_CALC
     &          uVel, vVel, wVel, theta, salt,
./pkg/obcs/obcs_calc.F 49 OBCS_CALC
      _RL theta(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
./pkg/obcs/obcs_calc.F 96 OBCS_CALC
     &          uVel, vVel, wVel, theta, salt,
./pkg/obcs/obcs_calc_stevens.F 252 OBCS_CALC_STEVENS
            dTracSpace = (theta(I-1,J,K,bi,bj)-theta(I-2,J,K,bi,bj))
./pkg/obcs/obcs_calc_stevens.F 254 OBCS_CALC_STEVENS
            dTracTime  = (theta(I-1,J,K,bi,bj)-OBEtStevens(J,K,bi,bj))
./pkg/obcs/obcs_calc_stevens.F 265 OBCS_CALC_STEVENS
     &          *(theta(I,J,K,bi,bj)-theta(I-1,J,K,bi,bj))
./pkg/obcs/obcs_calc_stevens.F 268 OBCS_CALC_STEVENS
     &           * (theta(I,J,K,bi,bj)-OBEt(J,K,bi,bj)) )
./pkg/obcs/obcs_calc_stevens.F 390 OBCS_CALC_STEVENS
            dTracSpace = (theta(I+2,J,K,bi,bj)-theta(I+1,J,K,bi,bj))
./pkg/obcs/obcs_calc_stevens.F 392 OBCS_CALC_STEVENS
            dTracTime  = (theta(I+1,J,K,bi,bj)-OBWtStevens(J,K,bi,bj))
./pkg/obcs/obcs_calc_stevens.F 403 OBCS_CALC_STEVENS
     &        *(theta(I+1,J,K,bi,bj)-theta(I,J,K,bi,bj))
./pkg/obcs/obcs_calc_stevens.F 406 OBCS_CALC_STEVENS
     &           * (theta(I,J,K,bi,bj)-OBWt(J,K,bi,bj)) )
./pkg/obcs/obcs_calc_stevens.F 527 OBCS_CALC_STEVENS
            dTracSpace = (theta(I,J-1,K,bi,bj)-theta(I,J-2,K,bi,bj))
./pkg/obcs/obcs_calc_stevens.F 529 OBCS_CALC_STEVENS
            dTracTime  = (theta(I,J-1,K,bi,bj)-OBNtStevens(I,K,bi,bj))
./pkg/obcs/obcs_calc_stevens.F 540 OBCS_CALC_STEVENS
     &        *(theta(I,J,K,bi,bj)-theta(I,J-1,K,bi,bj))
./pkg/obcs/obcs_calc_stevens.F 543 OBCS_CALC_STEVENS
     &           * (theta(I,J,K,bi,bj)-OBNt(I,K,bi,bj)) )
./pkg/obcs/obcs_calc_stevens.F 665 OBCS_CALC_STEVENS
            dTracSpace = (theta(I,J+2,K,bi,bj)-theta(I,J+1,K,bi,bj))
./pkg/obcs/obcs_calc_stevens.F 667 OBCS_CALC_STEVENS
            dTracTime  = (theta(I,J+1,K,bi,bj)-OBStStevens(I,K,bi,bj))
./pkg/obcs/obcs_calc_stevens.F 678 OBCS_CALC_STEVENS
     &        *(theta(I,J+1,K,bi,bj)-theta(I,J,K,bi,bj))
./pkg/obcs/obcs_calc_stevens.F 681 OBCS_CALC_STEVENS
     &           * (theta(I,J,K,bi,bj)-OBSt(I,K,bi,bj)) )
./pkg/obcs/obcs_calc_stevens.F 815 OBCS_STEVENS_SAVE_TRACERS
            OBNtStevens(i,k,bi,bj) = theta(i,Jobc-1,k,bi,bj)
./pkg/obcs/obcs_calc_stevens.F 831 OBCS_STEVENS_SAVE_TRACERS
            OBStStevens(i,k,bi,bj) = theta(i,Jobc+1,k,bi,bj)
./pkg/obcs/obcs_calc_stevens.F 847 OBCS_STEVENS_SAVE_TRACERS
            OBEtStevens(j,k,bi,bj) = theta(Iobc-1,j,k,bi,bj)
./pkg/obcs/obcs_calc_stevens.F 863 OBCS_STEVENS_SAVE_TRACERS
            OBWtStevens(j,k,bi,bj) = theta(Iobc+1,j,k,bi,bj)
./model/src/convective_adjustment.F 150 CONVECTIVE_ADJUSTMENT
     I           theta(1-OLx,1-OLy,K-1,bi,bj),
./model/src/convective_adjustment.F 164 CONVECTIVE_ADJUSTMENT
     I           theta(1-OLx,1-OLy,K,bi,bj),
./model/src/convective_adjustment.F 190 CONVECTIVE_ADJUSTMENT
     U                              theta,
./pkg/obcs/obcs_init_variables.F 427 OBCS_INIT_VARIABLES
     &              uVel, vVel, wVel, theta, salt, myThid )
./pkg/obcs/obcs_init_variables.F 438 OBCS_INIT_VARIABLES
          CALL OBCS_APPLY_TS( bi, bj, 0, theta, salt, myThid )
./pkg/obcs/obcs_init_variables.F 450 OBCS_INIT_VARIABLES
        _EXCH_XYZ_RL( theta, myThid )
./pkg/obcs/obcs_monitor.F 132 OBCS_MONITOR
     I             Nr, Nr, 0, theta, hFacW, dyG, drF, maskInC,
./pkg/obcs/obcs_monitor.F 138 OBCS_MONITOR
     I             Nr, Nr, 0, theta, hFacS, dxG, drF, maskInC,
./model/src/convective_adjustment_ini.F 150 CONVECTIVE_ADJUSTMENT_INI
     I           theta(1-OLx,1-OLy,K-1,bi,bj),
./model/src/convective_adjustment_ini.F 164 CONVECTIVE_ADJUSTMENT_INI
     I           theta(1-OLx,1-OLy,K,bi,bj),
./model/src/convective_adjustment_ini.F 188 CONVECTIVE_ADJUSTMENT_INI
     U                              theta,
./pkg/obcs/obcs_sponge.F 475 OBCS_SPONGE_T
     &    + float(jsl)*theta(i,j,k,bi,bj) )
./pkg/obcs/obcs_sponge.F 490 OBCS_SPONGE_T
     &    * ( theta(i,j,k,bi,bj) - trelax )
./pkg/obcs/obcs_sponge.F 515 OBCS_SPONGE_T
     &    + float(jsl)*theta(i,j,k,bi,bj) )
./pkg/obcs/obcs_sponge.F 531 OBCS_SPONGE_T
     &   * ( theta(i,j,k,bi,bj) - trelax )
./pkg/obcs/obcs_sponge.F 556 OBCS_SPONGE_T
     &    + float(isl)*theta(i,j,k,bi,bj) )
./pkg/obcs/obcs_sponge.F 572 OBCS_SPONGE_T
     &    * ( theta(i,j,k,bi,bj) - trelax )
./pkg/obcs/obcs_sponge.F 599 OBCS_SPONGE_T
     &    + float(isl)*theta(i,j,k,bi,bj) )
./pkg/obcs/obcs_sponge.F 615 OBCS_SPONGE_T
     &    * ( theta(i,j,k,bi,bj) - trelax )
./pkg/obcs/orlanski_east.F 167 ORLANSKI_EAST
                  CL=-(theta(I_obc-1,J,K,bi,bj)-TE_STORE_1(J,K,bi,bj))/
./pkg/obcs/orlanski_east.F 184 ORLANSKI_EAST
               OBEt(J,K,bi,bj)=theta(I_obc,J,K,bi,bj)-
./pkg/obcs/orlanski_east.F 186 ORLANSKI_EAST
     &           (ab1*(theta(I_obc,J,K,bi,bj)-theta(I_obc-1,J,K,bi,bj))+
./pkg/obcs/orlanski_east.F 265 ORLANSKI_EAST
               TE_STORE_2(J,K,bi,bj)=theta(I_obc-1,J,K,bi,bj) -
./pkg/obcs/orlanski_east.F 266 ORLANSKI_EAST
     &         theta(I_obc-2,J,K,bi,bj)
./pkg/obcs/orlanski_east.F 267 ORLANSKI_EAST
               TE_STORE_1(J,K,bi,bj)=theta(I_obc-1,J,K,bi,bj)
./pkg/obcs/orlanski_east.F 268 ORLANSKI_EAST
               TE_STORE_4(J,K,bi,bj)=theta(I_obc,J,K,bi,bj)
./pkg/obcs/orlanski_east.F 85 ORLANSKI_EAST
      _RL theta(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
./pkg/obcs/orlanski_east.F 8 ORLANSKI_EAST
     I                      uVel, vVel, wVel, theta, salt,
./pkg/obcs/orlanski_north.F 146 ORLANSKI_NORTH
                  CL=-(theta(I,J_obc-1,K,bi,bj)-TN_STORE_1(I,K,bi,bj))/
./pkg/obcs/orlanski_north.F 157 ORLANSKI_NORTH
               OBNt(I,K,bi,bj)=theta(I,J_obc,K,bi,bj)-
./pkg/obcs/orlanski_north.F 159 ORLANSKI_NORTH
     &          (ab1*(theta(I,J_obc,K,bi,bj)-theta(I,J_obc-1,K,bi,bj))+
./pkg/obcs/orlanski_north.F 226 ORLANSKI_NORTH
               TN_STORE_2(I,K,bi,bj)=theta(I,J_obc-1,K,bi,bj) -
./pkg/obcs/orlanski_north.F 227 ORLANSKI_NORTH
     &         theta(I,J_obc-2,K,bi,bj)
./pkg/obcs/orlanski_north.F 228 ORLANSKI_NORTH
               TN_STORE_1(I,K,bi,bj)=theta(I,J_obc-1,K,bi,bj)
./pkg/obcs/orlanski_north.F 229 ORLANSKI_NORTH
               TN_STORE_4(I,K,bi,bj)=theta(I,J_obc,K,bi,bj)
./pkg/obcs/orlanski_north.F 76 ORLANSKI_NORTH
      _RL theta(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
./pkg/obcs/orlanski_north.F 7 ORLANSKI_NORTH
     I                      uVel, vVel, wVel, theta, salt,
./pkg/obcs/orlanski_south.F 146 ORLANSKI_SOUTH
                  CL=(theta(I,J_obc+1,K,bi,bj)-TS_STORE_1(I,K,bi,bj))/
./pkg/obcs/orlanski_south.F 157 ORLANSKI_SOUTH
               OBSt(I,K,bi,bj)=theta(I,J_obc,K,bi,bj)+
./pkg/obcs/orlanski_south.F 159 ORLANSKI_SOUTH
     &          (ab1*(theta(I,J_obc+1,K,bi,bj)-theta(I,J_obc,K,bi,bj))+
./pkg/obcs/orlanski_south.F 226 ORLANSKI_SOUTH
               TS_STORE_2(I,K,bi,bj)=theta(I,J_obc+2,K,bi,bj) -
./pkg/obcs/orlanski_south.F 227 ORLANSKI_SOUTH
     &         theta(I,J_obc+1,K,bi,bj)
./pkg/obcs/orlanski_south.F 228 ORLANSKI_SOUTH
               TS_STORE_1(I,K,bi,bj)=theta(I,J_obc+1,K,bi,bj)
./pkg/obcs/orlanski_south.F 229 ORLANSKI_SOUTH
               TS_STORE_4(I,K,bi,bj)=theta(I,J_obc,K,bi,bj)
./pkg/obcs/orlanski_south.F 76 ORLANSKI_SOUTH
      _RL theta(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
./pkg/obcs/orlanski_south.F 7 ORLANSKI_SOUTH
     I                      uVel, vVel, wVel, theta, salt,
./pkg/obcs/orlanski_west.F 169 ORLANSKI_WEST
                  CL=(theta(I_obc+1,J,K,bi,bj)-TW_STORE_1(J,K,bi,bj))/
./pkg/obcs/orlanski_west.F 186 ORLANSKI_WEST
               OBWt(J,K,bi,bj)=theta(I_obc,J,K,bi,bj)+
./pkg/obcs/orlanski_west.F 188 ORLANSKI_WEST
     &          (ab1*(theta(I_obc+1,J,K,bi,bj)-theta(I_obc,J,K,bi,bj))+
./pkg/obcs/orlanski_west.F 267 ORLANSKI_WEST
               TW_STORE_2(J,K,bi,bj)=theta(I_obc+2,J,K,bi,bj) -
./pkg/obcs/orlanski_west.F 268 ORLANSKI_WEST
     &         theta(I_obc+1,J,K,bi,bj)
./pkg/obcs/orlanski_west.F 269 ORLANSKI_WEST
               TW_STORE_1(J,K,bi,bj)=theta(I_obc+1,J,K,bi,bj)
./pkg/obcs/orlanski_west.F 270 ORLANSKI_WEST
               TW_STORE_4(J,K,bi,bj)=theta(I_obc,J,K,bi,bj)
./pkg/obcs/orlanski_west.F 87 ORLANSKI_WEST
      _RL theta(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
./pkg/obcs/orlanski_west.F 8 ORLANSKI_WEST
     I                      uVel, vVel, wVel, theta, salt,
./pkg/ocn_compon_interf/ocn_store_my_data.F 61 OCN_STORE_MY_DATA
          SSTocn2cpl(i,j,bi,bj) = theta(i,j,1,bi,bj)
./pkg/offline/offline_fields_load.F 366 OFFLINE_FIELDS_LOAD
           theta(i,j,k,bi,bj)= bWght*tave0(i,j,k,bi,bj)
./pkg/openad/openad_dumpAdjoint.F 103 OPENAD_DUMPADJOINT
         foo5=thetad
./pkg/opps/opps_interface.F 108 OPPS_INTERFACE
           tracerLoc(K,1) = theta(I,J,K,bi,bj)
./pkg/opps/opps_interface.F 114 OPPS_INTERFACE
           tracerLoc(K,1) = theta(I,J,K,bi,bj)
./pkg/opps/opps_interface.F 172 OPPS_INTERFACE
          theta(I,J,K,bi,bj) = tracerLoc(K,1)
./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
./pkg/pp81/pp81_ri_number.F 100 PP81_RI_NUMBER
     I     theta(1-OLx,1-OLy,K,bi,bj), salt(1-OLx,1-OLy,K,bi,bj),
./pkg/pp81/pp81_ri_number.F 95 PP81_RI_NUMBER
     I     theta(1-OLx,1-OLy,Km1,bi,bj), salt(1-OLx,1-OLy,Km1,bi,bj),
./pkg/profiles/profiles_interp.F 150 PROFILES_INTERP
               tab_coeffs1(1,1)=theta(prof_i,prof_j,k,bi,bj) !SO
./pkg/profiles/profiles_interp.F 151 PROFILES_INTERP
               tab_coeffs1(1,2)=theta(prof_i+1,prof_j,k,bi,bj) !SE
./pkg/profiles/profiles_interp.F 152 PROFILES_INTERP
               tab_coeffs1(2,1)=theta(prof_i,prof_j+1,k,bi,bj) !NO
./pkg/profiles/profiles_interp.F 153 PROFILES_INTERP
               tab_coeffs1(2,2)=theta(prof_i+1,prof_j+1,k,bi,bj) !NZ              
./pkg/profiles/profiles_interp_gg.F 74 PROFILES_INTERP_GG
               tab_coeffs1(q)=theta(i_cur(q),j_cur(q),k,bi,bj)
./model/src/do_atmospheric_phys.F 79 DO_ATMOSPHERIC_PHYS
     &         ( theta(i,j,k,bi,bj)
./model/src/do_fields_blocking_exchanges.F 63 DO_FIELDS_BLOCKING_EXCHANGES
       _EXCH_XYZ_RL( theta, myThid )
./model/src/do_oceanic_phys.F 189 DO_OCEANIC_PHYS
     I                 uVel, vVel, wVel, theta, salt, myThid )
./model/src/do_oceanic_phys.F 668 DO_OCEANIC_PHYS
     I                  theta, salt,
./model/src/do_oceanic_phys.F 681 DO_OCEANIC_PHYS
     I                  theta, salt,
./model/src/do_oceanic_phys.F 692 DO_OCEANIC_PHYS
     I                theta(1-OLx,1-OLy,k,bi,bj),
./model/src/do_oceanic_phys.F 712 DO_OCEANIC_PHYS
     &         ( theta(i,j,k,bi,bj)
./model/src/do_oceanic_phys.F 760 DO_OCEANIC_PHYS
     I                 theta(1-OLx,1-OLy,k-1,bi,bj),
./model/src/do_oceanic_phys.F 846 DO_OCEANIC_PHYS
     I              theta, 0,
./model/src/do_stagger_fields_exchanges.F 52 DO_STAGGER_FIELDS_EXCHANGES
        CALL EXCH_SM_3D_RL( theta, .FALSE., Nr, myThid )
./pkg/rbcs/rbcs_add_tendency.F 97 RBCS_ADD_TENDENCY
     &        *( theta(i,j,k,bi,bj)- RBCtemp(i,j,k,bi,bj) )
./model/src/do_the_model_io.F 62 DO_THE_MODEL_IO
          CALL PLOT_FIELD_XYZRL( theta, 'Current theta ',
./model/src/dynamics.F 483 DYNAMICS
     I        theta, salt,
./model/src/dynamics.F 705 DYNAMICS
       CALL DEBUG_STATS_RL(Nr,theta,'Theta (DYNAMICS)',myThid)
./pkg/salt_plume/salt_plume_calc_depth.F 113 SALT_PLUME_CALC_DEPTH
     I       theta(1-OLx,1-OLy,K,bi,bj), salt(1-OLx,1-OLy,K,bi,bj),
./pkg/salt_plume/salt_plume_forcing_surf.F 61 SALT_PLUME_FORCING_SURF
     &   *theta(i,j,1,bi,bj)*rhoConst
./model/src/external_forcing.F 393 EXTERNAL_FORCING_T
     &          *( temp_addMass - theta(i,j,kLev,bi,bj) )
./model/src/external_forcing.F 463 EXTERNAL_FORCING_T
     &                  *( theta(i,j,km,bi,bj)*salt(i,j,km,bi,bj)
./model/src/external_forcing.F 464 EXTERNAL_FORCING_T
     &                   + theta(i,j,kc,bi,bj)*salt(i,j,kc,bi,bj)
./model/src/external_forcing.F 476 EXTERNAL_FORCING_T
     &                  *( theta(i,j,kc,bi,bj)*salt(i,j,kc,bi,bj)
./model/src/external_forcing.F 477 EXTERNAL_FORCING_T
     &                   + theta(i,j,kp,bi,bj)*salt(i,j,kp,bi,bj)
./model/src/external_forcing_surf.F 259 EXTERNAL_FORCING_SURF
     &          *( temp_EvPrRn - theta(i,j,ks,bi,bj) )
./model/src/external_forcing_surf.F 294 EXTERNAL_FORCING_SURF
     &           *( theta(i,j,ks,bi,bj) - temp_EvPrRn )
./model/src/find_alpha.F 126 FIND_ALPHA
     I        theta(1-OLx,1-OLy,k,bi,bj), salt(1-OLx,1-OLy,k,bi,bj),
./model/src/find_alpha.F 132 FIND_ALPHA
     I        theta(1-OLx,1-OLy,k,bi,bj), salt(1-OLx,1-OLy,k,bi,bj),
./model/src/find_alpha.F 140 FIND_ALPHA
               t1 = theta(i,j,k,bi,bj)
./model/src/find_alpha.F 233 FIND_ALPHA
     I        theta(1-OLx,1-OLy,k,bi,bj), salt(1-OLx,1-OLy,k,bi,bj),
./model/src/find_alpha.F 239 FIND_ALPHA
     I        theta(1-OLx,1-OLy,k,bi,bj), salt(1-OLx,1-OLy,k,bi,bj),
./model/src/find_alpha.F 245 FIND_ALPHA
               t1  = theta(i,j,k,bi,bj)
./model/src/find_alpha.F 291 FIND_ALPHA
     I      theta(1-OLx,1-OLy,k,bi,bj), salt(1-OLx,1-OLy,k,bi,bj),
./model/src/find_alpha.F 297 FIND_ALPHA
         ct      = theta(i,j,k,bi,bj)
./model/src/find_alpha.F 423 FIND_BETA
               tP=theta(i,j,k,bi,bj)-refTemp
./model/src/find_alpha.F 454 FIND_BETA
     I        theta(1-OLx,1-OLy,k,bi,bj), salt(1-OLx,1-OLy,k,bi,bj),
./model/src/find_alpha.F 460 FIND_BETA
     I        theta(1-OLx,1-OLy,k,bi,bj), salt(1-OLx,1-OLy,k,bi,bj),
./model/src/find_alpha.F 468 FIND_BETA
               t1 = theta(i,j,k,bi,bj)
./model/src/find_alpha.F 547 FIND_BETA
     I        theta(1-OLx,1-OLy,k,bi,bj), salt(1-OLx,1-OLy,k,bi,bj),
./model/src/find_alpha.F 553 FIND_BETA
     I        theta(1-OLx,1-OLy,k,bi,bj), salt(1-OLx,1-OLy,k,bi,bj),
./model/src/find_alpha.F 559 FIND_BETA
               t1  = theta(i,j,k,bi,bj)
./model/src/find_alpha.F 598 FIND_BETA
     I      theta(1-OLx,1-OLy,k,bi,bj), salt(1-OLx,1-OLy,k,bi,bj),
./model/src/find_alpha.F 604 FIND_BETA
         ct      = theta(i,j,k,bi,bj)
./model/src/find_alpha.F 93 FIND_ALPHA
               tP=theta(i,j,k,bi,bj)-refTemp
./pkg/seaice/cost_ice_test.F 92 COST_ICE_TEST
      print *, 'ph-ice B ', myiter, theta(4,4,kSrf,1,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 85 FORCING_SURF_RELAX
     &         *(theta(i,j,ks,bi,bj)-SST(i,j,bi,bj))
./pkg/seaice/seaice_fake.F 74 SEAICE_FAKE
            theta(i,j,1,bi,bj)=theta(i,j,1,bi,bj) + fac *
./pkg/seaice/seaice_fake.F 75 SEAICE_FAKE
     &        ( tempFrz-theta(i,j,1,bi,bj) ) *
./pkg/seaice/seaice_growth.F 1024 SEAICE_GROWTH
          IF ( theta(I,J,kSurface,bi,bj) .GE. tempFrz ) THEN
./pkg/seaice/seaice_growth.F 1042 SEAICE_GROWTH
     &         * (theta(I,J,kSurface,bi,bj)-tempFrz)
./pkg/seaice/seaice_growth.F 1061 SEAICE_GROWTH
           tmpscal1=(theta(I,J,kSurface,bi,bj)-tempFrz)
./pkg/seaice/seaice_growth.F 2239 SEAICE_GROWTH
     &         HeatCapacity_Cp * theta(I,J,kSurface,bi,bj)
./pkg/seaice/seaice_growth.F 2242 SEAICE_GROWTH
     &       ( temp_EvPrRn - theta(I,J,kSurface,bi,bj) )
./pkg/seaice/seaice_growth.F 2307 SEAICE_GROWTH
     &        ( ZERO + HeatCapacity_Cp * theta(I,J,kSurface,bi,bj) )
./pkg/seaice/seaice_growth.F 2310 SEAICE_GROWTH
     &        ( temp_EvPrRn - theta(I,J,kSurface,bi,bj) )
./pkg/seaice/seaice_growth.F 2449 SEAICE_GROWTH
     &       HeatCapacity_Cp * theta(I,J,kSurface,bi,bj)
./pkg/seaice/seaice_growth.F 2609 SEAICE_GROWTH
     &       ( ZERO + HeatCapacity_Cp * theta(I,J,kSurface,bi,bj) )
./pkg/seaice/seaice_growth.F 2612 SEAICE_GROWTH
     &       HeatCapacity_Cp*(temp_EvPrRn - theta(I,J,kSurface,bi,bj))
./pkg/seaice/seaice_growth.F 716 SEAICE_GROWTH
          TmixLoc(i,j) = theta(i,j,kSurface,bi,bj)+celsius2K
./model/src/freeze_surface.F 58 FREEZE_SURFACE
          IF (theta(i,j,k,bi,bj) .LT. Tfreezing) THEN
./model/src/freeze_surface.F 60 FREEZE_SURFACE
     &            ( Tfreezing - theta(i,j,k,bi,bj) )
./model/src/freeze_surface.F 62 FREEZE_SURFACE
             theta(i,j,k,bi,bj) = Tfreezing
./pkg/shelfice/shelfice_thermodynamics.F 240 SHELFICE_THERMODYNAMICS
          tLoc(I,J) = theta(I,J,K,bi,bj)
./pkg/shelfice/shelfice_thermodynamics.F 264 SHELFICE_THERMODYNAMICS
     &           + theta(I,J,Kp1,bi,bj) *drKp1 )
./model/src/ini_dynvars.F 58 INI_DYNVARS
           theta(i,j,k,bi,bj) = 0. _d 0
./pkg/thsice/thsice_get_ocean.F 96 THSICE_GET_OCEAN
          tOceMxL(i,j,bi,bj) = theta(i,j,ks,bi,bj)
./pkg/timeave/timeave_statvars.F 115 TIMEAVE_STATVARS
      CALL TIMEAVE_CUMULATE(thetatave,theta, Nr, DDTT, bi, bj, myThid)
./pkg/timeave/timeave_statvars.F 124 TIMEAVE_STATVARS
      CALL TIMEAVE_CUMUL_2V(TTtave, theta,theta, Nr, 0,
./pkg/timeave/timeave_statvars.F 134 TIMEAVE_STATVARS
      CALL TIMEAVE_CUMUL_2V(WTtave, theta, wVel, Nr, 3,
./pkg/timeave/timeave_statvars.F 154 TIMEAVE_STATVARS
      CALL TIMEAVE_CUMUL_2VFC(UTtave, theta, uVel,  hFacW, Nr, 1,
./pkg/timeave/timeave_statvars.F 156 TIMEAVE_STATVARS
      CALL TIMEAVE_CUMUL_2VFC(VTtave, theta, vVel,  hFacS, Nr, 2,
./pkg/timeave/timeave_surf_flux.F 81 TIMEAVE_SURF_FLUX
     &      + PmEpR(i,j,bi,bj)*theta(i,j,k,bi,bj)*HeatCapacity_Cp
./model/src/ini_pressure.F 116 INI_PRESSURE
     I           theta, salt,
./verification/advect_xy/code/ini_theta.F 102 INI_THETA
     &      .AND. theta(i,j,k,bi,bj).EQ.0. ) THEN
./verification/advect_xy/code/ini_theta.F 139 INI_THETA
            IF (theta(i,j,k,bi,bj) .LT. Tfreezing) THEN
./verification/advect_xy/code/ini_theta.F 140 INI_THETA
               theta(i,j,k,bi,bj) = Tfreezing
./verification/advect_xy/code/ini_theta.F 150 INI_THETA
        CALL PLOT_FIELD_XYZRL( theta, 'Initial Temperature',
./verification/advect_xy/code/ini_theta.F 59 INI_THETA
           theta(i,j,k,bi,bj) = tRef(k)
./verification/advect_xy/code/ini_theta.F 70 INI_THETA
            theta(i,j,k,bi,bj) = EXP( -0.5 _d 0*( rD/20. _d 3 )**2 )
./verification/advect_xy/code/ini_theta.F 78 INI_THETA
      _EXCH_XYZ_RL( theta, myThid )
./verification/advect_xy/code/ini_theta.F 82 INI_THETA
        CALL READ_FLD_XYZ_RL( hydrogThetaFile, ' ', theta, 0, myThid )
./verification/advect_xy/code/ini_theta.F 83 INI_THETA
        _EXCH_XYZ_RL(theta,myThid)
./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/ini_theta.F 103 INI_THETA
     &      .AND. theta(i,j,k,bi,bj).EQ.0. ) THEN
./model/src/ini_theta.F 140 INI_THETA
            IF (theta(i,j,k,bi,bj) .LT. Tfreezing) THEN
./model/src/ini_theta.F 141 INI_THETA
               theta(i,j,k,bi,bj) = Tfreezing
./model/src/ini_theta.F 153 INI_THETA
        CALL PLOT_FIELD_XYZRL( theta, 'Initial Temperature',
./model/src/ini_theta.F 62 INI_THETA
           theta(i,j,k,bi,bj) = tRef(k)
./model/src/ini_theta.F 76 INI_THETA
          CALL MNC_CW_RL_R('D',hydrogThetaFile,0,0,'Temp',theta,myThid)
./model/src/ini_theta.F 80 INI_THETA
          CALL READ_FLD_XYZ_RL( hydrogThetaFile, ' ', theta, 0, myThid )
./model/src/ini_theta.F 84 INI_THETA
        _EXCH_XYZ_RL(theta,myThid)
./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/dome/code/obcs_calc.F 131 OBCS_CALC
     &          uVel, vVel, wVel, theta, salt,
./verification/dome/code/obcs_calc.F 150 OBCS_CALC
     &          uVel, vVel, wVel, theta, salt,
./verification/dome/code/obcs_calc.F 211 OBCS_CALC
     &          uVel, vVel, wVel, theta, salt,
./verification/dome/code/obcs_calc.F 32 OBCS_CALC
      _RL theta(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
./verification/dome/code/obcs_calc.F 7 OBCS_CALC
     &                      uVel, vVel, wVel, theta, salt,
./verification/dome/code/obcs_calc.F 99 OBCS_CALC
     &          uVel, vVel, wVel, theta, salt,
./verification/fizhi-cs-aqualev20/code/ini_salt.F 76 INI_SALT
           temperature(I,J,K,bi,bj) = theta(I,J,K,bi,bj) * pkappa(K)
./verification/fizhi-gridalt-hs/code/fizhi_init_vars.F 162 FIZHI_INIT_VARS
        call DYN2PHYS(theta,pedyn,im1,im2,jm1,jm2,Nr,nSx,nSy,
./verification/global_ocean.cs32x15/code_ad/cost_test.F 74 COST_TEST
               objf_test(bi,bj) = theta(i,j,kLocOut,bi,bj)
./verification/global_ocean.cs32x15/code_ad/cost_test.F 77 COST_TEST
     &              theta(i,j,kLocOut,bi,bj)
./verification/global_ocean.cs32x15/code_ad/cost_test.F 84 COST_TEST
     &        (theta(i,j,k,bi,bj)-thetaLev(i,j,k,bi,bj))**2
./verification/halfpipe_streamice/code_ad/ctrl_map_ini_genarr.F 153 CTRL_MAP_INI_GENARR
                theta(i,j,k,bi,bj) = theta(i,j,k,bi,bj)
./verification/halfpipe_streamice/code_ad/ctrl_map_ini_genarr.F 171 CTRL_MAP_INI_GENARR
       _EXCH_XYZ_RL( theta, myThid )
./verification/halfpipe_streamice/code_oad/ctrl_map_ini_genarr.F 161 CTRL_MAP_INI_GENARR
                theta(i,j,k,bi,bj) = theta(i,j,k,bi,bj)
./verification/halfpipe_streamice/code_oad/ctrl_map_ini_genarr.F 179 CTRL_MAP_INI_GENARR
       _EXCH_XYZ_RL( theta, myThid )
./verification/hs94.128x64x5/code/external_forcing.F 206 EXTERNAL_FORCING_T
     &        - kT*( theta(i,j,kLev,bi,bj)-thetaEq )
./verification/hs94.1x64x5/code/apply_forcing.F 290 APPLY_FORCING_T
     &               - kT*( theta(i,j,k,bi,bj)-thetaEq )
./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 108 INI_THETA
     &      .AND. theta(I,J,K,bi,bj).EQ.0. ) THEN
./verification/hs94.1x64x5/code/ini_theta.F 125 INI_THETA
      _EXCH_XYZ_RL(theta , myThid )
./verification/hs94.1x64x5/code/ini_theta.F 128 INI_THETA
        CALL PLOT_FIELD_XYZRL( theta, 'Initial Temperature' ,
./verification/hs94.1x64x5/code/ini_theta.F 73 INI_THETA
            theta(I,J,K,bi,bj) = MAX( thetaLim, thetaEq )
./verification/hs94.1x64x5/code/ini_theta.F 84 INI_THETA
     U                       theta(1-OLx,1-OLy,1,bi,bj),
./verification/hs94.1x64x5/code/ini_theta.F 92 INI_THETA
       CALL READ_FLD_XYZ_RL( hydrogThetaFile, ' ', theta, 0, myThid )
./verification/hs94.1x64x5/code_ad/apply_forcing.F 290 APPLY_FORCING_T
     &               - kT*( theta(i,j,k,bi,bj)-thetaEq )
./verification/hs94.1x64x5/code_ad/cost_test.F 69 COST_TEST
     &              theta(i,j,kLocOut,bi,bj)
./verification/hs94.1x64x5/code_ad/ctrl_map_gentim2d.F 162 CTRL_MAP_GENTIM2D
             theta(I,J,1,bi,bj) = theta(I,J,1,bi,bj)
./verification/hs94.1x64x5/code_ad/ctrl_map_ini_genarr.F 104 CTRL_MAP_INI_GENARR
                theta(i,j,1,bi,bj) = theta(i,j,1,bi,bj)
./verification/hs94.1x64x5/code_ad/ctrl_map_ini_genarr.F 118 CTRL_MAP_INI_GENARR
                theta(i,j,1,bi,bj) = theta(i,j,1,bi,bj)
./verification/hs94.1x64x5/code_ad/ctrl_map_ini_genarr.F 136 CTRL_MAP_INI_GENARR
       _EXCH_XYZ_RL( theta, myThid )
./verification/hs94.1x64x5/code_ad/ctrl_map_ini_genarr.F 165 CTRL_MAP_INI_GENARR
                theta(i,j,k,bi,bj) = theta(i,j,k,bi,bj)
./verification/hs94.1x64x5/code_ad/ctrl_map_ini_genarr.F 173 CTRL_MAP_INI_GENARR
                theta(i,j,k,bi,bj) = theta(i,j,k,bi,bj)
./verification/hs94.1x64x5/code_ad/ctrl_map_ini_genarr.F 189 CTRL_MAP_INI_GENARR
       _EXCH_XYZ_RL( theta, myThid )
./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 108 INI_THETA
     &      .AND. theta(I,J,K,bi,bj).EQ.0. ) THEN
./verification/hs94.1x64x5/code_ad/ini_theta.F 125 INI_THETA
      _EXCH_XYZ_RL(theta , myThid )
./verification/hs94.1x64x5/code_ad/ini_theta.F 128 INI_THETA
        CALL PLOT_FIELD_XYZRL( theta, 'Initial Temperature' ,
./verification/hs94.1x64x5/code_ad/ini_theta.F 73 INI_THETA
            theta(I,J,K,bi,bj) = MAX( thetaLim, thetaEq )
./verification/hs94.1x64x5/code_ad/ini_theta.F 84 INI_THETA
     U                       theta(1-OLx,1-OLy,1,bi,bj),
./verification/hs94.1x64x5/code_ad/ini_theta.F 92 INI_THETA
       CALL READ_FLD_XYZ_RL( hydrogThetaFile, ' ', theta, 0, myThid )
./verification/hs94.1x64x5/code_oad/apply_forcing.F 290 APPLY_FORCING_T
     &               - kT*( theta(i,j,k,bi,bj)-thetaEq )
./verification/hs94.1x64x5/code_oad/cost_test.F 69 COST_TEST
     &              theta(i,j,kLocOut,bi,bj)
./verification/hs94.1x64x5/code_oad/ctrl_map_gentim2d.F 162 CTRL_MAP_GENTIM2D
             theta(I,J,1,bi,bj) = theta(I,J,1,bi,bj)
./verification/hs94.1x64x5/code_oad/ctrl_map_ini_genarr.F 103 CTRL_MAP_INI_GENARR
                theta(i,j,1,bi,bj) = theta(i,j,1,bi,bj)
./verification/hs94.1x64x5/code_oad/ctrl_map_ini_genarr.F 117 CTRL_MAP_INI_GENARR
                theta(i,j,1,bi,bj) = theta(i,j,1,bi,bj)
./verification/hs94.1x64x5/code_oad/ctrl_map_ini_genarr.F 135 CTRL_MAP_INI_GENARR
       _EXCH_XYZ_RL( theta, myThid )
./verification/hs94.1x64x5/code_oad/ctrl_map_ini_genarr.F 164 CTRL_MAP_INI_GENARR
                theta(i,j,k,bi,bj) = theta(i,j,k,bi,bj)
./verification/hs94.1x64x5/code_oad/ctrl_map_ini_genarr.F 172 CTRL_MAP_INI_GENARR
                theta(i,j,k,bi,bj) = theta(i,j,k,bi,bj)
./verification/hs94.1x64x5/code_oad/ctrl_map_ini_genarr.F 188 CTRL_MAP_INI_GENARR
       _EXCH_XYZ_RL( theta, myThid )
./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 108 INI_THETA
     &      .AND. theta(I,J,K,bi,bj).EQ.0. ) THEN
./verification/hs94.1x64x5/code_oad/ini_theta.F 125 INI_THETA
      _EXCH_XYZ_RL(theta , myThid )
./verification/hs94.1x64x5/code_oad/ini_theta.F 128 INI_THETA
        CALL PLOT_FIELD_XYZRL( theta, 'Initial Temperature' ,
./verification/hs94.1x64x5/code_oad/ini_theta.F 73 INI_THETA
            theta(I,J,K,bi,bj) = MAX( thetaLim, thetaEq )
./verification/hs94.1x64x5/code_oad/ini_theta.F 84 INI_THETA
     U                       theta(1-OLx,1-OLy,1,bi,bj),
./verification/hs94.1x64x5/code_oad/ini_theta.F 92 INI_THETA
       CALL READ_FLD_XYZ_RL( hydrogThetaFile, ' ', theta, 0, myThid )
./verification/hs94.cs-32x32x5/code/apply_forcing.F 290 APPLY_FORCING_T
     &               - kT*( theta(i,j,k,bi,bj)-thetaEq )
./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 108 INI_THETA
     &      .AND. theta(I,J,K,bi,bj).EQ.0. ) THEN
./verification/hs94.cs-32x32x5/code/ini_theta.F 125 INI_THETA
      _EXCH_XYZ_RL(theta , myThid )
./verification/hs94.cs-32x32x5/code/ini_theta.F 128 INI_THETA
        CALL PLOT_FIELD_XYZRL( theta, 'Initial Temperature' ,
./verification/hs94.cs-32x32x5/code/ini_theta.F 73 INI_THETA
            theta(I,J,K,bi,bj) = MAX( thetaLim, thetaEq )
./verification/hs94.cs-32x32x5/code/ini_theta.F 84 INI_THETA
     U                       theta(1-OLx,1-OLy,1,bi,bj),
./verification/hs94.cs-32x32x5/code/ini_theta.F 92 INI_THETA
       CALL READ_FLD_XYZ_RL( hydrogThetaFile, ' ', theta, 0, myThid )
./verification/internal_wave/code/obcs_calc.F 110 OBCS_CALC
     &          uVel, vVel, wVel, theta, salt,
./verification/internal_wave/code/obcs_calc.F 163 OBCS_CALC
     &          uVel, vVel, wVel, theta, salt,
./verification/internal_wave/code/obcs_calc.F 183 OBCS_CALC
     &          uVel, vVel, wVel, theta, salt,
./verification/internal_wave/code/obcs_calc.F 36 OBCS_CALC
      _RL theta(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
./verification/internal_wave/code/obcs_calc.F 7 OBCS_CALC
     &                      uVel, vVel, wVel, theta, salt,
./verification/internal_wave/code/obcs_calc.F 90 OBCS_CALC
     &          uVel, vVel, wVel, theta, salt,
./verification/isomip/code_ad/cost_test.F 74 COST_TEST
               objf_test(bi,bj) = theta(i,j,kLocOut,bi,bj)
./verification/isomip/code_ad/cost_test.F 77 COST_TEST
     &              theta(i,j,kLocOut,bi,bj)
./verification/isomip/code_ad/cost_test.F 84 COST_TEST
     &        ( theta(i,j,k,bi,bj) )**2
./verification/isomip/code_oad/cost_test.F 74 COST_TEST
               objf_test(bi,bj) = theta(i,j,kLocOut,bi,bj)
./verification/isomip/code_oad/cost_test.F 77 COST_TEST
     &              theta(i,j,kLocOut,bi,bj)
./verification/isomip/code_oad/cost_test.F 84 COST_TEST
     &        ( theta(i,j,k,bi,bj) )**2
./verification/rotating_tank/code/apply_forcing.F 523 APPLY_FORCING_T
     &          *( temp_addMass - theta(i,j,k,bi,bj) )
./verification/rotating_tank/code/apply_forcing.F 593 APPLY_FORCING_T
     &                  *( theta(i,j,km,bi,bj)*salt(i,j,km,bi,bj)
./verification/rotating_tank/code/apply_forcing.F 594 APPLY_FORCING_T
     &                   + theta(i,j,kc,bi,bj)*salt(i,j,kc,bi,bj)
./verification/rotating_tank/code/apply_forcing.F 606 APPLY_FORCING_T
     &                  *( theta(i,j,kc,bi,bj)*salt(i,j,kc,bi,bj)
./verification/rotating_tank/code/apply_forcing.F 607 APPLY_FORCING_T
     &                   + theta(i,j,kp,bi,bj)*salt(i,j,kp,bi,bj)
./verification/rotating_tank/code/apply_forcing.F 717 APPLY_FORCING_T
     &      -faceArea*kDiffCyl*(theta(i,j,k,bi,bj) - tCyl)
./verification/rotating_tank/code/apply_forcing.F 724 APPLY_FORCING_T
     &      -faceArea*kDiffCyl*(tCyl - theta(i,j,k,bi,bj))
./verification/rotating_tank/code/apply_forcing.F 731 APPLY_FORCING_T
     &      -faceArea*kDiffCyl*(tCyl-theta(i,j,k,bi,bj))
./verification/rotating_tank/code/apply_forcing.F 738 APPLY_FORCING_T
     &      -faceArea*kDiffCyl*(theta(i,j,k,bi,bj) - tCyl)
./model/src/read_pickup.F 175 READ_PICKUP
          CALL READ_REC_3D_RL( fn, fp, Nr, theta, 7, myIter,myThid )
./model/src/read_pickup.F 207 READ_PICKUP
          CALL READ_REC_3D_RL( fn, fp, Nr, theta,nj+1, myIter,myThid )
./model/src/read_pickup.F 224 READ_PICKUP
          CALL READ_REC_3D_RL( fn, fp, Nr, theta, 5, myIter,myThid )
./model/src/read_pickup.F 275 READ_PICKUP
          CALL READ_MFLDS_3D_RL( 'Theta   ', theta,
./model/src/read_pickup.F 465 READ_PICKUP
        CALL MNC_CW_RL_R('D',fn,0,0,'Temp',theta, myThid)
./model/src/read_pickup.F 510 READ_PICKUP
      CALL EXCH_3D_RL( theta, Nr, myThid )
./verification/tutorial_global_oce_optim/code_ad/external_forcing_surf.F 254 EXTERNAL_FORCING_SURF
     &          *( temp_EvPrRn - theta(i,j,ks,bi,bj) )
./verification/tutorial_global_oce_optim/code_ad/external_forcing_surf.F 289 EXTERNAL_FORCING_SURF
     &           *( theta(i,j,ks,bi,bj) - temp_EvPrRn )
./verification/tutorial_global_oce_optim/code_oad/external_forcing_surf.F 254 EXTERNAL_FORCING_SURF
     &          *( temp_EvPrRn - theta(i,j,ks,bi,bj) )
./verification/tutorial_global_oce_optim/code_oad/external_forcing_surf.F 289 EXTERNAL_FORCING_SURF
     &           *( theta(i,j,ks,bi,bj) - temp_EvPrRn )
./verification/tutorial_held_suarez_cs/code/apply_forcing.F 290 APPLY_FORCING_T
     &               - kT*( theta(i,j,k,bi,bj)-thetaEq )
./model/src/temp_integrate.F 192 TEMP_INTEGRATE
     I                           theta(1-OLx,1-OLy,1,bi,bj),
./model/src/temp_integrate.F 198 TEMP_INTEGRATE
     I                           theta(1-OLx,1-OLy,1,bi,bj),
./model/src/temp_integrate.F 266 TEMP_INTEGRATE
     I             dTtracerLev, uFld, vFld, wFld, theta,
./model/src/temp_integrate.F 280 TEMP_INTEGRATE
     I             dTtracerLev, uFld, vFld, wFld, theta,
./model/src/temp_integrate.F 347 TEMP_INTEGRATE
     I           theta(1-OLx,1-OLy,1,bi,bj),
./model/src/temp_integrate.F 362 TEMP_INTEGRATE
     I           theta(1-OLx,1-OLy,1,bi,bj),
./model/src/temp_integrate.F 451 TEMP_INTEGRATE
     I                  theta(1-OLx,1-OLy,1,bi,bj),
./model/src/temp_integrate.F 458 TEMP_INTEGRATE
     I                  theta(1-OLx,1-OLy,1,bi,bj),
./model/src/temp_integrate.F 469 TEMP_INTEGRATE
     I                  theta(1-OLx,1-OLy,1,bi,bj),
./model/src/temp_integrate.F 492 TEMP_INTEGRATE
     I         theta(1-OLx,1-OLy,1,bi,bj),
./model/src/temp_integrate.F 526 TEMP_INTEGRATE
     U             theta(1-OLx,1-OLy,1,bi,bj),
./model/src/temp_integrate.F 532 TEMP_INTEGRATE
     U             theta(1-OLx,1-OLy,1,bi,bj),
./model/src/temp_integrate.F 540 TEMP_INTEGRATE
     O             theta(1-OLx,1-OLy,1,bi,bj),
./model/src/thermodynamics.F 127 THERMODYNAMICS
       CALL CALC_WSURF_TR( theta, salt, wVel,
./model/src/thermodynamics.F 316 THERMODYNAMICS
           CALL OBCS_APPLY_TS( bi, bj, 0, theta, salt, myThid )
./model/src/thermodynamics.F 329 THERMODYNAMICS
       CALL DEBUG_STATS_RL(Nr,theta,'Theta (THERMODYNAMICS)',myThid)
./model/src/tracers_correction_step.F 70 TRACERS_CORRECTION_STEP
        CALL SHAP_FILT_APPLY_TS( theta,salt, myTime, myIter, myThid )
./model/src/tracers_correction_step.F 77 TRACERS_CORRECTION_STEP
        CALL ZONAL_FILT_APPLY_TS( theta, salt, myThid )
./model/src/tracers_iigw_correction.F 82 TRACERS_IIGW_CORRECTION
          theta(i,j,k,bi,bj) = theta(i,j,k,bi,bj)
./model/src/write_pickup.F 137 WRITE_PICKUP
        CALL WRITE_REC_3D_RL( fn, fp, Nr, theta,  -j, myIter, myThid )
./model/src/write_pickup.F 372 WRITE_PICKUP
        CALL MNC_CW_RL_W('D',fn,0,0,'Temp', theta, myThid)
./model/src/write_state.F 106 WRITE_STATE
          CALL WRITE_REC_XYZ_RL( 'T', theta,iRec,myIter,myThid)
./model/src/write_state.F 132 WRITE_STATE
          CALL WRITE_FLD_XYZ_RL( 'T.',suff,thetav,myIter,myThid)
./model/src/write_state.F 143 WRITE_STATE
          CALL WRITE_FLD_XYZ_RL( 'T.',suff,theta,myIter,myThid)
./model/src/write_state.F 190 WRITE_STATE
          CALL MNC_CW_RL_W(pf,'state',0,0,'Temp', theta, myThid)
./pkg/aim_v23/aim_dyn2aim.F 172 AIM_DYN2AIM
          TA(I2,Katm)  = theta(i,j,k,bi,bj)*conv_theta2T
./pkg/aim_v23/aim_dyn2aim.F 178 AIM_DYN2AIM
          ThA(I2,Katm) = theta(i,j,k,bi,bj)
./pkg/atm2d/get_ocnvars.F 72 GET_OCNVARS
          sstFromOcn(i,j) = theta(i,j,1,1,1)
./pkg/atm_phys/atm_phys_driver.F 211 ATM_PHYS_DRIVER
        t3d(:,:,k) = theta(1:sNx,1:sNy,kc,bi,bj)*conv_theta2T
./pkg/autodiff/autodiff_restore.F 145 AUTODIFF_RESTORE
           theta(I,J,K,bi,bj)    = StoreDynVars3D(I,J,K,bi,bj,5)
./pkg/autodiff/autodiff_store.F 155 AUTODIFF_STORE
           StoreDynVars3D(I,J,K,bi,bj,5)  = theta(I,J,K,bi,bj)
./pkg/bbl/bbl_calc_rhs.F 90 BBL_CALC_RHS
          tLoc(i,j)   = theta(i,j,kBot,bi,bj)
./pkg/bbl/bbl_init_varia.F 42 BBL_INIT_VARIA
          bbl_theta(i,j,bi,bj)     = theta(i,j,kBot,bi,bj)
./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 256 BULKF_FORCING
     &              + ( theta(i,j,ks,bi,bj) - temp_EvPrRn )
./pkg/bulk_force/bulkf_forcing.F 88 BULKF_FORCING
     I            theta(i,j,ks,bi,bj), flwdwn(i,j,bi,bj),
./pkg/cfc/cfc11_surfforcing.F 74 CFC11_SURFFORCING
     &          + sca_11_2 * theta(i,j,k,bi,bj)
./pkg/cfc/cfc11_surfforcing.F 75 CFC11_SURFFORCING
     &          + sca_11_3 * theta(i,j,k,bi,bj)*theta(i,j,k,bi,bj)
./pkg/cfc/cfc11_surfforcing.F 76 CFC11_SURFFORCING
     &          + sca_11_4 * theta(i,j,k,bi,bj)*theta(i,j,k,bi,bj)
./pkg/cfc/cfc11_surfforcing.F 77 CFC11_SURFFORCING
     &                *theta(i,j,k,bi,bj)
./pkg/cfc/cfc11_surfforcing.F 79 CFC11_SURFFORCING
             ttemp=( theta(i,j,k,bi,bj) + 273.16 _d 0)* 0.01 _d 0
./pkg/cfc/cfc12_surfforcing.F 74 CFC12_SURFFORCING
     &          + sca_12_2 * theta(i,j,k,bi,bj)
./pkg/cfc/cfc12_surfforcing.F 75 CFC12_SURFFORCING
     &          + sca_12_3 * theta(i,j,k,bi,bj)*theta(i,j,k,bi,bj)
./pkg/cfc/cfc12_surfforcing.F 76 CFC12_SURFFORCING
     &          + sca_12_4 * theta(i,j,k,bi,bj)*theta(i,j,k,bi,bj)
./pkg/cfc/cfc12_surfforcing.F 77 CFC12_SURFFORCING
     &                *theta(i,j,k,bi,bj)
./pkg/cfc/cfc12_surfforcing.F 79 CFC12_SURFFORCING
             ttemp=( theta(i,j,k,bi,bj) + 273.16 _d 0)* 0.01 _d 0
./pkg/cheapaml/cheapaml.F 158 CHEAPAML
                sst_htf(i,j) = theta(i,j,1,bi,bj) + celsius2K
./pkg/cheapaml/cheapaml.F 448 CHEAPAML
     I                      theta(1-OLx,1-OLy,1,bi,bj), windSq,
./pkg/cheapaml/cheapaml.F 486 CHEAPAML
     I                      theta(1-OLx,1-OLy,1,bi,bj), windSq,
./pkg/cheapaml/cheapaml.F 524 CHEAPAML
             xlwnet = 0.98 _d 0*stefan*(theta(i,j,1,bi,bj)+celsius2K)**4
./pkg/cheapaml/cheapaml.F 528 CHEAPAML
     &        + 4.0*0.98 _d 0*stefan*(theta(i,j,1,bi,bj)+celsius2K)**3
./pkg/cheapaml/cheapaml.F 529 CHEAPAML
     &          *(theta(i,j,1,bi,bj)-Tair(i,j,bi,bj))
./pkg/cheapaml/cheapaml_lanl_flux.F 127 CHEAPAML_LANL_FLUX
        to=theta(i,j,1,bi,bj)
./pkg/cost/cost_accumulate_mean.F 50 COST_ACCUMULATE_MEAN
     &                + theta(i,j,k,bi,bj)
./pkg/cost/cost_accumulate_mean.F 61 COST_ACCUMULATE_MEAN
     &               + (theta(i,j,k,bi,bj)+theta(i-1,j,k,bi,bj))
./pkg/cost/cost_accumulate_mean.F 67 COST_ACCUMULATE_MEAN
     &               + (theta(i,j,k,bi,bj)+theta(i,j-1,k,bi,bj))
./pkg/cost/cost_state_final.F 36 COST_STATE_FINAL
           objf_state_final(i,j,bi,bj,0*Nr+k) = theta(i,j,k,bi,bj)
./pkg/cost/cost_test.F 71 COST_TEST
               objf_test(bi,bj) = theta(i,j,kLocOut,bi,bj)
./pkg/cost/cost_test.F 74 COST_TEST
     &              theta(i,j,kLocOut,bi,bj)
./pkg/cost/cost_test.F 79 COST_TEST
     &                         theta(i,j,kLocOut,bi,bj)**2
./pkg/ctrl/ctrl_map_ini.F 120 CTRL_MAP_INI
                theta(i,j,k,bi,bj) = theta(i,j,k,bi,bj) +
./pkg/ctrl/ctrl_map_ini.F 124 CTRL_MAP_INI
                theta(i,j,k,bi,bj) = theta(i,j,k,bi,bj) +
./pkg/ctrl/ctrl_map_ini.F 128 CTRL_MAP_INI
                if(theta(i,j,k,bi,bj).lt.-2.0)
./pkg/ctrl/ctrl_map_ini.F 129 CTRL_MAP_INI
     &               theta(i,j,k,bi,bj)= -2.0
./pkg/ctrl/ctrl_map_ini.F 219 CTRL_MAP_INI
              theta(i,j,1,bi,bj) = theta(i,j,1,bi,bj)
./pkg/ctrl/ctrl_map_ini.F 583 CTRL_MAP_INI
      _EXCH_XYZ_RL( theta, mythid )
./pkg/ctrl/ctrl_map_ini_ecco.F 143 CTRL_MAP_INI_ECCO
     & xx_theta_dummy, boundsVec, theta, maskc, paramSmooth,
./pkg/ctrl/ctrl_map_ini_ecco.F 145 CTRL_MAP_INI_ECCO
      _EXCH_XYZ_RL( theta, mythid)
./pkg/ctrl/ctrl_map_ini_genarr.F 114 CTRL_MAP_INI_GENARR
       _EXCH_XYZ_RL( theta, myThid )
./pkg/ctrl/ctrl_map_ini_genarr.F 141 CTRL_MAP_INI_GENARR
                theta(i,j,k,bi,bj) = theta(i,j,k,bi,bj)
./pkg/ctrl/ctrl_map_ini_genarr.F 159 CTRL_MAP_INI_GENARR
       _EXCH_XYZ_RL( theta, myThid )
./pkg/ctrl/ctrl_map_ini_genarr.F 99 CTRL_MAP_INI_GENARR
                theta(i,j,1,bi,bj) = theta(i,j,1,bi,bj)
./pkg/diagnostics/diagnostics_fill_state.F 201 DIAGNOSTICS_FILL_STATE
     &                  *(theta(i,j,k,bi,bj)+theta(i-1,j,k,bi,bj))
./pkg/diagnostics/diagnostics_fill_state.F 217 DIAGNOSTICS_FILL_STATE
     &                  *(theta(i,j,k,bi,bj)+theta(i,j-1,k,bi,bj))
./pkg/diagnostics/diagnostics_fill_state.F 234 DIAGNOSTICS_FILL_STATE
     &                  *(theta(i,j,k,bi,bj)+theta(i,j,km1,bi,bj))
./pkg/diagnostics/diagnostics_fill_state.F 366 DIAGNOSTICS_FILL_STATE
               tmpMk(i,j,k,bi,bj) = tmpFac*theta(i,j,k,bi,bj)
./pkg/diagnostics/diagnostics_fill_state.F 399 DIAGNOSTICS_FILL_STATE
        CALL DIAGNOSTICS_FILL(theta,'THETA   ',0,Nr,0,1,1,myThid)
./pkg/diagnostics/diagnostics_fill_state.F 410 DIAGNOSTICS_FILL_STATE
            dummy1 = theta(i,j,k,bi,bj) * ((rC(k)/100.)/1000.)**kappa
./pkg/diagnostics/diagnostics_fill_state.F 424 DIAGNOSTICS_FILL_STATE
        CALL DIAGNOSTICS_SCALE_FILL( theta, oneRL, 2,
./pkg/diagnostics/diagnostics_fill_state.F 540 DIAGNOSTICS_FILL_STATE
     &                  *(theta(i,j,k,bi,bj)+theta(i-1,j,k,bi,bj))
./pkg/diagnostics/diagnostics_fill_state.F 557 DIAGNOSTICS_FILL_STATE
     &                  *(theta(i,j,k,bi,bj)+theta(i,j-1,k,bi,bj))
./pkg/diagnostics/diagnostics_fill_state.F 575 DIAGNOSTICS_FILL_STATE
     &                  *(theta(i,j,k,bi,bj)+theta(i,j,km1,bi,bj))
./pkg/diagnostics/diagnostics_fill_state.F 660 DIAGNOSTICS_FILL_STATE
               tmpMk(i,j,k,bi,bj) = tmpFac*theta(i,j,k,bi,bj)
./pkg/dic/calcite_saturation.F 109 CALCITE_SATURATION
     I          theta(i,j,k,bi,bj),salt(i,j,k,bi,bj),
./pkg/dic/calcite_saturation.F 76 CALCITE_SATURATION
     I                       theta,salt,
./pkg/dic/dic_surfforcing.F 113 DIC_SURFFORCING
     &                          min(50. _d 0, theta(i,j,kLev,bi,bj)))
./pkg/dic/dic_surfforcing.F 119 DIC_SURFFORCING
            surftemp(i,j) = theta(i,j,kLev,bi,bj)
./pkg/dic/dic_surfforcing.F 186 DIC_SURFFORCING
     &          + sca2 * theta(i,j,kLev,bi,bj)
./pkg/dic/dic_surfforcing.F 187 DIC_SURFFORCING
     &          + sca3 * theta(i,j,kLev,bi,bj)*theta(i,j,kLev,bi,bj)
./pkg/dic/dic_surfforcing.F 188 DIC_SURFFORCING
     &          + sca4 * theta(i,j,kLev,bi,bj)*theta(i,j,kLev,bi,bj)
./pkg/dic/dic_surfforcing.F 189 DIC_SURFFORCING
     &                *theta(i,j,kLev,bi,bj)
./pkg/dic/dic_surfforcing_init.F 195 DIC_SURFFORCING_INIT
     &                          min(50. _d 0, theta(i,j,kLev,bi,bj)))
./pkg/dic/dic_surfforcing_init.F 201 DIC_SURFFORCING_INIT
            surftemp(i,j) = theta(i,j,kLev,bi,bj)
./pkg/dic/dic_surfforcing_init.F 250 DIC_SURFFORCING_INIT
     &        theta(iprt,jprt,1,bi,bj), salt(iprt,jprt,1,bi,bj),
./pkg/dic/o2_surfforcing.F 73 O2_SURFFORCING
              ttemp = theta(i,j,k,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 85 EBM_AREA_T
     &                       + maskC(i,j,k,bi,bj)*theta(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 97 EBM_AREA_T
     &                       + maskC(i,j,k,bi,bj)*theta(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_atlantic.F 165 COST_ATLANTIC
                       tv = tv+p5*(theta(i,j,k,bi,bj)
./pkg/ecco/cost_atlantic.F 166 COST_ATLANTIC
     $                     + theta(i,j-1,k,bi,bj))*vVel(i,j,k,bi,bj)
./pkg/ecco/cost_atlantic.F 213 COST_ATLANTIC
     $                     +p5*(theta(i,j,k,bi,bj)
./pkg/ecco/cost_atlantic.F 214 COST_ATLANTIC
     $                     + theta(i,j-1,k,bi,bj))*vVel(i,j,k,bi,bj)
./pkg/ecco/cost_averagesfields.F 260 COST_AVERAGESFIELDS
     &     tbar, theta, xx_tbar_mean_dummy,
./pkg/ecco/cost_averagesfields.F 267 COST_AVERAGESFIELDS
     &     tbar, theta(1-Olx,1-Oly,1,1,1), xx_tbar_mean_dummy,
./pkg/ecco/cost_averagesfields.F 282 COST_AVERAGESFIELDS
     &     sstbar, theta(1-Olx,1-Oly,1,1,1), xx_sstbar_mean_dummy,
./pkg/ecco/cost_averagesfields.F 554 COST_AVERAGESFIELDS
     &     smrsstbar, theta(1-Olx,1-Oly,1,1,1),
./pkg/ecco/cost_gencost_customize.F 77 COST_GENCOST_CUSTOMIZE
     &        = theta(i,j,1,bi,bj)*maskC(i,j,1,bi,bj)
./model/inc/DYNVARS.h 38 NO PROCEDURE
     &                   uVel,vVel,wVel,theta,salt,
./model/inc/DYNVARS.h 44 NO PROCEDURE
     &                   uVel,vVel,wVel,theta,salt,
./model/inc/DYNVARS.h 52 NO PROCEDURE
      _RL  theta(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
./pkg/ecco/cost_theta_ini_fin.F 110 COST_THETA_INI_FIN
                  tmpx = tmpfld3d(i,j,k,bi,bj) - theta(i,j,k,bi,bj)
./pkg/ecco/ecco_phys.F 65 ECCO_PHYS
     I                theta(1-OLx,1-OLy,k,bi,bj),