Changeset 4881 for palm/trunk/SOURCE/chemistry_model_mod.f90

Timestamp:

Feb 19, 2021 10:05:08 PM (3 years ago)

Author:

forkel

Message:

removed unused parameters and write statements in chemistry

File:

• palm/trunk/SOURCE/chemistry_model_mod.f90 (modified) (11 diffs)

palm/trunk/SOURCE/chemistry_model_mod.f90

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! removed unnecessarty namelist parameters and commented output statements
+! and cs_surface_initial_change
+!
+!
+! 4860 2021-02-01 08:10:59Z raasch
 ! further re-numbering of message IDs
 !
@@ -297 +302 @@
 ! ------------
 !> Chemistry model for PALM-4U
-!> @todo Extend chem_species type by nspec and nvar as addititional elements (RF)
-!> @todo Check possibility to reduce dimension of chem_species%conc from nspec to nvar (RF)
 !> @todo Adjust chem_rrd_local to CASE structure of others modules. It is not allowed to use the
 !>       chemistry model in a precursor run and additionally not using it in a main run
@@ -308 +311 @@
 !> @todo slight differences in passive scalar and chem spcs when chem reactions turned off. Need to
 !>       be fixed. bK
-!> @todo test nesting for chem spcs, was implemented by suehring (kanani)
 !> @todo chemistry error messages
 !
@@ -1521 +1523 @@
                   (variable(char_len-2:) == '_xz')  .OR.                                           &
                   (variable(char_len-2:) == '_yz') ) )  THEN
-!
-!--   todo: remove or replace by "CALL message" mechanism (kanani)
-!                    IF(myid == 0)  WRITE(6,*) 'Output of species ' // TRIM( variable )  //       &
-!                                                             TRIM( chem_species(lsp)%name )
              IF (av == 0)  THEN
                 DO  i = nxl, nxr
@@ -1638 +1636 @@
       DO  lsp = 1, nspec
          IF ( TRIM( spec_name ) == TRIM( chem_species(lsp)%name) )  THEN
-!
-!--   todo: remove or replace by "CALL message" mechanism (kanani)
-!              IF(myid == 0 .AND. chem_debug0 )  WRITE(6,*) 'Output of species ' // TRIM( variable )  // &
-!                                                           TRIM( chem_species(lsp)%name )
             IF (av == 0)  THEN
                DO  i = nxl, nxr
@@ -1717 +1711 @@
     DO  lsp=1,nspec
        IF (TRIM( spec_name ) == TRIM( chem_species(lsp)%name) )  THEN
-!
-!-- todo: remove or replace by "CALL message" mechanism (kanani)
-!              IF(myid == 0 .AND. chem_debug0 )  WRITE(6,*) 'Output of species ' // TRIM( variable )  // &
-!                                                        TRIM( chem_species(lsp)%name )
           IF (av == 0)  THEN
              IF ( .NOT. mask_surface(mid) )  THEN
@@ -2272 +2262 @@
        ENDIF
 !
-!--    If required, change the surface chem spcs at the start of the 3D run
-       IF ( cs_surface_initial_change(1) /= 0.0_wp )  THEN
-          DO  lsp = 1, nspec
-             DO  i = nxlg, nxrg
-                DO  j = nysg, nyng
-                   flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_total_0(nzb,j,i), 0 ) )
-                   chem_species(lsp)%conc(nzb,j,i) = chem_species(lsp)%conc(nzb,j,i) +             &
-                                                     cs_surface_initial_change(lsp) * flag
-                ENDDO
-             ENDDO
-          ENDDO
-       ENDIF
-
     ENDIF
 !
@@ -2295 +2272 @@
     DO  lsp = 1, nphot
        phot_frequen(lsp)%name = phot_names(lsp)
-!
-!-- todo: remove or replace by "CALL message" mechanism (kanani)
-!--       IF( myid == 0 )  THEN
-!--          WRITE(6,'(a,i4,3x,a)')  'Photolysis: ',lsp,TRIM( phot_names(lsp) )
-!--       ENDIF
        phot_frequen(lsp)%freq(nzb:nzt+1,nysg:nyng,nxlg:nxrg)  =>  freq_1(:,:,:,lsp)
     ENDDO
@@ -2546 +2518 @@
          bc_cs_t,                                                                                  &
          call_chem_at_all_substeps,                                                                &
-         chem_debug0,                                                                              &
-         chem_debug1,                                                                              &
-         chem_debug2,                                                                              &
          chem_gasphase_on,                                                                         &
          chem_mechanism,                                                                           &
@@ -2555 +2524 @@
          cs_profile,                                                                               &
          cs_surface,                                                                               &
-         cs_surface_initial_change,                                                                &
-         cs_vertical_gradient_level,                                                               &
          daytype_mdh,                                                                              &
          deposition_dry,                                                                           &
@@ -2575 +2542 @@
          photolysis_scheme,                                                                        &
          wall_csflux,                                                                              &
-         cs_vertical_gradient,                                                                     &
-         top_csflux,                                                                               &
          surface_csflux,                                                                           &
          surface_csflux_name,                                                                      &