Changeset 96 for palm/trunk/SOURCE/sum_up_3d_data.f90
- Timestamp:
- Jun 4, 2007 8:07:41 AM (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/sum_up_3d_data.f90
r77 r96 4 4 ! Actual revisions: 5 5 ! ----------------- 6 ! 6 ! +sum-up of density and salinity 7 7 ! 8 8 ! Former revisions: … … 132 132 qv_av = 0.0 133 133 134 CASE ( 'rho' ) 135 IF ( .NOT. ALLOCATED( rho_av ) ) THEN 136 ALLOCATE( rho_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) ) 137 ENDIF 138 rho_av = 0.0 139 134 140 CASE ( 's' ) 135 141 IF ( .NOT. ALLOCATED( s_av ) ) THEN … … 137 143 ENDIF 138 144 s_av = 0.0 145 146 CASE ( 'sa' ) 147 IF ( .NOT. ALLOCATED( sa_av ) ) THEN 148 ALLOCATE( sa_av(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) ) 149 ENDIF 150 sa_av = 0.0 139 151 140 152 CASE ( 't*' ) … … 337 349 ENDDO 338 350 351 CASE ( 'rho' ) 352 DO i = nxl-1, nxr+1 353 DO j = nys-1, nyn+1 354 DO k = nzb, nzt+1 355 rho_av(k,j,i) = rho_av(k,j,i) + rho(k,j,i) 356 ENDDO 357 ENDDO 358 ENDDO 359 339 360 CASE ( 's' ) 340 361 DO i = nxl-1, nxr+1 … … 342 363 DO k = nzb, nzt+1 343 364 s_av(k,j,i) = s_av(k,j,i) + q(k,j,i) 365 ENDDO 366 ENDDO 367 ENDDO 368 369 CASE ( 'sa' ) 370 DO i = nxl-1, nxr+1 371 DO j = nys-1, nyn+1 372 DO k = nzb, nzt+1 373 sa_av(k,j,i) = sa_av(k,j,i) + sa(k,j,i) 344 374 ENDDO 345 375 ENDDO
Note: See TracChangeset
for help on using the changeset viewer.