Ignore:
Timestamp:
Dec 11, 2009 2:15:58 PM (14 years ago)
Author:
heinze
Message:

Large scale vertical motion (subsidence/ascent) can be applied to the prognostic equation for the potential temperature

File:
1 edited

Legend:

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

    r410 r411  
    44! Current revisions:
    55! -----------------
     6! +large scale vertical motion (subsidence/ascent)
    67!
    78!
     
    8485! end of the run, then writing additional information about CPU-usage on file
    8586! header.
    86 !------------------------------------------------------------------------------!
     87!-----------------------------------------------------------------------------!
    8788
    8889    USE arrays_3d
     
    9697    USE particle_attributes
    9798    USE pegrid
     99    USE subsidence_mod
    98100    USE spectrum
    99101
     
    323325          WRITE ( io, 152 )  dpdxy, dp_level_b, '.'
    324326       ENDIF
     327    ENDIF
     328    IF ( large_scale_subsidence )  THEN
     329        WRITE ( io, 153 )
     330        WRITE ( io, 154 )
    325331    ENDIF
    326332    WRITE ( io, 99 )
     
    13571363
    13581364       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 ), &
    13591396                          TRIM( gradients ), TRIM( slices )
    13601397    ENDIF
     
    15681605           /'     ',2(1X,E12.5),'Pa/m in x/y direction', &
    15691606           /'     starting from dp_level_b =', F8.3, 'm', A /)
     1607153 FORMAT (' --> Large-scale vertical motion is used in the ', &
     1608                  'prognostic equation for')
     1609154 FORMAT ('     the potential temperature')
    15701610200 FORMAT (//' Run time and time step information:'/ &
    15711611             ' ----------------------------------'/)
     
    18141854            '       Gradient:   ',A,'  psu/100m'/ &
    18151855            '       Gridpoint:  ',A)
     1856426 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)
    18161861450 FORMAT (//' LES / Turbulence quantities:'/ &
    18171862              ' ---------------------------'/)
Note: See TracChangeset for help on using the changeset viewer.