- Timestamp:
- Apr 21, 2016 11:11:40 AM (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/lpm_droplet_collision.f90
r1861 r1884 19 19 ! Current revisions: 20 20 ! ------------------ 21 ! 21 ! Conservation of should only be checked if collisions took place. 22 22 ! 23 23 ! Former revisions: … … 505 505 ENDIF 506 506 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 507 520 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 523 522 CALL cpu_log( log_point_s(43), 'lpm_droplet_coll', 'stop' ) 524 523
Note: See TracChangeset
for help on using the changeset viewer.