Ignore:
Timestamp:
Apr 8, 2019 8:25:23 PM (2 years ago)
Author:
knoop
Message:

Moved loop over chem_species into chem_boundary_conds_decycle

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/prognostic_equations.f90

    r3878 r3879  
    4141!
    4242! 3840 2019-03-29 10:35:52Z knoop
    43 ! added USE chem_gasphase_mod for nvar, nspec and spc_names
     43! added USE chem_gasphase_mod for nspec, nspec and spc_names
    4444!
    4545! 3820 2019-03-27 11:53:41Z forkel
     
    388388
    389389    USE chem_modules,                                                          &
    390         ONLY:  call_chem_at_all_substeps, chem_gasphase_on, cs_name,           &
    391                deposition_dry, chem_species
     390        ONLY:  chem_gasphase_on, deposition_dry, chem_species
    392391
    393392    USE chem_gasphase_mod,                                                     &
    394         ONLY:  nspec, nvar, spc_names
     393        ONLY:  nspec, spc_names
    395394
    396395    USE chemistry_model_mod,                                                   &
    397         ONLY:  chem_boundary_conds, chem_prognostic_equations
     396        ONLY:  chem_boundary_conds_decycle, chem_prognostic_equations
    398397
    399398    USE control_parameters,                                                    &
     
    520519    LOGICAL      ::  loop_start          !<
    521520    INTEGER(iwp) ::  lsp
    522     INTEGER(iwp) ::  lsp_usr             !< lsp running index for chem spcs
    523521
    524522
     
    538536    !$OMP END PARALLEL
    539537
    540     IF ( air_chemistry )  THEN
    541 !
    542 !--    Loop over chemical species       
    543        CALL cpu_log( log_point_s(84), 'chem.exch-horiz', 'start' )
    544        DO  lsp = 1, nspec
    545           CALL exchange_horiz( chem_species(lsp)%conc, nbgp )   
    546           lsp_usr = 1 
    547           DO WHILE ( TRIM( cs_name( lsp_usr ) ) /= 'novalue' )
    548              IF ( TRIM(chem_species(lsp)%name) == TRIM(cs_name(lsp_usr)) )  THEN
    549 
    550                 CALL chem_boundary_conds( chem_species(lsp)%conc_p,                                 &
    551                                           chem_species(lsp)%conc_pr_init )
    552              
    553              ENDIF
    554              lsp_usr = lsp_usr +1
    555           ENDDO
    556 
    557 
    558        ENDDO
    559        CALL cpu_log( log_point_s(84), 'chem.exch-horiz', 'stop' )
    560 
    561     ENDIF
     538    IF ( air_chemistry )  CALL chem_boundary_conds_decycle
    562539!
    563540!-- Run SALSA and aerosol dynamic processes. SALSA is run with a longer time
     
    11421119!
    11431120!--          Loop over chemical species
    1144              DO  lsp = 1, nvar                         
     1121             DO  lsp = 1, nspec
    11451122                CALL chem_prognostic_equations( chem_species(lsp)%conc_p,      &
    11461123                                     chem_species(lsp)%conc,                   &
     
    19241901!
    19251902!--    Loop over chemical species
    1926        DO  lsp = 1, nvar                         
     1903       DO  lsp = 1, nspec
    19271904          CALL chem_prognostic_equations( chem_species(lsp)%conc_p,            &
    19281905                                          chem_species(lsp)%conc,              &
Note: See TracChangeset for help on using the changeset viewer.