TEMPO Microphysics for NCAR MPAS

.gitignore

@@ -19,6 +19,7 @@ src/operators/*.f90
 # Executables
 *_model
 build_tables
+tempo_build_tables
 
 # NetCDF Files
 *.nc

src/core_atmosphere/Externals.cfg

@@ -12,5 +12,12 @@ repo_url = https://github.com/NOAA-GSL/UGWP.git
 tag = MPAS_20241223
 required = True
 
+[TEMPO_MP]
+local_path = ./physics_noaa/TEMPO
+protocol = git
+repo_url = https://github.com/NCAR/TEMPO.git
+tag = tempo_v3.1.0
+required = True
+
 [externals_description]
 schema_version = 1.0.0

src/core_atmosphere/Makefile

@@ -52,6 +52,7 @@ physcore: mpas_atm_dimensions.o
 	( cd ../..; ln -sf ./src/core_atmosphere/physics/physics_wrf/files/*TBL .)
 	( cd ../..; ln -sf ./src/core_atmosphere/physics/physics_wrf/files/*DATA* .)
 	( cd ../..; ln -sf ./src/core_atmosphere/physics/physics_noahmp/parameters/*TBL .)
+	( cd ../..; ln -sf ./src/core_atmosphere/physics/physics_noaa/TEMPO/tables/ccn_activate.bin .)
 
 chemcore:
 	( cd chemistry; $(MAKE) all CHEMISTRY="$(CHEMISTRY)" )
@@ -83,6 +84,7 @@ clean:
 	( cd utils; $(MAKE) clean )
 	( cd ../..; rm -f *TBL )
 	( cd ../..; rm -f *DATA* )
+	( cd ../..; rm -f ccn_activate.bin )
 	$(RM) -r libphys
 	$(RM) *.o *.mod *.f90 libdycore.a
 	$(RM) Registry_processed.xml
@@ -95,7 +97,7 @@ clean:
 	$(RM) $@ $*.mod
 ifeq "$(GEN_F90)" "true"
 	$(CPP) $(CPPFLAGS) $(PHYSICS) $(CHEMISTRY) $(CPPINCLUDES) -I./inc -I../framework $< > $*.f90
-	$(FC) $(FFLAGS) -c $*.f90 $(FCINCLUDES) -I../framework -I../operators -I./physics -I./dynamics -I./diagnostics -I./physics/physics_wrf -I./physics/physics_mmm -I./physics/physics_noaa/UGWP $(MPAS_ESMF_INC) -I./chemistry
+	$(FC) $(FFLAGS) -c $*.f90 $(FCINCLUDES) -I../framework -I../operators -I./physics -I./dynamics -I./diagnostics -I./physics/physics_wrf -I./physics/physics_mmm -I./physics/physics_noaa/UGWP -I./physics/physics_noaa/TEMPO $(MPAS_ESMF_INC) -I./chemistry
 else
-	$(FC) $(CPPFLAGS) $(PHYSICS) $(CHEMISTRY) $(FFLAGS) -c $*.F $(CPPINCLUDES) $(FCINCLUDES) -I./inc -I../framework -I../operators -I./physics -I./dynamics -I./diagnostics -I./physics/physics_wrf -I./physics/physics_mmm -I./physics/physics_noaa/UGWP $(MPAS_ESMF_INC) -I./chemistry
+	$(FC) $(CPPFLAGS) $(PHYSICS) $(CHEMISTRY) $(FFLAGS) -c $*.F $(CPPINCLUDES) $(FCINCLUDES) -I./inc -I../framework -I../operators -I./physics -I./dynamics -I./diagnostics -I./physics/physics_wrf -I./physics/physics_mmm -I./physics/physics_noaa/UGWP -I./physics/physics_noaa/TEMPO $(MPAS_ESMF_INC) -I./chemistry
 endif

src/core_atmosphere/physics/Makefile

@@ -6,7 +6,7 @@ endif
 
 all:
 	./../tools/manage_externals/checkout_externals --externals ./../Externals.cfg
-	$(MAKE) lookup_tables core_physics_init core_physics_mmm core_UGWP core_physics_wrf core_physics_noahmp core_physics
+	$(MAKE) lookup_tables core_physics_init core_physics_mmm core_UGWP core_TEMPO core_physics_wrf core_physics_noahmp core_physics
 
 dummy:
 	echo "****** compiling physics ******"
@@ -63,6 +63,9 @@ core_physics_mmm: core_physics_init
 core_UGWP: core_physics_init
 	(cd physics_noaa/UGWP; $(MAKE) all)
 
+core_TEMPO: core_physics_init
+	(cd physics_noaa/TEMPO; $(MAKE) -f Makefile.mpas)
+
 core_physics_wrf: core_physics_init core_physics_mmm core_UGWP
 	(cd physics_wrf; $(MAKE) all COREDEF="$(COREDEF)")
 
@@ -254,6 +257,7 @@ clean:
 	( cd physics_noahmp/src; $(MAKE) clean )
 	( cd physics_noahmp/utility; $(MAKE) clean )
 	( if [ -d physics_noaa/UGWP ]; then cd physics_noaa/UGWP; $(MAKE) clean; fi )
+	( if [ -d physics_noaa/TEMPO ]; then cd physics_noaa/TEMPO; $(MAKE) -f Makefile.mpas clean; fi )
 	@# Certain systems with intel compilers generate *.i files
 	@# This removes them during the clean process
 	$(RM) *.i
@@ -262,7 +266,7 @@ clean:
 	$(RM) $@ $*.mod
 ifeq "$(GEN_F90)" "true"
 	$(CPP) $(CPPFLAGS) $(COREDEF) $(HYDROSTATIC) $(CPPINCLUDES) -I../../framework $< > $*.f90
-	$(FC) $(FFLAGS) -c $*.f90 $(FCINCLUDES) -I./physics_mmm -I./physics_wrf -I./physics_noahmp -I./physics_noahmp/utility -I./physics_noahmp/drivers/mpas -I./physics_noahmp/src -I./physics_noaa/UGWP -I.. -I../../framework $(MPAS_ESMF_INC)
+	$(FC) $(FFLAGS) -c $*.f90 $(FCINCLUDES) -I./physics_mmm -I./physics_wrf -I./physics_noaa/TEMPO -I./physics_noahmp -I./physics_noahmp/utility -I./physics_noahmp/drivers/mpas -I./physics_noahmp/src -I./physics_noaa/UGWP -I.. -I../../framework $(MPAS_ESMF_INC)
 else
-	$(FC) $(CPPFLAGS) $(COREDEF) $(HYDROSATIC) $(FFLAGS) -c $*.F $(CPPINCLUDES) $(FCINCLUDES) -I./physics_mmm -I./physics_wrf -I./physics_noahmp -I./physics_noahmp/utility -I./physics_noahmp/drivers/mpas -I./physics_noahmp/src -I./physics_noaa/UGWP -I.. -I../../framework $(MPAS_ESMF_INC)
+	$(FC) $(CPPFLAGS) $(COREDEF) $(HYDROSATIC) $(FFLAGS) -c $*.F $(CPPINCLUDES) $(FCINCLUDES) -I./physics_mmm -I./physics_wrf -I./physics_noaa/TEMPO -I./physics_noahmp -I./physics_noahmp/utility -I./physics_noahmp/drivers/mpas -I./physics_noahmp/src -I./physics_noaa/UGWP -I.. -I../../framework $(MPAS_ESMF_INC)
 endif

src/core_atmosphere/physics/Registry_tempo.xml

@@ -0,0 +1,17 @@
+<!-- ========================================================================================================= -->
+<!-- CONFIGURATION OPTIONS AND ARRAYS SPECIFIC TO TEMPO:                                                       -->
+<!-- ========================================================================================================= -->
+
+  <nml_record name="physics_mp_tempo" in_defaults="true">
+    <nml_option name="config_tempo_aerosolaware" type="logical" default_value="true" in_defaults="false"
+      units="-"
+      description="Logical flag to turn on/off prognostic cloud droplet and aerosol number concentrations"
+      possible_values=".true. or .false."/>
+
+    <nml_option name="config_tempo_hailaware" type="logical" default_value="true" in_defaults="false"
+      units="-"
+      description="Logical flag to turn on/off prognostic graupel number concentration and rime density"
+      possible_values=".true. or .false."/>
+  </nml_record>
+
+<!-- ========================================================================================================= -->

src/core_atmosphere/physics/mpas_atmphys_control.F

@@ -177,6 +177,7 @@ subroutine physics_namelist_check(configs)
            config_microp_scheme .eq. 'mp_kessler'           .or. &
            config_microp_scheme .eq. 'mp_thompson'          .or. &
            config_microp_scheme .eq. 'mp_thompson_aerosols' .or. &
+           config_microp_scheme .eq. 'mp_tempo'             .or. &
            config_microp_scheme .eq. 'mp_wsm6')) then
 
     write(mpas_err_message,'(A,A20)') 'illegal value for config_microp_scheme:', &