Ignore:
Timestamp:
Apr 21, 2016 11:11:40 AM (5 years ago)
Author:
hoffmann
Message:

conservation of mass only chekced if collisions take place

File:
1 edited

Legend:

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

    r1861 r1884  
    1919! Current revisions:
    2020! ------------------
    21 !
     21! Conservation of should only be checked if collisions took place.
    2222!
    2323! Former revisions:
     
    505505       ENDIF
    506506 
     507!
     508!--    Check if LWC is conserved during collision process
     509       IF ( ql_v(k,j,i) /= 0.0_wp )  THEN
     510          IF ( ql_vp(k,j,i) / ql_v(k,j,i) >= 1.0001_wp  .OR.                      &
     511               ql_vp(k,j,i) / ql_v(k,j,i) <= 0.9999_wp )  THEN
     512             WRITE( message_string, * ) ' LWC is not conserved during',           &
     513                                        ' collision! ',                           &
     514                                        ' LWC after condensation: ', ql_v(k,j,i), &
     515                                        ' LWC after collision: ', ql_vp(k,j,i)
     516             CALL message( 'lpm_droplet_collision', 'PA0040', 2, 2, -1, 6, 1 )
     517          ENDIF
     518       ENDIF
     519
    507520    ENDIF
    508    
    509 
    510 !
    511 !-- Check if LWC is conserved during collision process
    512     IF ( ql_v(k,j,i) /= 0.0_wp )  THEN
    513        IF ( ql_vp(k,j,i) / ql_v(k,j,i) >= 1.0001_wp  .OR.                      &
    514             ql_vp(k,j,i) / ql_v(k,j,i) <= 0.9999_wp )  THEN
    515           WRITE( message_string, * ) ' LWC is not conserved during',           &
    516                                      ' collision! ',                           &
    517                                      ' LWC after condensation: ', ql_v(k,j,i), &
    518                                      ' LWC after collision: ', ql_vp(k,j,i)
    519           CALL message( 'lpm_droplet_collision', 'PA0040', 2, 2, -1, 6, 1 )
    520        ENDIF
    521     ENDIF
    522 
     521 
    523522    CALL cpu_log( log_point_s(43), 'lpm_droplet_coll', 'stop' )
    524523
Note: See TracChangeset for help on using the changeset viewer.