Merge pull request #838 from NCAR/WRF_THM

Check for THM USE_THETA_M attribute

Author: hkershaw-brown

CHANGELOG.rst

@@ -22,6 +22,11 @@ individual files.
 
 The changes are now listed with the most recent at the top.
 
+**March 14 2025 :: WRF use_theta_m check. Tag v11.10.5**
+
+- Assert USE_THETA_M=0 in WRF when initializing wrf model_mod.
+- fixed documentation broken links.
+
 **February 6 2025 :: WRF tutorial and prepbufr documentation. Tag v11.10.4**
 
 - Observation section of WRF tutorial updated to match provided obs.   

conf.py

@@ -21,7 +21,7 @@
 author = 'Data Assimilation Research Section'
 
 # The full version, including alpha/beta/rc tags
-release = '11.10.4'
+release = '11.10.5'
 root_doc = 'index'
 
 # -- General configuration ---------------------------------------------------

models/wrf/model_mod.f90

@@ -314,6 +314,7 @@ module model_mod
    integer  :: domain_size
    integer  :: localization_coord
    integer  :: hybrid_opt
+   integer  :: theta_m
    real(r8), dimension(:),     pointer :: znu, dn, dnw, zs, znw, c1h, c2h
    real(r8), dimension(:,:),   pointer :: mub, hgt
    real(r8), dimension(:,:),   pointer :: latitude, latitude_u, latitude_v
@@ -708,7 +709,16 @@ subroutine static_init_model()
                         ' for WRFv4 and later'
       call error_handler(E_ERR,'static_init_model', errstring, source, revision, revdate)
 
-   endif        
+   endif
+
+   ! DART requires that THM is the 'perturbed dry adiabatic potential temperature'
+   ! WRF USE_THETA_M namelist option must be 0
+   if (get_type_ind_from_type_string(id,'THM') >=0 .and. wrf%dom(id)%theta_m /= 0) then
+
+      write(errstring,*)'WRF namelist option USE_THETA_M must be 0, where THM = perturbed dry adiabatic potential temperature' 
+      call error_handler(E_ERR,'static_init_model', errstring, source, revision, revdate)
+
+   endif   
 
    wrf%dom(id)%type_u      = get_type_ind_from_type_string(id,'U')
    wrf%dom(id)%type_v      = get_type_ind_from_type_string(id,'V')
@@ -7159,6 +7169,12 @@ subroutine read_wrf_file_attributes(ncid,id)
                      'static_init_model', 'get_att STAND_LON')
    if(debug) write(*,*) ' stdlon is ',stdlon
 
+   call nc_check( nf90_get_att(ncid, nf90_global, 'USE_THETA_M', wrf%dom(id)%theta_m), &
+                     'static_init_model', 'get_att USE_THETA_M')
+   if(debug) write(*,*) ' theta_m is ',wrf%dom(id)%theta_m
+
+
+
    ! this attribute is not present in older wrf files, which means it is not
    ! using a hybrid vertical coordinate.  not having this attribute should
    ! not cause an error.