Ignore:
Timestamp:
Jul 12, 2016 4:34:24 PM (5 years ago)
Author:
suehring
Message:

Separate balance equations for humidity and passive_scalar

File:
1 edited

Legend:

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

    r1874 r1960  
    1919! Current revisions:
    2020! -----------------
    21 !
     21! Separate humidity and passive scalar
    2222!
    2323! Former revisions:
     
    112112
    113113    USE arrays_3d,                                                             &
    114         ONLY:  p, pt, q, ql, ts, u, us, v, w, zu
     114        ONLY:  p, pt, q, ql, s, ts, u, us, v, w, zu
    115115       
    116116    USE cloud_parameters,                                                      &
     
    336336
    337337                CASE ( 'q', 'q_xy', 'q_xz', 'q_yz' )
    338                    IF ( humidity  .OR.  passive_scalar )  THEN
     338                   IF ( humidity )  THEN
    339339                      DO  i = nxl_dvrp, nxr_dvrp+1
    340340                         DO  j = nys_dvrp, nyn_dvrp+1
     
    345345                      ENDDO           
    346346                   ELSE                   
    347                       message_string = 'if humidity/passive_scalar = '    //   &
     347                      message_string = 'if humidity = '    //                  &
    348348                            'FALSE output of ' // TRIM( output_variable ) //   &
    349349                            'is not provided'
     
    366366                                  'is not provided'
    367367                      CALL message( 'data_output_dvrp', 'PA0184',&
     368                                                                 0, 0, 0, 6, 0 )
     369                   ENDIF
     370
     371                CASE ( 's', 's_xy', 's_xz', 's_yz' )
     372                   IF ( passive_scalar )  THEN
     373                      DO  i = nxl_dvrp, nxr_dvrp+1
     374                         DO  j = nys_dvrp, nyn_dvrp+1
     375                            DO  k = nzb, nz_do3d
     376                               local_pf(i,j,k) = s(k,j,i)
     377                            ENDDO
     378                         ENDDO
     379                      ENDDO           
     380                   ELSE                   
     381                      message_string = 'if passive_scalar = '    //            &
     382                            'FALSE output of ' // TRIM( output_variable ) //   &
     383                            'is not provided'
     384                      CALL message( 'data_output_dvrp', 'PA0183',&
    368385                                                                 0, 0, 0, 6, 0 )
    369386                   ENDIF
Note: See TracChangeset for help on using the changeset viewer.