THM is mandatory WRF temp variable

models/wrf/readme.rst

@@ -24,7 +24,10 @@ Some important WRF-DART updates include:
   operator calculations.
 
 - Version 11.5.0: Improves compatibility with WRFv4+ versions where
-  the prognostic 3D temperature variable is THM.
+  the prognostic 3D temperature variable is THM.  It is now mandatory to
+  include THM instead of T in the ``wrf_state_variables`` namelist.
+
+
 
 It is always recommended that you update your DART version to the 
 `latest release <https://github.com/NCAR/DART/releases>`__ before beginning new research.
@@ -336,6 +339,14 @@ For example:
                            'PSFC','QTY_PRESSURE','TYPE_PS','UPDATE','999',
 
 
+.. Important::
+
+   It is mandatory to include ``THM`` instead of ``T`` as the ``TYPE_T`` WRF 
+   temperature variable. This is because ``THM`` is the prognostic temperature variable
+   that will impact the forecast when updated.  Also, including ``T`` can give 
+   boundscheck errors as described in `issue #728. <https://github.com/NCAR/DART/issues/728>`__  
+
+
 - polar, periodic_x
 
 The ``Polar`` and ``periodic_x`` namelist values are used in global WRF simulations.

models/wrf/tutorial/README.rst

@@ -19,12 +19,16 @@ either WRF or DART.
   versions 11.4.0 and earlier because those older versions do not account
   for different coordinate systems including the sigma hybrid coordinates as 
   described in `DART Issue #650 <https://github.com/NCAR/DART/pull/650>`__.
+  
   Furthermore, older versions do not account for the prognostic temperature variable
   switch from ``T`` (perturbation potential temperature) to ``THM``, (either perturbation
   potential temperature or perturbation moist potential temperature) as described in
   `DART issue #661 <https://github.com/NCAR/DART/issues/661>`__. The current implementation
   of the code sets ``T=THM`` because within &dynamics section of ``namelist.input``
-  ``use_theta_m=0``.
+  ``use_theta_m=0``.  For this reason, It is mandatory to include ``THM`` instead of 
+  ``T`` as the ``TYPE_T`` within the wrf_state_variables namelist. Furthermore, including
+  ``T`` in the namelist can lead to errors in the WRF model_mod bounds check as 
+  described in `Issue #728. <https://github.com/NCAR/DART/issues/728>`__
 
   Earlier version of WRF (v3.9) may run without errors with more recent versions of
   DART (later than 11.4.0), but the assimilation performance will be deprecated.