Ignore:
Timestamp:
Oct 24, 2018 6:39:32 PM (5 years ago)
Author:
gronemeier
Message:

new surface-data output; renamed output variables (pt to theta, rho_air to rho, rho_ocean to rho_sea_water)

File:
1 edited

Legend:

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

    r3294 r3421  
    2525! -----------------
    2626! $Id$
     27! Renamed output variables
     28!
     29! 3419 2018-10-24 17:27:31Z gronemeier
    2730! ocean renamed ocean_mode
    2831!
     
    316319                unit = 'm2/s2'
    317320
    318              CASE ( 'lpt' )
     321             CASE ( 'thetal' )
    319322                IF ( .NOT. bulk_cloud_model )  THEN
    320323                   WRITE ( message_string, * ) 'output of "', TRIM( var ),     &
     
    358361                IF ( TRIM( var ) == 'pr' )  unit = 'm'
    359362
    360              CASE ( 'q', 'vpt' )
     363             CASE ( 'q', 'thetav' )
    361364                IF ( .NOT. humidity )  THEN
    362365                   WRITE ( message_string, * ) 'output of "', TRIM( var ),     &
     
    365368                ENDIF
    366369                IF ( TRIM( var ) == 'q'   )  unit = 'kg/kg'
    367                 IF ( TRIM( var ) == 'vpt' )  unit = 'K'
     370                IF ( TRIM( var ) == 'thetav' )  unit = 'K'
    368371
    369372             CASE ( 'qc' )
     
    414417                unit = 'kg/kg'
    415418
    416              CASE ( 'rho_ocean' )
     419             CASE ( 'rho_sea_water' )
    417420                IF ( .NOT. ocean_mode )  THEN
    418421                   WRITE ( message_string, * ) 'output of "', TRIM( var ),     &
     
    438441                unit = 'psu'
    439442
    440              CASE ( 'u*', 't*', 'lwp*', 'pra*', 'prr*', 'z0*', 'z0h*' )
     443             CASE ( 'us*', 't*', 'lwp*', 'pra*', 'prr*', 'z0*', 'z0h*' )
    441444                WRITE ( message_string, * ) 'illegal value for data_',         &
    442445                     'output: "', TRIM( var ), '" is only allowed',            &
     
    444447                CALL message( 'init_masks', 'PA0111', 1, 2, 0, 6, 0 )
    445448
    446              CASE ( 'p', 'pt', 'u', 'v', 'w' )
     449             CASE ( 'p', 'theta', 'u', 'v', 'w' )
    447450                IF ( TRIM( var ) == 'p'  )  unit = 'Pa'
    448                 IF ( TRIM( var ) == 'pt' )  unit = 'K'
     451                IF ( TRIM( var ) == 'theta' )  unit = 'K'
    449452                IF ( TRIM( var ) == 'u'  )  unit = 'm/s'
    450453                IF ( TRIM( var ) == 'v'  )  unit = 'm/s'
Note: See TracChangeset for help on using the changeset viewer.