Changeset 3933 for palm/trunk/SOURCE/land_surface_model_mod.f90

Timestamp:

Apr 25, 2019 12:33:20 PM (6 years ago)

Author:

kanani

Message:

Bugfixes and clean-up for output quantity theta_2m*

File:

• palm/trunk/SOURCE/land_surface_model_mod.f90 (modified) (5 diffs)

palm/trunk/SOURCE/land_surface_model_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Remove unused subroutine and allocation of pt_2m, this is done in surface_mod
+! now (surfaces%pt_2m)
+! 
+! 
 ! Changes related to global restructuring of location messages and introduction 
 ! of additional debug messages
@@ -1017 +1021 @@
 !-- Public functions
     PUBLIC lsm_boundary_condition, lsm_check_data_output,                      &
-           lsm_check_data_output_pr, lsm_calc_pt_near_surface,                 &
+           lsm_check_data_output_pr,                                           &
            lsm_check_parameters, lsm_define_netcdf_grid, lsm_3d_data_averaging,& 
            lsm_data_output_2d, lsm_data_output_3d, lsm_energy_balance,         &
@@ -1038 +1042 @@
     END INTERFACE lsm_boundary_condition
 
-    INTERFACE lsm_calc_pt_near_surface
-       MODULE PROCEDURE lsm_calc_pt_near_surface
-    END INTERFACE lsm_calc_pt_near_surface
-    
     INTERFACE lsm_check_data_output
        MODULE PROCEDURE lsm_check_data_output
@@ -5057 +5057 @@
        ALLOCATE ( surf_lsm_h%r_s(1:surf_lsm_h%ns)                 )
        ALLOCATE ( surf_lsm_h%r_canopy_min(1:surf_lsm_h%ns)        )
-       ALLOCATE ( surf_lsm_h%pt_2m(1:surf_lsm_h%ns)               )
        ALLOCATE ( surf_lsm_h%vegetation_surface(1:surf_lsm_h%ns)  )
        ALLOCATE ( surf_lsm_h%water_surface(1:surf_lsm_h%ns)       )
@@ -7322 +7321 @@
     END SUBROUTINE calc_z0_water_surface
 
-    
-!------------------------------------------------------------------------------!
-! Description:
-! ------------
-!> Calculates 2m temperature for data output at coarse resolution
-!------------------------------------------------------------------------------!
-    SUBROUTINE lsm_calc_pt_near_surface
-
-       IMPLICIT NONE
-
-       INTEGER(iwp)                          :: i, j, k, m   !< running indices
-
-
-       DO  m = 1, surf_lsm_h%ns
-
-          i = surf_lsm_h%i(m)
-          j = surf_lsm_h%j(m)
-          k = surf_lsm_h%k(m)
-
-          surf_lsm_h%pt_2m(m) = surf_lsm_h%pt_surface(m) + surf_lsm_h%ts(m) / kappa &
-                             * ( log( 2.0_wp /  surf_lsm_h%z0h(m) )                 &
-                               - psi_h( 2.0_wp / surf_lsm_h%ol(m) )                 &
-                               + psi_h( surf_lsm_h%z0h(m) / surf_lsm_h%ol(m) ) )
-
-       ENDDO
-
-    END SUBROUTINE lsm_calc_pt_near_surface
-
-    
-    
 !
 !-- Integrated stability function for heat and moisture