/

srcF/sfx_data.f90

@@ -40,13 +40,21 @@ module sfx_data
     !> @brief meteorological input for surface flux calculation
     !> &details using arrays as input
     type, public :: meteoDataVecType
-
+#if defined(INCLUDE_CXX)
         real, pointer :: h(:)       !< constant flux layer height [m]
         real, pointer :: U(:)       !< abs(wind speed) at 'h' [m/s]
         real, pointer :: dT(:)      !< difference between potential temperature at 'h' and at surface [K]
         real, pointer :: Tsemi(:)   !< semi-sum of potential temperature at 'h' and at surface [K]
         real, pointer :: dQ(:)      !< difference between humidity at 'h' and at surface [g/g]
         real, pointer :: z0_m(:)    !< surface aerodynamic roughness (should be < 0 for water bodies surface)
+#else
+        real, allocatable :: h(:)       !< constant flux layer height [m]
+        real, allocatable :: U(:)       !< abs(wind speed) at 'h' [m/s]
+        real, allocatable :: dT(:)      !< difference between potential temperature at 'h' and at surface [K]
+        real, allocatable :: Tsemi(:)   !< semi-sum of potential temperature at 'h' and at surface [K]
+        real, allocatable :: dQ(:)      !< difference between humidity at 'h' and at surface [g/g]
+        real, allocatable :: z0_m(:)    !< surface aerodynamic roughness (should be < 0 for water bodies surface)
+#endif
     end type
 
 #if defined(INCLUDE_CXX)
@@ -82,7 +90,7 @@ module sfx_data
     !> @brief surface flux output data
     !> &details using arrays as output
     type, public :: sfxDataVecType
-
+#if defined(INCLUDE_CXX)
         real, pointer :: zeta(:)            !< = z/L [n/d]
         real, pointer :: Rib(:)             !< bulk Richardson number [n/d]
         real, pointer :: Re(:)              !< Reynolds number [n/d]
@@ -94,6 +102,19 @@ module sfx_data
         real, pointer :: Ct(:)              !< transfer coefficient for heat [n/d]
         real, pointer :: Km(:)              !< eddy viscosity coeff. at h [m^2/s]
         real, pointer :: Pr_t_inv(:)        !< inverse turbulent Prandtl number at h [n/d]
+#else
+        real, allocatable :: zeta(:)            !< = z/L [n/d]
+        real, allocatable :: Rib(:)             !< bulk Richardson number [n/d]
+        real, allocatable :: Re(:)              !< Reynolds number [n/d]
+        real, allocatable :: B(:)               !< = log(z0_m / z0_h) [n/d]
+        real, allocatable :: z0_m(:)            !< aerodynamic roughness [m]
+        real, allocatable :: z0_t(:)            !< thermal roughness [m]
+        real, allocatable :: Rib_conv_lim(:)    !< Ri-bulk convection critical value [n/d]
+        real, allocatable :: Cm(:)              !< transfer coefficient for momentum [n/d]
+        real, allocatable :: Ct(:)              !< transfer coefficient for heat [n/d]
+        real, allocatable :: Km(:)              !< eddy viscosity coeff. at h [m^2/s]
+        real, allocatable :: Pr_t_inv(:)        !< inverse turbulent Prandtl number at h [n/d]
+#endif
     end type
 
 #if defined(INCLUDE_CXX)