Changeset 3395 for palm/trunk/UTIL/inifor/src/util.f90

Timestamp:

Oct 22, 2018 5:32:49 PM (7 years ago)

Author:

eckhard

Message:

inifor: Added computation of geostrophic winds from COSMO input

File:

• palm/trunk/UTIL/inifor/src/util.f90 (modified) (4 diffs)

palm/trunk/UTIL/inifor/src/util.f90

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! New routines for computing potential temperature and moist air density
+! Increased number of digits in real-to-str conversion
+! 
+! 3183 2018-07-27 14:25:55Z suehring
 ! Improved real-to-string conversion
 ! 
@@ -48 +52 @@
     USE defs,                                                                  &
         ONLY :  dp, PI, DATE, SNAME
+    USE types,                                                                 &
+        ONLY :  grid_definition
 
     IMPLICIT NONE
@@ -281 +287 @@
 
 
+!------------------------------------------------------------------------------!
+! Description:
+! ------------
+!> Converts the absolute temperature to the potential temperature in place using
+!> the identity a^b = e^(b ln(a)).
+!>
+!>     theta = T * (p_ref/p)^(R/c_p) = T * e^( R/c_p * ln(p_ref/p) )
+!------------------------------------------------------------------------------!
+    SUBROUTINE potential_temperature(t, p, p_ref, r, cp)
+       REAL(dp), DIMENSION(:,:,:), INTENT(INOUT) ::  t
+       REAL(dp), DIMENSION(:,:,:), INTENT(IN)    ::  p
+       REAL(dp), INTENT(IN)                      ::  p_ref, r, cp
+       REAL(dp)                                  ::  rcp
+
+       rcp = r/cp
+       t(:,:,:) =  t(:,:,:) * EXP( rcp * LOG(p_ref / p(:,:,:)) )
+
+    END SUBROUTINE potential_temperature
+
+
+!------------------------------------------------------------------------------!
+! Description:
+! ------------
+!> Compute the density in place of the given temperature (t_rho).
+!------------------------------------------------------------------------------!
+   SUBROUTINE moist_density(t_rho, p, qv, rd, rv)
+       REAL(dp), DIMENSION(:,:,:), INTENT(INOUT) ::  t_rho
+       REAL(dp), DIMENSION(:,:,:), INTENT(IN)    ::  p, qv
+       REAL(dp), INTENT(IN)                      ::  rd, rv
+
+       t_rho(:,:,:) = p(:,:,:) / (                                             &
+          (rv * qv(:,:,:) + rd * (1.0_dp - qv(:,:,:))) * t_rho(:,:,:)          &
+       )
+
+    END SUBROUTINE moist_density
+
+
     ! Convert a real number to a string in scientific notation
     ! showing four significant digits.
@@ -298 +341 @@
 
 
-    CHARACTER(LEN=12) FUNCTION real_to_str_f(val)
+    CHARACTER(LEN=16) FUNCTION real_to_str_f(val)
 
         REAL(dp), INTENT(IN) ::  val
 
-        WRITE(real_to_str_f, '(F12.4)') val
+        WRITE(real_to_str_f, '(F16.8)') val
         real_to_str_f = ADJUSTL(real_to_str_f)