Symbol Description Uses
theta sea surface temperature (K) 350



File Line number Procedure Code
./pkg/ebm/ebm_forcing_surf.F 78 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/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/exf/exf_radiation.F 68 EXF_RADIATION
     &         ((theta(i,j,k,bi,bj)+cen2kel)**4)
./pkg/exf/exf_radiation.F 73 EXF_RADIATION
     &         ((theta(i,j,k,bi,bj)+cen2kel)**4)
./pkg/exf/exf_bulkformulae.F 242 EXF_BULKFORMULAE
             ssttmp = theta(i,j,k,bi,bj)
./pkg/shelfice/shelfice_thermodynamics.F 136 SHELFICE_THERMODYNAMICS
            tLoc = theta(I,J,K,bi,bj)
./pkg/shelfice/shelfice_thermodynamics.F 86 SHELFICE_THERMODYNAMICS
            tLoc = theta(I,J,K,bi,bj)
./pkg/pp81/pp81_ri_number.F 102 PP81_RI_NUMBER
     I     theta, salt,
./pkg/pp81/pp81_ri_number.F 97 PP81_RI_NUMBER
     I     theta, salt,
./pkg/gmredi/gmredi_calc_tensor.F 446 GMREDI_CALC_TENSOR
     &           (theta(i-1,j,km1,bi,bj)-theta(i-1,j,k,bi,bj))
./pkg/gmredi/gmredi_calc_tensor.F 448 GMREDI_CALC_TENSOR
     &           (theta( i ,j,km1,bi,bj)-theta( i ,j,k,bi,bj))
./pkg/gmredi/gmredi_calc_tensor.F 452 GMREDI_CALC_TENSOR
     &           (theta(i-1,j,k,bi,bj)-theta(i-1,j,kp1,bi,bj))
./pkg/gmredi/gmredi_calc_tensor.F 454 GMREDI_CALC_TENSOR
     &           (theta( i ,j,k,bi,bj)-theta( i ,j,kp1,bi,bj))
./pkg/gmredi/gmredi_calc_tensor.F 565 GMREDI_CALC_TENSOR
     &           (theta(i,j-1,km1,bi,bj)-theta(i,j-1,k,bi,bj))
./pkg/gmredi/gmredi_calc_tensor.F 567 GMREDI_CALC_TENSOR
     &           (theta(i, j ,km1,bi,bj)-theta(i, j ,k,bi,bj))
./pkg/gmredi/gmredi_calc_tensor.F 571 GMREDI_CALC_TENSOR
     &           (theta(i,j-1,k,bi,bj)-theta(i,j-1,kp1,bi,bj))
./pkg/gmredi/gmredi_calc_tensor.F 573 GMREDI_CALC_TENSOR
     &           (theta(i, j ,k,bi,bj)-theta(i, j ,kp1,bi,bj))
./model/inc/DYNVARS.h 36 NO PROCEDURE
     &                   uVel,vVel,wVel,theta,salt,
./model/inc/DYNVARS.h 42 NO PROCEDURE
     &                   uVel,vVel,wVel,theta,salt,
./model/inc/DYNVARS.h 50 NO PROCEDURE
      _RL  theta(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
./pkg/diagnostics/diagnostics_fill_state.F 183 DIAGNOSTICS_FILL_STATE
     &                  *(theta(i,j,K,bi,bj)+theta(i-1,j,K,bi,bj))
./pkg/diagnostics/diagnostics_fill_state.F 199 DIAGNOSTICS_FILL_STATE
     &                  *(theta(i,j,K,bi,bj)+theta(i,j-1,K,bi,bj))
./pkg/diagnostics/diagnostics_fill_state.F 216 DIAGNOSTICS_FILL_STATE
     &                  *(theta(i,j,K,bi,bj)+theta(i,j,km1,bi,bj))
./pkg/diagnostics/diagnostics_fill_state.F 371 DIAGNOSTICS_FILL_STATE
               tmpMk(i,j,K,bi,bj) = -theta(i,j,K,bi,bj)
./pkg/diagnostics/diagnostics_fill_state.F 405 DIAGNOSTICS_FILL_STATE
        CALL DIAGNOSTICS_FILL(theta,'THETA   ',0,Nr,0,1,1,myThid)
./pkg/diagnostics/diagnostics_fill_state.F 416 DIAGNOSTICS_FILL_STATE
           dummy1 = theta(i,j,k,bi,bj) * ((rc(k)/100.)/1000.)**kappa
./pkg/diagnostics/diagnostics_fill_state.F 430 DIAGNOSTICS_FILL_STATE
        CALL DIAGNOSTICS_SCALE_FILL(theta,tmpFac,2,
./pkg/diagnostics/diagnostics_fill_state.F 523 DIAGNOSTICS_FILL_STATE
     &                  *(theta(i,j,K,bi,bj)+theta(i-1,j,K,bi,bj))
./pkg/diagnostics/diagnostics_fill_state.F 540 DIAGNOSTICS_FILL_STATE
     &                  *(theta(i,j,K,bi,bj)+theta(i,j-1,K,bi,bj))
./pkg/diagnostics/diagnostics_fill_state.F 558 DIAGNOSTICS_FILL_STATE
     &                  *(theta(i,j,K,bi,bj)+theta(i,j,km1,bi,bj))
./pkg/diagnostics/diagnostics_fill_state.F 664 DIAGNOSTICS_FILL_STATE
               tmpMk(i,j,K,bi,bj) = theta(i,j,K,bi,bj)
./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/flt/flt_traj.F 140 FLT_TRAJ
               call FLT_BILINEAR  (xx,yy,tt,kp,theta, 1,bi,bj)
./pkg/flt/flt_up.F 177 FLT_UP
               call FLT_BILINEAR  (xx,yy,tt,k,theta, 1,bi,bj)
./model/src/find_alpha.F 117 FIND_ALPHA
     I        theta, salt,
./model/src/find_alpha.F 123 FIND_ALPHA
     I        locPres, theta, salt,
./model/src/find_alpha.F 131 FIND_ALPHA
               t1 = theta(i,j,k,bi,bj) 
./model/src/find_alpha.F 223 FIND_ALPHA
     &      locPres, theta, salt, rhoLoc, myThid )
./model/src/find_alpha.F 226 FIND_ALPHA
     &        locPres, theta, salt, rhoDen, myThid )
./model/src/find_alpha.F 230 FIND_ALPHA
               t1  = theta(i,j,k,bi,bj)
./model/src/find_alpha.F 339 FIND_BETA
               tP=theta(i,j,k,bi,bj)-refTemp
./model/src/find_alpha.F 370 FIND_BETA
     I        theta, salt,
./model/src/find_alpha.F 376 FIND_BETA
     I        locPres, theta, salt,
./model/src/find_alpha.F 384 FIND_BETA
               t1 = theta(i,j,k,bi,bj) 
./model/src/find_alpha.F 462 FIND_BETA
     &      locPres, theta, salt, rhoLoc, myThid )
./model/src/find_alpha.F 465 FIND_BETA
     &        locPres, theta, salt, rhoDen, myThid )
./model/src/find_alpha.F 469 FIND_BETA
               t1  = theta(i,j,k,bi,bj)
./model/src/find_alpha.F 84 FIND_ALPHA
               tP=theta(i,j,k,bi,bj)-refTemp
./model/src/do_the_model_io.F 62 DO_THE_MODEL_IO
          CALL PLOT_FIELD_XYZRL( theta, 'Current theta ' , 
./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),
./model/src/do_fields_blocking_exchanges.F 59 DO_FIELDS_BLOCKING_EXCHANGES
       _EXCH_XYZ_R8( theta, myThid )
./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/dic/o2_surfforcing.F 81 O2_SURFFORCING
     &          + sox2 * theta(i,j,k,bi,bj)
./pkg/dic/o2_surfforcing.F 82 O2_SURFFORCING
     &          + sox3 * theta(i,j,k,bi,bj)*theta(i,j,k,bi,bj)  
./pkg/dic/o2_surfforcing.F 83 O2_SURFFORCING
     &          + sox4 * theta(i,j,k,bi,bj)*theta(i,j,k,bi,bj) 
./pkg/dic/o2_surfforcing.F 84 O2_SURFFORCING
     &                *theta(i,j,k,bi,bj)
./pkg/dic/o2_surfforcing.F 94 O2_SURFFORCING
              ttemp = theta(i,j,k,bi,bj)
./model/src/tracers_correction_step.F 104 TRACERS_CORRECTION_STEP
        CALL SHAP_FILT_APPLY_TS( theta,salt, myTime, myIter, myThid )
./model/src/tracers_correction_step.F 111 TRACERS_CORRECTION_STEP
        CALL ZONAL_FILT_APPLY_TS( theta, salt, myThid )
./model/src/tracers_correction_step.F 69 TRACERS_CORRECTION_STEP
     U           theta,gT,
./pkg/dic/dic_surfforcing.F 110 DIC_SURFFORCING
     I        theta(i,j,kLev,bi,bj),salt(i,j,kLev,bi,bj),
./pkg/dic/dic_surfforcing.F 132 DIC_SURFFORCING
     &          + sca2 * theta(i,j,kLev,bi,bj)
./pkg/dic/dic_surfforcing.F 133 DIC_SURFFORCING
     &          + sca3 * theta(i,j,kLev,bi,bj)*theta(i,j,kLev,bi,bj)  
./pkg/dic/dic_surfforcing.F 134 DIC_SURFFORCING
     &          + sca4 * theta(i,j,kLev,bi,bj)*theta(i,j,kLev,bi,bj) 
./pkg/dic/dic_surfforcing.F 135 DIC_SURFFORCING
     &                *theta(i,j,kLev,bi,bj)
./pkg/dic/dic_surfforcing.F 98 DIC_SURFFORCING
     I                       theta,salt,
./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 75 CALCITE_SATURATION
     I                       theta,salt,
./pkg/dic/dic_surfforcing_init.F 150 DIC_SURFFORCING_INIT
     I                       theta,salt,
./pkg/dic/dic_surfforcing_init.F 174 DIC_SURFFORCING_INIT
     I        theta(i,j,kLev,bi,bj),salt(i,j,kLev,bi,bj),
./pkg/dic/dic_surfforcing_init.F 190 DIC_SURFFORCING_INIT
        print*,'QQ first guess pH', pH(20,20,1,1), theta(20,20,1,1,1),
./model/src/ini_dynvars.F 54 INI_DYNVARS
           theta(I,J,K,bi,bj)=0. _d 0
./model/src/tracers_iigw_correction.F 82 TRACERS_IIGW_CORRECTION
          theta(i,j,k,bi,bj) = theta(i,j,k,bi,bj)
./pkg/offline/offline_fields_load.F 397 OFFLINE_FIELDS_LOAD
          theta(i,j,k,bi,bj)    = bWght*tave0(i,j,k,bi,bj)
./pkg/fizhi/fizhi_init_vars.F 162 FIZHI_INIT_VARS
        call DYN2PHYS(theta,pedyn,im1,im2,jm1,jm2,Nr,Nsx,Nsy,
./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/fizhi/fizhi_wrapper.F 126 FIZHI_WRAPPER
         thdynbef(i,j,L,bi,bj) = theta(i,j,L,bi,bj)
./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/advect_xy/code/ini_theta.F 121 INI_THETA
            IF (theta(I,J,k,bi,bj) .LT. Tfreezing) THEN
./verification/advect_xy/code/ini_theta.F 122 INI_THETA
               theta(I,J,K,bi,bj) = Tfreezing
./verification/advect_xy/code/ini_theta.F 131 INI_THETA
      _EXCH_XYZ_R8(theta,myThid)
./verification/advect_xy/code/ini_theta.F 133 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 68 INI_THETA
            theta(i,j,k,bi,bj)=exp(-0.5*( rD/20.E3 )**2 )
./verification/advect_xy/code/ini_theta.F 77 INI_THETA
       CALL READ_FLD_XYZ_RL( hydrogThetaFile, ' ', theta, 0, myThid )
./verification/advect_xy/code/ini_theta.F 79 INI_THETA
       _EXCH_XYZ_R8(theta,myThid)
./verification/advect_xy/code/ini_theta.F 89 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 96 INI_THETA
     &      .AND. theta(I,J,K,bi,bj).EQ.0. ) THEN
./verification/hs94.128x64x5/code/external_forcing.F 206 EXTERNAL_FORCING_T
     &        - kT*( theta(i,j,kLev,bi,bj)-thetaEq )
./verification/rotating_tank/code/external_forcing.F 387 EXTERNAL_FORCING_T
     &      -faceArea*kDiffCyl*(theta(i,j,kLev,bi,bj) - tCyl)
./verification/rotating_tank/code/external_forcing.F 394 EXTERNAL_FORCING_T
     &      -faceArea*kDiffCyl*(tCyl - theta(i,j,kLev,bi,bj))
./verification/rotating_tank/code/external_forcing.F 401 EXTERNAL_FORCING_T
     &      -faceArea*kDiffCyl*(tCyl-theta(i,j,kLev,bi,bj))
./verification/rotating_tank/code/external_forcing.F 408 EXTERNAL_FORCING_T
     &      -faceArea*kDiffCyl*(theta(i,j,kLev,bi,bj) - tCyl)
./verification/rotating_tank/code/checkpoint.F 119 READ_CHECKPOINT
        CALL MDSREADFIELD(fn,prec,'RL',Nr,theta,  8,myThid)
./verification/rotating_tank/code/checkpoint.F 158 READ_CHECKPOINT
      _EXCH_XYZ_R8(theta , myThid )
./verification/rotating_tank/code/checkpoint.F 285 WRITE_CHECKPOINT
        CALL WRITE_REC_XYZ_RL( 'theta',  theta,  1,myIter, myThid)
./verification/rotating_tank/code/checkpoint.F 319 WRITE_CHECKPOINT
        CALL MDSWRITEFIELD(fn,prec,lgf,'RL',Nr,theta, 8,myIter,myThid)
./verification/rotating_tank/code/checkpoint.F 90 READ_CHECKPOINT
        CALL READ_REC_XYZ_RL( 'theta',  theta, 1,myIter, myThid)
./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/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 104 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 111 INI_THETA
     &      .AND. theta(I,J,K,bi,bj).EQ.0. ) THEN
./verification/hs94.1x64x5/code/ini_theta.F 128 INI_THETA
      _EXCH_XYZ_R8(theta , myThid )
./verification/hs94.1x64x5/code/ini_theta.F 130 INI_THETA
      CALL PLOT_FIELD_XYZRL( theta, 'Initial Temperature' ,
./verification/hs94.1x64x5/code/ini_theta.F 75 INI_THETA
            theta(I,J,K,bi,bj) = MAX( thetaLim, thetaEq )
./verification/hs94.1x64x5/code/ini_theta.F 86 INI_THETA
     U      theta, hFacC,
./verification/hs94.1x64x5/code/ini_theta.F 94 INI_THETA
       CALL READ_FLD_XYZ_RL( hydrogThetaFile, ' ', theta, 0, myThid )
./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/cost_test.F 74 COST_TEST
     &              theta(i,j,kLocOut,bi,bj)
./verification/hs94.1x64x5/code_ad/ini_theta.F 104 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 111 INI_THETA
     &      .AND. theta(I,J,K,bi,bj).EQ.0. ) THEN
./verification/hs94.1x64x5/code_ad/ini_theta.F 128 INI_THETA
      _EXCH_XYZ_R8(theta , myThid )
./verification/hs94.1x64x5/code_ad/ini_theta.F 130 INI_THETA
      CALL PLOT_FIELD_XYZRL( theta, 'Initial Temperature' ,
./verification/hs94.1x64x5/code_ad/ini_theta.F 75 INI_THETA
            theta(I,J,K,bi,bj) = MAX( thetaLim, thetaEq )
./verification/hs94.1x64x5/code_ad/ini_theta.F 86 INI_THETA
     U      theta, hFacC,
./verification/hs94.1x64x5/code_ad/ini_theta.F 94 INI_THETA
       CALL READ_FLD_XYZ_RL( hydrogThetaFile, ' ', theta, 0, myThid )
./verification/hs94.cs-32x32x5/code/external_forcing.F 206 EXTERNAL_FORCING_T
     &        - kT*( theta(i,j,kLev,bi,bj)-thetaEq )
./verification/hs94.cs-32x32x5/code/ini_theta.F 104 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 111 INI_THETA
     &      .AND. theta(I,J,K,bi,bj).EQ.0. ) THEN
./verification/hs94.cs-32x32x5/code/ini_theta.F 128 INI_THETA
      _EXCH_XYZ_R8(theta , myThid )
./verification/hs94.cs-32x32x5/code/ini_theta.F 130 INI_THETA
      CALL PLOT_FIELD_XYZRL( theta, 'Initial Temperature' ,
./verification/hs94.cs-32x32x5/code/ini_theta.F 75 INI_THETA
            theta(I,J,K,bi,bj) = MAX( thetaLim, thetaEq )
./verification/hs94.cs-32x32x5/code/ini_theta.F 86 INI_THETA
     U      theta, hFacC,
./verification/hs94.cs-32x32x5/code/ini_theta.F 94 INI_THETA
       CALL READ_FLD_XYZ_RL( hydrogThetaFile, ' ', theta, 0, myThid )
./model/src/dynamics.F 467 DYNAMICS
     I        theta, salt,
./model/src/dynamics.F 663 DYNAMICS
       CALL DEBUG_STATS_RL(Nr,theta,'Theta (DYNAMICS)',myThid)
./verification/internal_wave/code/obcs_calc.F 152 OBCS_CALC
     &          uVel, vVel, wVel, theta, salt, 
./verification/internal_wave/code/obcs_calc.F 172 OBCS_CALC
     &          uVel, vVel, wVel, theta, salt, 
./verification/internal_wave/code/obcs_calc.F 32 OBCS_CALC
      _RL theta(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
./verification/internal_wave/code/obcs_calc.F 79 OBCS_CALC
     &          uVel, vVel, wVel, theta, salt, 
./verification/internal_wave/code/obcs_calc.F 7 OBCS_CALC
     &                      uVel, vVel, wVel, theta, salt, 
./verification/internal_wave/code/obcs_calc.F 99 OBCS_CALC
     &          uVel, vVel, wVel, theta, salt, 
./verification/dome/code/obcs_calc.F 123 OBCS_CALC
     &          uVel, vVel, wVel, theta, salt, 
./verification/dome/code/obcs_calc.F 142 OBCS_CALC
     &          uVel, vVel, wVel, theta, salt, 
./verification/dome/code/obcs_calc.F 204 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 91 OBCS_CALC
     &          uVel, vVel, wVel, theta, salt, 
./verification/tutorial_held_suarez_cs/code/external_forcing.F 206 EXTERNAL_FORCING_T
     &        - kT*( theta(i,j,kLev,bi,bj)-thetaEq )
./model/src/checkpoint.F 121 READ_CHECKPOINT
        CALL READ_REC_XYZ_RL( 'theta',  theta, 1,myIter,myThid)
./model/src/checkpoint.F 155 READ_CHECKPOINT
          CALL MDSREADFIELD(fn,prec,'RL',Nr,theta,  7,myThid)
./model/src/checkpoint.F 185 READ_CHECKPOINT
          CALL MDSREADFIELD(fn,prec,'RL',Nr,theta,nj+1,myThid)
./model/src/checkpoint.F 202 READ_CHECKPOINT
          CALL MDSREADFIELD(fn,prec,'RL',Nr,theta,  5,myThid)
./model/src/checkpoint.F 246 READ_CHECKPOINT
        CALL MNC_CW_RL_R('D',fn,0,0,'Temp',theta, myThid)
./model/src/checkpoint.F 285 READ_CHECKPOINT
      _EXCH_XYZ_R8(theta , myThid )
./model/src/checkpoint.F 492 WRITE_CHECKPOINT_NOW
        CALL WRITE_REC_XYZ_RL( 'theta',  theta, 1,myIter,myThid)
./model/src/checkpoint.F 527 WRITE_CHECKPOINT_NOW
        CALL MDSWRITEFIELD(fn,prec,lgf,'RL',Nr,theta, 7,myIter,myThid)
./model/src/checkpoint.F 543 WRITE_CHECKPOINT_NOW
        CALL MDSWRITEFIELD(fn,prec,lgf,'RL',Nr,theta,5,myIter,myThid)
./model/src/checkpoint.F 608 WRITE_CHECKPOINT_NOW
        CALL MNC_CW_RL_W('D',fn,0,0,'Temp', theta, myThid)
./verification/lab_sea/code_ad_seaice_only/exf_getclim.F 104 EXF_GETCLIM
           theta(i,j,1,bi,bj) = climsst(i,j,bi,bj)
./verification/lab_sea/code_ad_seaice_only/exf_getclim.F 107 EXF_GETCLIM
     &             theta(i,j,1,bi,bj) .EQ. 0. ) then
./verification/lab_sea/code_ad_seaice_only/exf_getclim.F 108 EXF_GETCLIM
                print *, 'ph-warn-exf-clim ', i, j, theta(i,j,1,bi,bj)
./model/src/thermodynamics.F 291 THERMODYNAMICS
     U                                  theta, gtNm,
./model/src/thermodynamics.F 344 THERMODYNAMICS
     I             uVel, vVel, wVel, theta,
./model/src/thermodynamics.F 542 THERMODYNAMICS
     I         theta,
./model/src/thermodynamics.F 654 THERMODYNAMICS
     I         kappaRk, wVel, theta,
./model/src/thermodynamics.F 748 THERMODYNAMICS
       CALL DEBUG_STATS_RL(Nr,theta,'Theta (THERMODYNAMICS)',myThid)
./model/src/calc_buoyancy.F 41 CALC_BUOYANCY
         buoy(i,j)=(theta(i,j,k,bi,bj)-tRef(K))/tRef(K)
./model/src/external_forcing_surf.F 227 EXTERNAL_FORCING_SURF
     &         *( temp_EvPrRn - theta(i,j,ks,bi,bj) )
./model/src/external_forcing_surf.F 79 EXTERNAL_FORCING_SURF
     &         *(theta(i,j,ks,bi,bj)-SST(i,j,bi,bj))
./model/src/external_forcing_surf.F 90 EXTERNAL_FORCING_SURF
     &         *(theta(i,j,ks,bi,bj)-SST(i,j,bi,bj))
./model/src/convective_adjustment.F 120 CONVECTIVE_ADJUSTMENT
     I           theta,salt,
./model/src/convective_adjustment.F 131 CONVECTIVE_ADJUSTMENT
     I           theta,salt,
./model/src/convective_adjustment.F 154 CONVECTIVE_ADJUSTMENT
     U                              theta,
./model/src/calc_gt.F 151 CALC_GT
     I           gtNm(1-Olx,1-Oly,1,1,1,m2), theta,
./model/src/calc_gt.F 161 CALC_GT
     I           diffKhT, diffK4T, KappaRT, gtNm1, theta,
./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
./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)
./model/src/ini_theta.F 111 INI_THETA
            IF (theta(I,J,k,bi,bj) .LT. Tfreezing) THEN
./model/src/ini_theta.F 112 INI_THETA
               theta(I,J,K,bi,bj) = Tfreezing
./model/src/ini_theta.F 121 INI_THETA
      CALL PLOT_FIELD_XYZRL( theta, 'Initial Temperature' ,
./model/src/ini_theta.F 57 INI_THETA
           theta(I,J,K,bi,bj) = tRef(K)
./model/src/ini_theta.F 67 INI_THETA
       CALL READ_FLD_XYZ_RL( hydrogThetaFile, ' ', theta, 0, myThid )
./model/src/ini_theta.F 69 INI_THETA
       _EXCH_XYZ_R8(theta,myThid)
./model/src/ini_theta.F 79 INI_THETA
           IF (maskC(I,J,K,bi,bj).EQ.0.) theta(I,J,K,bi,bj) = 0.
./model/src/ini_theta.F 86 INI_THETA
     &      .AND. theta(I,J,K,bi,bj).EQ.0. ) THEN
./model/src/the_correction_step.F 101 THE_CORRECTION_STEP
     U           theta,gT,
./model/src/do_oceanic_phys.F 343 DO_OCEANIC_PHYS
     I        theta, salt,
./model/src/do_oceanic_phys.F 354 DO_OCEANIC_PHYS
     I        theta, salt,
./model/src/do_oceanic_phys.F 416 DO_OCEANIC_PHYS
     I            uVel, vVel, wVel, theta, salt,
./model/src/convective_adjustment_ini.F 127 CONVECTIVE_ADJUSTMENT_INI
     I           theta,salt,
./model/src/convective_adjustment_ini.F 140 CONVECTIVE_ADJUSTMENT_INI
     I           theta,salt,
./model/src/convective_adjustment_ini.F 163 CONVECTIVE_ADJUSTMENT_INI
     U                              theta,
./model/src/ini_pressure.F 108 INI_PRESSURE
     I          theta, salt, 
./model/src/diags_oceanic_surf_flux.F 88 DIAGS_OCEANIC_SURF_FLUX
     &                         *theta(i,j,ks,bi,bj)*HeatCapacity_Cp
./pkg/ggl90/ggl90_calc.F 152 GGL90_CALC
     I      theta, salt,
./pkg/ggl90/ggl90_calc.F 157 GGL90_CALC
     I      theta, salt,
./pkg/cost/cost_accumulate_mean.F 50 COST_ACCUMULATE_MEAN
     &                                    + theta(i,j,k,bi,bj)
./pkg/cost/cost_state_final.F 34 COST_STATE_FINAL
           objf_state_final(i,j,bi,bj,0*Nr+k) = theta(i,j,k,bi,bj)
./pkg/cost/cost_test.F 72 COST_TEST
               objf_test(bi,bj) = theta(i,j,kLocOut,bi,bj)
./pkg/cost/cost_test.F 75 COST_TEST
     &              theta(i,j,kLocOut,bi,bj)
./pkg/kpp/kpp_routines.F 1458 STATEKPP
     I     theta, salt,
./pkg/kpp/kpp_routines.F 1508 STATEKPP
     I        theta, salt,
./pkg/kpp/kpp_routines.F 1518 STATEKPP
     I        theta, salt,
./pkg/kpp/kpp_routines.F 1528 STATEKPP
     I        theta, salt,
./pkg/kpp/kpp_routines.F 1563 STATEKPP
     I           bi, bj, 1, sNx, 1, sNy, K-1, 1, theta, salt,
./pkg/kpp/kpp_routines.F 1567 STATEKPP
     I           bi, bj, 1, sNx, 1, sNy, K  , 1, theta, salt,
./pkg/opps/opps_interface.F 117 OPPS_INTERFACE
           tracerLoc(K,1) = theta(I,J,K,bi,bj)
./pkg/opps/opps_interface.F 123 OPPS_INTERFACE
           tracerLoc(K,1) = theta(I,J,K,bi,bj)
./pkg/opps/opps_interface.F 181 OPPS_INTERFACE
          theta(I,J,K,bi,bj) = tracerLoc(K,1)
./pkg/cfc/cfc12_surfforcing.F 54 CFC12_SURFFORCING
     &          + sca_12_2 * theta(i,j,k,bi,bj)
./pkg/cfc/cfc12_surfforcing.F 55 CFC12_SURFFORCING
     &          + sca_12_3 * theta(i,j,k,bi,bj)*theta(i,j,k,bi,bj)  
./pkg/cfc/cfc12_surfforcing.F 56 CFC12_SURFFORCING
     &          + sca_12_4 * theta(i,j,k,bi,bj)*theta(i,j,k,bi,bj) 
./pkg/cfc/cfc12_surfforcing.F 57 CFC12_SURFFORCING
     &                *theta(i,j,k,bi,bj)
./pkg/cfc/cfc12_surfforcing.F 59 CFC12_SURFFORCING
              ttemp=( theta(i,j,k,bi,bj) + 273.16)* 0.01
./pkg/cfc/cfc11_surfforcing.F 54 CFC11_SURFFORCING
     &          + sca_11_2 * theta(i,j,k,bi,bj)
./pkg/cfc/cfc11_surfforcing.F 55 CFC11_SURFFORCING
     &          + sca_11_3 * theta(i,j,k,bi,bj)*theta(i,j,k,bi,bj)  
./pkg/cfc/cfc11_surfforcing.F 56 CFC11_SURFFORCING
     &          + sca_11_4 * theta(i,j,k,bi,bj)*theta(i,j,k,bi,bj) 
./pkg/cfc/cfc11_surfforcing.F 57 CFC11_SURFFORCING
     &                *theta(i,j,k,bi,bj)
./pkg/cfc/cfc11_surfforcing.F 59 CFC11_SURFFORCING
              ttemp=( theta(i,j,k,bi,bj) + 273.16)* 0.01
./pkg/timeave/timeave_statvars.F 106 TIMEAVE_STATVARS
      CALL TIMEAVE_CUMULATE(thetatave,theta, Nr, DDTT, bi, bj, myThid)
./pkg/timeave/timeave_statvars.F 115 TIMEAVE_STATVARS
      CALL TIMEAVE_CUMUL_2V(