Compute temperature and spechum at output times via diagnostics

Author: jihyeonjang

src/core_atmosphere/Registry.xml

@@ -1767,8 +1767,7 @@
                      description="Relative humidity"/>
 
                 <var name="spechum" type="real" dimensions="nVertLevels nCells Time" units="kg kg^{-1}"
-                     description="Specific humidity"
-                     packages="jedi_da"/>
+                     description="Specific humidity"/>
 
                 <var name="v" type="real" dimensions="nVertLevels nEdges Time" units="m s^{-1}"
                      description="Horizontal tangential velocity at edges"/>

src/core_atmosphere/diagnostics/Makefile

@@ -10,6 +10,7 @@ DIAGNOSTIC_MODULES = \
 	mpas_convective_diagnostics.o \
 	mpas_pv_diagnostics.o \
 	mpas_soundings.o \
+	mpas_modellevel_diagnostics.o \
 
 mpas_isobaric_diagnostics.o: mpas_atm_diagnostics_utils.o
 
@@ -21,6 +22,7 @@ mpas_pv_diagnostics.o: mpas_atm_diagnostics_utils.o
 
 mpas_soundings.o:
 
+mpas_modellevel_diagnostics.o: mpas_atm_diagnostics_utils.o
 
 ################### Generally no need to modify below here ###################
 

src/core_atmosphere/diagnostics/Registry_diagnostics.xml

@@ -19,6 +19,9 @@
 <!-- soundings -->
 #include "Registry_soundings.xml"
 
+<!-- additional model-level diagnostics -->
+#include "Registry_modellevel.xml"
+
 <!-- ******************************* -->
 <!-- End includes from diagnostics -->
 <!-- ******************************* -->

src/core_atmosphere/diagnostics/mpas_atm_diagnostics_manager.F

@@ -38,6 +38,7 @@ subroutine mpas_atm_diag_setup(stream_mgr, configs, structs, clock, dminfo)
         use mpas_convective_diagnostics, only : convective_diagnostics_setup
         use mpas_pv_diagnostics, only : pv_diagnostics_setup
         use mpas_soundings, only : soundings_setup
+        use mpas_modellevel_diagnostics, only : modellevel_diagnostics_setup
 
         implicit none
 
@@ -59,6 +60,7 @@ subroutine mpas_atm_diag_setup(stream_mgr, configs, structs, clock, dminfo)
         call convective_diagnostics_setup(structs, clock)
         call pv_diagnostics_setup(structs, clock)
         call soundings_setup(configs, structs, clock, dminfo)
+        call modellevel_diagnostics_setup(structs, clock)
 
     end subroutine mpas_atm_diag_setup
 
@@ -105,6 +107,7 @@ subroutine mpas_atm_diag_compute()
         use mpas_convective_diagnostics, only : convective_diagnostics_compute
         use mpas_pv_diagnostics, only : pv_diagnostics_compute
         use mpas_soundings, only : soundings_compute
+        use mpas_modellevel_diagnostics, only : modellevel_diagnostics_compute
 
         implicit none
 
@@ -115,6 +118,7 @@ subroutine mpas_atm_diag_compute()
         call convective_diagnostics_compute()
         call pv_diagnostics_compute()
         call soundings_compute()
+        call modellevel_diagnostics_compute()
 
     end subroutine mpas_atm_diag_compute
 

src/core_atmosphere/mpas_atm_core.F

@@ -918,7 +918,6 @@ subroutine atm_compute_output_diagnostics(state, time_lev, diag, mesh)
       integer, pointer :: nCells, nVertLevels, index_qv
       real (kind=RKIND), dimension(:,:), pointer :: theta, rho, theta_m, rho_zz, zz
       real (kind=RKIND), dimension(:,:), pointer :: pressure_base, pressure_p, pressure
-      real (kind=RKIND), dimension(:,:), pointer :: exner, temperature
       real (kind=RKIND), dimension(:,:,:), pointer :: scalars
 
       call mpas_pool_get_dimension(mesh, 'nCells', nCells)
@@ -934,8 +933,6 @@ subroutine atm_compute_output_diagnostics(state, time_lev, diag, mesh)
       call mpas_pool_get_array(diag, 'pressure_p', pressure_p)
       call mpas_pool_get_array(diag, 'pressure_base', pressure_base)
       call mpas_pool_get_array(diag, 'pressure', pressure)
-      call mpas_pool_get_array(diag, 'exner', exner)
-      call mpas_pool_get_array(diag, 'temperature', temperature)
 
       call mpas_pool_get_array(mesh, 'zz', zz)
 
@@ -944,7 +941,6 @@ subroutine atm_compute_output_diagnostics(state, time_lev, diag, mesh)
             theta(k,iCell) = theta_m(k,iCell) / (1._RKIND + rvord * scalars(index_qv,k,iCell))
             rho(k,iCell) = rho_zz(k,iCell) * zz(k,iCell)
             pressure(k,iCell) = pressure_base(k,iCell) + pressure_p(k,iCell)
-            temperature(k,iCell) = theta(k,iCell) * exner(k,iCell)
          end do
       end do