Ignore:
Timestamp:
Apr 25, 2019 12:33:20 PM (2 years ago)
Author:
kanani
Message:

Bugfixes and clean-up for output quantity theta_2m*

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/land_surface_model_mod.f90

    r3885 r3933  
    2525! -----------------
    2626! $Id$
     27! Remove unused subroutine and allocation of pt_2m, this is done in surface_mod
     28! now (surfaces%pt_2m)
     29!
     30!
    2731! Changes related to global restructuring of location messages and introduction
    2832! of additional debug messages
     
    10171021!-- Public functions
    10181022    PUBLIC lsm_boundary_condition, lsm_check_data_output,                      &
    1019            lsm_check_data_output_pr, lsm_calc_pt_near_surface,                 &
     1023           lsm_check_data_output_pr,                                           &
    10201024           lsm_check_parameters, lsm_define_netcdf_grid, lsm_3d_data_averaging,&
    10211025           lsm_data_output_2d, lsm_data_output_3d, lsm_energy_balance,         &
     
    10381042    END INTERFACE lsm_boundary_condition
    10391043
    1040     INTERFACE lsm_calc_pt_near_surface
    1041        MODULE PROCEDURE lsm_calc_pt_near_surface
    1042     END INTERFACE lsm_calc_pt_near_surface
    1043    
    10441044    INTERFACE lsm_check_data_output
    10451045       MODULE PROCEDURE lsm_check_data_output
     
    50575057       ALLOCATE ( surf_lsm_h%r_s(1:surf_lsm_h%ns)                 )
    50585058       ALLOCATE ( surf_lsm_h%r_canopy_min(1:surf_lsm_h%ns)        )
    5059        ALLOCATE ( surf_lsm_h%pt_2m(1:surf_lsm_h%ns)               )
    50605059       ALLOCATE ( surf_lsm_h%vegetation_surface(1:surf_lsm_h%ns)  )
    50615060       ALLOCATE ( surf_lsm_h%water_surface(1:surf_lsm_h%ns)       )
     
    73227321    END SUBROUTINE calc_z0_water_surface
    73237322
    7324    
    7325 !------------------------------------------------------------------------------!
    7326 ! Description:
    7327 ! ------------
    7328 !> Calculates 2m temperature for data output at coarse resolution
    7329 !------------------------------------------------------------------------------!
    7330     SUBROUTINE lsm_calc_pt_near_surface
    7331 
    7332        IMPLICIT NONE
    7333 
    7334        INTEGER(iwp)                          :: i, j, k, m   !< running indices
    7335 
    7336 
    7337        DO  m = 1, surf_lsm_h%ns
    7338 
    7339           i = surf_lsm_h%i(m)
    7340           j = surf_lsm_h%j(m)
    7341           k = surf_lsm_h%k(m)
    7342 
    7343           surf_lsm_h%pt_2m(m) = surf_lsm_h%pt_surface(m) + surf_lsm_h%ts(m) / kappa &
    7344                              * ( log( 2.0_wp /  surf_lsm_h%z0h(m) )                 &
    7345                                - psi_h( 2.0_wp / surf_lsm_h%ol(m) )                 &
    7346                                + psi_h( surf_lsm_h%z0h(m) / surf_lsm_h%ol(m) ) )
    7347 
    7348        ENDDO
    7349 
    7350     END SUBROUTINE lsm_calc_pt_near_surface
    7351 
    7352    
    7353    
    73547323!
    73557324!-- Integrated stability function for heat and moisture
Note: See TracChangeset for help on using the changeset viewer.