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



File Line number Procedure Code
./pkg/ecco/ecco_phys.F 203 ECCO_PHYS
     &                *(theta(i,j,k,bi,bj)+theta(i-1,j,k,bi,bj))*halfRL
./pkg/ecco/ecco_phys.F 206 ECCO_PHYS
     &                *(theta(i,j,k,bi,bj)+theta(i,j-1,k,bi,bj))*halfRL
./pkg/ecco/ecco_phys.F 264 ECCO_PHYS
              tmpfld=theta(i,j,k,bi,bj)
./pkg/ecco/ecco_phys.F 79 ECCO_PHYS
     I                theta(1-OLx,1-OLy,k,bi,bj),
./pkg/exf/exf_bulkformulae.F 290 EXF_BULKFORMULAE
           Tsf = theta(i,j,ksrf,bi,bj) + cen2kel
./pkg/exf/exf_bulkformulae.F 293 EXF_BULKFORMULAE
     &              *( theta(i,j,ksrf,bi,bj)-theta(i,j,ksrfp1,bi,bj) )
./pkg/exf/exf_getclim.F 176 EXF_GETCLIM
     &          theta(i,j,1,bi,bj) = climsst(i,j,bi,bj)
./pkg/exf/exf_getclim.F 193 EXF_GETCLIM
     &             theta(i,j,1,bi,bj) .EQ. 0. ) THEN
./pkg/exf/exf_getclim.F 194 EXF_GETCLIM
                print *, 'ph-warn-exf-clim ', i, j, theta(i,j,1,bi,bj)
./pkg/exf/exf_init_varia.F 462 EXF_INIT_VARIA
     &          theta(i,j,1,bi,bj) = climsst(i,j,bi,bj)
./pkg/exf/exf_init_varia.F 478 EXF_INIT_VARIA
     &             theta(i,j,1,bi,bj) .EQ. 0. ) then
./pkg/exf/exf_init_varia.F 479 EXF_INIT_VARIA
                print *, 'ph-warn-exf-init ', i, j, theta(i,j,1,bi,bj)
./pkg/exf/exf_mapfields.F 190 EXF_MAPFIELDS
     &               *( theta(i,j,ks,bi,bj) - temp_EvPrRn )
./pkg/exf/exf_mapfields.F 200 EXF_MAPFIELDS
     &               *( theta(i,j,ks,bi,bj) - temp_EvPrRn )
./pkg/exf/exf_mapfields.F 215 EXF_MAPFIELDS
     &              *( theta(i,j,ks,bi,bj) - runoftemp(i,j,bi,bj) )
./pkg/exf/exf_radiation.F 100 EXF_RADIATION
     &          ((theta(i,j,ks,bi,bj)+cen2kel)**4)
./pkg/exf/exf_radiation.F 75 EXF_RADIATION
            Tsf = theta(i,j,ks,bi,bj) + cen2kel
./pkg/exf/exf_radiation.F 77 EXF_RADIATION
     &             *( theta(i,j,ks,bi,bj)-theta(i,j,kl,bi,bj) )
./pkg/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 1045 GMREDI_CALC_TENSOR
     &           (theta(i,j-1,km1,bi,bj)-theta(i,j-1,k,bi,bj))
./pkg/gmredi/gmredi_calc_tensor.F 1047 GMREDI_CALC_TENSOR
     &           (theta(i, j ,km1,bi,bj)-theta(i, j ,k,bi,bj))
./pkg/gmredi/gmredi_calc_tensor.F 1051 GMREDI_CALC_TENSOR
     &           (theta(i,j-1,k,bi,bj)-theta(i,j-1,kp1,bi,bj))
./pkg/gmredi/gmredi_calc_tensor.F 1053 GMREDI_CALC_TENSOR
     &           (theta(i, j ,k,bi,bj)-theta(i, j ,kp1,bi,bj))
./pkg/gmredi/gmredi_calc_tensor.F 827 GMREDI_CALC_TENSOR
     &           (theta(i-1,j,km1,bi,bj)-theta(i-1,j,k,bi,bj))
./pkg/gmredi/gmredi_calc_tensor.F 829 GMREDI_CALC_TENSOR
     &           (theta( i ,j,km1,bi,bj)-theta( i ,j,k,bi,bj))
./pkg/gmredi/gmredi_calc_tensor.F 833 GMREDI_CALC_TENSOR
     &           (theta(i-1,j,k,bi,bj)-theta(i-1,j,kp1,bi,bj))
./pkg/gmredi/gmredi_calc_tensor.F 835 GMREDI_CALC_TENSOR
     &           (theta( i ,j,k,bi,bj)-theta( i ,j,kp1,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
./model/src/apply_forcing.F 521 APPLY_FORCING_T
     &          *( temp_addMass - theta(i,j,k,bi,bj) )
./model/src/apply_forcing.F 590 APPLY_FORCING_T
     &                  *( theta(i,j,km,bi,bj)*salt(i,j,km,bi,bj)
./model/src/apply_forcing.F 591 APPLY_FORCING_T
     &                   + theta(i,j,kc,bi,bj)*salt(i,j,kc,bi,bj)
./model/src/apply_forcing.F 603 APPLY_FORCING_T
     &                  *( theta(i,j,kc,bi,bj)*salt(i,j,kc,bi,bj)
./model/src/apply_forcing.F 604 APPLY_FORCING_T
     &                   + theta(i,j,kp,bi,bj)*salt(i,j,kp,bi,bj)
./pkg/kpp/kpp_routines.F 1713 STATEKPP
     I     theta(1-OLx,1-OLy,1,bi,bj), salt(1-OLx,1-OLy,1,bi,bj),
./pkg/kpp/kpp_routines.F 1755 STATEKPP
     I        theta(1-OLx,1-OLy,k,bi,bj), salt(1-OLx,1-OLy,k,bi,bj),
./pkg/kpp/kpp_routines.F 1767 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 1779 STATEKPP
     I        theta(1-OLx,1-OLy,1,bi,bj), salt(1-OLx,1-OLy,1,bi,bj),
./pkg/kpp/kpp_routines.F 1926 KPP_DOUBLEDIFF
         alphaDT(I,J) = ( theta(I,J,Km1,bi,bj)-theta(I,J,K,bi,bj) )
./pkg/layers/layers_calc.F 125 LAYERS_CALC
        CALL LAYERS_FLUXCALC( uVel,vVel,theta,iLa,
./pkg/layers/layers_calc.F 132 LAYERS_CALC
        CALL LAYERS_DIAPYCNAL( theta,iLa,
./pkg/layers/layers_calc.F 182 LAYERS_CALC
     &                       theta(1-OLx,1-OLy,k,bi,bj),
./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),
./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,
./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_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 1054 OBCS_STEVENS_SAVE_TRACERS
            OBNtStevens(i,k,bi,bj) = theta(i,Jobc-1,k,bi,bj)
./pkg/obcs/obcs_calc_stevens.F 1070 OBCS_STEVENS_SAVE_TRACERS
            OBStStevens(i,k,bi,bj) = theta(i,Jobc+1,k,bi,bj)
./pkg/obcs/obcs_calc_stevens.F 1086 OBCS_STEVENS_SAVE_TRACERS
            OBEtStevens(j,k,bi,bj) = theta(Iobc-1,j,k,bi,bj)
./pkg/obcs/obcs_calc_stevens.F 1102 OBCS_STEVENS_SAVE_TRACERS
            OBWtStevens(j,k,bi,bj) = theta(Iobc+1,j,k,bi,bj)
./pkg/obcs/obcs_calc_stevens.F 252 OBCS_CALC_STEVENS
     I        OBEtStevens, theta, gammat, 
./pkg/obcs/obcs_calc_stevens.F 379 OBCS_CALC_STEVENS
     I        OBWtStevens, theta, gammat, 
./pkg/obcs/obcs_calc_stevens.F 490 OBCS_CALC_STEVENS
     I        OBNtStevens, theta, gammat,
./pkg/obcs/obcs_calc_stevens.F 602 OBCS_CALC_STEVENS
     I        OBStStevens, theta, gammat,
./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,
./pkg/obcs/obcs_sponge.F 483 OBCS_SPONGE_T
     &             + float(jsl)*theta(i,j,k,bi,bj) )
./pkg/obcs/obcs_sponge.F 497 OBCS_SPONGE_T
     &            * ( theta(i,j,k,bi,bj) - trelax )
./pkg/obcs/obcs_sponge.F 523 OBCS_SPONGE_T
     &             + float(jsl)*theta(i,j,k,bi,bj) )
./pkg/obcs/obcs_sponge.F 537 OBCS_SPONGE_T
     &            * ( theta(i,j,k,bi,bj) - trelax )
./pkg/obcs/obcs_sponge.F 563 OBCS_SPONGE_T
     &             + float(isl)*theta(i,j,k,bi,bj) )
./pkg/obcs/obcs_sponge.F 577 OBCS_SPONGE_T
     &            * ( theta(i,j,k,bi,bj) - trelax )
./pkg/obcs/obcs_sponge.F 603 OBCS_SPONGE_T
     &             + float(isl)*theta(i,j,k,bi,bj) )
./pkg/obcs/obcs_sponge.F 617 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,
./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/ocn_compon_interf/ocn_store_my_data.F 67 OCN_STORE_MY_DATA
            SSTocn2cpl(i,j,bi,bj) = theta(i,j,1,bi,bj)
./pkg/offline/offline_fields_load.F 375 OFFLINE_FIELDS_LOAD
             theta(i,j,k,bi,bj)= bWght*tave0(i,j,k,bi,bj)
./pkg/openad/externalDummies.F 21 DUMMY_IN_STEPPING
      theta(1,1,1,1,1)=2*theta(1,1,1,1,1)
./model/src/do_atmospheric_phys.F 93 DO_ATMOSPHERIC_PHYS
     &         ( theta(i,j,k,bi,bj)
./pkg/openad/openad_dumpAdjoint.F 106 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/do_fields_blocking_exchanges.F 64 DO_FIELDS_BLOCKING_EXCHANGES
       _EXCH_XYZ_RL( theta, myThid )
./model/src/do_oceanic_phys.F 274 DO_OCEANIC_PHYS
     I                 uVel, vVel, wVel, theta, salt, myThid )
./model/src/do_oceanic_phys.F 770 DO_OCEANIC_PHYS
     I                  theta, salt,
./model/src/do_oceanic_phys.F 783 DO_OCEANIC_PHYS
     I                  theta, salt,
./model/src/do_oceanic_phys.F 794 DO_OCEANIC_PHYS
     I                theta(1-OLx,1-OLy,k,bi,bj),
./model/src/do_oceanic_phys.F 814 DO_OCEANIC_PHYS
     &         ( theta(i,j,k,bi,bj)
./model/src/do_oceanic_phys.F 866 DO_OCEANIC_PHYS
     I                 theta(1-OLx,1-OLy,k-1,bi,bj),
./model/src/do_oceanic_phys.F 952 DO_OCEANIC_PHYS
     I              theta, 0,
./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/do_stagger_fields_exchanges.F 52 DO_STAGGER_FIELDS_EXCHANGES
        CALL EXCH_SM_3D_RL( theta, .FALSE., Nr, myThid )
./pkg/profiles/profiles_interp.F 86 PROFILES_INTERP
               tab_coeffs1(q)=theta(i_cur(q),j_cur(q),k,bi,bj)
./model/src/do_the_model_io.F 117 DO_THE_MODEL_IO
          CALL PLOT_FIELD_XYZRL( theta, 'Current theta ',
./model/src/dynamics.F 483 DYNAMICS
     I        theta, salt,
./model/src/dynamics.F 708 DYNAMICS
       CALL DEBUG_STATS_RL(Nr,theta,'Theta (DYNAMICS)',myThid)
./pkg/rbcs/rbcs_add_tendency.F 98 RBCS_ADD_TENDENCY
     &        *( theta(i,j,k,bi,bj)- RBCtemp(i,j,k,bi,bj) )
./model/src/external_forcing.F 393 EXTERNAL_FORCING_T
     &          *( temp_addMass - theta(i,j,kLev,bi,bj) )
./model/src/external_forcing.F 462 EXTERNAL_FORCING_T
     &                  *( theta(i,j,km,bi,bj)*salt(i,j,km,bi,bj)
./model/src/external_forcing.F 463 EXTERNAL_FORCING_T
     &                   + theta(i,j,kc,bi,bj)*salt(i,j,kc,bi,bj)
./model/src/external_forcing.F 475 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,kp,bi,bj)*salt(i,j,kp,bi,bj)
./model/src/external_forcing_surf.F 262 EXTERNAL_FORCING_SURF
     &          *( temp_EvPrRn - theta(i,j,ks,bi,bj) )
./model/src/external_forcing_surf.F 297 EXTERNAL_FORCING_SURF
     &           *( theta(i,j,ks,bi,bj) - temp_EvPrRn )
./model/src/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/salt_plume/salt_plume_calc_depth.F 116 SALT_PLUME_CALC_DEPTH
     I       theta(1-OLx,1-OLy,K,bi,bj), salt(1-OLx,1-OLy,K,bi,bj),
./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/salt_plume/salt_plume_forcing_surf.F 61 SALT_PLUME_FORCING_SURF
     &   *theta(i,j,1,bi,bj)*rhoConst
./pkg/seaice/cost_ice_test.F 93 COST_ICE_TEST
     &        theta(4,4,kSrf,1,1), area(4,4,1,1), heff(4,4,1,1)
./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/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 1042 SEAICE_GROWTH
          IF ( theta(I,J,kSurface,bi,bj) .GE. tempFrz ) THEN
./pkg/seaice/seaice_growth.F 1060 SEAICE_GROWTH
     &         * (theta(I,J,kSurface,bi,bj)-tempFrz)
./pkg/seaice/seaice_growth.F 1079 SEAICE_GROWTH
           tmpscal1=(theta(I,J,kSurface,bi,bj)-tempFrz)
./pkg/seaice/seaice_growth.F 2277 SEAICE_GROWTH
     &         HeatCapacity_Cp * theta(I,J,kSurface,bi,bj)
./pkg/seaice/seaice_growth.F 2280 SEAICE_GROWTH
     &       ( temp_EvPrRn - theta(I,J,kSurface,bi,bj) )
./pkg/seaice/seaice_growth.F 2345 SEAICE_GROWTH
     &        ( ZERO + HeatCapacity_Cp * theta(I,J,kSurface,bi,bj) )
./pkg/seaice/seaice_growth.F 2348 SEAICE_GROWTH
     &        ( temp_EvPrRn - theta(I,J,kSurface,bi,bj) )
./pkg/seaice/seaice_growth.F 2489 SEAICE_GROWTH
     &       HeatCapacity_Cp * theta(I,J,kSurface,bi,bj)
./pkg/seaice/seaice_growth.F 2649 SEAICE_GROWTH
     &       ( ZERO + HeatCapacity_Cp * theta(I,J,kSurface,bi,bj) )
./pkg/seaice/seaice_growth.F 2652 SEAICE_GROWTH
     &       HeatCapacity_Cp*(temp_EvPrRn - theta(I,J,kSurface,bi,bj))
./pkg/seaice/seaice_growth.F 734 SEAICE_GROWTH
          TmixLoc(i,j) = theta(i,j,kSurface,bi,bj)+celsius2K
./model/src/ini_dynvars.F 57 INI_DYNVARS
           theta(i,j,k,bi,bj) = 0. _d 0
./pkg/shelfice/shelfice_thermodynamics.F 263 SHELFICE_THERMODYNAMICS
          tLoc(I,J) = theta(I,J,K,bi,bj)
./pkg/shelfice/shelfice_thermodynamics.F 319 SHELFICE_THERMODYNAMICS
     &           + theta(I,J,Kp1,bi,bj) *drKp1 )
./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,
./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/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.
./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 157 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 175 CTRL_MAP_INI_GENARR
       _EXCH_XYZ_RL( theta, myThid )
./verification/halfpipe_streamice/code_oad/ctrl_map_ini_genarr.F 157 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 175 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 163 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 110 CTRL_MAP_INI_GENARR
       _EXCH_XYZ_RL( theta, myThid )
./verification/hs94.1x64x5/code_ad/ctrl_map_ini_genarr.F 97 CTRL_MAP_INI_GENARR
                theta(i,j,1,bi,bj) = theta(i,j,1,bi,bj)
./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 163 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 110 CTRL_MAP_INI_GENARR
       _EXCH_XYZ_RL( theta, myThid )
./verification/hs94.1x64x5/code_oad/ctrl_map_ini_genarr.F 97 CTRL_MAP_INI_GENARR
                theta(i,j,1,bi,bj) = theta(i,j,1,bi,bj)
./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
./model/src/read_pickup.F 174 READ_PICKUP
          CALL READ_REC_3D_RL( fn, fp, Nr, theta, 7, myIter,myThid )
./model/src/read_pickup.F 206 READ_PICKUP
          CALL READ_REC_3D_RL( fn, fp, Nr, theta,nj+1, myIter,myThid )
./model/src/read_pickup.F 223 READ_PICKUP
          CALL READ_REC_3D_RL( fn, fp, Nr, theta, 5, myIter,myThid )
./model/src/read_pickup.F 274 READ_PICKUP
          CALL READ_MFLDS_3D_RL( 'Theta   ', theta,
./model/src/read_pickup.F 471 READ_PICKUP
        CALL MNC_CW_RL_R('D',fn,0,0,'Temp',theta, myThid)
./model/src/read_pickup.F 516 READ_PICKUP
      CALL EXCH_3D_RL( theta, Nr, myThid )
./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)
./verification/tutorial_global_oce_optim/code_ad/external_forcing_surf.F 265 EXTERNAL_FORCING_SURF
     &          *( temp_EvPrRn - theta(i,j,ks,bi,bj) )
./verification/tutorial_global_oce_optim/code_ad/external_forcing_surf.F 300 EXTERNAL_FORCING_SURF
     &           *( theta(i,j,ks,bi,bj) - temp_EvPrRn )
./verification/tutorial_global_oce_optim/code_oad/external_forcing_surf.F 265 EXTERNAL_FORCING_SURF
     &          *( temp_EvPrRn - theta(i,j,ks,bi,bj) )
./verification/tutorial_global_oce_optim/code_oad/external_forcing_surf.F 300 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 200 TEMP_INTEGRATE
     I                           theta(1-OLx,1-OLy,1,bi,bj),
./model/src/temp_integrate.F 206 TEMP_INTEGRATE
     I                           theta(1-OLx,1-OLy,1,bi,bj),
./model/src/temp_integrate.F 274 TEMP_INTEGRATE
     I             dTtracerLev, uFld, vFld, wFld, theta,
./model/src/temp_integrate.F 288 TEMP_INTEGRATE
     I             dTtracerLev, uFld, vFld, wFld, theta,
./model/src/temp_integrate.F 355 TEMP_INTEGRATE
     I           theta(1-OLx,1-OLy,1,bi,bj),
./model/src/temp_integrate.F 370 TEMP_INTEGRATE
     I           theta(1-OLx,1-OLy,1,bi,bj),
./model/src/temp_integrate.F 459 TEMP_INTEGRATE
     I                  theta(1-OLx,1-OLy,1,bi,bj),
./model/src/temp_integrate.F 466 TEMP_INTEGRATE
     I                  theta(1-OLx,1-OLy,1,bi,bj),
./model/src/temp_integrate.F 477 TEMP_INTEGRATE
     I                  theta(1-OLx,1-OLy,1,bi,bj),
./model/src/temp_integrate.F 497 TEMP_INTEGRATE
     I         theta(1-OLx,1-OLy,1,bi,bj),
./model/src/temp_integrate.F 531 TEMP_INTEGRATE
     U             theta(1-OLx,1-OLy,1,bi,bj),
./model/src/temp_integrate.F 537 TEMP_INTEGRATE
     U             theta(1-OLx,1-OLy,1,bi,bj),
./model/src/temp_integrate.F 545 TEMP_INTEGRATE
     O             theta(1-OLx,1-OLy,1,bi,bj),
./model/src/thermodynamics.F 143 THERMODYNAMICS
       CALL CALC_WSURF_TR( theta, salt, wVel,
./model/src/thermodynamics.F 148 THERMODYNAMICS
        CALL LAYERS_WSURF_TR( theta, salt, wVel,
./model/src/thermodynamics.F 357 THERMODYNAMICS
          CALL OBCS_APPLY_TS( bi, bj, 0, theta, salt, myThid )
./model/src/thermodynamics.F 396 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 132 WRITE_PICKUP
        CALL WRITE_REC_3D_RL( fn, fp, Nr, theta,  -j, myIter, myThid )
./model/src/write_pickup.F 376 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 136 WRITE_STATE
          CALL WRITE_FLD_XYZ_RL( 'T.',suff,thetav,myIter,myThid)
./model/src/write_state.F 147 WRITE_STATE
          CALL WRITE_FLD_XYZ_RL( 'T.',suff,theta,myIter,myThid)
./model/src/write_state.F 194 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_dyn2phys.F 120 ATM_PHYS_DYN2PHYS
        t3d(:,:,k) = theta(1:sNx,1:sNy,kc,bi,bj)*conv_theta2T
./pkg/autodiff/addummy_in_stepping.F 455 ADDUMMY_IN_STEPPING
          foo3D = thetad
./pkg/autodiff/addummy_in_stepping.F 535 ADDUMMY_IN_STEPPING
          foo3D = thetad
./pkg/autodiff/addummy_in_stepping.F 623 ADDUMMY_IN_STEPPING
          CALL COPY_ADVAR_OUTP( dumRS, thetad,var3Du, Nr, 12, myThid )
./pkg/autodiff/autodiff_restore.F 147 AUTODIFF_RESTORE
           theta(I,J,K,bi,bj)    = StoreDynVars3D(I,J,K,bi,bj,5)
./pkg/autodiff/autodiff_store.F 157 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/bling/bling_airseaflux.F 124 BLING_AIRSEAFLUX
             surftemp(i,j) = theta(i,j,1,bi,bj)
./pkg/bling/bling_airseaflux.F 190 BLING_AIRSEAFLUX
     &          + sca2 * theta(i,j,1,bi,bj)
./pkg/bling/bling_airseaflux.F 191 BLING_AIRSEAFLUX
     &          + sca3 * theta(i,j,1,bi,bj)*theta(i,j,1,bi,bj)
./pkg/bling/bling_airseaflux.F 192 BLING_AIRSEAFLUX
     &          + sca4 * theta(i,j,1,bi,bj)*theta(i,j,1,bi,bj)
./pkg/bling/bling_airseaflux.F 193 BLING_AIRSEAFLUX
     &                *theta(i,j,1,bi,bj)
./pkg/bling/bling_airseaflux.F 250 BLING_AIRSEAFLUX
              ttemp = theta(i,j,1,bi,bj)
./pkg/bling/bling_carbonate_init.F 181 BLING_CARBONATE_INIT
             thetatmp(i,j) = theta(i,j,k,bi,bj)
./pkg/bling/bling_carbonate_sys.F 92 BLING_CARBONATE_SYS
             ttmp(i,j) = theta(i,j,k,bi,bj)
./pkg/bling/bling_dvm.F 150 BLING_DVM
            temp_upper = temp_upper + theta(i,j,k,bi,bj)*drf(k)
./pkg/bling/bling_dvm.F 156 BLING_DVM
            temp_lower = temp_lower + theta(i,j,k,bi,bj)*drf(k)
./pkg/bling/bling_mixedlayer.F 68 BLING_MIXEDLAYER
     I     theta(1-OLx,1-OLy,1,bi,bj), salt(1-OLx,1-OLy,1,bi,bj),
./pkg/bling/bling_mixedlayer.F 89 BLING_MIXEDLAYER
     I        theta(1-OLx,1-OLy,k,bi,bj), salt(1-OLx,1-OLy,k,bi,bj),
./pkg/bling/bling_production.F 421 BLING_PROD
          expkT(i,j,k) = exp(kappa_eppley * theta(i,j,k,bi,bj))   
./pkg/bling/bling_production.F 431 BLING_PROD
     &           * exp(kappa_eppley_diaz * theta(i,j,k,bi,bj))
./pkg/bling/bling_production.F 487 BLING_PROD
          IF (theta(i,j,k,bi,bj) .lt. 14) THEN
./pkg/bling/bling_production.F 621 BLING_PROD
     &           exp(kappa_remin * 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 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 361 CHEAPAML
     I                      theta(1-OLx,1-OLy,1,bi,bj), windSq,
./pkg/cheapaml/cheapaml.F 399 CHEAPAML
     I                      theta(1-OLx,1-OLy,1,bi,bj), windSq,
./pkg/cheapaml/cheapaml.F 437 CHEAPAML
             xlwnet = 0.98 _d 0*stefan*(theta(i,j,1,bi,bj)+celsius2K)**4
./pkg/cheapaml/cheapaml.F 442 CHEAPAML
     &        + 4.0*0.98 _d 0*stefan*(theta(i,j,1,bi,bj)+celsius2K)**3
./pkg/cheapaml/cheapaml.F 443 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 127 CTRL_MAP_INI
                theta(i,j,k,bi,bj) = theta(i,j,k,bi,bj) +
./pkg/ctrl/ctrl_map_ini.F 130 CTRL_MAP_INI
                if(theta(i,j,k,bi,bj).lt.-2.0)
./pkg/ctrl/ctrl_map_ini.F 131 CTRL_MAP_INI
     &               theta(i,j,k,bi,bj)= -2.0
./pkg/ctrl/ctrl_map_ini.F 215 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 142 CTRL_MAP_INI_ECCO
     & xx_theta_dummy, boundsVec, theta, maskc, paramSmooth,
./pkg/ctrl/ctrl_map_ini_ecco.F 144 CTRL_MAP_INI_ECCO
      _EXCH_XYZ_RL( theta, mythid)
./pkg/ctrl/ctrl_map_ini_genarr.F 128 CTRL_MAP_INI_GENARR
     &   call CTRL_MAP_GENARR3D(theta,igen_theta0,myThid)
./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 409 DIAGNOSTICS_FILL_STATE
        CALL DIAGNOSTICS_FILL(theta,'THETA   ',0,Nr,0,1,1,myThid)
./pkg/diagnostics/diagnostics_fill_state.F 420 DIAGNOSTICS_FILL_STATE
            dummy1 = theta(i,j,k,bi,bj) * ((rC(k)/100.)/1000.)**kappa
./pkg/diagnostics/diagnostics_fill_state.F 434 DIAGNOSTICS_FILL_STATE
        CALL DIAGNOSTICS_SCALE_FILL( theta, oneRL, 2,
./pkg/diagnostics/diagnostics_fill_state.F 550 DIAGNOSTICS_FILL_STATE
     &                  *(theta(i,j,k,bi,bj)+theta(i-1,j,k,bi,bj))
./pkg/diagnostics/diagnostics_fill_state.F 567 DIAGNOSTICS_FILL_STATE
     &                  *(theta(i,j,k,bi,bj)+theta(i,j-1,k,bi,bj))
./pkg/diagnostics/diagnostics_fill_state.F 585 DIAGNOSTICS_FILL_STATE
     &                  *(theta(i,j,k,bi,bj)+theta(i,j,km1,bi,bj))
./pkg/diagnostics/diagnostics_fill_state.F 670 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 112 DIC_SURFFORCING
     &                          min(50. _d 0, theta(i,j,kLev,bi,bj)))
./pkg/dic/dic_surfforcing.F 118 DIC_SURFFORCING
            surftemp(i,j) = theta(i,j,kLev,bi,bj)
./pkg/dic/dic_surfforcing.F 185 DIC_SURFFORCING
     &          + sca2 * theta(i,j,kLev,bi,bj)
./pkg/dic/dic_surfforcing.F 186 DIC_SURFFORCING
     &          + sca3 * theta(i,j,kLev,bi,bj)*theta(i,j,kLev,bi,bj)
./pkg/dic/dic_surfforcing.F 187 DIC_SURFFORCING
     &          + sca4 * theta(i,j,kLev,bi,bj)*theta(i,j,kLev,bi,bj)
./pkg/dic/dic_surfforcing.F 188 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_averagesfields.F 274 COST_AVERAGESFIELDS
     &     tbar, theta, xx_tbar_mean_dummy,
./pkg/ecco/cost_averagesfields.F 281 COST_AVERAGESFIELDS
     &     tbar, theta(1-Olx,1-Oly,1,1,1), xx_tbar_mean_dummy,
./pkg/ecco/cost_averagesfields.F 296 COST_AVERAGESFIELDS
     &     sstbar, theta(1-Olx,1-Oly,1,1,1), xx_sstbar_mean_dummy,
./pkg/ecco/cost_averagesfields.F 568 COST_AVERAGESFIELDS
     &     smrsstbar, theta(1-Olx,1-Oly,1,1,1),
./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_gencost_customize.F 203 COST_GENCOST_CUSTOMIZE
     &      theta(i,j,1,bi,bj)*maskC(i,j,1,bi,bj)
./pkg/ecco/cost_gencost_customize.F 217 COST_GENCOST_CUSTOMIZE
     &       theta(i,j,k2,bi,bj)*maskC(i,j,k2,bi,bj)