diff --git a/constants/Makefile.am b/constants/Makefile.am
index b8e141a05d..429dd6c329 100644
--- a/constants/Makefile.am
+++ b/constants/Makefile.am
@@ -32,12 +32,12 @@ noinst_LTLIBRARIES = libconstants.la
# The convenience library depends on its source.
libconstants_la_SOURCES = \
fmsconstants.F90 \
- gfdl_constants.h \
- gfs_constants.h \
- geos_constants.h \
+ gfdl_constants.fh \
+ gfs_constants.fh \
+ geos_constants.fh \
constants.F90
-FMSconstants.$(FC_MODEXT): gfdl_constants.h gfs_constants.h geos_constants.h
+FMSconstants.$(FC_MODEXT): gfdl_constants.fh gfs_constants.fh geos_constants.fh
constants_mod.$(FC_MODEXT): fmsconstants.$(FC_MODEXT)
# Mod files are built and then installed as headers
diff --git a/constants/fmsconstants.F90 b/constants/fmsconstants.F90
index 3abc5da70a..6734f88166 100644
--- a/constants/fmsconstants.F90
+++ b/constants/fmsconstants.F90
@@ -28,9 +28,9 @@
!! statement.
!!
!! The currently support contant systems are:
-!! GFDL constants (gfdl_constants.h)
-!! GEOS constants (geos_constants.h)
-!! GFS constants (gfs_constants.h)
+!! GFDL constants (gfdl_constants.fh)
+!! GEOS constants (geos_constants.fh)
+!! GFS constants (gfs_constants.fh)
!!
!!
!! The name given to a particular constant may be changed.
@@ -71,13 +71,13 @@ module FMSconstants
!--- perform error checking and include the correct system of constants
#if defined(GFDL_CONSTANTS) && !defined(GFS_CONSTANTS) && !defined(GEOS_CONSTANTS)
#warning "Using GFDL constants"
-#include
+#include
#elif !defined(GFDL_CONSTANTS) && defined(GFS_CONSTANTS) && !defined(GEOS_CONSTANTS)
#warning "Using GFS constants"
-#include
+#include
#elif !defined(GFDL_CONSTANTS) && !defined(GFS_CONSTANTS) && defined(GEOS_CONSTANTS)
#warning "Using GEOS constants"
-#include
+#include
#else
#error FATAL FMSConstants error - multiple constants macros are defined for FMS
#endif
diff --git a/constants/geos_constants.h b/constants/geos_constants.fh
similarity index 73%
rename from constants/geos_constants.h
rename to constants/geos_constants.fh
index edc298b41e..ce2939985e 100644
--- a/constants/geos_constants.h
+++ b/constants/geos_constants.fh
@@ -23,11 +23,14 @@ character(len=18), public, parameter :: constants_version = 'FMSConstants: GEOS'
real(kind=RKIND), public, parameter :: small_fac = 1._r8_kind
!--- Spherical coordinate conversion constants
-real(kind=r8_kind), public, parameter :: PI_8 = 3.14159265358979323846_r8_kind !< Ratio of circle circumference to diameter [N/A]
-real(kind=RKIND), public, parameter :: PI = PI_8 !< Ratio of circle circumference to diameter [N/A]
+real(kind=r8_kind), public, parameter :: PI_8 = 3.14159265358979323846_r8_kind !< Ratio of circle circumference
+ !! to diameter [N/A]
+real(kind=RKIND), public, parameter :: PI = PI_8 !< Ratio of circle circumference
+ !! to diameter [N/A]
real(kind=RKIND), public, parameter :: RAD_TO_DEG = 180._r8_kind/PI_8 !< Degrees per radian [deg/rad]
real(kind=RKIND), public, parameter :: DEG_TO_RAD = PI_8/180._r8_kind !< Radians per degree [rad/deg]
-real(kind=RKIND), public, parameter :: RADIAN = RAD_TO_DEG !< Equal to RAD_TO_DEG for backward compatability. [rad/deg]
+real(kind=RKIND), public, parameter :: RADIAN = RAD_TO_DEG !< Equal to RAD_TO_DEG for backward
+ !!compatability. [rad/deg]
!--- Earth physical constants
real(kind=RKIND), public, parameter :: RADIUS = 6371.0E3_r8_kind !< Radius of the Earth [m]
@@ -39,7 +42,8 @@ real(kind=RKIND), public, parameter :: SECONDS_PER_MINUTE = 60._r8_kind !
!--- Various gas constants
real(kind=RKIND), public, parameter :: RDGAS = 8314.47 /28.965 !< Gas constant for dry air [J/kg/deg]
-real(kind=RKIND), public, parameter :: RVGAS = 8314.47 /18.015 !< Gas constant for water vapor [J/kg/deg]
+real(kind=RKIND), public, parameter :: RVGAS = 8314.47 /18.015 !< Gas constant for water vapor
+ !! [J/kg/deg]
real(kind=RKIND), public, parameter :: HLV = 2.4665E6_r8_kind !< Latent heat of evaporation [J/kg]
real(kind=RKIND), public, parameter :: HLF = 3.3370E5_r8_kind !< Latent heat of fusion [J/kg]
real(kind=RKIND), public, parameter :: HLS = HLV + HLF !< Latent heat of sublimation [J/kg]
@@ -48,30 +52,38 @@ real(kind=RKIND), public, parameter :: CP_AIR = RDGAS/KAPPA !< S
!! at constant pressure [J/kg/deg]
real(kind=RKIND), public, parameter :: CP_VAPOR = 4.0_r8_kind*RVGAS !< Specific heat capacity of water vapor
!! at constant pressure [J/kg/deg]
-real(kind=RKIND), public, parameter :: CP_OCEAN = 3989.24495292815_r8_kind !< Specific heat capacity taken from McDougall (2002)
+real(kind=RKIND), public, parameter :: CP_OCEAN = 3989.24495292815_r8_kind !< Specific heat capacity taken
+ !! from McDougall (2002)
!! "Potential Enthalpy ..." [J/kg/deg]
real(kind=RKIND), public, parameter :: DENS_H2O = 1000._r8_kind !< Density of liquid water [kg/m^3]
real(kind=RKIND), public, parameter :: RHOAIR = 1.292269_r8_kind !< Reference atmospheric density [kg/m^3]
real(kind=RKIND), public, parameter :: RHO0 = 1.035E3_r8_kind !< Average density of sea water [kg/m^3]
-real(kind=RKIND), public, parameter :: RHO0R = 1.0_r8_kind/RHO0 !< Reciprocal of average density of sea water [m^3/kg]
+real(kind=RKIND), public, parameter :: RHO0R = 1.0_r8_kind/RHO0 !< Reciprocal of average density
+ !! of sea water [m^3/kg]
real(kind=RKIND), public, parameter :: RHO_CP = RHO0*CP_OCEAN !< (kg/m^3)*(cal/kg/deg C)(joules/cal) =
!! (joules/m^3/deg C) [J/m^3/deg]
-real(kind=RKIND), public, parameter :: O2MIXRAT = 2.0953E-01_r8_kind !< Mixing ratio of molecular oxygen in air [dimensionless]
+real(kind=RKIND), public, parameter :: O2MIXRAT = 2.0953E-01_r8_kind !< Mixing ratio of molecular oxygen
+ !! in air [dimensionless]
real(kind=RKIND), public, parameter :: WTMAIR = 2.896440E+01_r8_kind !< Molecular weight of air [AMU]
real(kind=RKIND), public, parameter :: WTMH2O = WTMAIR*(RDGAS/RVGAS) !< Molecular weight of water [AMU]
real(kind=RKIND), public, parameter :: WTMOZONE = 47.99820_r8_kind !< Molecular weight of ozone [AMU]
real(kind=RKIND), public, parameter :: WTMC = 12.00000_r8_kind !< Molecular weight of carbon [AMU]
-real(kind=RKIND), public, parameter :: WTMCO2 = 44.00995_r8_kind !< Molecular weight of carbon dioxide [AMU]
+real(kind=RKIND), public, parameter :: WTMCO2 = 44.00995_r8_kind !< Molecular weight of carbon dioxide
+ !! [AMU]
real(kind=RKIND), public, parameter :: WTMCH4 = 16.0425_r8_kind !< Molecular weight of methane [AMU]
-real(kind=RKIND), public, parameter :: WTMO2 = 31.9988_r8_kind !< Molecular weight of molecular oxygen [AMU]
-real(kind=RKIND), public, parameter :: WTMCFC11 = 137.3681_r8_kind !< Molecular weight of CFC-11 (CCl3F) [AMU]
-real(kind=RKIND), public, parameter :: WTMCFC12 = 120.9135_r8_kind !< Molecular weight of CFC-21 (CCl2F2) [AMU]
+real(kind=RKIND), public, parameter :: WTMO2 = 31.9988_r8_kind !< Molecular weight of molecular oxygen
+ !! [AMU]
+real(kind=RKIND), public, parameter :: WTMCFC11 = 137.3681_r8_kind !< Molecular weight of CFC-11
+ !! (CCl3F) [AMU]
+real(kind=RKIND), public, parameter :: WTMCFC12 = 120.9135_r8_kind !< Molecular weight of CFC-21
+ !! (CCl2F2) [AMU]
real(kind=RKIND), public, parameter :: WTMN = 14.0067_r8_kind !< Molecular weight of Nitrogen [AMU]
real(kind=RKIND), public, parameter :: DIFFAC = 1.660_r8_kind !< Diffusivity factor [dimensionless]
real(kind=RKIND), public, parameter :: ES0 = 1.0_r8_kind !< Humidity factor [dimensionless]
- !! Controls the humidity content of the atmosphere through
- !! the Saturation Vapour Pressure expression
- !! when using DO_SIMPLE
+ !! Controls the humidity content of
+ !! the atmosphere through
+ !! the Saturation Vapour Pressure
+ !! expression when using DO_SIMPLE
!--- Pressure and Temperature constants
real(kind=RKIND), public, parameter :: PSTD = 1.013250E+06_r8_kind !< Mean sea level pressure [dynes/cm^2]
@@ -87,11 +99,15 @@ real(kind=RKIND), public, parameter :: AVOGNO = 6.023000E+23_r8_kind !< A
real(kind=RKIND), public, parameter :: VONKARM = 0.40_r8_kind !< Von Karman constant [dimensionless]
!--- Miscellaneous constants
-real(kind=RKIND), public, parameter :: ALOGMIN = -50.0_r8_kind !< Minimum value allowed as argument to log function [N/A]
-real(kind=RKIND), public, parameter :: EPSLN = 1.0E-40_r8_kind !< A small number to prevent divide by zero exceptions [N/A]
-real(kind=RKIND), public, parameter :: RADCON = ((1.0E+02*GRAV)/(1.0D+04*CP_AIR))*SECONDS_PER_DAY !< Factor to convert flux divergence
- !! to heating rate in degrees per day
- !! [deg sec/(cm day)]
-real(kind=RKIND), public, parameter :: RADCON_MKS = (GRAV/CP_AIR)*SECONDS_PER_DAY !< Factor to convert flux divergence
- !! to heating rate in degrees per day
- !! [deg sec/(m day)]
+real(kind=RKIND), public, parameter :: ALOGMIN = -50.0_r8_kind !< Minimum value allowed as argument
+ !! to log function [N/A]
+real(kind=RKIND), public, parameter :: EPSLN = 1.0E-40_r8_kind !< A small number to prevent divide
+ !! by zero exceptions [N/A]
+real(kind=RKIND), public, parameter :: RADCON = ((1.0E+02*GRAV)/(1.0D+04*CP_AIR))*SECONDS_PER_DAY !< Factor to
+ !! convert flux divergence
+ !! to heating rate in degrees per day
+ !! [deg sec/(cm day)]
+real(kind=RKIND), public, parameter :: RADCON_MKS = (GRAV/CP_AIR)*SECONDS_PER_DAY !< Factor to
+ !! convert flux divergence
+ !! to heating rate in degrees per day
+ !! [deg sec/(m day)]
diff --git a/constants/gfdl_constants.h b/constants/gfdl_constants.fh
similarity index 73%
rename from constants/gfdl_constants.h
rename to constants/gfdl_constants.fh
index 9e0b8812a1..333cf44051 100644
--- a/constants/gfdl_constants.h
+++ b/constants/gfdl_constants.fh
@@ -23,11 +23,14 @@ character(len=18), public, parameter :: constants_version = 'FMSConstants: GFDL'
real(kind=RKIND), public, parameter :: small_fac = 1._r8_kind
!--- Spherical coordinate conversion constants
-real(kind=r8_kind), public, parameter :: PI_8 = 3.14159265358979323846_r8_kind !< Ratio of circle circumference to diameter [N/A]
-real(kind=RKIND), public, parameter :: PI = PI_8 !< Ratio of circle circumference to diameter [N/A]
+real(kind=r8_kind), public, parameter :: PI_8 = 3.14159265358979323846_r8_kind !< Ratio of circle circumference to
+ !! diameter [N/A]
+real(kind=RKIND), public, parameter :: PI = PI_8 !< Ratio of circle circumference to
+ !! diameter [N/A]
real(kind=RKIND), public, parameter :: RAD_TO_DEG = 180._r8_kind/PI_8 !< Degrees per radian [deg/rad]
real(kind=RKIND), public, parameter :: DEG_TO_RAD = PI_8/180._r8_kind !< Radians per degree [rad/deg]
-real(kind=RKIND), public, parameter :: RADIAN = RAD_TO_DEG !< Equal to RAD_TO_DEG for backward compatability. [rad/deg]
+real(kind=RKIND), public, parameter :: RADIAN = RAD_TO_DEG !< Equal to RAD_TO_DEG for backward
+ !! compatability. [rad/deg]
!--- Earth physical constants
real(kind=RKIND), public, parameter :: RADIUS = 6371.0E+3_r8_kind !< Radius of the Earth [m]
@@ -39,7 +42,8 @@ real(kind=RKIND), public, parameter :: SECONDS_PER_MINUTE = 60._r8_kind !
!--- Various gas constants
real(kind=RKIND), public, parameter :: RDGAS = 287.04_r8_kind !< Gas constant for dry air [J/kg/deg]
-real(kind=RKIND), public, parameter :: RVGAS = 461.50_r8_kind !< Gas constant for water vapor [J/kg/deg]
+real(kind=RKIND), public, parameter :: RVGAS = 461.50_r8_kind !< Gas constant for water vapor
+ !! [J/kg/deg]
real(kind=RKIND), public, parameter :: HLV = 2.500E6_r8_kind !< Latent heat of evaporation [J/kg]
real(kind=RKIND), public, parameter :: HLF = 3.34E5_r8_kind !< Latent heat of fusion [J/kg]
real(kind=RKIND), public, parameter :: HLS = HLV + HLF !< Latent heat of sublimation [J/kg]
@@ -48,30 +52,38 @@ real(kind=RKIND), public, parameter :: CP_AIR = RDGAS/KAPPA !< S
!! at constant pressure [J/kg/deg]
real(kind=RKIND), public, parameter :: CP_VAPOR = 4.0_r8_kind*RVGAS !< Specific heat capacity of water vapor
!! at constant pressure [J/kg/deg]
-real(kind=RKIND), public, parameter :: CP_OCEAN = 3989.24495292815_r8_kind !< Specific heat capacity taken from McDougall (2002)
+real(kind=RKIND), public, parameter :: CP_OCEAN = 3989.24495292815_r8_kind !< Specific heat capacity taken
+ !! from McDougall (2002)
!! "Potential Enthalpy ..." [J/kg/deg]
real(kind=RKIND), public, parameter :: DENS_H2O = 1000._r8_kind !< Density of liquid water [kg/m^3]
real(kind=RKIND), public, parameter :: RHOAIR = 1.292269_r8_kind !< Reference atmospheric density [kg/m^3]
real(kind=RKIND), public, parameter :: RHO0 = 1.035E3_r8_kind !< Average density of sea water [kg/m^3]
-real(kind=RKIND), public, parameter :: RHO0R = 1.0_r8_kind/RHO0 !< Reciprocal of average density of sea water [m^3/kg]
+real(kind=RKIND), public, parameter :: RHO0R = 1.0_r8_kind/RHO0 !< Reciprocal of average density of
+ !! sea water [m^3/kg]
real(kind=RKIND), public, parameter :: RHO_CP = RHO0*CP_OCEAN !< (kg/m^3)*(cal/kg/deg C)(joules/cal) =
!! (joules/m^3/deg C) [J/m^3/deg]
-real(kind=RKIND), public, parameter :: O2MIXRAT = 2.0953E-01_r8_kind !< Mixing ratio of molecular oxygen in air [dimensionless]
+real(kind=RKIND), public, parameter :: O2MIXRAT = 2.0953E-01_r8_kind !< Mixing ratio of molecular oxygen
+ !! in air [dimensionless]
real(kind=RKIND), public, parameter :: WTMAIR = 2.896440E+01_r8_kind !< Molecular weight of air [AMU]
real(kind=RKIND), public, parameter :: WTMH2O = WTMAIR*(RDGAS/RVGAS) !< Molecular weight of water [AMU]
real(kind=RKIND), public, parameter :: WTMOZONE = 47.99820_r8_kind !< Molecular weight of ozone [AMU]
real(kind=RKIND), public, parameter :: WTMC = 12.00000_r8_kind !< Molecular weight of carbon [AMU]
-real(kind=RKIND), public, parameter :: WTMCO2 = 44.00995_r8_kind !< Molecular weight of carbon dioxide [AMU]
+real(kind=RKIND), public, parameter :: WTMCO2 = 44.00995_r8_kind !< Molecular weight of carbon dioxide
+ !! [AMU]
real(kind=RKIND), public, parameter :: WTMCH4 = 16.0425_r8_kind !< Molecular weight of methane [AMU]
-real(kind=RKIND), public, parameter :: WTMO2 = 31.9988_r8_kind !< Molecular weight of molecular oxygen [AMU]
-real(kind=RKIND), public, parameter :: WTMCFC11 = 137.3681_r8_kind !< Molecular weight of CFC-11 (CCl3F) [AMU]
-real(kind=RKIND), public, parameter :: WTMCFC12 = 120.9135_r8_kind !< Molecular weight of CFC-21 (CCl2F2) [AMU]
+real(kind=RKIND), public, parameter :: WTMO2 = 31.9988_r8_kind !< Molecular weight of molecular
+ !! oxygen [AMU]
+real(kind=RKIND), public, parameter :: WTMCFC11 = 137.3681_r8_kind !< Molecular weight of CFC-11
+ !! (CCl3F) [AMU]
+real(kind=RKIND), public, parameter :: WTMCFC12 = 120.9135_r8_kind !< Molecular weight of CFC-21
+ !! (CCl2F2) [AMU]
real(kind=RKIND), public, parameter :: WTMN = 14.0067_r8_kind !< Molecular weight of Nitrogen [AMU]
real(kind=RKIND), public, parameter :: DIFFAC = 1.660_r8_kind !< Diffusivity factor [dimensionless]
real(kind=RKIND), public, parameter :: ES0 = 1.0_r8_kind !< Humidity factor [dimensionless]
- !! Controls the humidity content of the atmosphere through
- !! the Saturation Vapour Pressure expression
- !! when using DO_SIMPLE
+ !! Controls the humidity content of
+ !! the atmosphere through
+ !! the Saturation Vapour Pressure
+ !! expression when using DO_SIMPLE
!--- Pressure and Temperature constants
real(kind=RKIND), public, parameter :: PSTD = 1.013250E+06_r8_kind !< Mean sea level pressure [dynes/cm^2]
@@ -87,11 +99,15 @@ real(kind=RKIND), public, parameter :: AVOGNO = 6.023000E+23_r8_kind !< A
real(kind=RKIND), public, parameter :: VONKARM = 0.40_r8_kind !< Von Karman constant [dimensionless]
!--- Miscellaneous constants
-real(kind=RKIND), public, parameter :: ALOGMIN = -50.0_r8_kind !< Minimum value allowed as argument to log function [N/A]
-real(kind=RKIND), public, parameter :: EPSLN = 1.0E-40_r8_kind !< A small number to prevent divide by zero exceptions [N/A]
-real(kind=RKIND), public, parameter :: RADCON = ((1.0E+02*GRAV)/(1.0D+04*CP_AIR))*SECONDS_PER_DAY !< Factor to convert flux divergence
- !! to heating rate in degrees per day
- !! [deg sec/(cm day)]
-real(kind=RKIND), public, parameter :: RADCON_MKS = (GRAV/CP_AIR)*SECONDS_PER_DAY !< Factor to convert flux divergence
- !! to heating rate in degrees per day
- !! [deg sec/(m day)]
+real(kind=RKIND), public, parameter :: ALOGMIN = -50.0_r8_kind !< Minimum value allowed as argument
+ !! to log function [N/A]
+real(kind=RKIND), public, parameter :: EPSLN = 1.0E-40_r8_kind !< A small number to prevent
+ !! divide by zero exceptions [N/A]
+real(kind=RKIND), public, parameter :: RADCON = ((1.0E+02*GRAV)/(1.0D+04*CP_AIR))*SECONDS_PER_DAY !< Factor to
+ !! convert flux divergence
+ !! to heating rate in degrees per day
+ !! [deg sec/(cm day)]
+real(kind=RKIND), public, parameter :: RADCON_MKS = (GRAV/CP_AIR)*SECONDS_PER_DAY !< Factor to
+ !! convert flux divergence
+ !! to heating rate in degrees per day
+ !! [deg sec/(m day)]
diff --git a/constants/gfs_constants.h b/constants/gfs_constants.fh
similarity index 73%
rename from constants/gfs_constants.h
rename to constants/gfs_constants.fh
index b732f6d43a..d5eb711739 100644
--- a/constants/gfs_constants.h
+++ b/constants/gfs_constants.fh
@@ -23,17 +23,22 @@ character(len=18), public, parameter :: constants_version = 'FMSConstants: GFS '
real(kind=RKIND), public, parameter :: small_fac = 1._r8_kind
!--- Spherical coordinate conversion constants
-real(kind=r8_kind), public, parameter :: PI_8 = 3.1415926535897931_r8_kind !< Ratio of circle circumference to diameter [N/A]
-real(kind=RKIND), public, parameter :: PI = PI_8 !< Ratio of circle circumference to diameter [N/A]
+real(kind=r8_kind), public, parameter :: PI_8 = 3.1415926535897931_r8_kind !< Ratio of circle circumference
+ !! to diameter [N/A]
+real(kind=RKIND), public, parameter :: PI = PI_8 !< Ratio of circle circumference
+ !! to diameter [N/A]
real(kind=RKIND), public, parameter :: RAD_TO_DEG = 180._r8_kind/PI_8 !< Degrees per radian [deg/rad]
real(kind=RKIND), public, parameter :: DEG_TO_RAD = PI_8/180._r8_kind !< Radians per degree [rad/deg]
-real(kind=RKIND), public, parameter :: RADIAN = RAD_TO_DEG !< Equal to RAD_TO_DEG for backward compatability. [rad/deg]
+real(kind=RKIND), public, parameter :: RADIAN = RAD_TO_DEG !< Equal to RAD_TO_DEG for backward
+ !! compatability. [rad/deg]
!--- Earth physical constants
real(kind=RKIND), public, parameter :: RADIUS = 6.3712E+6_r8_kind !< Radius of the Earth [m]
real(kind=RKIND), public, parameter :: OMEGA = 7.2921E-5_r8_kind !< Rotation rate of the Earth [1/s]
-real(kind=r8_kind), public, parameter :: GRAV_8 = 9.80665_r8_kind !< Acceleration due to gravity [m/s^2] (REAL(KIND=8))
-real(kind=RKIND), public, parameter :: GRAV = GRAV_8 !< Acceleration due to gravity [m/s^2]
+real(kind=r8_kind), public, parameter :: GRAV_8 = 9.80665_r8_kind !< Acceleration due to gravity
+ !! [m/s^2] (REAL(KIND=8))
+real(kind=RKIND), public, parameter :: GRAV = GRAV_8 !< Acceleration due to gravity
+ !! [m/s^2]
real(kind=RKIND), public, parameter :: SECONDS_PER_DAY = 86400._r8_kind !< Seconds in a day [s]
real(kind=RKIND), public, parameter :: SECONDS_PER_HOUR = 3600._r8_kind !< Seconds in an hour [s]
real(kind=RKIND), public, parameter :: SECONDS_PER_MINUTE = 60._r8_kind !< Seconds in a minute [s]
@@ -48,31 +53,39 @@ real(kind=RKIND), public, parameter :: CP_AIR = 1004.6_r8_kind !< S
!! at constant pressure [J/kg/deg]
real(kind=RKIND), public, parameter :: CP_VAPOR = 4.0_r8_kind*RVGAS !< Specific heat capacity of water vapor
!! at constant pressure [J/kg/deg]
-real(kind=RKIND), public, parameter :: CP_OCEAN = 3989.24495292815_r8_kind !< Specific heat capacity taken from McDougall (2002)
+real(kind=RKIND), public, parameter :: CP_OCEAN = 3989.24495292815_r8_kind !< Specific heat capacity taken
+ !! from McDougall (2002)
!! "Potential Enthalpy ..." [J/kg/deg]
real(kind=RKIND), public, parameter :: KAPPA = RDGAS/CP_AIR !< RDGAS / CP_AIR [dimensionless]
real(kind=RKIND), public, parameter :: DENS_H2O = 1000._r8_kind !< Density of liquid water [kg/m^3]
real(kind=RKIND), public, parameter :: RHOAIR = 1.292269_r8_kind !< Reference atmospheric density [kg/m^3]
real(kind=RKIND), public, parameter :: RHO0 = 1.035E3_r8_kind !< Average density of sea water [kg/m^3]
-real(kind=RKIND), public, parameter :: RHO0R = 1.0_r8_kind/RHO0 !< Reciprocal of average density of sea water [m^3/kg]
+real(kind=RKIND), public, parameter :: RHO0R = 1.0_r8_kind/RHO0 !< Reciprocal of average density of
+ !! sea water [m^3/kg]
real(kind=RKIND), public, parameter :: RHO_CP = RHO0*CP_OCEAN !< (kg/m^3)*(cal/kg/deg C)(joules/cal) =
!! (joules/m^3/deg C) [J/m^3/deg]
-real(kind=RKIND), public, parameter :: O2MIXRAT = 2.0953E-01_r8_kind !< Mixing ratio of molecular oxygen in air [dimensionless]
+real(kind=RKIND), public, parameter :: O2MIXRAT = 2.0953E-01_r8_kind !< Mixing ratio of molecular oxygen
+ !! in air [dimensionless]
real(kind=RKIND), public, parameter :: WTMAIR = 2.896440E+01_r8_kind !< Molecular weight of air [AMU]
real(kind=RKIND), public, parameter :: WTMH2O = WTMAIR*(RDGAS/RVGAS) !< Molecular weight of water [AMU]
real(kind=RKIND), public, parameter :: WTMOZONE = 47.99820_r8_kind !< Molecular weight of ozone [AMU]
real(kind=RKIND), public, parameter :: WTMC = 12.00000_r8_kind !< Molecular weight of carbon [AMU]
-real(kind=RKIND), public, parameter :: WTMCO2 = 44.00995_r8_kind !< Molecular weight of carbon dioxide [AMU]
+real(kind=RKIND), public, parameter :: WTMCO2 = 44.00995_r8_kind !< Molecular weight of carbon dioxide
+ !! [AMU]
real(kind=RKIND), public, parameter :: WTMCH4 = 16.0425_r8_kind !< Molecular weight of methane [AMU]
-real(kind=RKIND), public, parameter :: WTMO2 = 31.9988_r8_kind !< Molecular weight of molecular oxygen [AMU]
-real(kind=RKIND), public, parameter :: WTMCFC11 = 137.3681_r8_kind !< Molecular weight of CFC-11 (CCl3F) [AMU]
-real(kind=RKIND), public, parameter :: WTMCFC12 = 120.9135_r8_kind !< Molecular weight of CFC-21 (CCl2F2) [AMU]
+real(kind=RKIND), public, parameter :: WTMO2 = 31.9988_r8_kind !< Molecular weight of molecular
+ !! oxygen [AMU]
+real(kind=RKIND), public, parameter :: WTMCFC11 = 137.3681_r8_kind !< Molecular weight of CFC-11
+ !! (CCl3F) [AMU]
+real(kind=RKIND), public, parameter :: WTMCFC12 = 120.9135_r8_kind !< Molecular weight of CFC-21
+ !! (CCl2F2) [AMU]
real(kind=RKIND), public, parameter :: WTMN = 14.0067_r8_kind !< Molecular weight of Nitrogen [AMU]
real(kind=RKIND), public, parameter :: DIFFAC = 1.660_r8_kind !< Diffusivity factor [dimensionless]
real(kind=RKIND), public, parameter :: ES0 = 1.0_r8_kind !< Humidity factor [dimensionless]
- !! Controls the humidity content of the atmosphere through
- !! the Saturation Vapour Pressure expression
- !! when using DO_SIMPLE
+ !! Controls the humidity content of
+ !! the atmosphere through
+ !! the Saturation Vapour Pressure
+ !! expression when using DO_SIMPLE
real(kind=RKIND), public, parameter :: CON_CLIQ = 4.1855E+3_r8_kind !< Specific heat H2O liq [J/kg/K]
real(kind=RKIND), public, parameter :: CON_CSOL = 2.1060E+3_r8_kind !< Specific heat H2O ice [J/kg/K]
@@ -90,11 +103,15 @@ real(kind=RKIND), public, parameter :: AVOGNO = 6.023000E+23_r8_kind !< A
real(kind=RKIND), public, parameter :: VONKARM = 0.40_r8_kind !< Von Karman constant [dimensionless]
!--- Miscellaneous constants
-real(kind=RKIND), public, parameter :: ALOGMIN = -50.0_r8_kind !< Minimum value allowed as argument to log function [N/A]
-real(kind=RKIND), public, parameter :: EPSLN = 1.0E-40_r8_kind !< A small number to prevent divide by zero exceptions [N/A]
-real(kind=RKIND), public, parameter :: RADCON = ((1.0D+02*GRAV)/(1.0D+04*CP_AIR))*SECONDS_PER_DAY !< Factor to convert flux divergence
- !! to heating rate in degrees per day
- !! [deg sec/(cm day)]
-real(kind=RKIND), public, parameter :: RADCON_MKS = (GRAV/CP_AIR)*SECONDS_PER_DAY !< Factor to convert flux divergence
- !! to heating rate in degrees per day
- !! [deg sec/(m day)]
+real(kind=RKIND), public, parameter :: ALOGMIN = -50.0_r8_kind !< Minimum value allowed as argument
+ !! to log function [N/A]
+real(kind=RKIND), public, parameter :: EPSLN = 1.0E-40_r8_kind !< A small number to prevent divide by
+ !! zero exceptions [N/A]
+real(kind=RKIND), public, parameter :: RADCON = ((1.0D+02*GRAV)/(1.0D+04*CP_AIR))*SECONDS_PER_DAY !< Factor to
+ !! convert flux divergence
+ !! to heating rate in degrees per day
+ !! [deg sec/(cm day)]
+real(kind=RKIND), public, parameter :: RADCON_MKS = (GRAV/CP_AIR)*SECONDS_PER_DAY !< Factor to
+ !! convert flux divergence
+ !! to heating rate in degrees per day
+ !! [deg sec/(m day)]
diff --git a/constants4/Makefile.am b/constants4/Makefile.am
index 73a4726675..9bba1fddf1 100644
--- a/constants4/Makefile.am
+++ b/constants4/Makefile.am
@@ -32,12 +32,12 @@ noinst_LTLIBRARIES = libconstants4.la
# The convenience library depends on its source.
libconstants4_la_SOURCES = \
fmsconstantsr4.F90 \
- gfdl_constantsR4.h \
- gfs_constantsR4.h \
- geos_constantsR4.h \
+ gfdl_constantsR4.fh \
+ gfs_constantsR4.fh \
+ geos_constantsR4.fh \
constantsr4.F90
-fmsconstantsr4.$(FC_MODEXT): gfdl_constantsR4.h gfs_constantsR4.h geos_constantsR4.h
+fmsconstantsr4.$(FC_MODEXT): gfdl_constantsR4.fh gfs_constantsR4.fh geos_constantsR4.fh
constantsr4_mod.$(FC_MODEXT): fmsconstantsr4.$(FC_MODEXT)
# Mod files are built and then installed as headers
diff --git a/constants4/fmsconstantsr4.F90 b/constants4/fmsconstantsr4.F90
index 76267101d6..fd80cc062e 100644
--- a/constants4/fmsconstantsr4.F90
+++ b/constants4/fmsconstantsr4.F90
@@ -28,9 +28,9 @@
!! statement.
!!
!! The currently support contant systems are:
-!! GFDL constants (gfdl_constantsR4.h)
-!! GEOS constants (geos_constantsR4.h)
-!! GFS constants (gfs_constantsR4.h)
+!! GFDL constants (gfdl_constantsR4.fh)
+!! GEOS constants (geos_constantsR4.fh)
+!! GFS constants (gfs_constantsR4.fh)
!!
!!
!! The name given to a particular constant may be changed.
@@ -68,13 +68,13 @@ module FMSconstantsR4
!--- perform error checking and include the correct system of constants
#if defined(GFDL_CONSTANTS) && !defined(GFS_CONSTANTS) && !defined(GEOS_CONSTANTS)
#warning "Using GFDL constantsR4"
-#include
+#include
#elif !defined(GFDL_CONSTANTS) && defined(GFS_CONSTANTS) && !defined(GEOS_CONSTANTS)
#warning "Using GFS constantsR4"
-#include
+#include
#elif !defined(GFDL_CONSTANTS) && !defined(GFS_CONSTANTS) && defined(GEOS_CONSTANTS)
#warning "Using GEOS constantsR4"
-#include
+#include
#else
#error FATAL FMSConstantsR4 error - multiple constants macros are defined for FMS
#endif
diff --git a/constants4/geos_constantsR4.h b/constants4/geos_constantsR4.fh
similarity index 73%
rename from constants4/geos_constantsR4.h
rename to constants4/geos_constantsR4.fh
index a2b719b29e..d2840e5a6b 100644
--- a/constants4/geos_constantsR4.h
+++ b/constants4/geos_constantsR4.fh
@@ -23,11 +23,14 @@ character(len=20), public, parameter :: constantsR4_version = 'FMSConstantsR4: G
real(kind=RKIND), public, parameter :: small_fac = 1._r8_kind
!--- Spherical coordinate conversion constants
-real(kind=r8_kind), public, parameter :: PI_8 = 3.14159265358979323846_r8_kind !< Ratio of circle circumference to diameter [N/A]
-real(kind=RKIND), public, parameter :: PI = PI_8 !< Ratio of circle circumference to diameter [N/A]
+real(kind=r8_kind), public, parameter :: PI_8 = 3.14159265358979323846_r8_kind !< Ratio of circle circumference
+ !! to diameter [N/A]
+real(kind=RKIND), public, parameter :: PI = PI_8 !< Ratio of circle circumference
+ !! to diameter [N/A]
real(kind=RKIND), public, parameter :: RAD_TO_DEG = 180._r8_kind/PI_8 !< Degrees per radian [deg/rad]
real(kind=RKIND), public, parameter :: DEG_TO_RAD = PI_8/180._r8_kind !< Radians per degree [rad/deg]
-real(kind=RKIND), public, parameter :: RADIAN = RAD_TO_DEG !< Equal to RAD_TO_DEG for backward compatability. [rad/deg]
+real(kind=RKIND), public, parameter :: RADIAN = RAD_TO_DEG !< Equal to RAD_TO_DEG for backward
+ !! compatability.[rad/deg]
!--- Earth physical constants
real(kind=RKIND), public, parameter :: RADIUS = 6371.0E3_r8_kind !< Radius of the Earth [m]
@@ -39,7 +42,8 @@ real(kind=RKIND), public, parameter :: SECONDS_PER_MINUTE = 60._r8_kind !
!--- Various gas constants
real(kind=RKIND), public, parameter :: RDGAS = 8314.47 /28.965 !< Gas constant for dry air [J/kg/deg]
-real(kind=RKIND), public, parameter :: RVGAS = 8314.47 /18.015 !< Gas constant for water vapor [J/kg/deg]
+real(kind=RKIND), public, parameter :: RVGAS = 8314.47 /18.015 !< Gas constant for water vapor
+ !! [J/kg/deg]
real(kind=RKIND), public, parameter :: HLV = 2.4665E6_r8_kind !< Latent heat of evaporation [J/kg]
real(kind=RKIND), public, parameter :: HLF = 3.3370E5_r8_kind !< Latent heat of fusion [J/kg]
real(kind=RKIND), public, parameter :: HLS = HLV + HLF !< Latent heat of sublimation [J/kg]
@@ -48,30 +52,38 @@ real(kind=RKIND), public, parameter :: CP_AIR = RDGAS/KAPPA !< S
!! at constant pressure [J/kg/deg]
real(kind=RKIND), public, parameter :: CP_VAPOR = 4.0_r8_kind*RVGAS !< Specific heat capacity of water vapor
!! at constant pressure [J/kg/deg]
-real(kind=RKIND), public, parameter :: CP_OCEAN = 3989.24495292815_r8_kind !< Specific heat capacity taken from McDougall (2002)
+real(kind=RKIND), public, parameter :: CP_OCEAN = 3989.24495292815_r8_kind !< Specific heat capacity taken
+ !! from McDougall (2002)
!! "Potential Enthalpy ..." [J/kg/deg]
real(kind=RKIND), public, parameter :: DENS_H2O = 1000._r8_kind !< Density of liquid water [kg/m^3]
real(kind=RKIND), public, parameter :: RHOAIR = 1.292269_r8_kind !< Reference atmospheric density [kg/m^3]
real(kind=RKIND), public, parameter :: RHO0 = 1.035E3_r8_kind !< Average density of sea water [kg/m^3]
-real(kind=RKIND), public, parameter :: RHO0R = 1.0_r8_kind/RHO0 !< Reciprocal of average density of sea water [m^3/kg]
+real(kind=RKIND), public, parameter :: RHO0R = 1.0_r8_kind/RHO0 !< Reciprocal of average density of sea
+ !! water [m^3/kg]
real(kind=RKIND), public, parameter :: RHO_CP = RHO0*CP_OCEAN !< (kg/m^3)*(cal/kg/deg C)(joules/cal) =
!! (joules/m^3/deg C) [J/m^3/deg]
-real(kind=RKIND), public, parameter :: O2MIXRAT = 2.0953E-01_r8_kind !< Mixing ratio of molecular oxygen in air [dimensionless]
+real(kind=RKIND), public, parameter :: O2MIXRAT = 2.0953E-01_r8_kind !< Mixing ratio of molecular oxygen
+ !! in air [dimensionless]
real(kind=RKIND), public, parameter :: WTMAIR = 2.896440E+01_r8_kind !< Molecular weight of air [AMU]
real(kind=RKIND), public, parameter :: WTMH2O = WTMAIR*(RDGAS/RVGAS) !< Molecular weight of water [AMU]
real(kind=RKIND), public, parameter :: WTMOZONE = 47.99820_r8_kind !< Molecular weight of ozone [AMU]
real(kind=RKIND), public, parameter :: WTMC = 12.00000_r8_kind !< Molecular weight of carbon [AMU]
-real(kind=RKIND), public, parameter :: WTMCO2 = 44.00995_r8_kind !< Molecular weight of carbon dioxide [AMU]
+real(kind=RKIND), public, parameter :: WTMCO2 = 44.00995_r8_kind !< Molecular weight of carbon dioxide
+ !! [AMU]
real(kind=RKIND), public, parameter :: WTMCH4 = 16.0425_r8_kind !< Molecular weight of methane [AMU]
-real(kind=RKIND), public, parameter :: WTMO2 = 31.9988_r8_kind !< Molecular weight of molecular oxygen [AMU]
-real(kind=RKIND), public, parameter :: WTMCFC11 = 137.3681_r8_kind !< Molecular weight of CFC-11 (CCl3F) [AMU]
-real(kind=RKIND), public, parameter :: WTMCFC12 = 120.9135_r8_kind !< Molecular weight of CFC-21 (CCl2F2) [AMU]
+real(kind=RKIND), public, parameter :: WTMO2 = 31.9988_r8_kind !< Molecular weight of molecular oxygen
+ !! [AMU]
+real(kind=RKIND), public, parameter :: WTMCFC11 = 137.3681_r8_kind !< Molecular weight of CFC-11 (CCl3F)
+ !! [AMU]
+real(kind=RKIND), public, parameter :: WTMCFC12 = 120.9135_r8_kind !< Molecular weight of CFC-21 (CCl2F2)
+ !! [AMU]
real(kind=RKIND), public, parameter :: WTMN = 14.0067_r8_kind !< Molecular weight of Nitrogen [AMU]
real(kind=RKIND), public, parameter :: DIFFAC = 1.660_r8_kind !< Diffusivity factor [dimensionless]
real(kind=RKIND), public, parameter :: ES0 = 1.0_r8_kind !< Humidity factor [dimensionless]
- !! Controls the humidity content of the atmosphere through
- !! the Saturation Vapour Pressure expression
- !! when using DO_SIMPLE
+ !! Controls the humidity content of
+ !! the atmosphere through
+ !! the Saturation Vapour Pressure
+ !! expression when using DO_SIMPLE
!--- Pressure and Temperature constants
real(kind=RKIND), public, parameter :: PSTD = 1.013250E+06_r8_kind !< Mean sea level pressure [dynes/cm^2]
@@ -87,11 +99,16 @@ real(kind=RKIND), public, parameter :: AVOGNO = 6.023000E+23_r8_kind !< A
real(kind=RKIND), public, parameter :: VONKARM = 0.40_r8_kind !< Von Karman constant [dimensionless]
!--- Miscellaneous constants
-real(kind=RKIND), public, parameter :: ALOGMIN = -50.0_r8_kind !< Minimum value allowed as argument to log function [N/A]
-real(kind=RKIND), public, parameter :: EPSLN = 1.0E-40_r8_kind !< A small number to prevent divide by zero exceptions [N/A]
-real(kind=RKIND), public, parameter :: RADCON = ((1.0E+02*GRAV)/(1.0D+04*CP_AIR))*SECONDS_PER_DAY !< Factor to convert flux divergence
- !! to heating rate in degrees per day
- !! [deg sec/(cm day)]
-real(kind=RKIND), public, parameter :: RADCON_MKS = (GRAV/CP_AIR)*SECONDS_PER_DAY !< Factor to convert flux divergence
- !! to heating rate in degrees per day
- !! [deg sec/(m day)]
+real(kind=RKIND), public, parameter :: ALOGMIN = -50.0_r8_kind !< Minimum value allowed as argument to
+ !! log function [N/A]
+real(kind=RKIND), public, parameter :: EPSLN = 1.0E-40_r8_kind !< A small number to prevent divide by zero
+ !! exceptions [N/A]
+real(kind=RKIND), public, parameter :: RADCON = ((1.0E+02*GRAV)/(1.0D+04*CP_AIR))*SECONDS_PER_DAY !< Factor to
+ !! convert flux divergence
+ !! to heating rate in degrees per day
+ !! [deg sec/(cm day)]
+
+real(kind=RKIND), public, parameter :: RADCON_MKS = (GRAV/CP_AIR)*SECONDS_PER_DAY !< Factor to
+ !! convert flux divergence
+ !! to heating rate in degrees per day
+ !! [deg sec/(m day)]
diff --git a/constants4/gfdl_constantsR4.h b/constants4/gfdl_constantsR4.fh
similarity index 74%
rename from constants4/gfdl_constantsR4.h
rename to constants4/gfdl_constantsR4.fh
index e0bd9573ba..935445aa56 100644
--- a/constants4/gfdl_constantsR4.h
+++ b/constants4/gfdl_constantsR4.fh
@@ -23,11 +23,14 @@ character(len=20), public, parameter :: constantsR4_version = 'FMSConstantsR4: G
real(kind=RKIND), public, parameter :: small_fac = 1._r8_kind
!--- Spherical coordinate conversion constants
-real(kind=r8_kind), public, parameter :: PI_8 = 3.14159265358979323846_r8_kind !< Ratio of circle circumference to diameter [N/A]
-real(kind=RKIND), public, parameter :: PI = PI_8 !< Ratio of circle circumference to diameter [N/A]
+real(kind=r8_kind), public, parameter :: PI_8 = 3.14159265358979323846_r8_kind !< Ratio of circle circumference
+ !! to diameter [N/A]
+real(kind=RKIND), public, parameter :: PI = PI_8 !< Ratio of circle circumferen
+ !! to diameter [N/A]
real(kind=RKIND), public, parameter :: RAD_TO_DEG = 180._r8_kind/PI_8 !< Degrees per radian [deg/rad]
real(kind=RKIND), public, parameter :: DEG_TO_RAD = PI_8/180._r8_kind !< Radians per degree [rad/deg]
-real(kind=RKIND), public, parameter :: RADIAN = RAD_TO_DEG !< Equal to RAD_TO_DEG for backward compatability. [rad/deg]
+real(kind=RKIND), public, parameter :: RADIAN = RAD_TO_DEG !< Equal to RAD_TO_DEG for backward
+ !! compatability. [rad/deg]
!--- Earth physical constants
real(kind=RKIND), public, parameter :: RADIUS = 6371.0E+3_r8_kind !< Radius of the Earth [m]
@@ -39,7 +42,8 @@ real(kind=RKIND), public, parameter :: SECONDS_PER_MINUTE = 60._r8_kind !
!--- Various gas constants
real(kind=RKIND), public, parameter :: RDGAS = 287.04_r8_kind !< Gas constant for dry air [J/kg/deg]
-real(kind=RKIND), public, parameter :: RVGAS = 461.50_r8_kind !< Gas constant for water vapor [J/kg/deg]
+real(kind=RKIND), public, parameter :: RVGAS = 461.50_r8_kind !< Gas constant for water vapor
+ !! [J/kg/deg]
real(kind=RKIND), public, parameter :: HLV = 2.500E6_r8_kind !< Latent heat of evaporation [J/kg]
real(kind=RKIND), public, parameter :: HLF = 3.34E5_r8_kind !< Latent heat of fusion [J/kg]
real(kind=RKIND), public, parameter :: HLS = HLV + HLF !< Latent heat of sublimation [J/kg]
@@ -48,29 +52,37 @@ real(kind=RKIND), public, parameter :: CP_AIR = RDGAS/KAPPA !< S
!! at constant pressure [J/kg/deg]
real(kind=RKIND), public, parameter :: CP_VAPOR = 4.0_r8_kind*RVGAS !< Specific heat capacity of water vapor
!! at constant pressure [J/kg/deg]
-real(kind=RKIND), public, parameter :: CP_OCEAN = 3989.24495292815_r8_kind !< Specific heat capacity taken from McDougall (2002)
+real(kind=RKIND), public, parameter :: CP_OCEAN = 3989.24495292815_r8_kind !< Specific heat capacity taken
+ !! from McDougall (2002)
!! "Potential Enthalpy ..." [J/kg/deg]
real(kind=RKIND), public, parameter :: DENS_H2O = 1000._r8_kind !< Density of liquid water [kg/m^3]
real(kind=RKIND), public, parameter :: RHOAIR = 1.292269_r8_kind !< Reference atmospheric density [kg/m^3]
real(kind=RKIND), public, parameter :: RHO0 = 1.035E3_r8_kind !< Average density of sea water [kg/m^3]
-real(kind=RKIND), public, parameter :: RHO0R = 1.0_r8_kind/RHO0 !< Reciprocal of average density of sea water [m^3/kg]
+real(kind=RKIND), public, parameter :: RHO0R = 1.0_r8_kind/RHO0 !< Reciprocal of average density of
+ !! sea water [m^3/kg]
real(kind=RKIND), public, parameter :: RHO_CP = RHO0*CP_OCEAN !< (kg/m^3)*(cal/kg/deg C)(joules/cal) =
!! (joules/m^3/deg C) [J/m^3/deg]
-real(kind=RKIND), public, parameter :: O2MIXRAT = 2.0953E-01_r8_kind !< Mixing ratio of molecular oxygen in air [dimensionless]
+real(kind=RKIND), public, parameter :: O2MIXRAT = 2.0953E-01_r8_kind !< Mixing ratio of molecular oxygen
+ !! in air [dimensionless]
real(kind=RKIND), public, parameter :: WTMAIR = 2.896440E+01_r8_kind !< Molecular weight of air [AMU]
real(kind=RKIND), public, parameter :: WTMH2O = WTMAIR*(RDGAS/RVGAS) !< Molecular weight of water [AMU]
real(kind=RKIND), public, parameter :: WTMOZONE = 47.99820_r8_kind !< Molecular weight of ozone [AMU]
real(kind=RKIND), public, parameter :: WTMC = 12.00000_r8_kind !< Molecular weight of carbon [AMU]
-real(kind=RKIND), public, parameter :: WTMCO2 = 44.00995_r8_kind !< Molecular weight of carbon dioxide [AMU]
+real(kind=RKIND), public, parameter :: WTMCO2 = 44.00995_r8_kind !< Molecular weight of carbon dioxide
+ !! [AMU]
real(kind=RKIND), public, parameter :: WTMCH4 = 16.0425_r8_kind !< Molecular weight of methane [AMU]
-real(kind=RKIND), public, parameter :: WTMO2 = 31.9988_r8_kind !< Molecular weight of molecular oxygen [AMU]
-real(kind=RKIND), public, parameter :: WTMCFC11 = 137.3681_r8_kind !< Molecular weight of CFC-11 (CCl3F) [AMU]
-real(kind=RKIND), public, parameter :: WTMCFC12 = 120.9135_r8_kind !< Molecular weight of CFC-21 (CCl2F2) [AMU]
+real(kind=RKIND), public, parameter :: WTMO2 = 31.9988_r8_kind !< Molecular weight of molecular oxygen
+ !! [AMU]
+real(kind=RKIND), public, parameter :: WTMCFC11 = 137.3681_r8_kind !< Molecular weight of CFC-11
+ !! (CCl3F) [AMU]
+real(kind=RKIND), public, parameter :: WTMCFC12 = 120.9135_r8_kind !< Molecular weight of CFC-21
+ !! (CCl2F2) [AMU]
real(kind=RKIND), public, parameter :: WTMN = 14.0067_r8_kind !< Molecular weight of Nitrogen [AMU]
real(kind=RKIND), public, parameter :: DIFFAC = 1.660_r8_kind !< Diffusivity factor [dimensionless]
real(kind=RKIND), public, parameter :: ES0 = 1.0_r8_kind !< Humidity factor [dimensionless]
- !! Controls the humidity content of the atmosphere through
- !! the Saturation Vapour Pressure expression
+ !! Controls the humidity content of
+ !! the atmosphere through the
+ !! Saturation Vapour Pressure expression
!! when using DO_SIMPLE
!--- Pressure and Temperature constants
@@ -87,11 +99,15 @@ real(kind=RKIND), public, parameter :: AVOGNO = 6.023000E+23_r8_kind !< A
real(kind=RKIND), public, parameter :: VONKARM = 0.40_r8_kind !< Von Karman constant [dimensionless]
!--- Miscellaneous constants
-real(kind=RKIND), public, parameter :: ALOGMIN = -50.0_r8_kind !< Minimum value allowed as argument to log function [N/A]
-real(kind=RKIND), public, parameter :: EPSLN = 1.0E-40_r8_kind !< A small number to prevent divide by zero exceptions [N/A]
-real(kind=RKIND), public, parameter :: RADCON = ((1.0E+02*GRAV)/(1.0D+04*CP_AIR))*SECONDS_PER_DAY !< Factor to convert flux divergence
- !! to heating rate in degrees per day
- !! [deg sec/(cm day)]
-real(kind=RKIND), public, parameter :: RADCON_MKS = (GRAV/CP_AIR)*SECONDS_PER_DAY !< Factor to convert flux divergence
- !! to heating rate in degrees per day
- !! [deg sec/(m day)]
+real(kind=RKIND), public, parameter :: ALOGMIN = -50.0_r8_kind !< Minimum value allowed as argument
+ !! to log function [N/A]
+real(kind=RKIND), public, parameter :: EPSLN = 1.0E-40_r8_kind !< A small number to prevent divide by zero
+ !! exceptions [N/A]
+real(kind=RKIND), public, parameter :: RADCON = ((1.0E+02*GRAV)/(1.0D+04*CP_AIR))*SECONDS_PER_DAY !< Factor to
+ !! convert flux divergence
+ !! to heating rate in degrees per day
+ !! [deg sec/(cm day)]
+real(kind=RKIND), public, parameter :: RADCON_MKS = (GRAV/CP_AIR)*SECONDS_PER_DAY !< Factor to
+ !! convert flux divergence
+ !! to heating rate in degrees per day
+ !! [deg sec/(m day)]
diff --git a/constants4/gfs_constantsR4.h b/constants4/gfs_constantsR4.fh
similarity index 73%
rename from constants4/gfs_constantsR4.h
rename to constants4/gfs_constantsR4.fh
index 04da7b0014..0a7178b8e6 100644
--- a/constants4/gfs_constantsR4.h
+++ b/constants4/gfs_constantsR4.fh
@@ -23,17 +23,22 @@ character(len=20), public, parameter :: constantsR4_version = 'FMSConstantsR4: G
real(kind=RKIND), public, parameter :: small_fac = 1._r8_kind
!--- Spherical coordinate conversion constants
-real(kind=r8_kind), public, parameter :: PI_8 = 3.1415926535897931_r8_kind !< Ratio of circle circumference to diameter [N/A]
-real(kind=RKIND), public, parameter :: PI = PI_8 !< Ratio of circle circumference to diameter [N/A]
+real(kind=r8_kind), public, parameter :: PI_8 = 3.1415926535897931_r8_kind !< Ratio of circle circumference
+ !! to diameter [N/A]
+real(kind=RKIND), public, parameter :: PI = PI_8 !< Ratio of circle circumference
+ !! to diameter [N/A]
real(kind=RKIND), public, parameter :: RAD_TO_DEG = 180._r8_kind/PI_8 !< Degrees per radian [deg/rad]
real(kind=RKIND), public, parameter :: DEG_TO_RAD = PI_8/180._r8_kind !< Radians per degree [rad/deg]
-real(kind=RKIND), public, parameter :: RADIAN = RAD_TO_DEG !< Equal to RAD_TO_DEG for backward compatability. [rad/deg]
+real(kind=RKIND), public, parameter :: RADIAN = RAD_TO_DEG !< Equal to RAD_TO_DEG for
+ !! backward compatability. [rad/deg]
!--- Earth physical constants
real(kind=RKIND), public, parameter :: RADIUS = 6.3712E+6_r8_kind !< Radius of the Earth [m]
real(kind=RKIND), public, parameter :: OMEGA = 7.2921E-5_r8_kind !< Rotation rate of the Earth [1/s]
-real(kind=r8_kind), public, parameter :: GRAV_8 = 9.80665_r8_kind !< Acceleration due to gravity [m/s^2] (REAL(KIND=8))
-real(kind=RKIND), public, parameter :: GRAV = GRAV_8 !< Acceleration due to gravity [m/s^2]
+real(kind=r8_kind), public, parameter :: GRAV_8 = 9.80665_r8_kind !< Acceleration due to gravity
+ !! [m/s^2] (REAL(KIND=8))
+real(kind=RKIND), public, parameter :: GRAV = GRAV_8 !< Acceleration due to gravity
+ !! [m/s^2]
real(kind=RKIND), public, parameter :: SECONDS_PER_DAY = 86400._r8_kind !< Seconds in a day [s]
real(kind=RKIND), public, parameter :: SECONDS_PER_HOUR = 3600._r8_kind !< Seconds in an hour [s]
real(kind=RKIND), public, parameter :: SECONDS_PER_MINUTE = 60._r8_kind !< Seconds in a minute [s]
@@ -48,31 +53,39 @@ real(kind=RKIND), public, parameter :: CP_AIR = 1004.6_r8_kind !< S
!! at constant pressure [J/kg/deg]
real(kind=RKIND), public, parameter :: CP_VAPOR = 4.0_r8_kind*RVGAS !< Specific heat capacity of water vapor
!! at constant pressure [J/kg/deg]
-real(kind=RKIND), public, parameter :: CP_OCEAN = 3989.24495292815_r8_kind !< Specific heat capacity taken from McDougall (2002)
+real(kind=RKIND), public, parameter :: CP_OCEAN = 3989.24495292815_r8_kind !< Specific heat capacity taken
+ !! from McDougall (2002)
!! "Potential Enthalpy ..." [J/kg/deg]
real(kind=RKIND), public, parameter :: KAPPA = RDGAS/CP_AIR !< RDGAS / CP_AIR [dimensionless]
real(kind=RKIND), public, parameter :: DENS_H2O = 1000._r8_kind !< Density of liquid water [kg/m^3]
real(kind=RKIND), public, parameter :: RHOAIR = 1.292269_r8_kind !< Reference atmospheric density [kg/m^3]
real(kind=RKIND), public, parameter :: RHO0 = 1.035E3_r8_kind !< Average density of sea water [kg/m^3]
-real(kind=RKIND), public, parameter :: RHO0R = 1.0_r8_kind/RHO0 !< Reciprocal of average density of sea water [m^3/kg]
+real(kind=RKIND), public, parameter :: RHO0R = 1.0_r8_kind/RHO0 !< Reciprocal of average density
+ !! of sea water [m^3/kg]
real(kind=RKIND), public, parameter :: RHO_CP = RHO0*CP_OCEAN !< (kg/m^3)*(cal/kg/deg C)(joules/cal) =
!! (joules/m^3/deg C) [J/m^3/deg]
-real(kind=RKIND), public, parameter :: O2MIXRAT = 2.0953E-01_r8_kind !< Mixing ratio of molecular oxygen in air [dimensionless]
+real(kind=RKIND), public, parameter :: O2MIXRAT = 2.0953E-01_r8_kind !< Mixing ratio of molecular oxygen
+ !! in air [dimensionless]
real(kind=RKIND), public, parameter :: WTMAIR = 2.896440E+01_r8_kind !< Molecular weight of air [AMU]
real(kind=RKIND), public, parameter :: WTMH2O = WTMAIR*(RDGAS/RVGAS) !< Molecular weight of water [AMU]
real(kind=RKIND), public, parameter :: WTMOZONE = 47.99820_r8_kind !< Molecular weight of ozone [AMU]
real(kind=RKIND), public, parameter :: WTMC = 12.00000_r8_kind !< Molecular weight of carbon [AMU]
-real(kind=RKIND), public, parameter :: WTMCO2 = 44.00995_r8_kind !< Molecular weight of carbon dioxide [AMU]
+real(kind=RKIND), public, parameter :: WTMCO2 = 44.00995_r8_kind !< Molecular weight of carbon dioxide
+ !! [AMU]
real(kind=RKIND), public, parameter :: WTMCH4 = 16.0425_r8_kind !< Molecular weight of methane [AMU]
-real(kind=RKIND), public, parameter :: WTMO2 = 31.9988_r8_kind !< Molecular weight of molecular oxygen [AMU]
-real(kind=RKIND), public, parameter :: WTMCFC11 = 137.3681_r8_kind !< Molecular weight of CFC-11 (CCl3F) [AMU]
-real(kind=RKIND), public, parameter :: WTMCFC12 = 120.9135_r8_kind !< Molecular weight of CFC-21 (CCl2F2) [AMU]
+real(kind=RKIND), public, parameter :: WTMO2 = 31.9988_r8_kind !< Molecular weight of molecular oxygen
+ !! [AMU]
+real(kind=RKIND), public, parameter :: WTMCFC11 = 137.3681_r8_kind !< Molecular weight of CFC-11
+ !! (CCl3F) [AMU]
+real(kind=RKIND), public, parameter :: WTMCFC12 = 120.9135_r8_kind !< Molecular weight of CFC-21
+ !! (CCl2F2) [AMU]
real(kind=RKIND), public, parameter :: WTMN = 14.0067_r8_kind !< Molecular weight of Nitrogen [AMU]
real(kind=RKIND), public, parameter :: DIFFAC = 1.660_r8_kind !< Diffusivity factor [dimensionless]
real(kind=RKIND), public, parameter :: ES0 = 1.0_r8_kind !< Humidity factor [dimensionless]
- !! Controls the humidity content of the atmosphere through
- !! the Saturation Vapour Pressure expression
- !! when using DO_SIMPLE
+ !! Controls the humidity content of the
+ !! atmosphere through
+ !! the Saturation Vapour Pressure
+ !! expression when using DO_SIMPLE
real(kind=RKIND), public, parameter :: CON_CLIQ = 4.1855E+3_r8_kind !< Specific heat H2O liq [J/kg/K]
real(kind=RKIND), public, parameter :: CON_CSOL = 2.1060E+3_r8_kind !< Specific heat H2O ice [J/kg/K]
@@ -90,11 +103,16 @@ real(kind=RKIND), public, parameter :: AVOGNO = 6.023000E+23_r8_kind !< A
real(kind=RKIND), public, parameter :: VONKARM = 0.40_r8_kind !< Von Karman constant [dimensionless]
!--- Miscellaneous constants
-real(kind=RKIND), public, parameter :: ALOGMIN = -50.0_r8_kind !< Minimum value allowed as argument to log function [N/A]
-real(kind=RKIND), public, parameter :: EPSLN = 1.0E-40_r8_kind !< A small number to prevent divide by zero exceptions [N/A]
-real(kind=RKIND), public, parameter :: RADCON = ((1.0D+02*GRAV)/(1.0D+04*CP_AIR))*SECONDS_PER_DAY !< Factor to convert flux divergence
- !! to heating rate in degrees per day
- !! [deg sec/(cm day)]
-real(kind=RKIND), public, parameter :: RADCON_MKS = (GRAV/CP_AIR)*SECONDS_PER_DAY !< Factor to convert flux divergence
- !! to heating rate in degrees per day
- !! [deg sec/(m day)]
+real(kind=RKIND), public, parameter :: ALOGMIN = -50.0_r8_kind !< Minimum value allowed as argument to
+ !! log function [N/A]
+real(kind=RKIND), public, parameter :: EPSLN = 1.0E-40_r8_kind !< A small number to prevent divide by
+ !! zero exceptions [N/A]
+real(kind=RKIND), public, parameter :: RADCON = ((1.0D+02*GRAV)/(1.0D+04*CP_AIR))*SECONDS_PER_DAY !< Factor to
+ !! convert flux divergence
+ !! to heating rate in degrees per day
+ !! [deg sec/(cm day)]
+
+real(kind=RKIND), public, parameter :: RADCON_MKS = (GRAV/CP_AIR)*SECONDS_PER_DAY !< Factor to
+ !! convert flux divergence
+ !! to heating rate in degrees per day
+ !! [deg sec/(m day)]
diff --git a/drifters/Makefile.am b/drifters/Makefile.am
index 521e18907e..f365f8e1f2 100644
--- a/drifters/Makefile.am
+++ b/drifters/Makefile.am
@@ -41,9 +41,9 @@ libdrifters_la_SOURCES = drifters.F90 \
drifters_io.F90 \
cloud_interpolator.F90 \
quicksort.F90 \
- drifters_compute_k.h \
- drifters_push.h \
- drifters_set_field.h \
+ drifters_compute_k.fh \
+ drifters_push.fh \
+ drifters_set_field.fh \
fms_switches.h
drifters_mod.$(FC_MODEXT): drifters_core_mod.$(FC_MODEXT) drifters_input_mod.$(FC_MODEXT) drifters_io_mod.$(FC_MODEXT) \
diff --git a/drifters/drifters.F90 b/drifters/drifters.F90
index ab1bcdc395..eac1d6cbd8 100644
--- a/drifters/drifters.F90
+++ b/drifters/drifters.F90
@@ -89,8 +89,9 @@ module drifters_mod
drifters_io_set_position_names, drifters_io_set_position_units, &
drifters_io_set_field_names, drifters_io_set_field_units, drifters_io_write
- use drifters_comm_mod, only: drifters_comm_type,drifters_comm_new,drifters_comm_del,drifters_comm_set_pe_neighbors,&
- drifters_comm_set_domain, drifters_comm_gather, drifters_comm_update
+ use drifters_comm_mod, only: drifters_comm_type,drifters_comm_new,drifters_comm_del, &
+ drifters_comm_set_pe_neighbors, drifters_comm_set_domain, &
+ drifters_comm_gather, drifters_comm_update
use cloud_interpolator_mod, only: cld_ntrp_linear_cell_interp, cld_ntrp_locate_cell, cld_ntrp_get_cell_values
implicit none
@@ -462,14 +463,14 @@ end subroutine drifters_set_pe_neighbors
!============================================================================
#define _DIMS 2
#define drifters_push_XXX drifters_push_2
-#include "drifters_push.h"
+#include "drifters_push.fh"
#undef _DIMS
#undef drifters_push_XXX
!============================================================================
#define _DIMS 3
#define drifters_push_XXX drifters_push_3
-#include "drifters_push.h"
+#include "drifters_push.fh"
#undef _DIMS
#undef drifters_push_XXX
@@ -506,14 +507,14 @@ end subroutine drifters_modulo
!============================================================================
#define _DIMS 2
#define drifters_set_field_XXX drifters_set_field_2d
-#include "drifters_set_field.h"
+#include "drifters_set_field.fh"
#undef _DIMS
#undef drifters_set_field_XXX
!============================================================================
#define _DIMS 3
#define drifters_set_field_XXX drifters_set_field_3d
-#include "drifters_set_field.h"
+#include "drifters_set_field.fh"
#undef _DIMS
#undef drifters_set_field_XXX
!============================================================================
@@ -644,14 +645,14 @@ end subroutine drifters_write_restart
!============================================================================
#define _DIMS 2
#define drifters_compute_k_XXX drifters_computek2d
-#include "drifters_compute_k.h"
+#include "drifters_compute_k.fh"
#undef _DIMS
#undef drifters_compute_k_XXX
!============================================================================
#define _DIMS 3
#define drifters_compute_k_XXX drifters_computek3d
-#include "drifters_compute_k.h"
+#include "drifters_compute_k.fh"
#undef _DIMS
#undef drifters_compute_k_XXX
diff --git a/drifters/drifters_compute_k.h b/drifters/drifters_compute_k.fh
similarity index 100%
rename from drifters/drifters_compute_k.h
rename to drifters/drifters_compute_k.fh
diff --git a/drifters/drifters_push.h b/drifters/drifters_push.fh
similarity index 100%
rename from drifters/drifters_push.h
rename to drifters/drifters_push.fh
diff --git a/drifters/drifters_set_field.h b/drifters/drifters_set_field.fh
similarity index 93%
rename from drifters/drifters_set_field.h
rename to drifters/drifters_set_field.fh
index ab8a4b583f..9ca2a2acf5 100644
--- a/drifters/drifters_set_field.h
+++ b/drifters/drifters_set_field.fh
@@ -76,7 +76,8 @@ subroutine drifters_set_field_XXX(self, index_field, x, y, &
ij(1) = i
#ifdef _DEBUG
if(i<1) then
- print *,'drifters_set_field_XXX::OUT OF BOUND ERROR. xmin, x, xmax=', minval(x), self%core%positions(1,ip), maxval(x)
+ print *,'drifters_set_field_XXX::OUT OF BOUND ERROR. xmin, x, xmax=', &
+ & minval(x), self%core%positions(1,ip), maxval(x)
endif
#endif
ts(1) = (self%core%positions(1,ip) - x(i))/(x(i+1) - x(i))
@@ -85,7 +86,8 @@ subroutine drifters_set_field_XXX(self, index_field, x, y, &
ij(2) = j
#ifdef _DEBUG
if(j<1) then
- print *,'drifters_set_field_XXX::OUT OF BOUND ERROR. ymin, y, ymax=', minval(y), self%core%positions(2,ip), maxval(y)
+ print *,'drifters_set_field_XXX::OUT OF BOUND ERROR. ymin, y, ymax=', &
+ & minval(y), self%core%positions(2,ip), maxval(y)
endif
#endif
ts(2) = (self%core%positions(2,ip) - y(j))/(y(j+1) - y(j))
@@ -95,7 +97,8 @@ subroutine drifters_set_field_XXX(self, index_field, x, y, &
ij(3) = j
#ifdef _DEBUG
if(j<1) then
- print *,'drifters_set_field_XXX::OUT OF BOUND ERROR. ymin, y, ymax=', minval(y), self%core%positions(2,ip), maxval(y)
+ print *,'drifters_set_field_XXX::OUT OF BOUND ERROR. ymin, y, ymax=', &
+ & minval(y), self%core%positions(2,ip), maxval(y)
endif
#endif
ts(3) = (self%core%positions(3,ip) - z(j))/(z(j+1) - z(j))
diff --git a/mpp/Makefile.am b/mpp/Makefile.am
index 150127c315..c7e482ea13 100644
--- a/mpp/Makefile.am
+++ b/mpp/Makefile.am
@@ -22,7 +22,7 @@
# Ed Hartnett 2/22/19
-# Include .h and .mod files.
+# Include .fh and .mod files.
AM_CPPFLAGS = \
-I$(top_srcdir)/include \
-I$(top_srcdir)/mpp/include
@@ -43,51 +43,51 @@ libmpp_la_SOURCES = \
mpp_pset.F90 \
include/group_update_pack.inc \
include/group_update_unpack.inc \
- include/mpp_alltoall_mpi.h \
- include/mpp_alltoall_nocomm.h \
- include/mpp_chksum.h \
- include/mpp_chksum_int.h \
- include/mpp_chksum_scalar.h \
+ include/mpp_alltoall_mpi.fh \
+ include/mpp_alltoall_nocomm.fh \
+ include/mpp_chksum.fh \
+ include/mpp_chksum_int.fh \
+ include/mpp_chksum_scalar.fh \
include/mpp_comm.inc \
include/mpp_comm_mpi.inc \
include/mpp_comm_nocomm.inc \
include/mpp_data_mpi.inc \
include/mpp_data_nocomm.inc \
include/mpp_define_nest_domains.inc \
- include/mpp_do_check.h \
- include/mpp_do_checkV.h \
- include/mpp_do_get_boundary.h \
- include/mpp_do_get_boundary_ad.h \
- include/mpp_do_global_field.h \
- include/mpp_do_global_field_ad.h \
- include/mpp_do_redistribute.h \
- include/mpp_do_update.h \
- include/mpp_do_updateV.h \
- include/mpp_do_updateV_ad.h \
- include/mpp_do_updateV_nonblock.h \
- include/mpp_do_update_ad.h \
- include/mpp_do_update_nest.h \
- include/mpp_do_update_nonblock.h \
+ include/mpp_do_check.fh \
+ include/mpp_do_checkV.fh \
+ include/mpp_do_get_boundary.fh \
+ include/mpp_do_get_boundary_ad.fh \
+ include/mpp_do_global_field.fh \
+ include/mpp_do_global_field_ad.fh \
+ include/mpp_do_redistribute.fh \
+ include/mpp_do_update.fh \
+ include/mpp_do_updateV.fh \
+ include/mpp_do_updateV_ad.fh \
+ include/mpp_do_updateV_nonblock.fh \
+ include/mpp_do_update_ad.fh \
+ include/mpp_do_update_nest.fh \
+ include/mpp_do_update_nonblock.fh \
include/mpp_domains_comm.inc \
include/mpp_domains_define.inc \
include/mpp_domains_misc.inc \
include/mpp_domains_reduce.inc \
include/mpp_domains_util.inc \
- include/mpp_error_a_a.h \
- include/mpp_error_a_s.h \
- include/mpp_error_s_a.h \
- include/mpp_error_s_s.h \
- include/mpp_gather.h \
- include/mpp_get_boundary.h \
- include/mpp_get_boundary_ad.h \
- include/mpp_global_field.h \
- include/mpp_global_field_ad.h \
- include/mpp_global_field_ug.h \
- include/mpp_global_reduce.h \
- include/mpp_global_sum.h \
- include/mpp_global_sum_ad.h \
- include/mpp_global_sum_tl.h \
- include/mpp_group_update.h \
+ include/mpp_error_a_a.fh \
+ include/mpp_error_a_s.fh \
+ include/mpp_error_s_a.fh \
+ include/mpp_error_s_s.fh \
+ include/mpp_gather.fh \
+ include/mpp_get_boundary.fh \
+ include/mpp_get_boundary_ad.fh \
+ include/mpp_global_field.fh \
+ include/mpp_global_field_ad.fh \
+ include/mpp_global_field_ug.fh \
+ include/mpp_global_reduce.fh \
+ include/mpp_global_sum.fh \
+ include/mpp_global_sum_ad.fh \
+ include/mpp_global_sum_tl.fh \
+ include/mpp_group_update.fh \
include/mpp_io_connect.inc \
include/mpp_io_misc.inc \
include/mpp_io_read.inc \
@@ -95,73 +95,73 @@ libmpp_la_SOURCES = \
include/mpp_io_unstructured_write.inc \
include/mpp_io_util.inc \
include/mpp_io_write.inc \
- include/mpp_read_2Ddecomp.h \
- include/mpp_read_compressed.h \
- include/mpp_read_distributed_ascii.h \
+ include/mpp_read_2Ddecomp.fh \
+ include/mpp_read_compressed.fh \
+ include/mpp_read_distributed_ascii.fh \
include/mpp_read_distributed_ascii.inc \
- include/mpp_reduce_mpi.h \
- include/mpp_reduce_nocomm.h \
- include/mpp_scatter.h \
+ include/mpp_reduce_mpi.fh \
+ include/mpp_reduce_nocomm.fh \
+ include/mpp_scatter.fh \
include/mpp_sum.inc \
include/mpp_sum_ad.inc \
- include/mpp_sum_mpi.h \
- include/mpp_sum_mpi_ad.h \
- include/mpp_sum_nocomm.h \
- include/mpp_sum_nocomm_ad.h \
+ include/mpp_sum_mpi.fh \
+ include/mpp_sum_mpi_ad.fh \
+ include/mpp_sum_nocomm.fh \
+ include/mpp_sum_nocomm_ad.fh \
include/mpp_transmit.inc \
- include/mpp_transmit_mpi.h \
- include/mpp_transmit_nocomm.h \
- include/mpp_type_mpi.h \
- include/mpp_type_nocomm.h \
+ include/mpp_transmit_mpi.fh \
+ include/mpp_transmit_nocomm.fh \
+ include/mpp_type_mpi.fh \
+ include/mpp_type_nocomm.fh \
include/mpp_unstruct_domain.inc \
- include/mpp_unstruct_pass_data.h \
- include/mpp_update_domains2D.h \
- include/mpp_update_domains2D_ad.h \
- include/mpp_update_domains2D_nonblock.h \
- include/mpp_update_nest_domains.h \
+ include/mpp_unstruct_pass_data.fh \
+ include/mpp_update_domains2D.fh \
+ include/mpp_update_domains2D_ad.fh \
+ include/mpp_update_domains2D_nonblock.fh \
+ include/mpp_update_nest_domains.fh \
include/mpp_util.inc \
include/mpp_util_mpi.inc \
include/mpp_util_nocomm.inc \
- include/mpp_write.h \
- include/mpp_write_2Ddecomp.h \
- include/mpp_write_compressed.h \
- include/mpp_write_unlimited_axis.h \
- include/system_clock.h
+ include/mpp_write.fh \
+ include/mpp_write_2Ddecomp.fh \
+ include/mpp_write_compressed.fh \
+ include/mpp_write_unlimited_axis.fh \
+ include/system_clock.fh
mpp_mod.$(FC_MODEXT): \
mpp_parameter_mod.$(FC_MODEXT) \
mpp_data_mod.$(FC_MODEXT) \
include/mpp_util.inc \
- include/mpp_error_a_a.h \
- include/mpp_error_a_s.h \
- include/mpp_error_s_a.h \
- include/mpp_error_s_s.h \
+ include/mpp_error_a_a.fh \
+ include/mpp_error_a_s.fh \
+ include/mpp_error_s_a.fh \
+ include/mpp_error_s_s.fh \
include/mpp_util_mpi.inc \
include/mpp_util_nocomm.inc \
include/mpp_comm.inc \
- include/mpp_chksum.h \
- include/mpp_chksum_int.h \
- include/mpp_chksum_scalar.h \
+ include/mpp_chksum.fh \
+ include/mpp_chksum_int.fh \
+ include/mpp_chksum_scalar.fh \
include/mpp_comm_mpi.inc \
- include/mpp_alltoall_mpi.h \
- include/mpp_reduce_mpi.h \
- include/mpp_sum_mpi.h \
+ include/mpp_alltoall_mpi.fh \
+ include/mpp_reduce_mpi.fh \
+ include/mpp_sum_mpi.fh \
include/mpp_sum.inc \
- include/mpp_sum_mpi_ad.h \
+ include/mpp_sum_mpi_ad.fh \
include/mpp_sum_ad.inc \
- include/mpp_transmit_mpi.h \
+ include/mpp_transmit_mpi.fh \
include/mpp_transmit.inc \
- include/mpp_type_mpi.h \
+ include/mpp_type_mpi.fh \
include/mpp_comm_nocomm.inc \
- include/mpp_alltoall_nocomm.h \
- include/mpp_reduce_nocomm.h \
- include/mpp_sum_nocomm.h \
- include/mpp_sum_nocomm_ad.h \
- include/mpp_transmit_nocomm.h \
- include/mpp_type_nocomm.h \
- include/mpp_gather.h \
- include/mpp_scatter.h \
- include/system_clock.h
+ include/mpp_alltoall_nocomm.fh \
+ include/mpp_reduce_nocomm.fh \
+ include/mpp_sum_nocomm.fh \
+ include/mpp_sum_nocomm_ad.fh \
+ include/mpp_transmit_nocomm.fh \
+ include/mpp_type_nocomm.fh \
+ include/mpp_gather.fh \
+ include/mpp_scatter.fh \
+ include/system_clock.fh
mpp_data_mod.$(FC_MODEXT): \
mpp_parameter_mod.$(FC_MODEXT) \
include/mpp_data_mpi.inc \
@@ -180,39 +180,39 @@ mpp_domains_mod.$(FC_MODEXT): \
include/mpp_domains_comm.inc \
include/mpp_domains_define.inc \
include/mpp_domains_misc.inc \
- include/mpp_do_check.h \
- include/mpp_do_checkV.h \
- include/mpp_do_get_boundary.h \
- include/mpp_do_get_boundary_ad.h \
- include/mpp_do_redistribute.h \
- include/mpp_do_update.h \
- include/mpp_do_updateV.h \
- include/mpp_do_updateV_ad.h \
- include/mpp_do_updateV_nonblock.h \
- include/mpp_do_update_ad.h \
- include/mpp_do_update_nest.h \
- include/mpp_do_update_nonblock.h \
- include/mpp_get_boundary.h \
- include/mpp_get_boundary_ad.h \
- include/mpp_group_update.h \
+ include/mpp_do_check.fh \
+ include/mpp_do_checkV.fh \
+ include/mpp_do_get_boundary.fh \
+ include/mpp_do_get_boundary_ad.fh \
+ include/mpp_do_redistribute.fh \
+ include/mpp_do_update.fh \
+ include/mpp_do_updateV.fh \
+ include/mpp_do_updateV_ad.fh \
+ include/mpp_do_updateV_nonblock.fh \
+ include/mpp_do_update_ad.fh \
+ include/mpp_do_update_nest.fh \
+ include/mpp_do_update_nonblock.fh \
+ include/mpp_get_boundary.fh \
+ include/mpp_get_boundary_ad.fh \
+ include/mpp_group_update.fh \
include/group_update_pack.inc \
include/group_update_unpack.inc \
- include/mpp_update_domains2D.h \
- include/mpp_update_domains2D_ad.h \
- include/mpp_update_domains2D_nonblock.h \
- include/mpp_update_nest_domains.h \
+ include/mpp_update_domains2D.fh \
+ include/mpp_update_domains2D_ad.fh \
+ include/mpp_update_domains2D_nonblock.fh \
+ include/mpp_update_nest_domains.fh \
include/mpp_domains_reduce.inc \
- include/mpp_do_global_field.h \
- include/mpp_do_global_field_ad.h \
- include/mpp_global_field.h \
- include/mpp_global_field_ad.h \
- include/mpp_global_reduce.h \
- include/mpp_global_sum.h \
- include/mpp_global_sum_ad.h \
- include/mpp_global_sum_tl.h \
+ include/mpp_do_global_field.fh \
+ include/mpp_do_global_field_ad.fh \
+ include/mpp_global_field.fh \
+ include/mpp_global_field_ad.fh \
+ include/mpp_global_reduce.fh \
+ include/mpp_global_sum.fh \
+ include/mpp_global_sum_ad.fh \
+ include/mpp_global_sum_tl.fh \
include/mpp_unstruct_domain.inc \
- include/mpp_global_field_ug.h \
- include/mpp_unstruct_pass_data.h
+ include/mpp_global_field_ug.fh \
+ include/mpp_unstruct_pass_data.fh
mpp_io_mod.$(FC_MODEXT): \
mpp_parameter_mod.$(FC_MODEXT) \
mpp_mod.$(FC_MODEXT) \
@@ -221,15 +221,15 @@ mpp_io_mod.$(FC_MODEXT): \
include/mpp_io_misc.inc \
include/mpp_io_connect.inc \
include/mpp_io_read.inc \
- include/mpp_read_2Ddecomp.h \
- include/mpp_read_compressed.h \
+ include/mpp_read_2Ddecomp.fh \
+ include/mpp_read_compressed.fh \
include/mpp_read_distributed_ascii.inc \
- include/mpp_read_distributed_ascii.h \
+ include/mpp_read_distributed_ascii.fh \
include/mpp_io_write.inc \
- include/mpp_write.h \
- include/mpp_write_2Ddecomp.h \
- include/mpp_write_compressed.h \
- include/mpp_write_unlimited_axis.h \
+ include/mpp_write.fh \
+ include/mpp_write_2Ddecomp.fh \
+ include/mpp_write_compressed.fh \
+ include/mpp_write_unlimited_axis.fh \
include/mpp_io_unstructured_write.inc \
include/mpp_io_unstructured_read.inc
mpp_efp_mod.$(FC_MODEXT): mpp_parameter_mod.$(FC_MODEXT) mpp_mod.$(FC_MODEXT)
diff --git a/mpp/include/mpp_alltoall_mpi.h b/mpp/include/mpp_alltoall_mpi.fh
similarity index 100%
rename from mpp/include/mpp_alltoall_mpi.h
rename to mpp/include/mpp_alltoall_mpi.fh
diff --git a/mpp/include/mpp_alltoall_nocomm.h b/mpp/include/mpp_alltoall_nocomm.fh
similarity index 100%
rename from mpp/include/mpp_alltoall_nocomm.h
rename to mpp/include/mpp_alltoall_nocomm.fh
diff --git a/mpp/include/mpp_chksum.h b/mpp/include/mpp_chksum.fh
similarity index 100%
rename from mpp/include/mpp_chksum.h
rename to mpp/include/mpp_chksum.fh
diff --git a/mpp/include/mpp_chksum_int.h b/mpp/include/mpp_chksum_int.fh
similarity index 100%
rename from mpp/include/mpp_chksum_int.h
rename to mpp/include/mpp_chksum_int.fh
diff --git a/mpp/include/mpp_chksum_scalar.h b/mpp/include/mpp_chksum_scalar.fh
similarity index 100%
rename from mpp/include/mpp_chksum_scalar.h
rename to mpp/include/mpp_chksum_scalar.fh
diff --git a/mpp/include/mpp_comm.inc b/mpp/include/mpp_comm.inc
index 22489fbb37..2355102ea9 100644
--- a/mpp/include/mpp_comm.inc
+++ b/mpp/include/mpp_comm.inc
@@ -39,7 +39,7 @@
#define MPP_TYPE_ integer(i8_kind)
#undef MPP_RANK_
#define MPP_RANK_ (:)
-#include
+#include
#undef MPP_CHKSUM_INT_
#define MPP_CHKSUM_INT_ mpp_chksum_i8_2d
@@ -49,7 +49,7 @@
#define MPP_TYPE_ integer(i8_kind)
#undef MPP_RANK_
#define MPP_RANK_ (:,:)
-#include
+#include
#undef MPP_CHKSUM_INT_
#define MPP_CHKSUM_INT_ mpp_chksum_i8_3d
@@ -59,7 +59,7 @@
#define MPP_TYPE_ integer(i8_kind)
#undef MPP_RANK_
#define MPP_RANK_ (:,:,:)
-#include
+#include
#undef MPP_CHKSUM_INT_
#define MPP_CHKSUM_INT_ mpp_chksum_i8_4d
@@ -69,7 +69,7 @@
#define MPP_TYPE_ integer(i8_kind)
#undef MPP_RANK_
#define MPP_RANK_ (:,:,:,:)
-#include
+#include
#undef MPP_CHKSUM_INT_
#define MPP_CHKSUM_INT_ mpp_chksum_i8_5d
@@ -79,7 +79,7 @@
#define MPP_TYPE_ integer(i8_kind)
#undef MPP_RANK_
#define MPP_RANK_ (:,:,:,:,:)
-#include
+#include
#undef MPP_CHKSUM_INT_
#define MPP_CHKSUM_INT_ mpp_chksum_i4_1d
@@ -89,7 +89,7 @@
#define MPP_TYPE_ integer(i4_kind)
#undef MPP_RANK_
#define MPP_RANK_ (:)
-#include
+#include
#undef MPP_CHKSUM_INT_
#define MPP_CHKSUM_INT_ mpp_chksum_i4_2d
@@ -99,7 +99,7 @@
#define MPP_TYPE_ integer(i4_kind)
#undef MPP_RANK_
#define MPP_RANK_ (:,:)
-#include
+#include
#undef MPP_CHKSUM_INT_
#define MPP_CHKSUM_INT_ mpp_chksum_i4_3d
@@ -109,7 +109,7 @@
#define MPP_TYPE_ integer(i4_kind)
#undef MPP_RANK_
#define MPP_RANK_ (:,:,:)
-#include
+#include
#undef MPP_CHKSUM_INT_
#define MPP_CHKSUM_INT_ mpp_chksum_i4_4d
@@ -119,7 +119,7 @@
#define MPP_TYPE_ integer(i4_kind)
#undef MPP_RANK_
#define MPP_RANK_ (:,:,:,:)
-#include
+#include
#undef MPP_CHKSUM_INT_
#define MPP_CHKSUM_INT_ mpp_chksum_i4_5d
@@ -129,7 +129,7 @@
#define MPP_TYPE_ integer(i4_kind)
#undef MPP_RANK_
#define MPP_RANK_ (:,:,:,:,:)
-#include
+#include
#undef MPP_CHKSUM_
#define MPP_CHKSUM_ mpp_chksum_r8_0d
@@ -137,7 +137,7 @@
#define MPP_TYPE_ real(r8_kind)
#undef MPP_RANK_
#define MPP_RANK_ !
-#include
+#include
#undef MPP_CHKSUM_
#define MPP_CHKSUM_ mpp_chksum_r8_1d
@@ -147,7 +147,7 @@
#define MPP_TYPE_ real(r8_kind)
#undef MPP_RANK_
#define MPP_RANK_ (:)
-#include
+#include
#undef MPP_CHKSUM_
#define MPP_CHKSUM_ mpp_chksum_r8_2d
@@ -157,7 +157,7 @@
#define MPP_TYPE_ real(r8_kind)
#undef MPP_RANK_
#define MPP_RANK_ (:,:)
-#include
+#include
#undef MPP_CHKSUM_
#define MPP_CHKSUM_ mpp_chksum_r8_3d
@@ -167,7 +167,7 @@
#define MPP_TYPE_ real(r8_kind)
#undef MPP_RANK_
#define MPP_RANK_ (:,:,:)
-#include
+#include
#undef MPP_CHKSUM_
#define MPP_CHKSUM_ mpp_chksum_r8_4d
@@ -177,7 +177,7 @@
#define MPP_TYPE_ real(r8_kind)
#undef MPP_RANK_
#define MPP_RANK_ (:,:,:,:)
-#include
+#include
#undef MPP_CHKSUM_
#define MPP_CHKSUM_ mpp_chksum_r8_5d
@@ -187,7 +187,7 @@
#define MPP_TYPE_ real(r8_kind)
#undef MPP_RANK_
#define MPP_RANK_ (:,:,:,:,:)
-#include
+#include
#ifdef OVERLOAD_C8
#undef MPP_CHKSUM_
@@ -198,7 +198,7 @@
#define MPP_TYPE_ complex(c8_kind)
#undef MPP_RANK_
#define MPP_RANK_ !
-#include
+#include
#undef MPP_CHKSUM_
#define MPP_CHKSUM_ mpp_chksum_c8_1d
@@ -208,7 +208,7 @@
#define MPP_TYPE_ complex(c8_kind)
#undef MPP_RANK_
#define MPP_RANK_ (:)
-#include
+#include
#undef MPP_CHKSUM_
#define MPP_CHKSUM_ mpp_chksum_c8_2d
@@ -218,7 +218,7 @@
#define MPP_TYPE_ complex(c8_kind)
#undef MPP_RANK_
#define MPP_RANK_ (:,:)
-#include
+#include
#undef MPP_CHKSUM_
#define MPP_CHKSUM_ mpp_chksum_c8_3d
@@ -228,7 +228,7 @@
#define MPP_TYPE_ complex(c8_kind)
#undef MPP_RANK_
#define MPP_RANK_ (:,:,:)
-#include
+#include
#undef MPP_CHKSUM_
#define MPP_CHKSUM_ mpp_chksum_c8_4d
@@ -238,7 +238,7 @@
#define MPP_TYPE_ complex(c8_kind)
#undef MPP_RANK_
#define MPP_RANK_ (:,:,:,:)
-#include
+#include
#undef MPP_CHKSUM_
#define MPP_CHKSUM_ mpp_chksum_c8_5d
@@ -248,7 +248,7 @@
#define MPP_TYPE_ complex(c8_kind)
#undef MPP_RANK_
#define MPP_RANK_ (:,:,:,:,:)
-#include
+#include
#endif
#undef MPP_CHKSUM_
@@ -259,7 +259,7 @@
#define MPP_TYPE_ real(r4_kind)
#undef MPP_RANK_
#define MPP_RANK_ !
-#include
+#include
#undef MPP_CHKSUM_
#define MPP_CHKSUM_ mpp_chksum_r4_1d
@@ -269,7 +269,7 @@
#define MPP_TYPE_ real(r4_kind)
#undef MPP_RANK_
#define MPP_RANK_ (:)
-#include
+#include
#undef MPP_CHKSUM_
#define MPP_CHKSUM_ mpp_chksum_r4_2d
@@ -279,7 +279,7 @@
#define MPP_TYPE_ real(r4_kind)
#undef MPP_RANK_
#define MPP_RANK_ (:,:)
-#include
+#include
#undef MPP_CHKSUM_
#define MPP_CHKSUM_ mpp_chksum_r4_3d
@@ -289,7 +289,7 @@
#define MPP_TYPE_ real(r4_kind)
#undef MPP_RANK_
#define MPP_RANK_ (:,:,:)
-#include
+#include
#undef MPP_CHKSUM_
#define MPP_CHKSUM_ mpp_chksum_r4_4d
@@ -299,7 +299,7 @@
#define MPP_TYPE_ real(r4_kind)
#undef MPP_RANK_
#define MPP_RANK_ (:,:,:,:)
-#include
+#include
#undef MPP_CHKSUM_
#define MPP_CHKSUM_ mpp_chksum_r4_5d
@@ -309,7 +309,7 @@
#define MPP_TYPE_ real(r4_kind)
#undef MPP_RANK_
#define MPP_RANK_ (:,:,:,:,:)
-#include
+#include
#ifdef OVERLOAD_C4
#undef MPP_CHKSUM_
@@ -320,7 +320,7 @@
#define MPP_TYPE_ complex(c4_kind)
#undef MPP_RANK_
#define MPP_RANK_ !
-#include
+#include
#undef MPP_CHKSUM_
#define MPP_CHKSUM_ mpp_chksum_c4_1d
@@ -330,7 +330,7 @@
#define MPP_TYPE_ complex(c4_kind)
#undef MPP_RANK_
#define MPP_RANK_ (:)
-#include
+#include
#undef MPP_CHKSUM_
#define MPP_CHKSUM_ mpp_chksum_c4_2d
@@ -340,7 +340,7 @@
#define MPP_TYPE_ complex(c4_kind)
#undef MPP_RANK_
#define MPP_RANK_ (:,:)
-#include
+#include
#undef MPP_CHKSUM_
#define MPP_CHKSUM_ mpp_chksum_c4_3d
@@ -350,7 +350,7 @@
#define MPP_TYPE_ complex(c4_kind)
#undef MPP_RANK_
#define MPP_RANK_ (:,:,:)
-#include
+#include
#undef MPP_CHKSUM_
#define MPP_CHKSUM_ mpp_chksum_c4_4d
@@ -360,7 +360,7 @@
#define MPP_TYPE_ complex(c4_kind)
#undef MPP_RANK_
#define MPP_RANK_ (:,:,:,:)
-#include
+#include
#undef MPP_CHKSUM_
#define MPP_CHKSUM_ mpp_chksum_c4_5d
@@ -370,7 +370,7 @@
#define MPP_TYPE_ complex(c4_kind)
#undef MPP_RANK_
#define MPP_RANK_ (:,:,:,:,:)
-#include
+#include
#endif
!#################################################
@@ -384,7 +384,7 @@
#undef MPP_GATHER_PELIST_3D_
#define MPP_GATHER_PELIST_2D_ mpp_gather_pelist_logical_2d
#define MPP_GATHER_PELIST_3D_ mpp_gather_pelist_logical_3d
-#include
+#include
#undef MPP_GATHER_1D_
#undef MPP_GATHER_1DV_
@@ -396,7 +396,7 @@
#undef MPP_GATHER_PELIST_3D_
#define MPP_GATHER_PELIST_2D_ mpp_gather_pelist_int4_2d
#define MPP_GATHER_PELIST_3D_ mpp_gather_pelist_int4_3d
-#include
+#include
#undef MPP_GATHER_1D_
#undef MPP_GATHER_1DV_
@@ -408,7 +408,7 @@
#undef MPP_GATHER_PELIST_3D_
#define MPP_GATHER_PELIST_2D_ mpp_gather_pelist_real4_2d
#define MPP_GATHER_PELIST_3D_ mpp_gather_pelist_real4_3d
-#include
+#include
#undef MPP_GATHER_1D_
#undef MPP_GATHER_1DV_
@@ -420,7 +420,7 @@
#undef MPP_GATHER_PELIST_3D_
#define MPP_GATHER_PELIST_2D_ mpp_gather_pelist_real8_2d
#define MPP_GATHER_PELIST_3D_ mpp_gather_pelist_real8_3d
-#include
+#include
!#################################################
#undef MPP_SCATTER_PELIST_2D_
@@ -429,7 +429,7 @@
#define MPP_TYPE_ integer(i4_kind)
#define MPP_SCATTER_PELIST_2D_ mpp_scatter_pelist_int4_2d
#define MPP_SCATTER_PELIST_3D_ mpp_scatter_pelist_int4_3d
-#include
+#include
#undef MPP_SCATTER_PELIST_2D_
#undef MPP_SCATTER_PELIST_3D_
@@ -437,7 +437,7 @@
#define MPP_TYPE_ real(r4_kind)
#define MPP_SCATTER_PELIST_2D_ mpp_scatter_pelist_real4_2d
#define MPP_SCATTER_PELIST_3D_ mpp_scatter_pelist_real4_3d
-#include
+#include
#undef MPP_SCATTER_PELIST_2D_
#undef MPP_SCATTER_PELIST_3D_
@@ -445,5 +445,5 @@
#define MPP_TYPE_ real(r8_kind)
#define MPP_SCATTER_PELIST_2D_ mpp_scatter_pelist_real8_2d
#define MPP_SCATTER_PELIST_3D_ mpp_scatter_pelist_real8_3d
-#include
+#include
!> @}
diff --git a/mpp/include/mpp_comm_mpi.inc b/mpp/include/mpp_comm_mpi.inc
index ac6df35894..928f9fcb92 100644
--- a/mpp/include/mpp_comm_mpi.inc
+++ b/mpp/include/mpp_comm_mpi.inc
@@ -414,7 +414,7 @@ end subroutine mpp_exit
#define MPP_TYPE_BYTELEN_ 8
#undef MPI_TYPE_
#define MPI_TYPE_ MPI_REAL8
-#include
+#include
#ifdef OVERLOAD_C8
#undef MPP_TRANSMIT_
@@ -471,7 +471,7 @@ end subroutine mpp_exit
#define MPP_TYPE_BYTELEN_ 16
#undef MPI_TYPE_
#define MPI_TYPE_ MPI_DOUBLE_COMPLEX
-#include
+#include
#endif
#undef MPP_TRANSMIT_
@@ -528,7 +528,7 @@ end subroutine mpp_exit
#define MPP_TYPE_BYTELEN_ 4
#undef MPI_TYPE_
#define MPI_TYPE_ MPI_REAL4
-#include
+#include
#ifdef OVERLOAD_C4
#undef MPP_TRANSMIT_
@@ -585,7 +585,7 @@ end subroutine mpp_exit
#define MPP_TYPE_BYTELEN_ 8
#undef MPI_TYPE_
#define MPI_TYPE_ MPI_COMPLEX
-#include
+#include
#endif
#undef MPP_TYPE_INIT_VALUE
@@ -644,7 +644,7 @@ end subroutine mpp_exit
#define MPP_TYPE_BYTELEN_ 8
#undef MPI_TYPE_
#define MPI_TYPE_ MPI_INTEGER8
-#include
+#include
#undef MPP_TRANSMIT_
#define MPP_TRANSMIT_ mpp_transmit_int4
@@ -700,7 +700,7 @@ end subroutine mpp_exit
#define MPP_TYPE_BYTELEN_ 4
#undef MPI_TYPE_
#define MPI_TYPE_ MPI_INTEGER4
-#include
+#include
#undef MPP_TYPE_INIT_VALUE
#define MPP_TYPE_INIT_VALUE .false.
@@ -758,7 +758,7 @@ end subroutine mpp_exit
#define MPP_TYPE_BYTELEN_ 8
#undef MPI_TYPE_
#define MPI_TYPE_ MPI_INTEGER8
-#include
+#include
#undef MPP_TRANSMIT_
#define MPP_TRANSMIT_ mpp_transmit_logical4
@@ -814,7 +814,7 @@ end subroutine mpp_exit
#define MPP_TYPE_BYTELEN_ 4
#undef MPI_TYPE_
#define MPI_TYPE_ MPI_INTEGER4
-#include
+#include
#undef MPP_TYPE_INIT_VALUE
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
@@ -834,7 +834,7 @@ end subroutine mpp_exit
#define MPI_TYPE_ MPI_REAL8
#undef MPI_REDUCE_
#define MPI_REDUCE_ MPI_MAX
-#include
+#include
#undef MPP_REDUCE_0D_
#define MPP_REDUCE_0D_ mpp_max_real4_0d
@@ -848,7 +848,7 @@ end subroutine mpp_exit
#define MPI_TYPE_ MPI_REAL4
#undef MPI_REDUCE_
#define MPI_REDUCE_ MPI_MAX
-#include
+#include
#undef MPP_REDUCE_0D_
#define MPP_REDUCE_0D_ mpp_max_int8_0d
@@ -862,7 +862,7 @@ end subroutine mpp_exit
#define MPI_TYPE_ MPI_INTEGER8
#undef MPI_REDUCE_
#define MPI_REDUCE_ MPI_MAX
-#include
+#include
#undef MPP_REDUCE_0D_
#define MPP_REDUCE_0D_ mpp_max_int4_0d
@@ -876,7 +876,7 @@ end subroutine mpp_exit
#define MPI_TYPE_ MPI_INTEGER4
#undef MPI_REDUCE_
#define MPI_REDUCE_ MPI_MAX
-#include
+#include
#undef MPP_REDUCE_0D_
#define MPP_REDUCE_0D_ mpp_min_real8_0d
@@ -890,7 +890,7 @@ end subroutine mpp_exit
#define MPI_TYPE_ MPI_REAL8
#undef MPI_REDUCE_
#define MPI_REDUCE_ MPI_MIN
-#include
+#include
#undef MPP_REDUCE_0D_
#define MPP_REDUCE_0D_ mpp_min_real4_0d
@@ -904,7 +904,7 @@ end subroutine mpp_exit
#define MPI_TYPE_ MPI_REAL4
#undef MPI_REDUCE_
#define MPI_REDUCE_ MPI_MIN
-#include
+#include
#undef MPP_REDUCE_0D_
#define MPP_REDUCE_0D_ mpp_min_int8_0d
@@ -918,7 +918,7 @@ end subroutine mpp_exit
#define MPI_TYPE_ MPI_INTEGER8
#undef MPI_REDUCE_
#define MPI_REDUCE_ MPI_MIN
-#include
+#include
#undef MPP_REDUCE_0D_
#define MPP_REDUCE_0D_ mpp_min_int4_0d
@@ -932,7 +932,7 @@ end subroutine mpp_exit
#define MPI_TYPE_ MPI_INTEGER4
#undef MPI_REDUCE_
#define MPI_REDUCE_ MPI_MIN
-#include
+#include
#undef MPP_SUM_
#define MPP_SUM_ mpp_sum_real8
@@ -952,7 +952,7 @@ end subroutine mpp_exit
#define MPI_TYPE_ MPI_REAL8
#undef MPP_TYPE_BYTELEN_
#define MPP_TYPE_BYTELEN_ 8
-#include
+#include
#ifdef OVERLOAD_C8
#undef MPP_SUM_
@@ -973,7 +973,7 @@ end subroutine mpp_exit
#define MPI_TYPE_ MPI_DOUBLE_COMPLEX
#undef MPP_TYPE_BYTELEN_
#define MPP_TYPE_BYTELEN_ 16
-#include
+#include
#endif
#undef MPP_SUM_
@@ -994,7 +994,7 @@ end subroutine mpp_exit
#define MPI_TYPE_ MPI_REAL4
#undef MPP_TYPE_BYTELEN_
#define MPP_TYPE_BYTELEN_ 4
-#include
+#include
#ifdef OVERLOAD_C4
#undef MPP_SUM_
@@ -1015,7 +1015,7 @@ end subroutine mpp_exit
#define MPI_TYPE_ MPI_COMPLEX
#undef MPP_TYPE_BYTELEN_
#define MPP_TYPE_BYTELEN_ 8
-#include
+#include
#endif
#undef MPP_SUM_
@@ -1036,7 +1036,7 @@ end subroutine mpp_exit
#define MPI_TYPE_ MPI_INTEGER8
#undef MPP_TYPE_BYTELEN_
#define MPP_TYPE_BYTELEN_ 8
-#include
+#include
#undef MPP_SUM_
#define MPP_SUM_ mpp_sum_int4
@@ -1056,7 +1056,7 @@ end subroutine mpp_exit
#define MPI_TYPE_ MPI_INTEGER4
#undef MPP_TYPE_BYTELEN_
#define MPP_TYPE_BYTELEN_ 4
-#include
+#include
!--------------------------------
#undef MPP_SUM_AD_
#define MPP_SUM_AD_ mpp_sum_real8_ad
@@ -1076,7 +1076,7 @@ end subroutine mpp_exit
#define MPI_TYPE_ MPI_REAL8
#undef MPP_TYPE_BYTELEN_
#define MPP_TYPE_BYTELEN_ 8
-#include
+#include
#ifdef OVERLOAD_C8
#undef MPP_SUM_AD_
@@ -1097,7 +1097,7 @@ end subroutine mpp_exit
#define MPI_TYPE_ MPI_DOUBLE_COMPLEX
#undef MPP_TYPE_BYTELEN_
#define MPP_TYPE_BYTELEN_ 16
-#include
+#include
#endif
#undef MPP_SUM_AD_
@@ -1118,7 +1118,7 @@ end subroutine mpp_exit
#define MPI_TYPE_ MPI_REAL4
#undef MPP_TYPE_BYTELEN_
#define MPP_TYPE_BYTELEN_ 4
-#include
+#include
#ifdef OVERLOAD_C4
#undef MPP_SUM_AD_
@@ -1139,7 +1139,7 @@ end subroutine mpp_exit
#define MPI_TYPE_ MPI_COMPLEX
#undef MPP_TYPE_BYTELEN_
#define MPP_TYPE_BYTELEN_ 8
-#include
+#include
#endif
#undef MPP_SUM_AD_
@@ -1160,7 +1160,7 @@ end subroutine mpp_exit
#define MPI_TYPE_ MPI_INTEGER8
#undef MPP_TYPE_BYTELEN_
#define MPP_TYPE_BYTELEN_ 8
-#include
+#include
#undef MPP_SUM_AD_
#define MPP_SUM_AD_ mpp_sum_int4_ad
@@ -1180,7 +1180,7 @@ end subroutine mpp_exit
#define MPI_TYPE_ MPI_INTEGER4
#undef MPP_TYPE_BYTELEN_
#define MPP_TYPE_BYTELEN_ 4
-#include
+#include
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! !
@@ -1200,7 +1200,7 @@ end subroutine mpp_exit
#define MPP_TYPE_ integer(i4_kind)
#define MPP_TYPE_BYTELEN_ 4
#define MPI_TYPE_ MPI_INTEGER4
-#include
+#include
#undef MPP_ALLTOALL_
#undef MPP_ALLTOALLV_
@@ -1214,7 +1214,7 @@ end subroutine mpp_exit
#define MPP_TYPE_ integer(i8_kind)
#define MPP_TYPE_BYTELEN_ 8
#define MPI_TYPE_ MPI_INTEGER8
-#include
+#include
#undef MPP_ALLTOALL_
#undef MPP_ALLTOALLV_
@@ -1228,7 +1228,7 @@ end subroutine mpp_exit
#define MPP_TYPE_ real(r4_kind)
#define MPP_TYPE_BYTELEN_ 4
#define MPI_TYPE_ MPI_REAL4
-#include
+#include
#undef MPP_ALLTOALL_
#undef MPP_ALLTOALLV_
@@ -1242,7 +1242,7 @@ end subroutine mpp_exit
#define MPP_TYPE_ real(r8_kind)
#define MPP_TYPE_BYTELEN_ 8
#define MPI_TYPE_ MPI_REAL8
-#include
+#include
#ifdef OVERLOAD_C4
#undef MPP_ALLTOALL_
@@ -1257,7 +1257,7 @@ end subroutine mpp_exit
#define MPP_TYPE_ complex(c4_kind)
#define MPP_TYPE_BYTELEN_ 8
#define MPI_TYPE_ MPI_COMPLEX8
-#include
+#include
#endif
#ifdef OVERLOAD_C8
@@ -1273,7 +1273,7 @@ end subroutine mpp_exit
#define MPP_TYPE_ complex(c8_kind)
#define MPP_TYPE_BYTELEN_ 16
#define MPI_TYPE_ MPI_COMPLEX16
-#include
+#include
#endif
#undef MPP_ALLTOALL_
@@ -1288,7 +1288,7 @@ end subroutine mpp_exit
#define MPP_TYPE_ logical(l4_kind)
#define MPP_TYPE_BYTELEN_ 4
#define MPI_TYPE_ MPI_INTEGER4
-#include
+#include
#undef MPP_ALLTOALL_
#undef MPP_ALLTOALLV_
@@ -1302,7 +1302,7 @@ end subroutine mpp_exit
#define MPP_TYPE_ logical(l8_kind)
#define MPP_TYPE_BYTELEN_ 8
#define MPI_TYPE_ MPI_INTEGER8
-#include
+#include
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! !
@@ -1316,7 +1316,7 @@ end subroutine mpp_exit
#define MPP_TYPE_CREATE_ mpp_type_create_int4
#define MPP_TYPE_ integer(i4_kind)
#define MPI_TYPE_ MPI_INTEGER4
-#include
+#include
#undef MPP_TYPE_CREATE_
#undef MPP_TYPE_
@@ -1324,7 +1324,7 @@ end subroutine mpp_exit
#define MPP_TYPE_CREATE_ mpp_type_create_int8
#define MPP_TYPE_ integer(i8_kind)
#define MPI_TYPE_ MPI_INTEGER8
-#include
+#include
#undef MPP_TYPE_CREATE_
#undef MPP_TYPE_
@@ -1332,7 +1332,7 @@ end subroutine mpp_exit
#define MPP_TYPE_CREATE_ mpp_type_create_real4
#define MPP_TYPE_ real(r4_kind)
#define MPI_TYPE_ MPI_REAL4
-#include
+#include
#undef MPP_TYPE_CREATE_
#undef MPP_TYPE_
@@ -1340,7 +1340,7 @@ end subroutine mpp_exit
#define MPP_TYPE_CREATE_ mpp_type_create_real8
#define MPP_TYPE_ real(r8_kind)
#define MPI_TYPE_ MPI_REAL8
-#include
+#include
#undef MPP_TYPE_CREATE_
#undef MPP_TYPE_
@@ -1348,7 +1348,7 @@ end subroutine mpp_exit
#define MPP_TYPE_CREATE_ mpp_type_create_cmplx4
#define MPP_TYPE_ complex(c4_kind)
#define MPI_TYPE_ MPI_COMPLEX8
-#include
+#include
#undef MPP_TYPE_CREATE_
#undef MPP_TYPE_
@@ -1356,7 +1356,7 @@ end subroutine mpp_exit
#define MPP_TYPE_CREATE_ mpp_type_create_cmplx8
#define MPP_TYPE_ complex(c8_kind)
#define MPI_TYPE_ MPI_COMPLEX16
-#include
+#include
#undef MPP_TYPE_CREATE_
#undef MPP_TYPE_
@@ -1364,7 +1364,7 @@ end subroutine mpp_exit
#define MPP_TYPE_CREATE_ mpp_type_create_logical4
#define MPP_TYPE_ logical(l4_kind)
#define MPI_TYPE_ MPI_INTEGER4
-#include
+#include
#undef MPP_TYPE_CREATE_
#undef MPP_TYPE_
@@ -1372,7 +1372,7 @@ end subroutine mpp_exit
#define MPP_TYPE_CREATE_ mpp_type_create_logical8
#define MPP_TYPE_ logical(l8_kind)
#define MPI_TYPE_ MPI_INTEGER8
-#include
+#include
! Clear preprocessor flags
#undef MPI_TYPE_
diff --git a/mpp/include/mpp_comm_nocomm.inc b/mpp/include/mpp_comm_nocomm.inc
index efe25205e3..a1d849b831 100644
--- a/mpp/include/mpp_comm_nocomm.inc
+++ b/mpp/include/mpp_comm_nocomm.inc
@@ -317,7 +317,7 @@ end subroutine mpp_exit
#define MPP_TYPE_BYTELEN_ 8
#undef MPI_TYPE_
#define MPI_TYPE_ MPI_REAL8
-#include
+#include
#ifdef OVERLOAD_C8
#undef MPP_TRANSMIT_
@@ -374,7 +374,7 @@ end subroutine mpp_exit
#define MPP_TYPE_BYTELEN_ 16
#undef MPI_TYPE_
#define MPI_TYPE_ MPI_DOUBLE_COMPLEX
-#include
+#include
#endif
#undef MPP_TRANSMIT_
@@ -431,7 +431,7 @@ end subroutine mpp_exit
#define MPP_TYPE_BYTELEN_ 4
#undef MPI_TYPE_
#define MPI_TYPE_ MPI_REAL4
-#include
+#include
#ifdef OVERLOAD_C4
#undef MPP_TRANSMIT_
@@ -488,7 +488,7 @@ end subroutine mpp_exit
#define MPP_TYPE_BYTELEN_ 8
#undef MPI_TYPE_
#define MPI_TYPE_ MPI_COMPLEX
-#include
+#include
#endif
#undef MPP_TYPE_INIT_VALUE
@@ -547,7 +547,7 @@ end subroutine mpp_exit
#define MPP_TYPE_BYTELEN_ 8
#undef MPI_TYPE_
#define MPI_TYPE_ MPI_INTEGER8
-#include
+#include
#undef MPP_TRANSMIT_
#define MPP_TRANSMIT_ mpp_transmit_int4
@@ -603,7 +603,7 @@ end subroutine mpp_exit
#define MPP_TYPE_BYTELEN_ 4
#undef MPI_TYPE_
#define MPI_TYPE_ MPI_INTEGER4
-#include
+#include
#undef MPP_TYPE_INIT_VALUE
#define MPP_TYPE_INIT_VALUE .false.
@@ -661,7 +661,7 @@ end subroutine mpp_exit
#define MPP_TYPE_BYTELEN_ 8
#undef MPI_TYPE_
#define MPI_TYPE_ MPI_INTEGER8
-#include
+#include
#undef MPP_TRANSMIT_
#define MPP_TRANSMIT_ mpp_transmit_logical4
@@ -717,7 +717,7 @@ end subroutine mpp_exit
#define MPP_TYPE_BYTELEN_ 4
#undef MPI_TYPE_
#define MPI_TYPE_ MPI_INTEGER4
-#include
+#include
#undef MPP_TYPE_INIT_VALUE
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
@@ -737,7 +737,7 @@ end subroutine mpp_exit
#define MPI_TYPE_ MPI_REAL8
#undef MPI_REDUCE_
#define MPI_REDUCE_ MPI_MAX
-#include
+#include
#undef MPP_REDUCE_0D_
#define MPP_REDUCE_0D_ mpp_max_real4_0d
@@ -751,7 +751,7 @@ end subroutine mpp_exit
#define MPI_TYPE_ MPI_REAL4
#undef MPI_REDUCE_
#define MPI_REDUCE_ MPI_MAX
-#include
+#include
#undef MPP_REDUCE_0D_
#define MPP_REDUCE_0D_ mpp_max_int8_0d
@@ -765,7 +765,7 @@ end subroutine mpp_exit
#define MPI_TYPE_ MPI_INTEGER8
#undef MPI_REDUCE_
#define MPI_REDUCE_ MPI_MAX
-#include
+#include
#undef MPP_REDUCE_0D_
#define MPP_REDUCE_0D_ mpp_max_int4_0d
@@ -779,7 +779,7 @@ end subroutine mpp_exit
#define MPI_TYPE_ MPI_INTEGER4
#undef MPI_REDUCE_
#define MPI_REDUCE_ MPI_MAX
-#include
+#include
#undef MPP_REDUCE_0D_
#define MPP_REDUCE_0D_ mpp_min_real8_0d
@@ -793,7 +793,7 @@ end subroutine mpp_exit
#define MPI_TYPE_ MPI_REAL8
#undef MPI_REDUCE_
#define MPI_REDUCE_ MPI_MIN
-#include
+#include
#undef MPP_REDUCE_0D_
#define MPP_REDUCE_0D_ mpp_min_real4_0d
@@ -807,7 +807,7 @@ end subroutine mpp_exit
#define MPI_TYPE_ MPI_REAL4
#undef MPI_REDUCE_
#define MPI_REDUCE_ MPI_MIN
-#include
+#include
#undef MPP_REDUCE_0D_
#define MPP_REDUCE_0D_ mpp_min_int8_0d
@@ -821,7 +821,7 @@ end subroutine mpp_exit
#define MPI_TYPE_ MPI_INTEGER8
#undef MPI_REDUCE_
#define MPI_REDUCE_ MPI_MIN
-#include
+#include
#undef MPP_REDUCE_0D_
#define MPP_REDUCE_0D_ mpp_min_int4_0d
@@ -835,7 +835,7 @@ end subroutine mpp_exit
#define MPI_TYPE_ MPI_INTEGER4
#undef MPI_REDUCE_
#define MPI_REDUCE_ MPI_MIN
-#include
+#include
#undef MPP_SUM_
#define MPP_SUM_ mpp_sum_real8
@@ -855,7 +855,7 @@ end subroutine mpp_exit
#define MPI_TYPE_ MPI_REAL8
#undef MPP_TYPE_BYTELEN_
#define MPP_TYPE_BYTELEN_ 8
-#include
+#include
#ifdef OVERLOAD_C8
#undef MPP_SUM_
@@ -876,7 +876,7 @@ end subroutine mpp_exit
#define MPI_TYPE_ MPI_DOUBLE_COMPLEX
#undef MPP_TYPE_BYTELEN_
#define MPP_TYPE_BYTELEN_ 16
-#include
+#include
#endif
#undef MPP_SUM_
@@ -897,7 +897,7 @@ end subroutine mpp_exit
#define MPI_TYPE_ MPI_REAL4
#undef MPP_TYPE_BYTELEN_
#define MPP_TYPE_BYTELEN_ 4
-#include
+#include
#ifdef OVERLOAD_C4
#undef MPP_SUM_
@@ -918,7 +918,7 @@ end subroutine mpp_exit
#define MPI_TYPE_ MPI_COMPLEX
#undef MPP_TYPE_BYTELEN_
#define MPP_TYPE_BYTELEN_ 8
-#include
+#include
#endif
#undef MPP_SUM_
@@ -939,7 +939,7 @@ end subroutine mpp_exit
#define MPI_TYPE_ MPI_INTEGER8
#undef MPP_TYPE_BYTELEN_
#define MPP_TYPE_BYTELEN_ 8
-#include
+#include
#undef MPP_SUM_
#define MPP_SUM_ mpp_sum_int4
@@ -959,7 +959,7 @@ end subroutine mpp_exit
#define MPI_TYPE_ MPI_INTEGER4
#undef MPP_TYPE_BYTELEN_
#define MPP_TYPE_BYTELEN_ 4
-#include
+#include
!--------------------------------
#undef MPP_SUM_AD_
#define MPP_SUM_AD_ mpp_sum_real8_ad
@@ -979,7 +979,7 @@ end subroutine mpp_exit
#define MPI_TYPE_ MPI_REAL8
#undef MPP_TYPE_BYTELEN_
#define MPP_TYPE_BYTELEN_ 8
-#include
+#include
#ifdef OVERLOAD_C8
#undef MPP_SUM_AD_
@@ -1000,7 +1000,7 @@ end subroutine mpp_exit
#define MPI_TYPE_ MPI_DOUBLE_COMPLEX
#undef MPP_TYPE_BYTELEN_
#define MPP_TYPE_BYTELEN_ 16
-#include
+#include
#endif
#undef MPP_SUM_AD_
@@ -1021,7 +1021,7 @@ end subroutine mpp_exit
#define MPI_TYPE_ MPI_REAL4
#undef MPP_TYPE_BYTELEN_
#define MPP_TYPE_BYTELEN_ 4
-#include
+#include
#ifdef OVERLOAD_C4
#undef MPP_SUM_AD_
@@ -1042,7 +1042,7 @@ end subroutine mpp_exit
#define MPI_TYPE_ MPI_COMPLEX
#undef MPP_TYPE_BYTELEN_
#define MPP_TYPE_BYTELEN_ 8
-#include
+#include
#endif
#undef MPP_SUM_AD_
@@ -1063,7 +1063,7 @@ end subroutine mpp_exit
#define MPI_TYPE_ MPI_INTEGER8
#undef MPP_TYPE_BYTELEN_
#define MPP_TYPE_BYTELEN_ 8
-#include
+#include
#undef MPP_SUM_AD_
#define MPP_SUM_AD_ mpp_sum_int4_ad
@@ -1083,7 +1083,7 @@ end subroutine mpp_exit
#define MPI_TYPE_ MPI_INTEGER4
#undef MPP_TYPE_BYTELEN_
#define MPP_TYPE_BYTELEN_ 4
-#include
+#include
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! !
@@ -1103,7 +1103,7 @@ end subroutine mpp_exit
#define MPP_TYPE_ integer(i4_kind)
#define MPP_TYPE_BYTELEN_ 4
#define MPI_TYPE_ MPI_INTEGER4
-#include
+#include
#undef MPP_ALLTOALL_
#undef MPP_ALLTOALLV_
@@ -1117,7 +1117,7 @@ end subroutine mpp_exit
#define MPP_TYPE_ integer(i8_kind)
#define MPP_TYPE_BYTELEN_ 8
#define MPI_TYPE_ MPI_INTEGER8
-#include
+#include
#undef MPP_ALLTOALL_
#undef MPP_ALLTOALLV_
@@ -1131,7 +1131,7 @@ end subroutine mpp_exit
#define MPP_TYPE_ real(r4_kind)
#define MPP_TYPE_BYTELEN_ 4
#define MPI_TYPE_ MPI_REAL4
-#include
+#include
#undef MPP_ALLTOALL_
#undef MPP_ALLTOALLV_
@@ -1145,7 +1145,7 @@ end subroutine mpp_exit
#define MPP_TYPE_ real(r8_kind)
#define MPP_TYPE_BYTELEN_ 8
#define MPI_TYPE_ MPI_REAL8
-#include
+#include
#undef MPP_ALLTOALL_
#undef MPP_ALLTOALLV_
@@ -1159,7 +1159,7 @@ end subroutine mpp_exit
#define MPP_TYPE_ logical(l4_kind)
#define MPP_TYPE_BYTELEN_ 4
#define MPI_TYPE_ MPI_INTEGER4
-#include
+#include
#undef MPP_ALLTOALL_
#undef MPP_ALLTOALLV_
@@ -1173,7 +1173,7 @@ end subroutine mpp_exit
#define MPP_TYPE_ logical(l8_kind)
#define MPP_TYPE_BYTELEN_ 8
#define MPI_TYPE_ MPI_INTEGER8
-#include
+#include
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! !
! DATA TRANSFER TYPES: mpp_type_create, mpp_type_free !
@@ -1183,42 +1183,42 @@ end subroutine mpp_exit
#define MPP_TYPE_CREATE_ mpp_type_create_int4
#define MPP_TYPE_ integer(i4_kind)
#define MPI_TYPE_ MPI_INTEGER4
-#include
+#include
#define MPP_TYPE_CREATE_ mpp_type_create_int8
#define MPP_TYPE_ integer(i8_kind)
#define MPI_TYPE_ MPI_INTEGER8
-#include
+#include