Changeset 2932 for palm/trunk/SOURCE/chemistry_model_mod.f90
- Timestamp:
- Mar 26, 2018 9:39:22 AM (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/chemistry_model_mod.f90
r2894 r2932 27 27 ! ----------------- 28 28 ! $Id$ 29 ! renamed chemistry_par to chemistry_parameters 30 ! 31 ! 2894 2018-03-15 09:17:58Z Giersch 29 32 ! Calculations of the index range of the subdomain on file which overlaps with 30 33 ! the current subdomain are already done in read_restart_data_mod, … … 600 603 ! 601 604 !-- Set initial concentration of profiles prescribed by parameters cs_profile 602 !-- and cs_heights in the namelist &chemistry_par 605 !-- and cs_heights in the namelist &chemistry_parameters 603 606 !-- (todo (FK): chem_init_profiles not ready yet, has some bugs) 604 607 ! CALL chem_init_profiles … … 780 783 ! Description: 781 784 ! ------------ 782 !> Subroutine defining parin for &chemistry_par for chemistry model785 !> Subroutine defining parin for &chemistry_parameters for chemistry model 783 786 !------------------------------------------------------------------------------! 784 787 SUBROUTINE chem_parin … … 797 800 REAL(wp), DIMENSION(nmaxfixsteps) :: my_steps !< List of fixed timesteps my_step(1) = 0.0 automatic stepping 798 801 799 NAMELIST /chemistry_par / bc_cs_b,&800 bc_cs_t,&801 call_chem_at_all_substeps,&802 chem_debug0,&803 chem_debug1,&804 chem_debug2,&805 chem_gasphase_on,&806 cs_heights,&807 cs_name,&808 cs_profile,&809 cs_profile_name,&810 cs_surface,&811 emiss_factor_main,&812 emiss_factor_side,&813 icntrl,&814 main_street_id,&815 max_street_id,&816 my_steps,&817 rcntrl,&818 side_street_id,&819 photolysis_scheme,&820 wall_csflux,&821 cs_vertical_gradient,&822 top_csflux,&823 surface_csflux,&824 surface_csflux_name,&825 cs_surface_initial_change,&826 cs_vertical_gradient_level802 NAMELIST /chemistry_parameters/ bc_cs_b, & 803 bc_cs_t, & 804 call_chem_at_all_substeps, & 805 chem_debug0, & 806 chem_debug1, & 807 chem_debug2, & 808 chem_gasphase_on, & 809 cs_heights, & 810 cs_name, & 811 cs_profile, & 812 cs_profile_name, & 813 cs_surface, & 814 emiss_factor_main, & 815 emiss_factor_side, & 816 icntrl, & 817 main_street_id, & 818 max_street_id, & 819 my_steps, & 820 rcntrl, & 821 side_street_id, & 822 photolysis_scheme, & 823 wall_csflux, & 824 cs_vertical_gradient, & 825 top_csflux, & 826 surface_csflux, & 827 surface_csflux_name, & 828 cs_surface_initial_change, & 829 cs_vertical_gradient_level 827 830 828 831 !-- analogue to chem_names(nspj) we could invent chem_surfaceflux(nspj) and chem_topflux(nspj) 829 832 !-- so this way we could prescribe a specific flux value for each species 830 !> chemistry_par for initial profiles833 !> chemistry_parameters for initial profiles 831 834 !> cs_names = 'O3', 'NO2', 'NO', ... to set initial profiles) 832 835 !> cs_heights(1,:) = 0.0, 100.0, 500.0, 2000.0, .... (height levels where concs will be prescribed for O3) … … 839 842 !-- Read chem namelist 840 843 !-- (todo: initialize these parameters in declaration part, do this for 841 !-- all chemistry_par namelist parameters)844 !-- all chemistry_parameters namelist parameters) 842 845 icntrl = 0 843 846 rcntrl = 0.0_wp … … 851 854 REWIND ( 11 ) 852 855 line = ' ' 853 DO WHILE ( INDEX( line, '&chemistry_par ' ) == 0 )856 DO WHILE ( INDEX( line, '&chemistry_parameters' ) == 0 ) 854 857 READ ( 11, '(A)', END=10 ) line 855 858 ENDDO … … 857 860 ! 858 861 !-- Read chemistry namelist 859 READ ( 11, chemistry_par )862 READ ( 11, chemistry_parameters ) 860 863 ! 861 864 !-- Enable chemistry model … … 1773 1776 ! 1774 1777 ! ! 1775 ! !-- Writing out input parameters that are not part of chemistry_par namelist1776 ! !-- (namelist parameters are anyway read in again in case of restart)1778 ! !-- Writing out input parameters that are not part of chemistry_parameters 1779 ! !-- namelist (namelist parameters are anyway read in again in case of restart) 1777 1780 ! DO lsp = 1, nvar 1778 1781 ! CALL wrd_write_string( 'conc_pr_init_'//chem_species(lsp)%name )
Note: See TracChangeset
for help on using the changeset viewer.