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



File Line number Procedure Code
./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 152 FLT_TRAJ
              CALL FLT_BILINEAR  (ix,jy,tt,theta, kp,0,bi,bj,myThid)
./pkg/flt/flt_up.F 145 FLT_UP
               CALL FLT_BILINEAR  (ix,jy,tt,theta, k,0,bi,bj,myThid)
./pkg/ggl90/ggl90_calc.F 162 GGL90_CALC
     I      theta(1-OLx,1-OLy,km1,bi,bj), salt(1-OLx,1-OLy,km1,bi,bj),
./pkg/ggl90/ggl90_calc.F 168 GGL90_CALC
     I      theta(1-OLx,1-OLy,k,bi,bj), salt(1-OLx,1-OLy,k,bi,bj),
./pkg/gmredi/gmredi_calc_tensor.F 692 GMREDI_CALC_TENSOR
     &           (theta(i-1,j,km1,bi,bj)-theta(i-1,j,k,bi,bj))
./pkg/gmredi/gmredi_calc_tensor.F 694 GMREDI_CALC_TENSOR
     &           (theta( i ,j,km1,bi,bj)-theta( i ,j,k,bi,bj))
./pkg/gmredi/gmredi_calc_tensor.F 698 GMREDI_CALC_TENSOR
     &           (theta(i-1,j,k,bi,bj)-theta(i-1,j,kp1,bi,bj))
./pkg/gmredi/gmredi_calc_tensor.F 700 GMREDI_CALC_TENSOR
     &           (theta( i ,j,k,bi,bj)-theta( i ,j,kp1,bi,bj))
./pkg/gmredi/gmredi_calc_tensor.F 880 GMREDI_CALC_TENSOR
     &           (theta(i,j-1,km1,bi,bj)-theta(i,j-1,k,bi,bj))
./pkg/gmredi/gmredi_calc_tensor.F 882 GMREDI_CALC_TENSOR
     &           (theta(i, j ,km1,bi,bj)-theta(i, j ,k,bi,bj))
./pkg/gmredi/gmredi_calc_tensor.F 886 GMREDI_CALC_TENSOR
     &           (theta(i,j-1,k,bi,bj)-theta(i,j-1,kp1,bi,bj))
./pkg/gmredi/gmredi_calc_tensor.F 888 GMREDI_CALC_TENSOR
     &           (theta(i, j ,k,bi,bj)-theta(i, j ,kp1,bi,bj))
./pkg/icefront/icefront_thermodynamics.F 107 ICEFRONT_THERMODYNAMICS
            tLoc = theta(I,J,K,bi,bj)
./pkg/kpp/kpp_routines.F 1554 STATEKPP
     I     theta(1-OLx,1-OLy,1,bi,bj), salt(1-OLx,1-OLy,1,bi,bj),
./pkg/kpp/kpp_routines.F 1596 STATEKPP
     I        theta(1-OLx,1-OLy,k,bi,bj), salt(1-OLx,1-OLy,k,bi,bj),
./pkg/kpp/kpp_routines.F 1608 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 1620 STATEKPP
     I        theta(1-OLx,1-OLy,1,bi,bj), salt(1-OLx,1-OLy,1,bi,bj),
./pkg/kpp/kpp_routines.F 1769 KPP_DOUBLEDIFF
         alphaDT(I,J) = ( theta(I,J,Km1,bi,bj)-theta(I,J,K,bi,bj) )
./pkg/layers/layers_calc.F 118 LAYERS_CALC
     &    0.5 _d 0 * (theta(i-1,j,k,bi,bj)+theta(i,j,k,bi,bj)) +
./pkg/layers/layers_calc.F 120 LAYERS_CALC
     &    0.5 _d 0 * (theta(i-1,j,k+1,bi,bj)+theta(i,j,k+1,bi,bj))
./pkg/layers/layers_calc.F 190 LAYERS_CALC
     &    0.5 _d 0 * (theta(i,j-1,k,bi,bj)+theta(i,j,k,bi,bj)) +
./pkg/layers/layers_calc.F 192 LAYERS_CALC
     &    0.5 _d 0 * (theta(i,j-1,k+1,bi,bj)+theta(i,j,k+1,bi,bj))
./pkg/longstep/longstep_average.F 93 LONGSTEP_AVERAGE
      CALL LONGSTEP_FILL_3D(LS_thetaCount, LS_theta, theta, Nr, myThid)
./model/src/calc_buoyancy.F 41 CALC_BUOYANCY
         buoy(i,j)=(theta(i,j,k,bi,bj)-tRef(k))/tRef(k)
./model/src/calc_gt.F 166 CALC_GT
     I           gtNm(1-Olx,1-Oly,1,1,1,m2), theta, dTtracerLev,
./model/src/calc_gt.F 177 CALC_GT
     I           diffKhT, diffK4T, KappaRT, gtNm1, theta, dTtracerLev,
./pkg/monitor/mon_solution.F 45 MON_SOLUTION
     I                Nr, theta, hFacC, maskInC, rA, drF,
./model/src/calc_oce_mxlayer.F 111 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 157 CALC_OCE_MXLAYER
     I        theta(1-OLx,1-OLy,k,bi,bj), salt(1-OLx,1-OLy,k,bi,bj),
./pkg/monitor/mon_surfcor.F 128 MON_SURFCOR
     &           *(theta(i,j,k,bi,bj)+theta(i,j,k-1,bi,bj))*0.5 _d 0
./pkg/monitor/mon_surfcor.F 149 MON_SURFCOR
              tmpVal  = rStarDhCDt(i,j,bi,bj)*theta(i,j,k,bi,bj)
./pkg/monitor/mon_surfcor.F 87 MON_SURFCOR
     &          rA(i,j,bi,bj)*wVel(i,j,ks,bi,bj)*theta(i,j,ks,bi,bj)
./pkg/monitor/monitor.F 116 MONITOR
        CALL MON_WRITESTATS_RL( Nr, theta,'_theta',
./pkg/monitor/monitor.F 125 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 141 OBCS_CALC
     &          uVel, vVel, wVel, theta, salt,
./pkg/obcs/obcs_calc.F 193 OBCS_CALC
     &          uVel, vVel, wVel, theta, salt,
./pkg/obcs/obcs_calc.F 228 OBCS_CALC
     &          uVel, vVel, wVel, theta, salt,
./pkg/obcs/obcs_calc.F 40 OBCS_CALC
      _RL theta(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
./pkg/obcs/obcs_calc.F 7 OBCS_CALC
     &                      uVel, vVel, wVel, theta, salt,
./pkg/obcs/obcs_calc.F 89 OBCS_CALC
     &          uVel, vVel, wVel, theta, salt,
./pkg/obcs/obcs_calc_stevens.F 205 OBCS_CALC_STEVENS
            dTracSpace = (theta(I-1,J,K,bi,bj)-theta(I-2,J,K,bi,bj))
./pkg/obcs/obcs_calc_stevens.F 223 OBCS_CALC_STEVENS
            OBEt(J,K,bi,bj) = theta(I,J,K,bi,bj) + dTtracerLev(K)
./pkg/obcs/obcs_calc_stevens.F 227 OBCS_CALC_STEVENS
     &          *(theta(I,J,K,bi,bj)-theta(I-1,J,K,bi,bj))
./pkg/obcs/obcs_calc_stevens.F 229 OBCS_CALC_STEVENS
     &          - gFac * gammat * (theta(I,J,K,bi,bj)-OBEt(J,K,bi,bj)) )
./pkg/obcs/obcs_calc_stevens.F 352 OBCS_CALC_STEVENS
            dTracSpace = (theta(I+2,J,K,bi,bj)-theta(I+1,J,K,bi,bj))
./pkg/obcs/obcs_calc_stevens.F 370 OBCS_CALC_STEVENS
            OBWt(J,K,bi,bj) = theta(I,J,K,bi,bj) + dTtracerLev(K)
./pkg/obcs/obcs_calc_stevens.F 374 OBCS_CALC_STEVENS
     &        *(theta(I+1,J,K,bi,bj)-theta(I,J,K,bi,bj))
./pkg/obcs/obcs_calc_stevens.F 376 OBCS_CALC_STEVENS
     &        - gFac * gammat * (theta(I,J,K,bi,bj)-OBWt(J,K,bi,bj)) )
./pkg/obcs/obcs_init_variables.F 343 OBCS_INIT_VARIABLES
     &              uVel, vVel, wVel, theta, salt, myThid )
./pkg/obcs/obcs_init_variables.F 354 OBCS_INIT_VARIABLES
          CALL OBCS_APPLY_TS( bi, bj, 0, theta, salt, myThid )
./pkg/obcs/obcs_init_variables.F 366 OBCS_INIT_VARIABLES
        _EXCH_XYZ_RL( theta, myThid )
./pkg/obcs/obcs_init_variables.F 382 OBCS_INIT_VARIABLES
         CALL OBCS_COPY_TRACER( theta(1-Olx,1-Oly,1,bi,bj),
./pkg/obcs/obcs_sponge.F 458 OBCS_SPONGE_T
     &    + float(jsl)*theta(i,j,kLev,bi,bj) )
./pkg/obcs/obcs_sponge.F 473 OBCS_SPONGE_T
     &    * ( theta(i,j,kLev,bi,bj) - trelax )
./pkg/obcs/obcs_sponge.F 498 OBCS_SPONGE_T
     &    + float(jsl)*theta(i,j,kLev,bi,bj) )
./pkg/obcs/obcs_sponge.F 514 OBCS_SPONGE_T
     &   * ( theta(i,j,kLev,bi,bj) - trelax )
./pkg/obcs/obcs_sponge.F 539 OBCS_SPONGE_T
     &    + float(isl)*theta(i,j,kLev,bi,bj) )
./pkg/obcs/obcs_sponge.F 555 OBCS_SPONGE_T
     &    * ( theta(i,j,kLev,bi,bj) - trelax )
./pkg/obcs/obcs_sponge.F 582 OBCS_SPONGE_T
     &    + float(isl)*theta(i,j,kLev,bi,bj) )
./pkg/obcs/obcs_sponge.F 598 OBCS_SPONGE_T
     &    * ( theta(i,j,kLev,bi,bj) - trelax )
./pkg/obcs/orlanski_east.F 165 ORLANSKI_EAST
                  CL=-(theta(I_obc-1,J,K,bi,bj)-TE_STORE_1(J,K,bi,bj))/
./pkg/obcs/orlanski_east.F 182 ORLANSKI_EAST
               OBEt(J,K,bi,bj)=theta(I_obc,J,K,bi,bj)-
./pkg/obcs/orlanski_east.F 184 ORLANSKI_EAST
     &           (ab1*(theta(I_obc,J,K,bi,bj)-theta(I_obc-1,J,K,bi,bj))+
./pkg/obcs/orlanski_east.F 263 ORLANSKI_EAST
               TE_STORE_2(J,K,bi,bj)=theta(I_obc-1,J,K,bi,bj) -
./pkg/obcs/orlanski_east.F 264 ORLANSKI_EAST
     &         theta(I_obc-2,J,K,bi,bj)
./pkg/obcs/orlanski_east.F 265 ORLANSKI_EAST
               TE_STORE_1(J,K,bi,bj)=theta(I_obc-1,J,K,bi,bj)
./pkg/obcs/orlanski_east.F 266 ORLANSKI_EAST
               TE_STORE_4(J,K,bi,bj)=theta(I_obc,J,K,bi,bj)
./pkg/obcs/orlanski_east.F 83 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 144 ORLANSKI_NORTH
                  CL=-(theta(I,J_obc-1,K,bi,bj)-TN_STORE_1(I,K,bi,bj))/
./pkg/obcs/orlanski_north.F 155 ORLANSKI_NORTH
               OBNt(I,K,bi,bj)=theta(I,J_obc,K,bi,bj)-
./pkg/obcs/orlanski_north.F 157 ORLANSKI_NORTH
     &          (ab1*(theta(I,J_obc,K,bi,bj)-theta(I,J_obc-1,K,bi,bj))+
./pkg/obcs/orlanski_north.F 224 ORLANSKI_NORTH
               TN_STORE_2(I,K,bi,bj)=theta(I,J_obc-1,K,bi,bj) -
./pkg/obcs/orlanski_north.F 225 ORLANSKI_NORTH
     &         theta(I,J_obc-2,K,bi,bj)
./pkg/obcs/orlanski_north.F 226 ORLANSKI_NORTH
               TN_STORE_1(I,K,bi,bj)=theta(I,J_obc-1,K,bi,bj)
./pkg/obcs/orlanski_north.F 227 ORLANSKI_NORTH
               TN_STORE_4(I,K,bi,bj)=theta(I,J_obc,K,bi,bj)
./pkg/obcs/orlanski_north.F 74 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 144 ORLANSKI_SOUTH
                  CL=(theta(I,J_obc+1,K,bi,bj)-TS_STORE_1(I,K,bi,bj))/
./pkg/obcs/orlanski_south.F 155 ORLANSKI_SOUTH
               OBSt(I,K,bi,bj)=theta(I,J_obc,K,bi,bj)+
./pkg/obcs/orlanski_south.F 157 ORLANSKI_SOUTH
     &          (ab1*(theta(I,J_obc+1,K,bi,bj)-theta(I,J_obc,K,bi,bj))+
./pkg/obcs/orlanski_south.F 224 ORLANSKI_SOUTH
               TS_STORE_2(I,K,bi,bj)=theta(I,J_obc+2,K,bi,bj) -
./pkg/obcs/orlanski_south.F 225 ORLANSKI_SOUTH
     &         theta(I,J_obc+1,K,bi,bj)
./pkg/obcs/orlanski_south.F 226 ORLANSKI_SOUTH
               TS_STORE_1(I,K,bi,bj)=theta(I,J_obc+1,K,bi,bj)
./pkg/obcs/orlanski_south.F 227 ORLANSKI_SOUTH
               TS_STORE_4(I,K,bi,bj)=theta(I,J_obc,K,bi,bj)
./pkg/obcs/orlanski_south.F 74 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 167 ORLANSKI_WEST
                  CL=(theta(I_obc+1,J,K,bi,bj)-TW_STORE_1(J,K,bi,bj))/
./pkg/obcs/orlanski_west.F 184 ORLANSKI_WEST
               OBWt(J,K,bi,bj)=theta(I_obc,J,K,bi,bj)+
./pkg/obcs/orlanski_west.F 186 ORLANSKI_WEST
     &          (ab1*(theta(I_obc+1,J,K,bi,bj)-theta(I_obc,J,K,bi,bj))+
./pkg/obcs/orlanski_west.F 265 ORLANSKI_WEST
               TW_STORE_2(J,K,bi,bj)=theta(I_obc+2,J,K,bi,bj) -
./pkg/obcs/orlanski_west.F 266 ORLANSKI_WEST
     &         theta(I_obc+1,J,K,bi,bj)
./pkg/obcs/orlanski_west.F 267 ORLANSKI_WEST
               TW_STORE_1(J,K,bi,bj)=theta(I_obc+1,J,K,bi,bj)
./pkg/obcs/orlanski_west.F 268 ORLANSKI_WEST
               TW_STORE_4(J,K,bi,bj)=theta(I_obc,J,K,bi,bj)
./pkg/obcs/orlanski_west.F 85 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 53 OCN_STORE_MY_DATA
          SSTocn2cpl(i,j,bi,bj) = theta(i,j,1,bi,bj)
./pkg/offline/offline_fields_load.F 346 OFFLINE_FIELDS_LOAD
          theta(i,j,k,bi,bj)    = bWght*tave0(i,j,k,bi,bj)
./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)
./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),
./model/src/convective_adjustment.F 145 CONVECTIVE_ADJUSTMENT
     I           theta(1-OLx,1-OLy,K-1,bi,bj),
./model/src/convective_adjustment.F 159 CONVECTIVE_ADJUSTMENT
     I           theta(1-OLx,1-OLy,K,bi,bj),
./model/src/convective_adjustment.F 185 CONVECTIVE_ADJUSTMENT
     U                              theta,
./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/convective_adjustment_ini.F 145 CONVECTIVE_ADJUSTMENT_INI
     I           theta(1-OLx,1-OLy,K-1,bi,bj),
./model/src/convective_adjustment_ini.F 159 CONVECTIVE_ADJUSTMENT_INI
     I           theta(1-OLx,1-OLy,K,bi,bj),
./model/src/convective_adjustment_ini.F 183 CONVECTIVE_ADJUSTMENT_INI
     U                              theta,
./pkg/rbcs/rbcs_add_tendency.F 60 RBCS_ADD_TENDENCY
     &        *( theta(i,j,k,bi,bj)- RBCtemp(i,j,k,bi,bj) )
./model/src/diags_oceanic_surf_flux.F 130 DIAGS_OCEANIC_SURF_FLUX
     &       + PmEpR(i,j,bi,bj)*theta(i,j,ks,bi,bj)*HeatCapacity_Cp
./pkg/salt_plume/salt_plume_calc_depth.F 79 SALT_PLUME_CALC_DEPTH
     I       theta(1-OLx,1-OLy,K,bi,bj), salt(1-OLx,1-OLy,K,bi,bj),
./model/src/do_fields_blocking_exchanges.F 63 DO_FIELDS_BLOCKING_EXCHANGES
       _EXCH_XYZ_RL( theta, myThid )
./model/src/do_fields_blocking_exchanges.F 70 DO_FIELDS_BLOCKING_EXCHANGES
         CALL OBCS_COPY_TRACER( theta(1-Olx,1-Oly,1,bi,bj),
./pkg/sbo/sbo_calc.F 204 SBO_CALC
     I           theta(1-OLx,1-OLy,kn0,bi,bj),
./model/src/do_oceanic_phys.F 158 DO_OCEANIC_PHYS
     I                 uVel, vVel, wVel, theta, salt, myThid )
./model/src/do_oceanic_phys.F 473 DO_OCEANIC_PHYS
     I                  theta, salt,
./model/src/do_oceanic_phys.F 480 DO_OCEANIC_PHYS
     I                theta(1-OLx,1-OLy,k,bi,bj),
./model/src/do_oceanic_phys.F 513 DO_OCEANIC_PHYS
     I                 theta(1-OLx,1-OLy,k-1,bi,bj),
./pkg/seaice/cost_ice_test.F 91 COST_ICE_TEST
      print *, 'ph-ice B ', myiter, theta(4,4,kSrf,1,1),
./model/src/do_the_model_io.F 63 DO_THE_MODEL_IO
          CALL PLOT_FIELD_XYZRL( theta, 'Current theta ',
./pkg/seaice/seaice_growth.F 502 SEAICE_GROWTH
          TMIX(I,J,bi,bj)=theta(I,J,kSurface,bi,bj)+celsius2K
./pkg/seaice/seaice_growth.F 710 SEAICE_GROWTH
           IF ( theta(I,J,kSurface,bi,bj) .GE. TBC ) THEN
./pkg/seaice/seaice_growth.F 712 SEAICE_GROWTH
     &             * (theta(I,J,kSurface,bi,bj)-TBC) * dRf(kSurface)
./pkg/seaice/seaice_growth.F 717 SEAICE_GROWTH
     &             * (theta(I,J,kSurface,bi,bj)-TBC) * dRf(kSurface)
./pkg/seaice/seaice_growth_if.F 439 SEAICE_GROWTH_IF
     &                 theta(I,J,kSurface,bi,bj) + TMELT
./pkg/seaice/seaice_growth_if.F 622 SEAICE_GROWTH_IF
                  surf_theta = max(theta(I,J,kSurface,bi,bj), TBC)
./model/src/dynamics.F 492 DYNAMICS
     I        theta, salt,
./model/src/dynamics.F 697 DYNAMICS
       CALL DEBUG_STATS_RL(Nr,theta,'Theta (DYNAMICS)',myThid)
./model/src/external_forcing.F 298 EXTERNAL_FORCING_T
     &          *( temp_addMass - theta(i,j,kLev,bi,bj) )
./pkg/shelfice/shelfice_thermodynamics.F 122 SHELFICE_THERMODYNAMICS
          tLoc(I,J) = theta(I,J,K,bi,bj)
./pkg/shelfice/shelfice_thermodynamics.F 140 SHELFICE_THERMODYNAMICS
     &           + theta(I,J,Kp1,bi,bj) *drKp1 )
./model/src/external_forcing_surf.F 103 EXTERNAL_FORCING_SURF
     &         *(theta(i,j,ks,bi,bj)-SST(i,j,bi,bj))
./model/src/external_forcing_surf.F 280 EXTERNAL_FORCING_SURF
     &          *( temp_EvPrRn - theta(i,j,ks,bi,bj) )
./model/src/external_forcing_surf.F 315 EXTERNAL_FORCING_SURF
     &           *( theta(i,j,ks,bi,bj) - temp_EvPrRn )
./model/src/external_forcing_surf.F 87 EXTERNAL_FORCING_SURF
     &         *(theta(i,j,ks,bi,bj)-SST(i,j,bi,bj))
./model/src/find_alpha.F 121 FIND_ALPHA
     I        theta(1-OLx,1-OLy,k,bi,bj), salt(1-OLx,1-OLy,k,bi,bj),
./model/src/find_alpha.F 127 FIND_ALPHA
     I        theta(1-OLx,1-OLy,k,bi,bj), salt(1-OLx,1-OLy,k,bi,bj),
./model/src/find_alpha.F 135 FIND_ALPHA
               t1 = theta(i,j,k,bi,bj)
./model/src/find_alpha.F 225 FIND_ALPHA
     I        theta(1-OLx,1-OLy,k,bi,bj), salt(1-OLx,1-OLy,k,bi,bj),
./model/src/find_alpha.F 231 FIND_ALPHA
     I        theta(1-OLx,1-OLy,k,bi,bj), salt(1-OLx,1-OLy,k,bi,bj),
./model/src/find_alpha.F 237 FIND_ALPHA
               t1  = theta(i,j,k,bi,bj)
./model/src/find_alpha.F 350 FIND_BETA
               tP=theta(i,j,k,bi,bj)-refTemp
./model/src/find_alpha.F 381 FIND_BETA
     I        theta(1-OLx,1-OLy,k,bi,bj), salt(1-OLx,1-OLy,k,bi,bj),
./model/src/find_alpha.F 387 FIND_BETA
     I        theta(1-OLx,1-OLy,k,bi,bj), salt(1-OLx,1-OLy,k,bi,bj),
./model/src/find_alpha.F 395 FIND_BETA
               t1 = theta(i,j,k,bi,bj)
./model/src/find_alpha.F 471 FIND_BETA
     I        theta(1-OLx,1-OLy,k,bi,bj), salt(1-OLx,1-OLy,k,bi,bj),
./model/src/find_alpha.F 477 FIND_BETA
     I        theta(1-OLx,1-OLy,k,bi,bj), salt(1-OLx,1-OLy,k,bi,bj),
./model/src/find_alpha.F 483 FIND_BETA
               t1  = theta(i,j,k,bi,bj)
./model/src/find_alpha.F 88 FIND_ALPHA
               tP=theta(i,j,k,bi,bj)-refTemp
./pkg/thsice/thsice_main.F 156 THSICE_MAIN
          tOceMxL(i,j,bi,bj) = theta(i,j,1,bi,bj)
./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/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 82 TIMEAVE_SURF_FLUX
     &                        *theta(i,j,k,bi,bj)*HeatCapacity_Cp
./verification/OpenAD/code_oad_all/openad_dumpAdjoint.F 91 OPENAD_DUMPADJOINT
         foo5=thetad
./verification/OpenAD/code_oad_all/write_state.F 106 WRITE_STATE
          CALL WRITE_REC_XYZ_RL( 'T', theta,iRec,myIter,myThid)
./verification/OpenAD/code_oad_all/write_state.F 129 WRITE_STATE
          CALL WRITE_FLD_XYZ_RL( 'T.',suff,thetav,myIter,myThid)
./verification/OpenAD/code_oad_all/write_state.F 174 WRITE_STATE
          CALL MNC_CW_RL_W(pf,'state',0,0,'Temp', theta, myThid)
./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_dynvars.F 54 INI_DYNVARS
           theta(I,J,K,bi,bj)=0. _d 0
./verification/dome/code/obcs_calc.F 126 OBCS_CALC
     &          uVel, vVel, wVel, theta, salt,
./verification/dome/code/obcs_calc.F 145 OBCS_CALC
     &          uVel, vVel, wVel, theta, salt,
./verification/dome/code/obcs_calc.F 206 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 94 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,
./model/src/ini_pressure.F 114 INI_PRESSURE
     I          theta, salt,
./verification/global_ocean.cs32x15/code_ad/cost_test.F 77 COST_TEST
               objf_test(bi,bj) = theta(i,j,kLocOut,bi,bj)
./verification/global_ocean.cs32x15/code_ad/cost_test.F 80 COST_TEST
     &              theta(i,j,kLocOut,bi,bj)
./verification/global_ocean.cs32x15/code_ad/cost_test.F 87 COST_TEST
     &        (theta(i,j,k,bi,bj)-thetaLev(i,j,k,bi,bj))**2
./verification/hs94.128x64x5/code/external_forcing.F 206 EXTERNAL_FORCING_T
     &        - kT*( theta(i,j,kLev,bi,bj)-thetaEq )
./verification/hs94.1x64x5/code/external_forcing.F 206 EXTERNAL_FORCING_T
     &        - kT*( theta(i,j,kLev,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 )
./model/src/ini_theta.F 102 INI_THETA
     &      .AND. theta(i,j,k,bi,bj).EQ.0. ) THEN
./model/src/ini_theta.F 139 INI_THETA
            IF (theta(i,j,k,bi,bj) .LT. Tfreezing) THEN
./model/src/ini_theta.F 140 INI_THETA
               theta(i,j,k,bi,bj) = Tfreezing
./model/src/ini_theta.F 150 INI_THETA
        CALL PLOT_FIELD_XYZRL( theta, 'Initial Temperature',
./model/src/ini_theta.F 61 INI_THETA
           theta(i,j,k,bi,bj) = tRef(k)
./model/src/ini_theta.F 75 INI_THETA
          CALL MNC_CW_RL_R('D',hydrogThetaFile,0,0,'Temp',theta,myThid)
./model/src/ini_theta.F 79 INI_THETA
          CALL READ_FLD_XYZ_RL( hydrogThetaFile, ' ', theta, 0, myThid )
./model/src/ini_theta.F 83 INI_THETA
        _EXCH_XYZ_RL(theta,myThid)
./model/src/ini_theta.F 94 INI_THETA
            IF (maskC(i,j,k,bi,bj).EQ.0.) theta(i,j,k,bi,bj) = 0.
./verification/hs94.1x64x5/code_ad/cost_test.F 75 COST_TEST
     &              theta(i,j,kLocOut,bi,bj)
./verification/hs94.1x64x5/code_ad/external_forcing.F 206 EXTERNAL_FORCING_T
     &        - kT*( theta(i,j,kLev,bi,bj)-thetaEq )
./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.cs-32x32x5/code/external_forcing.F 254 EXTERNAL_FORCING_T
     &        - kT*( theta(i,j,kLev,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 101 OBCS_CALC
     &          uVel, vVel, wVel, theta, salt, 
./verification/internal_wave/code/obcs_calc.F 154 OBCS_CALC
     &          uVel, vVel, wVel, theta, salt, 
./verification/internal_wave/code/obcs_calc.F 174 OBCS_CALC
     &          uVel, vVel, wVel, theta, salt, 
./verification/internal_wave/code/obcs_calc.F 31 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 81 OBCS_CALC
     &          uVel, vVel, wVel, theta, salt, 
./verification/rotating_tank/code/external_forcing.F 386 EXTERNAL_FORCING_T
     &      -faceArea*kDiffCyl*(theta(i,j,kLev,bi,bj) - tCyl)
./verification/rotating_tank/code/external_forcing.F 393 EXTERNAL_FORCING_T
     &      -faceArea*kDiffCyl*(tCyl - theta(i,j,kLev,bi,bj))
./verification/rotating_tank/code/external_forcing.F 400 EXTERNAL_FORCING_T
     &      -faceArea*kDiffCyl*(tCyl-theta(i,j,kLev,bi,bj))
./verification/rotating_tank/code/external_forcing.F 407 EXTERNAL_FORCING_T
     &      -faceArea*kDiffCyl*(theta(i,j,kLev,bi,bj) - tCyl)
./verification/tutorial_global_oce_optim/code_ad/external_forcing_surf.F 103 EXTERNAL_FORCING_SURF
     &         *(theta(i,j,ks,bi,bj)-SST(i,j,bi,bj))
./verification/tutorial_global_oce_optim/code_ad/external_forcing_surf.F 266 EXTERNAL_FORCING_SURF
     &          *( temp_EvPrRn - theta(i,j,ks,bi,bj) )
./verification/tutorial_global_oce_optim/code_ad/external_forcing_surf.F 301 EXTERNAL_FORCING_SURF
     &           *( theta(i,j,ks,bi,bj) - temp_EvPrRn )
./verification/tutorial_global_oce_optim/code_ad/external_forcing_surf.F 87 EXTERNAL_FORCING_SURF
     &         *(theta(i,j,ks,bi,bj)-SST(i,j,bi,bj))
./verification/tutorial_held_suarez_cs/code/external_forcing.F 206 EXTERNAL_FORCING_T
     &        - kT*( theta(i,j,kLev,bi,bj)-thetaEq )
./model/src/read_pickup.F 169 READ_PICKUP
          CALL READ_REC_3D_RL( fn, fp, Nr, theta, 7, myIter,myThid )
./model/src/read_pickup.F 201 READ_PICKUP
          CALL READ_REC_3D_RL( fn, fp, Nr, theta,nj+1, myIter,myThid )
./model/src/read_pickup.F 218 READ_PICKUP
          CALL READ_REC_3D_RL( fn, fp, Nr, theta, 5, myIter,myThid )
./model/src/read_pickup.F 261 READ_PICKUP
          CALL READ_MFLDS_3D_RL( 'Theta   ', theta,
./model/src/read_pickup.F 439 READ_PICKUP
        CALL MNC_CW_RL_R('D',fn,0,0,'Temp',theta, myThid)
./model/src/read_pickup.F 479 READ_PICKUP
      CALL EXCH_3D_RL( theta, Nr, myThid )
./model/src/thermodynamics.F 222 THERMODYNAMICS
       CALL CALC_WSURF_TR(theta,salt,wVel,
./model/src/thermodynamics.F 326 THERMODYNAMICS
     U                                  theta, gtNm,
./model/src/thermodynamics.F 381 THERMODYNAMICS
     I             uVel, vVel, wVel, theta,
./model/src/thermodynamics.F 395 THERMODYNAMICS
     I             uVel, vVel, wVel, theta,
./model/src/thermodynamics.F 642 THERMODYNAMICS
     I         theta,
./model/src/thermodynamics.F 746 THERMODYNAMICS
     I                  theta,
./model/src/thermodynamics.F 754 THERMODYNAMICS
     I                  theta,
./model/src/thermodynamics.F 822 THERMODYNAMICS
     I         kappaRk, wVel, theta,
./model/src/thermodynamics.F 921 THERMODYNAMICS
       CALL DEBUG_STATS_RL(Nr,theta,'Theta (THERMODYNAMICS)',myThid)
./model/src/tracers_correction_step.F 106 TRACERS_CORRECTION_STEP
        CALL ZONAL_FILT_APPLY_TS( theta, salt, myThid )
./model/src/tracers_correction_step.F 65 TRACERS_CORRECTION_STEP
     U           theta,gT,
./model/src/tracers_correction_step.F 99 TRACERS_CORRECTION_STEP
        CALL SHAP_FILT_APPLY_TS( theta,salt, myTime, myIter, 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 119 WRITE_PICKUP
        CALL WRITE_REC_3D_RL( fn, fp, Nr, theta,  -j, myIter, myThid )
./model/src/write_pickup.F 341 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 129 WRITE_STATE
          CALL WRITE_FLD_XYZ_RL( 'T.',suff,theta,myIter,myThid)
./model/src/write_state.F 174 WRITE_STATE
          CALL MNC_CW_RL_W(pf,'state',0,0,'Temp', theta, myThid)
./pkg/aim_v23/aim_dyn2aim.F 169 AIM_DYN2AIM
          TA(I2,Katm)  = theta(I,J,K,bi,bj)*conv_theta2T
./pkg/aim_v23/aim_dyn2aim.F 175 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/autodiff/autodiff_restore.F 146 AUTODIFF_RESTORE
       theta(I,J,K,bi,bj)    = StoreDynVars3D(I,J,K,bi,bj,7)        
./pkg/autodiff/autodiff_store.F 156 AUTODIFF_STORE
       StoreDynVars3D(I,J,K,bi,bj,7) = theta(I,J,K,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 88 BULKF_FORCING
     I            theta(i,j,ks,bi,bj), flwdwn(i,j,bi,bj),
./pkg/cfc/cfc11_surfforcing.F 56 CFC11_SURFFORCING
     &          + sca_11_2 * theta(i,j,k,bi,bj)
./pkg/cfc/cfc11_surfforcing.F 57 CFC11_SURFFORCING
     &          + sca_11_3 * theta(i,j,k,bi,bj)*theta(i,j,k,bi,bj)  
./pkg/cfc/cfc11_surfforcing.F 58 CFC11_SURFFORCING
     &          + sca_11_4 * theta(i,j,k,bi,bj)*theta(i,j,k,bi,bj) 
./pkg/cfc/cfc11_surfforcing.F 59 CFC11_SURFFORCING
     &                *theta(i,j,k,bi,bj)
./pkg/cfc/cfc11_surfforcing.F 61 CFC11_SURFFORCING
              ttemp=( theta(i,j,k,bi,bj) + 273.16 _d 0)* 0.01 _d 0
./pkg/cfc/cfc12_surfforcing.F 56 CFC12_SURFFORCING
     &          + sca_12_2 * theta(i,j,k,bi,bj)
./pkg/cfc/cfc12_surfforcing.F 57 CFC12_SURFFORCING
     &          + sca_12_3 * theta(i,j,k,bi,bj)*theta(i,j,k,bi,bj)  
./pkg/cfc/cfc12_surfforcing.F 58 CFC12_SURFFORCING
     &          + sca_12_4 * theta(i,j,k,bi,bj)*theta(i,j,k,bi,bj) 
./pkg/cfc/cfc12_surfforcing.F 59 CFC12_SURFFORCING
     &                *theta(i,j,k,bi,bj)
./pkg/cfc/cfc12_surfforcing.F 61 CFC12_SURFFORCING
              ttemp=( theta(i,j,k,bi,bj) + 273.16 _d 0)* 0.01 _d 0
./pkg/cheapaml/cheapaml.F 284 CHEAPAML
        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 74 COST_TEST
               objf_test(bi,bj) = theta(i,j,kLocOut,bi,bj)
./pkg/cost/cost_test.F 77 COST_TEST
     &              theta(i,j,kLocOut,bi,bj)
./pkg/cost/cost_test.F 82 COST_TEST
     &                         theta(i,j,kLocOut,bi,bj)**2
./pkg/ctrl/ctrl_map_ini.F 114 CTRL_MAP_INI
                theta(i,j,k,bi,bj) = theta(i,j,k,bi,bj) +
./pkg/ctrl/ctrl_map_ini.F 118 CTRL_MAP_INI
                theta(i,j,k,bi,bj) = theta(i,j,k,bi,bj) +
./pkg/ctrl/ctrl_map_ini.F 122 CTRL_MAP_INI
                if(theta(i,j,k,bi,bj).lt.-2.0)
./pkg/ctrl/ctrl_map_ini.F 123 CTRL_MAP_INI
     &               theta(i,j,k,bi,bj)= -2.0
./pkg/ctrl/ctrl_map_ini.F 207 CTRL_MAP_INI
              theta(i,j,1,bi,bj) = theta(i,j,1,bi,bj)
./pkg/ctrl/ctrl_map_ini.F 565 CTRL_MAP_INI
      _EXCH_XYZ_RL( theta, mythid )
./pkg/ctrl/ctrl_map_ini_ecco.F 147 CTRL_MAP_INI_ECCO
     & xx_theta_dummy, boundsVec, theta, maskc, paramSmooth,
./pkg/ctrl/ctrl_map_ini_ecco.F 149 CTRL_MAP_INI_ECCO
      _EXCH_XYZ_RL( theta, mythid)
./pkg/diagnostics/diagnostics_fill_state.F 195 DIAGNOSTICS_FILL_STATE
     &                  *(theta(i,j,K,bi,bj)+theta(i-1,j,K,bi,bj))
./pkg/diagnostics/diagnostics_fill_state.F 211 DIAGNOSTICS_FILL_STATE
     &                  *(theta(i,j,K,bi,bj)+theta(i,j-1,K,bi,bj))
./pkg/diagnostics/diagnostics_fill_state.F 228 DIAGNOSTICS_FILL_STATE
     &                  *(theta(i,j,K,bi,bj)+theta(i,j,km1,bi,bj))
./pkg/diagnostics/diagnostics_fill_state.F 385 DIAGNOSTICS_FILL_STATE
               tmpMk(i,j,K,bi,bj) = -theta(i,j,K,bi,bj)
./pkg/diagnostics/diagnostics_fill_state.F 419 DIAGNOSTICS_FILL_STATE
        CALL DIAGNOSTICS_FILL(theta,'THETA   ',0,Nr,0,1,1,myThid)
./pkg/diagnostics/diagnostics_fill_state.F 430 DIAGNOSTICS_FILL_STATE
           dummy1 = theta(i,j,k,bi,bj) * ((rc(k)/100.)/1000.)**kappa
./pkg/diagnostics/diagnostics_fill_state.F 444 DIAGNOSTICS_FILL_STATE
        CALL DIAGNOSTICS_SCALE_FILL(theta,tmpFac,2,
./pkg/diagnostics/diagnostics_fill_state.F 537 DIAGNOSTICS_FILL_STATE
     &                  *(theta(i,j,K,bi,bj)+theta(i-1,j,K,bi,bj))
./pkg/diagnostics/diagnostics_fill_state.F 554 DIAGNOSTICS_FILL_STATE
     &                  *(theta(i,j,K,bi,bj)+theta(i,j-1,K,bi,bj))
./pkg/diagnostics/diagnostics_fill_state.F 572 DIAGNOSTICS_FILL_STATE
     &                  *(theta(i,j,K,bi,bj)+theta(i,j,km1,bi,bj))
./pkg/diagnostics/diagnostics_fill_state.F 678 DIAGNOSTICS_FILL_STATE
               tmpMk(i,j,K,bi,bj) = 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 102 DIC_SURFFORCING
     I                       theta,salt,
./pkg/dic/dic_surfforcing.F 136 DIC_SURFFORCING
     I        theta(i,j,kLev,bi,bj),salt(i,j,kLev,bi,bj),
./pkg/dic/dic_surfforcing.F 159 DIC_SURFFORCING
     &          + sca2 * theta(i,j,kLev,bi,bj)
./pkg/dic/dic_surfforcing.F 160 DIC_SURFFORCING
     &          + sca3 * theta(i,j,kLev,bi,bj)*theta(i,j,kLev,bi,bj)
./pkg/dic/dic_surfforcing.F 161 DIC_SURFFORCING
     &          + sca4 * theta(i,j,kLev,bi,bj)*theta(i,j,kLev,bi,bj)
./pkg/dic/dic_surfforcing.F 162 DIC_SURFFORCING
     &                *theta(i,j,kLev,bi,bj)
./pkg/dic/dic_surfforcing_init.F 197 DIC_SURFFORCING_INIT
     I                       theta,salt,
./pkg/dic/dic_surfforcing_init.F 217 DIC_SURFFORCING_INIT
     I        theta(i,j,kLev,bi,bj),salt(i,j,kLev,bi,bj),
./pkg/dic/dic_surfforcing_init.F 235 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 103 EBM_AREA_T
     &           theta(i,j,k,bi,bj)
./pkg/ebm/ebm_area_t.F 85 EBM_AREA_T
     &           theta(i,j,k,bi,bj)
./pkg/ebm/ebm_area_t.F 91 EBM_AREA_T
     &           theta(i,j,k,bi,bj)
./pkg/ebm/ebm_area_t.F 97 EBM_AREA_T
     &           theta(i,j,k,bi,bj)
./pkg/ebm/ebm_forcing_surf.F 76 EBM_FORCING_SURF
     &         (theta(i,j,kSurface,bi,bj)-ZonalMeanSST(j,bj))
./pkg/ebm/ebm_zonalmean.F 65 EBM_ZONALMEAN
     &           theta(i,j,k,bi,bj)
./pkg/ecco/cost_atlantic.F 160 COST_ATLANTIC
                       tv = tv+p5*(theta(i,j,k,bi,bj)
./pkg/ecco/cost_atlantic.F 161 COST_ATLANTIC
     $                     + theta(i,j-1,k,bi,bj))*vVel(i,j,k,bi,bj)
./pkg/ecco/cost_atlantic.F 204 COST_ATLANTIC
     $                     +p5*(theta(i,j,k,bi,bj)
./pkg/ecco/cost_atlantic.F 205 COST_ATLANTIC
     $                     + theta(i,j-1,k,bi,bj))*vVel(i,j,k,bi,bj)
./pkg/ecco/cost_averagesfields.F 176 COST_AVERAGESFIELDS
     &     tbar, theta, xx_tbar_mean_dummy,
./pkg/ecco/cost_averagesfields.F 183 COST_AVERAGESFIELDS
     &     tbar, theta(1-Olx,1-Oly,1,1,1), xx_tbar_mean_dummy,
./pkg/ecco/cost_averagesfields.F 197 COST_AVERAGESFIELDS
     &     sstbar, theta(1-Olx,1-Oly,1,1,1), xx_sstbar_mean_dummy,
./pkg/ecco/cost_averagesfields.F 466 COST_AVERAGESFIELDS
     &     smrsstbar, theta(1-Olx,1-Oly,1,1,1), 
./pkg/ecco/cost_theta_ini_fin.F 109 COST_THETA_INI_FIN
                  tmpx = tmpfld3d(i,j,k,bi,bj) - theta(i,j,k,bi,bj)
./model/inc/DYNVARS.h 37 NO PROCEDURE
     &                   uVel,vVel,wVel,theta,salt,
./model/inc/DYNVARS.h 43 NO PROCEDURE
     &                   uVel,vVel,wVel,theta,salt,
./model/inc/DYNVARS.h 51 NO PROCEDURE
      _RL  theta(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
./pkg/exf/exf_bulkformulae.F 282 EXF_BULKFORMULAE
           Tsf = theta(i,j,ksrf,bi,bj) + cen2kel
./pkg/exf/exf_bulkformulae.F 285 EXF_BULKFORMULAE
     &              *( theta(i,j,ksrf,bi,bj)-theta(i,j,ksrfp1,bi,bj) )
./pkg/exf/exf_getclim.F 154 EXF_GETCLIM
     &          theta(i,j,1,bi,bj) = climsst(i,j,bi,bj)
./pkg/exf/exf_getclim.F 171 EXF_GETCLIM
     &             theta(i,j,1,bi,bj) .EQ. 0. ) then
./pkg/exf/exf_getclim.F 172 EXF_GETCLIM
                print *, 'ph-warn-exf-clim ', i, j, theta(i,j,1,bi,bj)
./pkg/exf/exf_init.F 383 EXF_INIT
     &          theta(i,j,1,bi,bj) = climsst(i,j,bi,bj)
./pkg/exf/exf_init.F 399 EXF_INIT
     &             theta(i,j,1,bi,bj) .EQ. 0. ) then
./pkg/exf/exf_init.F 400 EXF_INIT
                print *, 'ph-warn-exf-init ', i, j, theta(i,j,1,bi,bj)
./pkg/exf/exf_radiation.F 68 EXF_RADIATION
            Tsf = theta(i,j,1,bi,bj) + cen2kel
./pkg/exf/exf_radiation.F 70 EXF_RADIATION
     &             *( theta(i,j,1,bi,bj)-theta(i,j,2,bi,bj) )
./pkg/exf/exf_radiation.F 84 EXF_RADIATION
     &          ((theta(i,j,k,bi,bj)+cen2kel)**4)