Changeset 3796 for palm/trunk/SOURCE
- Timestamp:
- Mar 15, 2019 11:15:38 AM (6 years ago)
- File:
-
- 1 edited
-
palm/trunk/SOURCE/chemistry_model_mod.f90 (modified) (14 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/chemistry_model_mod.f90
r3784 r3796 26 26 ! Former revisions: 27 27 ! ----------------- 28 ! $Id$ 28 ! $Id: chemistry_model_mod.f90 3784 2019-03-05 14:16:20Z banzhafs 29 ! Unused variables removed/taken care of 30 ! 31 ! 32 ! 3784 2019-03-05 14:16:20Z forkel 29 33 ! Replaced READ from unit 10 by CALL get_mechanismname also in chem_header 30 ! 31 ! 34 ! 35 ! 32 36 ! 3783 2019-03-05 13:23:50Z forkel 33 37 ! Removed forgotte write statements an some unused variables (did not touch the … … 578 582 END INTERFACE rc_rctot 579 583 580 INTERFACE rc_comp_point_rc_eff581 MODULE PROCEDURE rc_comp_point_rc_eff582 END INTERFACE rc_comp_point_rc_eff584 ! INTERFACE rc_comp_point_rc_eff 585 ! MODULE PROCEDURE rc_comp_point_rc_eff 586 ! END INTERFACE rc_comp_point_rc_eff 583 587 584 588 INTERFACE drydepo_aero_zhang_vd … … 737 741 738 742 USE control_parameters, & 739 ONLY: air_chemistry, bc_radiation_l, bc_radiation_n, bc_radiation_r, &740 bc_radiation_s 743 ONLY: bc_radiation_l, bc_radiation_n, bc_radiation_r, bc_radiation_s 744 741 745 USE indices, & 742 746 ONLY: nxl, nxr, nzt … … 1045 1049 CHARACTER(len=16) :: spec_name 1046 1050 1051 ! 1052 !-- Next statement is to avoid compiler warnings about unused variables 1053 IF ( ( i + ilen + k ) > 0 .OR. var(1:1) == ' ' ) CONTINUE 1054 1047 1055 unit = 'illegal' 1048 1056 … … 1278 1286 INTEGER(iwp) :: k !< grid index along z-direction 1279 1287 INTEGER(iwp) :: char_len !< length of a character string 1288 1289 ! 1290 !-- Next statement is to avoid compiler warnings about unused variables 1291 IF ( mode(1:1) == ' ' .OR. two_d ) CONTINUE 1292 1280 1293 found = .FALSE. 1281 1294 char_len = LEN_TRIM(variable) … … 1751 1764 INTEGER(iwp) :: j !< running index y dimension 1752 1765 INTEGER(iwp) :: n !< running index for chemical species 1766 1767 ! 1768 !-- Next statement is to avoid compiler warning about unused variables 1769 IF ( ( ilu_arable + ilu_coniferous_forest + ilu_deciduous_forest + ilu_mediterrean_scrub + & 1770 ilu_permanent_crops + ilu_savanna + ilu_semi_natural_veg + ilu_tropical_forest + & 1771 ilu_urban ) == 0 ) CONTINUE 1772 1773 1753 1774 1754 1775 IF ( do_emis ) CALL chem_emissions_init … … 2786 2807 INTEGER(iwp) :: k !< matching k to surface m at i,j 2787 2808 INTEGER(iwp) :: lsp !< running index for chem spcs. 2788 INTEGER(iwp) :: lu !< running index for landuse classes2809 ! INTEGER(iwp) :: lu !< running index for landuse classes 2789 2810 INTEGER(iwp) :: lu_palm !< index of PALM LSM vegetation_type at current surface element 2790 2811 INTEGER(iwp) :: lup_palm !< index of PALM LSM pavement_type at current surface element … … 4496 4517 REAL(wp), INTENT(OUT) :: ccomp_tot !< total compensation point (ug/m3) 4497 4518 !< [= 0 for species that don't have a compensation 4498 4499 4519 ! 4500 4520 !-- Local variables: … … 4524 4544 LOGICAL :: sai_present !< vegetation is present for current land use type 4525 4545 4526 REAL(wp) :: laimax !< maximum leaf area index (-)4546 ! REAL(wp) :: laimax !< maximum leaf area index (-) 4527 4547 REAL(wp) :: gw !< external leaf conductance (m/s) 4528 4548 REAL(wp) :: gstom !< stomatal conductance (m/s) … … 4533 4553 REAL(wp) :: csoil !< soil compensation point (ug/m3) 4534 4554 4555 4556 ! 4557 !-- Next statement is just to avoid compiler warning about unused variable 4558 IF ( day_of_year == 0 .OR. ( catm + lat + ra + rb ) > 0.0_wp ) CONTINUE 4535 4559 4536 4560 … … 4661 4685 !< = 1 -> constant Rc 4662 4686 !< = 2 -> temperature dependent Rc 4687 4688 ! 4689 !-- Next line is to avoid compiler warning about unused variable 4690 IF ( icmp == 0 ) CONTINUE 4663 4691 4664 4692 ! … … 4950 4978 REAL(wp), PARAMETER :: dO3 = 0.13e-4 !< diffusion coefficient of ozon (m2/s) 4951 4979 4980 ! 4981 !-- Next line is to avoid compiler warning about unused variables 4982 IF ( icmp == 0 ) CONTINUE 4952 4983 4953 4984 SELECT CASE( TRIM(compnam) ) … … 5526 5557 ! ------------------------------------------------------------------------------------------- 5527 5558 5528 SUBROUTINE rc_comp_point_rc_eff( ccomp_tot, catm, ra, rb, rc_tot, rc_eff )5529 5530 IMPLICIT NONE5531 5532 !5533 !-- Input/output variables:5534 REAL(wp), INTENT(IN) :: ccomp_tot !< total compensation point (weighed average of separate compensation points) (ug/m3)5535 REAL(wp), INTENT(IN) :: catm !< atmospheric concentration (ug/m3)5536 REAL(wp), INTENT(IN) :: ra !< aerodynamic resistance (s/m)5537 REAL(wp), INTENT(IN) :: rb !< boundary layer resistance (s/m)5538 REAL(wp), INTENT(IN) :: rc_tot !< total canopy resistance (s/m)5539 5540 REAL(wp), INTENT(OUT) :: rc_eff !< effective total canopy resistance (s/m)5541 5542 !5543 !-- Compute effective resistance:5544 IF ( ccomp_tot == 0.0_wp ) THEN5545 !5546 !-- trace with no compensiation point ( or compensation point equal to zero)5547 rc_eff = rc_tot5548 5549 ELSE IF ( ccomp_tot > 0.0_wp .AND. ( abs( catm - ccomp_tot ) < 1.e-8 ) ) THEN5550 !5551 !--surface concentration (almost) equal to atmospheric concentration5552 !-- no exchange between surface and atmosphere, infinite RC --> vd=05553 rc_eff = 9999999999.0_wp5554 5555 ELSE IF ( ccomp_tot > 0.0_wp ) THEN5556 !5557 !-- compensation point available, calculate effective resistance5558 rc_eff = ( ( ra + rb ) * ccomp_tot + rc_tot * catm ) / ( catm - ccomp_tot )5559 5560 ELSE5561 rc_eff = -999.0_wp5562 message_string = 'This should not be possible, check ccomp_tot'5563 CALL message( 'rc_comp_point_rc_eff', 'CM0461', 1, 2, 0, 6, 0 )5564 ENDIF5565 5566 RETURN5567 5568 END SUBROUTINE rc_comp_point_rc_eff5559 ! SUBROUTINE rc_comp_point_rc_eff( ccomp_tot, catm, ra, rb, rc_tot, rc_eff ) 5560 ! 5561 ! IMPLICIT NONE 5562 ! 5563 ! ! 5564 ! !-- Input/output variables: 5565 ! REAL(wp), INTENT(IN) :: ccomp_tot !< total compensation point (weighed average of separate compensation points) (ug/m3) 5566 ! REAL(wp), INTENT(IN) :: catm !< atmospheric concentration (ug/m3) 5567 ! REAL(wp), INTENT(IN) :: ra !< aerodynamic resistance (s/m) 5568 ! REAL(wp), INTENT(IN) :: rb !< boundary layer resistance (s/m) 5569 ! REAL(wp), INTENT(IN) :: rc_tot !< total canopy resistance (s/m) 5570 ! 5571 ! REAL(wp), INTENT(OUT) :: rc_eff !< effective total canopy resistance (s/m) 5572 ! 5573 ! ! 5574 ! !-- Compute effective resistance: 5575 ! IF ( ccomp_tot == 0.0_wp ) THEN 5576 ! ! 5577 ! !-- trace with no compensiation point ( or compensation point equal to zero) 5578 ! rc_eff = rc_tot 5579 ! 5580 ! ELSE IF ( ccomp_tot > 0.0_wp .AND. ( abs( catm - ccomp_tot ) < 1.e-8 ) ) THEN 5581 ! ! 5582 ! !--surface concentration (almost) equal to atmospheric concentration 5583 ! !-- no exchange between surface and atmosphere, infinite RC --> vd=0 5584 ! rc_eff = 9999999999.0_wp 5585 ! 5586 ! ELSE IF ( ccomp_tot > 0.0_wp ) THEN 5587 ! ! 5588 ! !-- compensation point available, calculate effective resistance 5589 ! rc_eff = ( ( ra + rb ) * ccomp_tot + rc_tot * catm ) / ( catm - ccomp_tot ) 5590 ! 5591 ! ELSE 5592 ! rc_eff = -999.0_wp 5593 ! message_string = 'This should not be possible, check ccomp_tot' 5594 ! CALL message( 'rc_comp_point_rc_eff', 'CM0461', 1, 2, 0, 6, 0 ) 5595 ! ENDIF 5596 ! 5597 ! RETURN 5598 ! 5599 ! END SUBROUTINE rc_comp_point_rc_eff 5569 5600 5570 5601 … … 5758 5789 REAL(wp), PARAMETER :: kappa_stab = 0.35 !< von Karman constant 5759 5790 5791 ! 5792 !-- Next line is to avoid compiler warning about unused variable 5793 IF ( is_water .OR. ( z0h + kappa_stab ) > 0.0_wp ) CONTINUE 5760 5794 5761 5795 !
Note: See TracChangeset
for help on using the changeset viewer.
