Changeset 3597 for palm/trunk/SOURCE/surface_layer_fluxes_mod.f90

Timestamp:

Dec 4, 2018 8:40:18 AM (5 years ago)

Author:

maronga

Message:

revised calculation of near surface air potential temperature

File:

• palm/trunk/SOURCE/surface_layer_fluxes_mod.f90 (modified) (9 diffs)

palm/trunk/SOURCE/surface_layer_fluxes_mod.f90

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! Added routine for calculating near surface air potential temperature (moved
+! from urban_surface_mod)
+! 
+! 3547 2018-11-21 13:21:24Z suehring
 ! variables documented
 ! 
@@ -250 +254 @@
         ONLY:  air_chemistry, cloud_droplets,                                  &
                constant_heatflux, constant_scalarflux,                         &
-               constant_waterflux, coupling_mode, humidity, ibc_e_b,           &
-               ibc_pt_b, initializing_actions,                                 &
+               constant_waterflux, coupling_mode, do_output_at_2m, humidity,   &
+               ibc_e_b, ibc_pt_b, indoor_model, initializing_actions,          &
                intermediate_timestep_count, intermediate_timestep_count_max,   &
                land_surface, large_scale_forcing, lsf_surf, message_string,    &
@@ -413 +417 @@
              CALL calc_us
              CALL calc_surface_fluxes
+             
+             IF ( do_output_at_2m )  THEN
+                CALL calc_pt_near_surface ( '2m' )
+             ENDIF
           ENDIF
 !
@@ -423 +431 @@
              CALL calc_us
              CALL calc_surface_fluxes
+             IF ( do_output_at_2m )  THEN
+                CALL calc_pt_near_surface ( '2m' )
+             ENDIF
           ENDIF
 !
@@ -433 +444 @@
              CALL calc_us
              CALL calc_surface_fluxes
+             IF ( do_output_at_2m )  THEN
+                CALL calc_pt_near_surface ( '2m' )
+             ENDIF
+             IF ( indoor_model )  THEN
+                CALL calc_pt_near_surface ( '10cm' )
+             ENDIF
           ENDIF
 !
@@ -447 +464 @@
              CALL calc_scaling_parameters
              CALL calc_surface_fluxes
+             IF ( do_output_at_2m )  THEN
+                CALL calc_pt_near_surface ( '2m' )
+             ENDIF
           ENDIF
 !
@@ -457 +477 @@
              CALL calc_scaling_parameters
              CALL calc_surface_fluxes
+             IF ( do_output_at_2m )  THEN
+                CALL calc_pt_near_surface ( '2m' )
+             ENDIF
           ENDIF
 !
@@ -467 +490 @@
              CALL calc_scaling_parameters
              CALL calc_surface_fluxes
+             IF ( do_output_at_2m )  THEN
+                CALL calc_pt_near_surface ( '2m' )
+             ENDIF
+             IF ( indoor_model )  THEN
+                CALL calc_pt_near_surface ( '10cm' )
+             ENDIF
           ENDIF
 
@@ -2305 +2334 @@
     END SUBROUTINE calc_surface_fluxes
 
+    
+!------------------------------------------------------------------------------!
+! Description:
+! ------------
+!> Calculates temperature near surface (10 cm) for indoor model or 2 m 
+!> temperature for output
+!------------------------------------------------------------------------------!
+    SUBROUTINE calc_pt_near_surface ( z_char )
+
+       IMPLICIT NONE
+
+       CHARACTER (LEN = *), INTENT(IN)       :: z_char          !< string identifier to identify z level
+       INTEGER(iwp)                          :: i, j, k, l, m   !< running indices
+
+       
+       SELECT CASE ( z_char)
+            
+       
+          CASE ( '10cm' )
+
+             DO  m = 1, surf%ns
+
+                i = surf%i(m)
+                j = surf%j(m)
+                k = surf%k(m)
+
+                surf%pt_10cm(m) = surf%pt_surface(m) + surf%ts(m) / kappa  &
+                                   * ( log( 0.1_wp /  surf%z0h(m) )                    &
+                                  - psi_h( 0.1_wp / surf%ol(m) )                    &
+                                     + psi_h( surf%z0h(m) / surf%ol(m) ) )
+
+             ENDDO
+
+
+          CASE ( '2m' )
+      
+             DO  m = 1, surf%ns
+
+                i = surf%i(m)
+                j = surf%j(m)
+                k = surf%k(m)
+
+                surf%pt_2m(m) = surf%pt_surface(m) + surf%ts(m) / kappa  &
+                                   * ( log( 2.0_wp /  surf%z0h(m) )                    &
+                                     - psi_h( 2.0_wp / surf%ol(m) )                    &
+                                     + psi_h( surf%z0h(m) / surf%ol(m) ) )
+
+             ENDDO
+          
+       
+       END SELECT
+
+    END SUBROUTINE calc_pt_near_surface
+    
 
 !