Index: palm/trunk/SOURCE/biometeorology_pt_mod.f90
===================================================================
--- palm/trunk/SOURCE/biometeorology_pt_mod.f90	(revision 3448)
+++ palm/trunk/SOURCE/biometeorology_pt_mod.f90	(revision 3464)
@@ -20,12 +20,17 @@
 !--------------------------------------------------------------------------------!
 !
-! Current revisions: 001
-! -----------------
+! Current revisions:
+! ------------------
 ! 
 ! 
-! Former revisions: 001
+! Former revisions:
 ! -----------------
 ! $Id$
-! Initial revision 001
+! From branch resler@3462, pavelkrc:
+! replace specific precision DLOG, DMOD intrinsics with generic precision for
+! compatibility with 'wp'
+! 
+! 
+! Initial revision
 !
 ! 
@@ -116,5 +121,5 @@
 ! Description:
 ! ------------
-!> PT_BASIC.F90   Version of perceived temperature (PT, °C) for
+!> PT_BASIC.F90   Version of perceived temperature (PT, ï¿œC) for
 !> - standard measured/predicted meteorological values and TMRT 
 !> as input;
@@ -131,12 +136,12 @@
 
 !-- Type of input of the argument list
-    REAL(wp), INTENT ( IN )  :: tt2m   !< Local air temperature (°C)
+    REAL(wp), INTENT ( IN )  :: tt2m   !< Local air temperature (ï¿œC)
     REAL(wp), INTENT ( IN )  :: el2m   !< Local vapour pressure (hPa)
-    REAL(wp), INTENT ( IN )  :: tmrt   !< Local mean radiant temperature (°C)
+    REAL(wp), INTENT ( IN )  :: tmrt   !< Local mean radiant temperature (ï¿œC)
     REAL(wp), INTENT ( IN )  :: vau1m  !< Local wind velocitry (m/s)
     REAL(wp), INTENT ( IN )  :: pb     !< Local barometric air pressure (hPa)
 
 !-- Type of output of the argument list
-    REAL(wp), INTENT ( OUT ) :: pt_basic  !< Perceived temperature (°C)
+    REAL(wp), INTENT ( OUT ) :: pt_basic  !< Perceived temperature (ï¿œC)
     REAL(wp), INTENT ( OUT ) :: clo       !< Clothing index (dimensionless)
 
@@ -301,5 +306,5 @@
 ! ------------
 !> The SUBROUTINE calculates the saturation water vapour pressure 
-!> (hPa = hecto Pascal) for a given temperature tt (°C).
+!> (hPa = hecto Pascal) for a given temperature tt (ï¿œC).
 !> For example, tt can be the air temperature or the dew point temperature.
 !------------------------------------------------------------------------------!
@@ -308,5 +313,5 @@
     IMPLICIT NONE
 
-    REAL(wp), INTENT ( IN )  ::  tt   !< ambient air temperature (°C)
+    REAL(wp), INTENT ( IN )  ::  tt   !< ambient air temperature (ï¿œC)
     REAL(wp), INTENT ( OUT ) ::  p_st !< saturation water vapour pressure (hPa)
 
@@ -316,9 +321,9 @@
 
     IF ( tt < 0._wp ) THEN
-!--    tt  < 0 (°C): saturation water vapour pressure over ice
+!--    tt  < 0 (ï¿œC): saturation water vapour pressure over ice
        b = 17.84362_wp
        c = 245.425_wp
     ELSE
-!--    tt >= 0 (°C): saturation water vapour pressure over water
+!--    tt >= 0 (ï¿œC): saturation water vapour pressure over water
        b = 17.08085_wp
        c = 234.175_wp
@@ -346,6 +351,6 @@
 
 !-- Input variables of argument list:
-    REAL(wp), INTENT ( IN )  :: tt2m     !< Ambient temperature (°C)
-    REAL(wp), INTENT ( IN )  :: tmrt     !< Mean radiant temperature (°C)
+    REAL(wp), INTENT ( IN )  :: tt2m     !< Ambient temperature (ï¿œC)
+    REAL(wp), INTENT ( IN )  :: tmrt     !< Mean radiant temperature (ï¿œC)
     REAL(wp), INTENT ( IN )  :: el2m     !< Water vapour pressure (hPa)
     REAL(wp), INTENT ( IN )  :: vau1m    !< Wind speed (m/s) 1 m above ground
@@ -362,5 +367,5 @@
 ! Output variables of argument list:
     REAL(wp), INTENT ( OUT ) :: pmva     !< 0 (set to zero, because clo is evaluated for comfort)
-    REAL(wp), INTENT ( OUT ) :: top      !< Operative temperature (°C) at found root of Fanger's PMV
+    REAL(wp), INTENT ( OUT ) :: top      !< Operative temperature (ï¿œC) at found root of Fanger's PMV
     REAL(wp), INTENT ( OUT ) :: clo_res  !< Resulting clothing insulation value (clo)
     INTEGER(iwp), INTENT ( OUT ) :: nerr !< Error status / quality flag
@@ -488,5 +493,5 @@
     REAL(wp), INTENT ( IN ) ::  pmv   !< Fangers predicted mean vote (dimensionless)
     REAL(wp), INTENT ( IN ) ::  clo   !< clothing insulation index (clo)
-    REAL(wp), INTENT ( OUT ) ::  pt   !< pt (°C) corresponding to given PMV / clo
+    REAL(wp), INTENT ( OUT ) ::  pt   !< pt (ï¿œC) corresponding to given PMV / clo
 
     IF ( pmv <= -0.11_wp ) THEN
@@ -517,6 +522,6 @@
 
 !-- Input variables of argument list:
-    REAL(wp), INTENT ( IN ) ::  tt       !< Ambient air temperature (°C)
-    REAL(wp), INTENT ( IN ) ::  tmrt     !< Mean radiant temperature (°C)
+    REAL(wp), INTENT ( IN ) ::  tt       !< Ambient air temperature (ï¿œC)
+    REAL(wp), INTENT ( IN ) ::  tmrt     !< Mean radiant temperature (ï¿œC)
     REAL(wp), INTENT ( IN ) ::  pa       !< Water vapour pressure (hPa)
     REAL(wp), INTENT ( IN ) ::  pb       !< Barometric pressure (hPa) at site
@@ -530,5 +535,5 @@
 !            to Fanger corresponding to meteorological (tt,tmrt,pa,vau,pb) 
 !            and individual variables (clo, actlev, eta)
-    REAL(wp), INTENT ( OUT ) ::  top      !< operative temperature (°C)
+    REAL(wp), INTENT ( OUT ) ::  top      !< operative temperature (ï¿œC)
 
 !-- Internal variables
@@ -638,8 +643,8 @@
 !-- Input variables of argument list:
     REAL(wp),     INTENT ( IN )  :: pmva     !< Actual Predicted Mean Vote (PMV) according to Fanger
-    REAL(wp),     INTENT ( IN )  :: tt2m     !< Ambient temperature (°C) at screen level
+    REAL(wp),     INTENT ( IN )  :: tt2m     !< Ambient temperature (ï¿œC) at screen level
     REAL(wp),     INTENT ( IN )  :: el2m     !< Water vapour pressure (hPa) at screen level
     REAL(wp),     INTENT ( IN )  :: svp_tt   !< Saturation water vapour pressure (hPa) at tt2m
-    REAL(wp),     INTENT ( IN )  :: tmrt     !< Mean radiant temperature (°C) at screen level
+    REAL(wp),     INTENT ( IN )  :: tmrt     !< Mean radiant temperature (ï¿œC) at screen level
     REAL(wp),     INTENT ( IN )  :: vau1m    !< Wind speed (m/s) 1 m above ground
 
@@ -724,8 +729,8 @@
        nerr = -3_iwp
        IF ( el2m < p10 ) THEN
-!--       Due to conditions of regressíon: r.H. >= 5 %
+!--       Due to conditions of regressï¿œon: r.H. >= 5 %
           pa = p10
        ELSE
-!--       Due to conditions of regressíon: r.H. <= 95 %
+!--       Due to conditions of regressï¿œon: r.H. <= 95 %
           pa = p95
        ENDIF
@@ -733,5 +738,5 @@
     IF ( pa > 0._wp ) THEN
 !--    Natural logarithm of pa
-       apa = DLOG ( pa )
+       apa = LOG ( pa )
     ELSE
        apa = -5._wp
@@ -741,5 +746,5 @@
     pa_p50   = 0.5_wp * svp_tt
     IF ( pa_p50 > 0._wp .AND. pa > 0._wp ) THEN
-       dapa   = apa - DLOG ( pa_p50 )
+       dapa   = apa - LOG ( pa_p50 )
        pa_p50 = pa / pa_p50
     ELSE
@@ -765,5 +770,5 @@
        RETURN
     ENDIF
-    gew = DMOD ( pmv, 1._wp )
+    gew = MOD ( pmv, 1._wp )
     IF ( gew < 0._wp ) gew = 0._wp
     IF ( nreg > 5_iwp ) THEN
@@ -834,5 +839,5 @@
 !> velocity:
 !> - defined only under warm conditions: gtc based on 0.5 (clo) and
-!>                                       gtc > 16.826 (°C)
+!>                                       gtc > 16.826 (ï¿œC)
 !> - undefined: sultr_res set at +99.
 !------------------------------------------------------------------------------!
@@ -904,7 +909,7 @@
 !-- Type of input arguments
     REAL(wp), INTENT ( IN ) :: pmva      !< Fanger's classical predicted mean vote
-    REAL(wp), INTENT ( IN ) :: tt2m      !< Air temperature (°C) 2 m above ground
+    REAL(wp), INTENT ( IN ) :: tt2m      !< Air temperature (ï¿œC) 2 m above ground
     REAL(wp), INTENT ( IN ) :: vau1m     !< Relative wind velocity 1 m above ground (m/s)
-    REAL(wp), INTENT ( IN ) :: tmrt      !< Mean radiant temperature (°C)
+    REAL(wp), INTENT ( IN ) :: tmrt      !< Mean radiant temperature (ï¿œC)
 
 !-- Type of output argument
@@ -1036,5 +1041,5 @@
 !> (neither body cooling nor body heating). It is related to the Klima- 
 !> Michel activity level (134.682 W/m2). IREQ_neutral is only defined 
-!> for gt < 10 (°C)
+!> for gt < 10 (ï¿œC)
 !------------------------------------------------------------------------------!
  REAL(wp) FUNCTION ireq_neutral( gt, ireq_minimal, nerr )
@@ -1059,5 +1064,5 @@
     ireq_neutral = 1.62_wp - 0.0564_wp * top02
 
-!-- Regression only defined for gt <= 10 (°C)
+!-- Regression only defined for gt <= 10 (ï¿œC)
     IF ( ireq_neutral < 0.5_wp ) THEN
        IF ( ireq_neutral < 0.48_wp ) THEN
