Changeset 411 for palm/trunk/SOURCE/header.f90
- Timestamp:
- Dec 11, 2009 2:15:58 PM (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/header.f90
r410 r411 4 4 ! Current revisions: 5 5 ! ----------------- 6 ! +large scale vertical motion (subsidence/ascent) 6 7 ! 7 8 ! … … 84 85 ! end of the run, then writing additional information about CPU-usage on file 85 86 ! header. 86 !----------------------------------------------------------------------------- -!87 !-----------------------------------------------------------------------------! 87 88 88 89 USE arrays_3d … … 96 97 USE particle_attributes 97 98 USE pegrid 99 USE subsidence_mod 98 100 USE spectrum 99 101 … … 323 325 WRITE ( io, 152 ) dpdxy, dp_level_b, '.' 324 326 ENDIF 327 ENDIF 328 IF ( large_scale_subsidence ) THEN 329 WRITE ( io, 153 ) 330 WRITE ( io, 154 ) 325 331 ENDIF 326 332 WRITE ( io, 99 ) … … 1357 1363 1358 1364 WRITE ( io, 425 ) TRIM( coordinates ), TRIM( temperatures ), & 1365 TRIM( gradients ), TRIM( slices ) 1366 ENDIF 1367 1368 ! 1369 !-- Profile for the large scale vertial velocity 1370 !-- Building output strings, starting with surface value 1371 IF ( large_scale_subsidence ) THEN 1372 temperatures = ' 0.0' 1373 gradients = '------' 1374 slices = ' 0' 1375 coordinates = ' 0.0' 1376 i = 1 1377 DO WHILE ( ws_vertical_gradient_level_ind(i) /= -9999 ) 1378 1379 WRITE (coor_chr,'(E10.2,7X)') & 1380 w_subs(ws_vertical_gradient_level_ind(i)) 1381 temperatures = TRIM( temperatures ) // ' ' // TRIM( coor_chr ) 1382 1383 WRITE (coor_chr,'(E10.2,7X)') ws_vertical_gradient(i) 1384 gradients = TRIM( gradients ) // ' ' // TRIM( coor_chr ) 1385 1386 WRITE (coor_chr,'(I10,7X)') ws_vertical_gradient_level_ind(i) 1387 slices = TRIM( slices ) // ' ' // TRIM( coor_chr ) 1388 1389 WRITE (coor_chr,'(F10.2,7X)') ws_vertical_gradient_level(i) 1390 coordinates = TRIM( coordinates ) // ' ' // TRIM( coor_chr ) 1391 1392 i = i + 1 1393 ENDDO 1394 1395 WRITE ( io, 426 ) TRIM( coordinates ), TRIM( temperatures ), & 1359 1396 TRIM( gradients ), TRIM( slices ) 1360 1397 ENDIF … … 1568 1605 /' ',2(1X,E12.5),'Pa/m in x/y direction', & 1569 1606 /' starting from dp_level_b =', F8.3, 'm', A /) 1607 153 FORMAT (' --> Large-scale vertical motion is used in the ', & 1608 'prognostic equation for') 1609 154 FORMAT (' the potential temperature') 1570 1610 200 FORMAT (//' Run time and time step information:'/ & 1571 1611 ' ----------------------------------'/) … … 1814 1854 ' Gradient: ',A,' psu/100m'/ & 1815 1855 ' Gridpoint: ',A) 1856 426 FORMAT (/' Characteristic levels of the subsidence/ascent profile:'// & 1857 ' Height: ',A,' m'/ & 1858 ' w_subs: ',A,' m/s'/ & 1859 ' Gradient: ',A,' (m/s)/100m'/ & 1860 ' Gridpoint: ',A) 1816 1861 450 FORMAT (//' LES / Turbulence quantities:'/ & 1817 1862 ' ---------------------------'/)
Note: See TracChangeset
for help on using the changeset viewer.