Changeset 96 for palm/trunk/SOURCE/buoyancy.f90
 Timestamp:
 Jun 4, 2007 8:07:41 AM (14 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

palm/trunk/SOURCE/buoyancy.f90
r77 r96 4 4 ! Actual revisions: 5 5 !  6 ! 6 ! calc_mean_pt_profile renamed calc_mean_profile 7 7 ! 8 8 ! Former revisions: … … 29 29 30 30 PRIVATE 31 PUBLIC buoyancy, calc_mean_p t_profile31 PUBLIC buoyancy, calc_mean_profile 32 32 33 33 INTERFACE buoyancy … … 36 36 END INTERFACE buoyancy 37 37 38 INTERFACE calc_mean_p t_profile39 MODULE PROCEDURE calc_mean_p t_profile40 END INTERFACE calc_mean_p t_profile38 INTERFACE calc_mean_profile 39 MODULE PROCEDURE calc_mean_profile 40 END INTERFACE calc_mean_profile 41 41 42 42 CONTAINS … … 207 207 208 208 209 SUBROUTINE calc_mean_p t_profile( theta, pr )209 SUBROUTINE calc_mean_profile( var, pr ) 210 210 211 211 !! … … 224 224 225 225 INTEGER :: i, j, k, omp_get_thread_num, pr, tn 226 REAL, DIMENSION(:,:,:), POINTER :: theta227 228 ! 229 ! Computation of the horizontally averaged temperature profile, unless226 REAL, DIMENSION(:,:,:), POINTER :: var 227 228 ! 229 ! Computation of the horizontally averaged profile of variable var, unless 230 230 ! already done by the relevant call from flow_statistics. The calculation 231 231 ! is done only for the first respective intermediate timestep in order to 232 232 ! spare communication time and to produce identical model results with jobs 233 233 ! which are calling flow_statistics at different time intervals. 234 ! Although this calculation is not required for model runs with a slope,235 ! it is nevertheless also computed.236 234 IF ( .NOT. flow_statistics_called .AND. & 237 235 intermediate_timestep_count == 1 ) THEN 238 236 239 237 ! 240 ! Horizontal average of the potential temperature238 ! Horizontal average of variable var 241 239 tn = 0 ! Default thread number in case of one thread 242 240 !$OMP PARALLEL PRIVATE( i, j, k, tn ) … … 247 245 DO j = nys, nyn 248 246 DO k = nzb_s_outer(j,i), nzt+1 249 sums_l(k,pr,tn) = sums_l(k,pr,tn) + theta(k,j,i)247 sums_l(k,pr,tn) = sums_l(k,pr,tn) + var(k,j,i) 250 248 ENDDO 251 249 ENDDO … … 272 270 ENDIF 273 271 274 END SUBROUTINE calc_mean_p t_profile272 END SUBROUTINE calc_mean_profile 275 273 276 274 END MODULE buoyancy_mod
Note: See TracChangeset
for help on using the changeset viewer.