- Timestamp:
- Mar 5, 2019 1:23:50 PM (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/chemistry_model_mod.f90
r3780 r3783 27 27 ! ----------------- 28 28 ! $Id$ 29 ! Removed forgotte write statements an some unused variables (did not touch the 30 ! parts related to deposition) 31 ! 32 ! 33 ! 3780 2019-03-05 11:19:45Z forkel 29 34 ! Removed READ from unit 10, added CALL get_mechanismname 30 35 ! … … 731 736 bc_radiation_s 732 737 USE indices, & 733 ONLY: nxl, nxr, nxlg, nxrg, nyng, nysg,nzt738 ONLY: nxl, nxr, nzt 734 739 735 740 USE arrays_3d, & … … 747 752 INTEGER(iwp) :: m !< running index surface elements. 748 753 INTEGER(iwp) :: lsp !< running index for chem spcs. 749 INTEGER(iwp) :: lph !< running index for photolysis frequencies750 754 751 755 … … 869 873 SUBROUTINE chem_boundary_conds_decycle( cs_3d, cs_pr_init ) 870 874 871 USE pegrid, &872 ONLY: myid873 874 875 IMPLICIT NONE 875 876 INTEGER(iwp) :: boundary !< … … 1027 1028 SUBROUTINE chem_check_data_output( var, unit, i, ilen, k ) 1028 1029 1029 1030 USE control_parameters, &1031 ONLY: data_output, message_string1032 1030 1033 1031 IMPLICIT NONE … … 1186 1184 !-- check for chemical mechanism used 1187 1185 CALL get_mechanismname 1188 WRITE(06,*) 'cs_mech ',cs_mech1189 WRITE(06,*) 'chem_mechanism ', chem_mechanism1190 1186 IF (chem_mechanism /= trim(cs_mech) ) THEN 1191 1187 message_string = 'Incorrect chemistry mechanism selected, check spelling in namelist and/or chem_gasphase_mod' … … 1257 1253 USE kinds 1258 1254 1259 USE pegrid, &1260 ONLY: myid, threads_per_task1261 1262 1255 IMPLICIT NONE 1263 1256 … … 1280 1273 INTEGER(iwp) :: j !< grid index along y-direction 1281 1274 INTEGER(iwp) :: k !< grid index along z-direction 1282 INTEGER(iwp) :: m !< running index surface elements1283 1275 INTEGER(iwp) :: char_len !< length of a character string 1284 1276 found = .FALSE. … … 1460 1452 USE indices 1461 1453 USE kinds 1462 USE pegrid, &1463 ONLY: myid, threads_per_task1464 1454 USE surface_mod, & 1465 1455 ONLY: get_topography_top_index_ji … … 1793 1783 SUBROUTINE chem_init_internal 1794 1784 1795 USE control_parameters, &1796 ONLY: message_string, io_blocks, io_group, turbulent_inflow1797 USE arrays_3d, &1798 ONLY: mean_inflow_profiles1799 1785 USE pegrid 1800 1786 … … 2063 2049 INTEGER(iwp) :: lsp !< running index for chem spcs. 2064 2050 INTEGER(iwp) :: lph !< running index for photolysis frequencies 2065 INTEGER :: k2066 INTEGER :: m2067 INTEGER :: istatf2068 2051 INTEGER, DIMENSION(20) :: istatus 2069 2052 REAL(kind=wp), DIMENSION(nzb+1:nzt,nspec) :: tmp_conc … … 2083 2066 ! REAL(wp), PARAMETER :: xm_air = 28.96_wp !< Mole mass of dry air 2084 2067 ! REAL(wp), PARAMETER :: xm_h2o = 18.01528_wp !< Mole mass of water vapor 2085 REAL(wp), PARAMETER :: pref_i = 1._wp / 100000.0_wp !< inverse reference pressure (1/Pa)2086 2068 REAL(wp), PARAMETER :: t_std = 273.15_wp !< standard pressure (Pa) 2087 2069 REAL(wp), PARAMETER :: p_std = 101325.0_wp !< standard pressure (Pa) … … 2182 2164 2183 2165 CHARACTER (LEN=80) :: line !< dummy string that contains the current line of the parameter file 2184 CHARACTER (LEN=3) :: cs_prefix2185 2166 2186 2167 REAL(wp), DIMENSION(nmaxfixsteps) :: my_steps !< List of fixed timesteps my_step(1) = 0.0 automatic stepping 2187 2168 INTEGER(iwp) :: i !< 2188 INTEGER(iwp) :: j !<2189 2169 INTEGER(iwp) :: max_pr_cs_tmp !< 2190 2170 … … 2242 2222 !-- Read chem namelist 2243 2223 2244 INTEGER :: ier2245 CHARACTER(LEN=64) :: text2246 2224 CHARACTER(LEN=8) :: solver_type 2247 2225 … … 2677 2655 INTEGER(iwp) :: sr !< statistical region 2678 2656 INTEGER(iwp) :: tn !< thread number 2679 INTEGER(iwp) :: n !<2680 INTEGER(iwp) :: m !<2681 2657 INTEGER(iwp) :: lpr !< running index chem spcs 2682 2658
Note: See TracChangeset
for help on using the changeset viewer.