Changeset 4090 for palm/trunk/SOURCE
- Timestamp:
- Jul 11, 2019 3:06:47 PM (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/init_3d_model.f90
r4088 r4090 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Unused variables removed 28 ! 29 ! 4088 2019-07-11 13:57:56Z Giersch 27 30 ! Pressure and density profile calculations revised using basic functions 28 31 ! … … 614 617 615 618 USE basic_constants_and_equations_mod, & 616 ONLY: c_p, g, l_v, pi, r_d, exner_function, exner_function_invers,&619 ONLY: c_p, g, l_v, pi, exner_function, exner_function_invers, & 617 620 ideal_gas_law_rho, ideal_gas_law_rho_pt, barometric_formula 618 621 … … 711 714 INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE :: ngp_2dh_outer_l !< number of horizontal non-wall bounded grid points on subdomain 712 715 INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE :: ngp_2dh_s_inner_l !< number of horizontal non-topography grid points on subdomain 713 714 REAL(wp) :: t_surface !< air temperature at the surface715 716 716 717 REAL(wp), DIMENSION(:), ALLOCATABLE :: init_l !< dummy array used for averaging 3D data to obtain inital profiles … … 852 853 !-- mainly for output purposes. This density do not need to be considered 853 854 !-- in the model's system of equations. 854 ! t_surface = pt_surface * ( surface_pressure / 1000.0_wp )**( r_d / c_p )855 855 IF ( TRIM( approximation ) == 'anelastic' ) THEN 856 856 DO k = nzb, nzt+1
Note: See TracChangeset
for help on using the changeset viewer.