Ignore:
Timestamp:
Oct 29, 2013 10:11:53 AM (8 years ago)
Author:
heinze
Message:

routines for nudging and large scale forcing from external file added

File:
1 edited

Legend:

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

    r1217 r1239  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! output for nudging + large scale forcing from external file
    2323!
    2424! Former revisions:
     
    428428        WRITE ( io, 154 )
    429429    ENDIF
     430    IF ( nudging )  THEN
     431        WRITE ( io, 155 )
     432        WRITE ( io, 156 )
     433    ENDIF
     434    IF ( large_scale_forcing )  THEN
     435        WRITE ( io, 157 )
     436        WRITE ( io, 158 )
     437    ENDIF
    430438    WRITE ( io, 99 )
    431439
     
    706714       WRITE ( io, 303 )
    707715       IF ( constant_heatflux )  THEN
    708           WRITE ( io, 306 )  surface_heatflux
     716          IF ( large_scale_forcing .AND. lsf_surf ) THEN
     717             WRITE ( io, 306 )  shf(0,0)
     718          ELSE
     719             WRITE ( io, 306 )  surface_heatflux
     720          ENDIF
    709721          IF ( random_heatflux )  WRITE ( io, 307 )
    710722       ENDIF
    711723       IF ( humidity  .AND.  constant_waterflux )  THEN
    712           WRITE ( io, 311 ) surface_waterflux
     724          IF ( large_scale_forcing .AND. lsf_surf ) THEN
     725             WRITE ( io, 311 ) qsws(0,0)
     726          ELSE
     727             WRITE ( io, 311 ) surface_waterflux
     728          ENDIF
    713729       ENDIF
    714730       IF ( passive_scalar  .AND.  constant_waterflux )  THEN
     
    12881304    ENDDO
    12891305
    1290     WRITE ( io, 423 )  TRIM( coordinates ), TRIM( ugcomponent ), &
    1291                        TRIM( gradients ), TRIM( slices )
     1306    IF ( .NOT. large_scale_forcing ) THEN
     1307       WRITE ( io, 423 )  TRIM( coordinates ), TRIM( ugcomponent ), &
     1308                          TRIM( gradients ), TRIM( slices )
     1309    ELSE
     1310       WRITE ( io, 429 )
     1311    ENDIF
    12921312
    12931313!-- Profile of the geostrophic wind (component vg)
     
    13201340    ENDDO
    13211341
    1322     WRITE ( io, 424 )  TRIM( coordinates ), TRIM( vgcomponent ), &
    1323                        TRIM( gradients ), TRIM( slices )
     1342    IF ( .NOT. large_scale_forcing ) THEN
     1343       WRITE ( io, 424 )  TRIM( coordinates ), TRIM( vgcomponent ), &
     1344                          TRIM( gradients ), TRIM( slices )
     1345    ENDIF
    13241346
    13251347!
     
    13571379    ENDDO
    13581380
    1359     WRITE ( io, 420 )  TRIM( coordinates ), TRIM( temperatures ), &
    1360                        TRIM( gradients ), TRIM( slices )
     1381    IF ( .NOT. nudging ) THEN
     1382       WRITE ( io, 420 )  TRIM( coordinates ), TRIM( temperatures ), &
     1383                          TRIM( gradients ), TRIM( slices )
     1384    ELSE
     1385       WRITE ( io, 428 )
     1386    ENDIF
    13611387
    13621388!
     
    13921418
    13931419       IF ( humidity )  THEN
    1394           WRITE ( io, 421 )  TRIM( coordinates ), TRIM( temperatures ), &
    1395                              TRIM( gradients ), TRIM( slices )
     1420          IF ( .NOT. nudging ) THEN
     1421             WRITE ( io, 421 )  TRIM( coordinates ), TRIM( temperatures ), &
     1422                                TRIM( gradients ), TRIM( slices )
     1423          ENDIF
    13961424       ELSE
    13971425          WRITE ( io, 422 )  TRIM( coordinates ), TRIM( temperatures ), &
     
    17051733                  'prognostic equation for')
    17061734154 FORMAT ('     the potential temperature')
     1735155 FORMAT (' --> Nudging is used - initial profiles for u, v, pt and q ',&
     1736                  'correspond to the')
     1737156 FORMAT ('     first profiles in NUDGING_DATA')
     1738157 FORMAT (' --> Large scale forcing from external file (LSF_DATA) is used: ')
     1739158 FORMAT ('     prescribed surfaces fluxes and geostrophic wind components')
    17071740200 FORMAT (//' Run time and time step information:'/ &
    17081741             ' ----------------------------------'/)
     
    19571990427 FORMAT (/'    Initial wind profiles (u,v) are interpolated from given'// &
    19581991                  ' profiles')
     1992428 FORMAT (/'    Initial profiles (u, v, pt, q) are taken from file '/ &
     1993             '    NUDGING_DATA')
     1994429 FORMAT (/'    Geostrophic wind profiles (ug, vg) are are taken from file '/ &
     1995             '    LSF_DATA_DATA')
    19591996430 FORMAT (//' Cloud physics quantities / methods:'/ &
    19601997              ' ----------------------------------'/)
Note: See TracChangeset for help on using the changeset viewer.