Changeset 2766 for palm/trunk


Ignore:
Timestamp:
Jan 22, 2018 5:17:47 PM (7 years ago)
Author:
kanani
Message:

Removal of chem directive, plus minor changes

Location:
palm/trunk/SOURCE
Files:
21 edited

Legend:

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

    r2742 r2766  
    2525! -----------------
    2626! $Id$
     27! Removed preprocessor directive __chem
     28!
     29! 2742 2018-01-12 14:59:47Z suehring
    2730! Enable output of surface temperature
    2831!
     
    126129    USE averaging
    127130
    128 #if defined( __chem )
    129131    USE chemistry_model_mod,                                                   &
    130132        ONLY:  chem_3d_data_averaging
    131 #endif
    132133
    133134    USE control_parameters,                                                    &
     
    514515!
    515516!--          Chemistry quantity
    516 #if defined( __chem )
    517517             IF ( air_chemistry )  THEN
    518518                CALL chem_3d_data_averaging( 'average', doav(ii) )
    519519             ENDIF
    520 #endif
    521520!
    522521!--          User-defined quantity
  • palm/trunk/SOURCE/boundary_conds.f90

    r2718 r2766  
    2525! -----------------
    2626! $Id$
     27! Removed preprocessor directive __chem
     28!
     29! 2718 2018-01-02 08:49:38Z maronga
    2730! Corrected "Former revisions" section
    2831!
     
    179182               w, w_p, w_m_l, w_m_n, w_m_r, w_m_s, pt_init
    180183
    181 #if defined( __chem )
    182184    USE chemistry_model_mod,                                                   &
    183185        ONLY:  chem_boundary_conds
    184 #endif
    185186             
    186187    USE control_parameters,                                                    &
     
    525526!
    526527!-- Top/bottom boundary conditions for chemical species
    527 #if defined( __chem )
    528528    IF ( air_chemistry )  CALL chem_boundary_conds( 'set_bc_bottomtop' )
    529 #endif
    530529!
    531530!-- In case of inflow or nest boundary at the south boundary the boundary for v
     
    630629!
    631630!-- Lateral boundary conditions for chemical species
    632 #if defined( __chem )
    633631    IF ( air_chemistry )  CALL chem_boundary_conds( 'set_bc_lateral' )   
    634 #endif
    635 
    636632
    637633!
  • palm/trunk/SOURCE/check_parameters.f90

    r2765 r2766  
    2525! -----------------
    2626! $Id$
     27! Removed preprocessor directive __chem
     28!
     29! 2765 2018-01-22 11:34:58Z maronga
    2730! Renamed simulation_time_since_reference to
    2831! time_to_be_simulated_from_reference_point
     
    578581
    579582    USE arrays_3d
    580 #if defined( __chem )
    581583    USE chemistry_model_mod,                                                   &
    582584        ONLY:  chem_boundary_conds, chem_check_data_output,                    &
    583585               chem_check_data_output_pr, chem_species                 
    584586    USE chem_modules
    585 #endif
    586587    USE cloud_parameters
    587588    USE constants
     
    14011402       IF ( ocean          )  sa_init = sa_surface
    14021403       IF ( passive_scalar )  s_init  = s_surface
    1403 #if defined(__chem )                             
    14041404       IF ( air_chemistry )  THEN
    14051405         DO lsp = 1, nvar
     
    14071407         ENDDO
    14081408       ENDIF
    1409 #endif
    14101409!
    14111410!--
     
    16191618!
    16201619!--    Compute initial chemistry profile using the given chemical species gradients
    1621 #if defined( __chem )
    16221620       IF ( air_chemistry )  THEN                                                         
    16231621         DO lsp = 1, nvar
     
    16291627         ENDDO
    16301628       ENDIF
    1631 #endif
    16321629!
    16331630!
     
    21132110!
    21142111!-- Boundary conditions for chemical species
    2115 #if defined( __chem )
    21162112    IF ( air_chemistry )  CALL chem_boundary_conds( 'init' )
    2117 #endif
    21182113
    21192114!
     
    30703065             ENDIF
    30713066
    3072 #if defined( __chem )
    30733067             IF ( unit == 'illegal' )  THEN                                       
    30743068                CALL chem_check_data_output_pr( data_output_pr(i), i,          &
    30753069                                                unit, dopr_unit(i) )
    30763070             ENDIF
    3077 #endif
    30783071
    30793072             IF ( unit == 'illegal' )  THEN
     
    33863379!
    33873380!--          Block of chemistry model outputs
    3388 #if defined( __chem )
    33893381             IF ( unit == 'illegal'  .AND.  air_chemistry                      &
    33903382                                     .AND.  var(1:3) == 'kc_' )  THEN
    33913383                CALL chem_check_data_output( var, unit, i, ilen, k )
    33923384             ENDIF
    3393 #endif
    33943385
    33953386!
  • palm/trunk/SOURCE/chem_gasphase_mod.f90

    r2718 r2766  
    11MODULE chem_gasphase_mod
    2 #if defined( __chem )
    32!------------------------------------------------------------------------------!
    43!
     
    3231! -----------------
    3332! $Id: module_header 2460 2017-09-13 14:47:48Z forkel $
     33! Removed preprocessor directive __chem
     34!
     35! 2460 2017-09-13 14:47:48Z forkel
    3436!
    3537!
     
    24332435 
    24342436  END   SUBROUTINE fill_temp
    2435 #endif
     2437
    24362438END MODULE chem_gasphase_mod
    24372439 
  • palm/trunk/SOURCE/chem_modules.f90

    r2718 r2766  
    2727! -----------------
    2828! $Id$
     29! Removed preprocessor directive __chem
     30!
     31! 2718 2018-01-02 08:49:38Z maronga
    2932! Initial revision
    3033!
     
    4548!
    4649 MODULE chem_modules
    47 #if defined( __chem ) 
     50
    4851    USE kinds
    4952    USE chem_gasphase_mod,                                                     &   
     
    128131
    129132    SAVE
    130 #endif
     133
    131134 END MODULE chem_modules
    132135
  • palm/trunk/SOURCE/chem_photolysis_mod.f90

    r2718 r2766  
    2626! -----------------
    2727! $Id$
     28! Removed preprocessor directive __chem
     29!
     30! 2718 2018-01-02 08:49:38Z maronga
    2831! Initial revision
    2932!
     
    4346!------------------------------------------------------------------------------!
    4447 MODULE chem_photolysis_mod
    45 #if defined( __chem )
     48
    4649    USE arrays_3d,                                                             &
    4750        ONLY:  dzw, hyp, pt, q, ql, zu, zw
     
    329332
    330333    END SUBROUTINE photolysis_simple
    331 #endif
     334
    332335 END MODULE chem_photolysis_mod
  • palm/trunk/SOURCE/chemistry_model_mod.f90

    r2756 r2766  
    2727! -----------------
    2828! $Id$
     29! Removed preprocessor directive __chem
     30!
     31! 2756 2018-01-16 18:11:14Z suehring
    2932! Fill values in 3D output introduced.
    3033!
     
    7477
    7578MODULE chemistry_model_mod
    76 #if defined( __chem )
     79
    7780   USE kinds,              ONLY: wp, iwp
    7881   USE indices,            ONLY: nz, nzb,nzt,nysg,nyng,nxlg,nxrg,nys,nyn
     
    18081811
    18091812
    1810 #endif
    18111813 END MODULE chemistry_model_mod
    18121814
  • palm/trunk/SOURCE/data_output_3d.f90

    r2756 r2766  
    2525! -----------------
    2626! $Id$
     27! Removed preprocessor directive __chem
     28!
     29! 2756 2018-01-16 18:11:14Z suehring
    2730! Fill values for 3D output of chemical species introduced.
    2831!
     
    185188    USE averaging
    186189       
    187 #if defined( __chem )
    188190    USE chemistry_model_mod,                                                   &
    189191        ONLY:  chem_data_output_3d
    190 #endif       
    191192
    192193    USE cloud_parameters,                                                      &
     
    680681!
    681682!--          Chemistry model output
    682 #if defined( __chem )
    683683             IF ( .NOT. found  .AND.  air_chemistry )  THEN
    684684                CALL chem_data_output_3d( av, do3d(av,if), found,              &
     
    686686                resorted = .TRUE.
    687687             ENDIF
    688 #endif
    689688
    690689!
  • palm/trunk/SOURCE/init_3d_model.f90

    r2758 r2766  
    2525! -----------------
    2626! $Id$
     27! Removed preprocessor directive __chem
     28!
     29! 2758 2018-01-17 12:55:21Z suehring
    2730! In case of spinup of land- and urban-surface model, do not mask wind velocity
    2831! at first computational grid level
     
    411414    USE arrays_3d
    412415
    413 #if defined( __chem )
    414416    USE chemistry_model_mod,                                                   &
    415417        ONLY:  chem_emissions
    416 #endif
    417418
    418419    USE cloud_parameters,                                                      &
     
    23022303!-- If required, set chemical emissions
    23032304!-- (todo(FK): This should later on be CALLed time-dependently in init_3d_model)
    2304 #if defined( __chem )
    23052305    IF ( air_chemistry ) THEN
    23062306       CALL chem_emissions
    23072307    ENDIF
    2308 #endif
    23092308
    23102309!
  • palm/trunk/SOURCE/netcdf_interface_mod.f90

    r2746 r2766  
    2525! -----------------
    2626! $Id$
     27! Removed preprocessor directive __chem
     28!
     29! 2746 2018-01-15 12:06:04Z suehring
    2730! Move flag plant canopy to modules
    2831!
     
    481484        ONLY:  zu, zw
    482485
    483 #if defined( __chem )
    484486    USE chemistry_model_mod,                                                   &
    485487        ONLY:  chem_define_netcdf_grid
    486 #endif
    487488
    488489    USE constants,                                                             &
     
    930931!
    931932!--                Check for chemistry quantities                   
    932 #if defined( __chem )
    933933                   IF ( .NOT. found  .AND.  air_chemistry )  THEN
    934934                      CALL chem_define_netcdf_grid( domask(mid,av,i),          &
     
    936936                                                    grid_z )
    937937                   ENDIF
    938 #endif
    939938
    940939!
     
    14741473!
    14751474!--                Check for chemistry quantities                   
    1476 #if defined( __chem )
    14771475                   IF ( .NOT. found  .AND.  air_chemistry )  THEN
    14781476                      CALL chem_define_netcdf_grid( do3d(av,i), found,         &
    14791477                                                    grid_x, grid_y, grid_z )
    14801478                   ENDIF
    1481 #endif
    14821479
    14831480!--                Check for user-defined quantities
     
    21292126!
    21302127!--                      Check for chemistry quantities
    2131 #if defined( __chem )
    21322128                         IF ( .NOT. found  .AND.  air_chemistry )  THEN
    21332129                            CALL chem_define_netcdf_grid( do2d(av,i), found,   &
     
    21352131                                                          grid_z )
    21362132                         ENDIF
    2137 #endif
    21382133
    21392134!
     
    28492844!
    28502845!--                   Check for chemistry quantities
    2851 #if defined( __chem )
    28522846                      IF ( .NOT. found  .AND.  air_chemistry )  THEN
    28532847                         CALL chem_define_netcdf_grid( do2d(av,i), found,      &
     
    28552849                                                       grid_z )
    28562850                      ENDIF
    2857 #endif
    28582851
    28592852!
     
    35183511!
    35193512!--                   Check for chemistry quantities
    3520 #if defined( __chem )
    35213513                      IF ( .NOT. found  .AND.  air_chemistry )  THEN
    35223514                         CALL chem_define_netcdf_grid( do2d(av,i), found,      &
     
    35243516                                                       grid_z )
    35253517                      ENDIF
    3526 #endif
    35273518!
    35283519!--                   Check for user-defined quantities
  • palm/trunk/SOURCE/palm.f90

    r2720 r2766  
    2525! -----------------
    2626! $Id$
     27! Removed preprocessor directive __chem
     28!
     29! 2720 2018-01-02 16:27:15Z kanani
    2730! Version update to 5.0
    2831!
     
    205208    USE arrays_3d
    206209
    207 #if defined( __chem )
    208210    USE chemistry_model_mod,                                                   &
    209211        ONLY:  chem_init, chem_last_actions
     212
    210213    USE chem_photolysis_mod,                                                   &
    211214        ONLY:  photolysis_init
    212 #endif
    213215
    214216    USE control_parameters,                                                    &
     
    374376!-- Initialize chemistry (called before check_parameters due to dependencies)
    375377!-- --> Needs to be moved!! What is the dependency about?
    376 #if defined( __chem )
    377378! IF (  TRIM( initializing_actions ) /= 'read_restart_data' )  THEN
    378379    IF ( air_chemistry )  THEN
     
    381382    ENDIF
    382383! END IF
    383 #endif
    384384!
    385385!-- Check control parameters and deduce further quantities
     
    518518             CALL radiation_last_actions
    519519          ENDIF
    520 #if defined( __chem )
    521520          IF ( air_chemistry )  THEN
    522521             CALL chem_last_actions
    523522          ENDIF
    524 #endif
    525523         
    526524          CALL user_last_actions
  • palm/trunk/SOURCE/parin.f90

    r2718 r2766  
    2525! -----------------
    2626! $Id$
     27! Removed preprocessor directive __chem
     28!
     29! 2718 2018-01-02 08:49:38Z maronga
    2730! Corrected "Former revisions" section
    2831!
     
    330333               ug, u_init, v_init, vg
    331334
    332 #if defined( __chem )
    333335    USE chemistry_model_mod,                                                   &
    334336        ONLY:  chem_parin
    335337       
    336338    USE chem_modules
    337 #endif
    338339
    339340    USE control_parameters
     
    743744!
    744745!--       Read chemistry variables
    745 #if defined( __chem )
    746746          CALL chem_parin
    747 #endif
    748747!
    749748!--       Check if uv exposure model is used and read &uvexposure_par
  • palm/trunk/SOURCE/plant_canopy_model_mod.f90

    r2746 r2766  
    2525! -----------------
    2626! $Id$
     27! Increased LEN of canopy mode to 30
     28!
     29! 2746 2018-01-15 12:06:04Z suehring
    2730! Move flag plant canopy to modules
    2831!
     
    170173
    171174
    172     CHARACTER (LEN=20)   ::  canopy_mode = 'block' !< canopy coverage
     175    CHARACTER (LEN=30)   ::  canopy_mode = 'block' !< canopy coverage
    173176
    174177    INTEGER(iwp) ::  pch_index = 0                               !< plant canopy height/top index
  • palm/trunk/SOURCE/prognostic_equations.f90

    r2746 r2766  
    2525! -----------------
    2626! $Id$
     27! Removed preprocessor directive __chem
     28!
     29! 2746 2018-01-15 12:06:04Z suehring
    2730! Move flag plant canopy to modules
    2831!
     
    266269               w, w_p
    267270
    268 #if defined( __chem )
    269271    USE chemistry_model_mod,                                                   &
    270272        ONLY:  chem_integrate, chem_species,                                   &
     
    276278    USE chem_modules,                                                          &
    277279        ONLY:  call_chem_at_all_substeps, chem_gasphase_on
    278 #endif
    279280
    280281    USE control_parameters,                                                    &
     
    444445!-- since exchange of ghost points is required after this update of the
    445446!-- concentrations of chemical species                                   
    446 #if defined( __chem )
    447447    IF ( air_chemistry )  THEN
    448448!
     
    478478
    479479    ENDIF       
    480 #endif     
    481480
    482481!
     
    12981297!
    12991298!--       If required, compute prognostic equation for chemical quantites
    1300 #if defined( __chem )
    13011299          IF ( air_chemistry )  THEN
    13021300             CALL cpu_log( log_point(83), '(chem advec+diff+prog)', 'start' )
     
    13171315             CALL cpu_log( log_point(83), '(chem advec+diff+prog)', 'stop' )             
    13181316          ENDIF   ! Chemicals equations                     
    1319 #endif
    13201317
    13211318       ENDDO
  • palm/trunk/SOURCE/read_3d_binary.f90

    r2718 r2766  
    2525! -----------------
    2626! $Id$
     27! Removed preprocessor directive __chem
     28!
     29! 2718 2018-01-02 08:49:38Z maronga
    2730! Corrected "Former revisions" section
    2831!
     
    155158    USE averaging
    156159
    157 #if defined( __chem )
    158160    USE chemistry_model_mod,                                                   &
    159161        ONLY:  chem_read_restart_data, chem_species, nspec                               
    160 #endif
    161162
    162163    USE control_parameters,                                                    &
     
    11141115!
    11151116!--    Read chemistry restart data
    1116 #if defined( __chem )
    11171117       IF ( air_chemistry )  THEN
    11181118          CALL chem_read_restart_data( i, nxlfa, nxl_on_file, nxrfa,           &
     
    11221122                                       tmp_2d, tmp_3d )
    11231123       ENDIF
    1124 #endif
    11251124
    11261125!
  • palm/trunk/SOURCE/sum_up_3d_data.f90

    r2743 r2766  
    2525! -----------------
    2626! $Id$
     27! Removed preprocessor directive __chem
     28!
     29! 2743 2018-01-12 16:03:39Z suehring
    2730! In case of natural- and urban-type surfaces output surfaces fluxes in W/m2.
    2831!
     
    172175               ts_av, tsurf_av, u_av, us_av, v_av, vpt_av, w_av, z0_av, z0h_av,&
    173176               z0q_av
    174 #if defined( __chem )
    175177    USE chemistry_model_mod,                                                   &
    176178        ONLY:  chem_3d_data_averaging, chem_integrate, chem_species, nspec                                   
    177 #endif
    178179
    179180    USE cloud_parameters,                                                      &
     
    10231024!
    10241025!--          Chemical quantity
    1025 #if defined( __chem )               
    10261026             IF ( air_chemistry  .AND.  trimvar(1:3) == 'kc_')  THEN
    10271027                CALL chem_3d_data_averaging( 'sum',doav(ii) )
    10281028             ENDIF
    1029 #endif
    10301029
    10311030!
  • palm/trunk/SOURCE/surface_layer_fluxes_mod.f90

    r2718 r2766  
    2626! -----------------
    2727! $Id$
     28! Removed preprocessor directive __chem
     29!
     30! 2718 2018-01-02 08:49:38Z maronga
    2831! Corrected "Former revisions" section
    2932!
     
    207210               drho_air_zw, rho_air_zw
    208211
    209 #if defined( __chem )
    210212    USE chem_modules,                                                          &
    211213        ONLY:  constant_csflux, nvar
    212 #endif
    213214
    214215    USE cloud_parameters,                                                      &
     
    18141815!
    18151816!--    If required compute cs* (chemical species)
    1816 #if defined( __chem )
    18171817       IF ( air_chemistry  )  THEN 
    18181818!
     
    18421842          ENDIF
    18431843       ENDIF
    1844 #endif
    18451844
    18461845!
     
    20502049!
    20512050!--       Compute the vertical chemical species' flux
    2052 #if defined( __chem )
    20532051          DO  lsp = 1, nvar
    20542052             IF (  .NOT.  constant_csflux(lsp)  .AND.  air_chemistry  .AND.    &
     
    20622060             ENDIF   
    20632061          ENDDO
    2064 #endif
    20652062
    20662063!
  • palm/trunk/SOURCE/surface_mod.f90

    r2759 r2766  
    2626! -----------------
    2727! $Id$
     28! Removed preprocessor directive __chem
     29!
     30! 2759 2018-01-17 16:24:59Z suehring
    2831! Bugfix, consider density in vertical fluxes of passive scalar as well as
    2932! chemical species.
     
    123126               rho_air, rho_air_zw, zu, zw, waterflux_input_conversion
    124127
    125 #if defined( __chem )
    126128    USE chem_modules
    127 #endif
    128129
    129130    USE control_parameters               
     
    968969!
    969970!--    Scaling parameter (cs*) and surface flux of chemical species
    970 #if defined ( __chem )
    971971       IF ( air_chemistry )  THEN
    972972          ALLOCATE ( surfaces%css(1:nvar,1:surfaces%ns)   )   
    973973          ALLOCATE ( surfaces%cssws(1:nvar,1:surfaces%ns) )
    974974       ENDIF
    975 #endif
    976975!
    977976!--    Arrays for storing potential temperature and
     
    10521051!
    10531052!--    Chemical species flux
    1054 #if defined( __chem )
    10551053       IF ( air_chemistry )  THEN
    10561054          ALLOCATE ( surfaces%cssws(1:nvar,1:surfaces%ns) )
    10571055       ENDIF
    1058 #endif
    10591056!
    10601057!--       
     
    11671164!
    11681165!--    Scaling parameter (cs*) and surface flux of chemical species
    1169 #if defined( __chem )
    11701166       IF ( air_chemistry )  THEN
    11711167             ALLOCATE ( surfaces%css(1:nvar,1:surfaces%ns)   )   
    11721168             ALLOCATE ( surfaces%cssws(1:nvar,1:surfaces%ns) )
    11731169       ENDIF
    1174 #endif
    11751170!
    11761171!--    Arrays for storing potential temperature and
     
    17371732             IF ( passive_scalar )  surf%ss(num_h) = 0.0_wp
    17381733
    1739 #if defined( __chem )
    17401734             DO  lsp = 1, nvar
    17411735                IF ( air_chemistry )  surf%css(lsp,num_h) = 0.0_wp
    17421736             ENDDO
    1743 #endif
    17441737!
    17451738!--          Set initial value for surface temperature
     
    18191812                ENDIF
    18201813
    1821 #if defined( __chem )
    18221814                IF ( air_chemistry )  THEN
    18231815                   lsp_pr = 1
     
    18491841                   ENDDO
    18501842                ENDIF
    1851 #endif
    18521843
    18531844                IF ( ocean )  THEN
     
    19201911!
    19211912!--          Prescribe top chemical species' flux
    1922 #if defined( __chem )
    19231913             DO  lsp = 1, nvar
    19241914                IF ( air_chemistry  .AND.  constant_top_csflux(lsp) )  THEN
     
    19261916                ENDIF
    19271917             ENDDO
    1928 #endif
    19291918!
    19301919!--          Prescribe top salinity flux
     
    20632052             ENDIF
    20642053
    2065 #if defined( __chem )
    20662054             IF ( air_chemistry )  THEN       
    20672055                DO  lsp = 1, nvar
     
    20702058                ENDDO
    20712059             ENDIF
    2072 #endif
    20732060
    20742061!
     
    22712258                   IF ( ALLOCATED( surf_def_h(l)%ssws ) )                      &
    22722259                      surf_h(l)%ssws(mm(l))    = surf_def_h(l)%ssws(m)
    2273 #if defined( __chem )
    22742260                   IF ( ALLOCATED( surf_def_h(l)%css ) )  THEN
    22752261                      DO  lsp = 1,nvar
     
    22822268                      ENDDO
    22832269                   ENDIF
    2284 #endif
    22852270                   IF ( ALLOCATED( surf_def_h(l)%ncsws ) )                     &
    22862271                      surf_h(l)%ncsws(mm(l))   = surf_def_h(l)%ncsws(m)
     
    23282313                      IF ( ALLOCATED( surf_lsm_h%ssws ) )                      &
    23292314                         surf_h(0)%ssws(mm(0))    = surf_lsm_h%ssws(m)
    2330 #if defined( __chem )
    23312315                      IF ( ALLOCATED( surf_lsm_h%css ) )  THEN                 
    23322316                         DO  lsp = 1, nvar
     
    23392323                         ENDDO
    23402324                      ENDIF
    2341 #endif
    23422325                      IF ( ALLOCATED( surf_lsm_h%ncsws ) )                     &
    23432326                         surf_h(0)%ncsws(mm(0))   = surf_lsm_h%ncsws(m)
     
    23852368                      IF ( ALLOCATED( surf_usm_h%ssws ) )                      &
    23862369                         surf_h(0)%ssws(mm(0))    = surf_usm_h%ssws(m)
    2387 #if defined( __chem )
    23882370                      IF ( ALLOCATED( surf_usm_h%css ) )  THEN             
    23892371                         DO lsp = 1, nvar
     
    23962378                         ENDDO
    23972379                      ENDIF
    2398 #endif
    23992380                      IF ( ALLOCATED( surf_usm_h%ncsws ) )                     &
    24002381                         surf_h(0)%ncsws(mm(0))   = surf_usm_h%ncsws(m)
     
    24742455                   IF ( ALLOCATED( surf_def_v(l)%ssws ) )                      &
    24752456                      surf_v(l)%ssws(mm(l))    = surf_def_v(l)%ssws(m)
    2476 #if defined( __chem )
    24772457                   IF ( ALLOCATED( surf_def_v(l)%css ) )  THEN               
    24782458                      DO  lsp = 1, nvar
     
    24852465                      ENDDO
    24862466                   ENDIF
    2487 #endif
    24882467                   IF ( ALLOCATED( surf_def_v(l)%ncsws ) )                     &
    24892468                      surf_v(l)%ncsws(mm(l))   = surf_def_v(l)%ncsws(m)
     
    25362515                   IF ( ALLOCATED( surf_lsm_v(l)%ssws ) )                      &
    25372516                      surf_v(l)%ssws(mm(l))    = surf_lsm_v(l)%ssws(m)
    2538 #if defined( __chem )
    25392517                   IF ( ALLOCATED( surf_lsm_v(l)%css ) )  THEN             
    25402518                      DO  lsp = 1, nvar
     
    25472525                      ENDDO
    25482526                   ENDIF
    2549 #endif
    25502527                   IF ( ALLOCATED( surf_lsm_v(l)%ncsws ) )                     &
    25512528                      surf_v(l)%ncsws(mm(l))   = surf_lsm_v(l)%ncsws(m)
     
    25982575                   IF ( ALLOCATED( surf_usm_v(l)%ssws ) )                      &
    25992576                      surf_v(l)%ssws(mm(l))    = surf_usm_v(l)%ssws(m)
    2600 #if defined( __chem )
    26012577                   IF ( ALLOCATED( surf_usm_v(l)%css ) )  THEN             
    26022578                      DO  lsp = 1, nvar
     
    26092585                      ENDDO
    26102586                   ENDIF
    2611 #endif
    26122587                   IF ( ALLOCATED( surf_usm_v(l)%ncsws ) )                     &
    26132588                      surf_v(l)%ncsws(mm(l))   = surf_usm_v(l)%ncsws(m)
     
    27172692             WRITE ( 14 )  surf_h(l)%ssws
    27182693          ENDIF
    2719 #if defined ( __chem )
    27202694          WRITE ( 14 )  'surf_h(' // dum // ')%css                 '
    27212695          IF ( ALLOCATED ( surf_h(l)%css ) )  THEN
     
    27262700             WRITE ( 14 )  surf_h(l)%cssws
    27272701          ENDIF
    2728 #endif
    27292702          WRITE ( 14 )  'surf_h(' // dum // ')%qcsws               ' 
    27302703          IF ( ALLOCATED ( surf_h(l)%qcsws ) )  THEN
     
    28142787             WRITE ( 14 )  surf_v(l)%ssws
    28152788          ENDIF
    2816 #if defined( __chem )
    28172789          WRITE ( 14 )  'surf_v(' // dum // ')%css                 ' 
    28182790          IF ( ALLOCATED ( surf_v(l)%css ) )  THEN
     
    28232795             WRITE ( 14 )  surf_v(l)%cssws
    28242796          ENDIF
    2825 #endif
    28262797          WRITE ( 14 )  'surf_v(' // dum // ')%qcsws               ' 
    28272798          IF ( ALLOCATED ( surf_v(l)%qcsws ) )  THEN
     
    30593030                      IF ( ALLOCATED( surf_h(0)%ssws )  .AND.  kk == 1 )       &
    30603031                         READ ( 13 )  surf_h(0)%ssws
    3061 #if defined( __chem )
    30623032                   CASE ( 'surf_h(0)%css' )
    30633033                      IF ( ALLOCATED( surf_h(0)%css )  .AND.  kk == 1 )        &
     
    30663036                      IF ( ALLOCATED( surf_h(0)%cssws )  .AND.  kk == 1 )      &
    30673037                         READ ( 13 )  surf_h(0)%cssws
    3068 #endif
    30693038                   CASE ( 'surf_h(0)%qcsws' )         
    30703039                      IF ( ALLOCATED( surf_h(0)%qcsws )  .AND.  kk == 1 )      &
     
    31383107                      IF ( ALLOCATED( surf_h(1)%ssws )  .AND.  kk == 1 )       &
    31393108                         READ ( 13 )  surf_h(1)%ssws
    3140 #if defined( __chem )
    31413109                   CASE ( 'surf_h(1)%css' )
    31423110                      IF ( ALLOCATED( surf_h(1)%css )  .AND.  kk == 1 )        &
     
    31453113                      IF ( ALLOCATED( surf_h(1)%cssws )  .AND.  kk == 1 )      &
    31463114                         READ ( 13 )  surf_h(1)%cssws
    3147 #endif
    31483115                   CASE ( 'surf_h(1)%qcsws' )         
    31493116                      IF ( ALLOCATED( surf_h(1)%qcsws )  .AND.  kk == 1 )      &
     
    32143181                      IF ( ALLOCATED( surf_h(2)%ssws )  .AND.  kk == 1 )       &
    32153182                         READ ( 13 )  surf_h(2)%ssws
    3216 #if defined( __chem )
    32173183                   CASE ( 'surf_h(2)%css' )
    32183184                      IF ( ALLOCATED( surf_h(2)%css )  .AND.  kk == 1 )        &
     
    32213187                      IF ( ALLOCATED( surf_h(2)%cssws )  .AND.  kk == 1 )      &
    32223188                         READ ( 13 )  surf_h(2)%cssws
    3223 #endif
    32243189                   CASE ( 'surf_h(2)%qcsws' )         
    32253190                      IF ( ALLOCATED( surf_h(2)%qcsws )  .AND.  kk == 1 )      &
     
    32893254                      IF ( ALLOCATED( surf_v(0)%ssws )  .AND.  kk == 1 )       &
    32903255                         READ ( 13 )  surf_v(0)%ssws
    3291 #if defined( __chem )
    32923256                   CASE ( 'surf_v(0)%css' )
    32933257                      IF ( ALLOCATED( surf_v(0)%css )  .AND.  kk == 1 )        &
     
    32963260                      IF ( ALLOCATED( surf_v(0)%cssws )  .AND.  kk == 1 )      &
    32973261                         READ ( 13 )  surf_v(0)%cssws
    3298 #endif
    32993262                   CASE ( 'surf_v(0)%qcsws' )         
    33003263                      IF ( ALLOCATED( surf_v(0)%qcsws )  .AND.  kk == 1 )      &
     
    33713334                      IF ( ALLOCATED( surf_v(1)%ssws )  .AND.  kk == 1 )       &
    33723335                         READ ( 13 )  surf_v(1)%ssws
    3373 #if defined( __chem )
    33743336                   CASE ( 'surf_v(1)%css' )
    33753337                      IF ( ALLOCATED( surf_v(1)%css )  .AND.  kk == 1 )        &
     
    33783340                      IF ( ALLOCATED( surf_v(1)%cssws )  .AND.  kk == 1 )      &
    33793341                         READ ( 13 )  surf_v(1)%cssws
    3380 #endif
    33813342                   CASE ( 'surf_v(1)%qcsws' )         
    33823343                      IF ( ALLOCATED( surf_v(1)%qcsws )  .AND.  kk == 1 )      &
     
    34533414                      IF ( ALLOCATED( surf_v(2)%ssws )  .AND.  kk == 1 )       &
    34543415                         READ ( 13 )  surf_v(2)%ssws
    3455 #if defined( __chem )
    34563416                   CASE ( 'surf_v(2)%css' )
    34573417                      IF ( ALLOCATED( surf_v(2)%css )  .AND.  kk == 1 )        &
     
    34603420                      IF ( ALLOCATED( surf_v(2)%cssws )  .AND.  kk == 1 )      &
    34613421                         READ ( 13 )  surf_v(2)%cssws
    3462 #endif
    34633422                   CASE ( 'surf_v(2)%qcsws' )         
    34643423                      IF ( ALLOCATED( surf_v(2)%qcsws )  .AND.  kk == 1 )      &
     
    35353494                      IF ( ALLOCATED( surf_v(3)%ssws )  .AND.  kk == 1 )       &
    35363495                         READ ( 13 )  surf_v(3)%ssws
    3537 #if defined( __chem )
    35383496                   CASE ( 'surf_v(3)%css' )
    35393497                      IF ( ALLOCATED( surf_v(3)%css )  .AND.  kk == 1 )        &
     
    35423500                      IF ( ALLOCATED( surf_v(3)%cssws )  .AND.  kk == 1 )      &
    35433501                         READ ( 13 )  surf_v(3)%cssws
    3544 #endif
    35453502                   CASE ( 'surf_v(3)%qcsws' )         
    35463503                      IF ( ALLOCATED( surf_v(3)%qcsws )  .AND.  kk == 1 )      &
     
    37623719             ENDIF
    37633720
    3764 #if defined( __chem )
    37653721             IF ( INDEX( TRIM( field_chr ), '%css' ) /= 0 )  THEN
    37663722                IF ( ALLOCATED( surf_target%css )  .AND.                     &
     
    37793735                ENDIF
    37803736             ENDIF
    3781 #endif
    37823737
    37833738             IF ( INDEX( TRIM( field_chr ), '%qcs' ) /= 0 )  THEN
  • palm/trunk/SOURCE/swap_timelevel.f90

    r2718 r2766  
    2525! -----------------
    2626! $Id$
     27! Removed preprocessor directive __chem
     28!
     29! 2718 2018-01-02 08:49:38Z maronga
    2730! Corrected "Former revisions" section
    2831!
     
    130133#endif
    131134
    132 #if defined( __chem )
    133135    USE chemistry_model_mod,                                                   &
    134136        ONLY:  chem_swap_timelevel
    135 #endif
    136137
    137138    USE land_surface_model_mod,                                                &
     
    246247          ENDIF
    247248
    248 #if defined( __chem )
    249249          IF ( air_chemistry )  CALL chem_swap_timelevel(0)
    250 #endif
    251250
    252251          swap_level = 1
     
    278277          ENDIF
    279278
    280 #if defined( __chem )
    281279          IF ( air_chemistry )  CALL chem_swap_timelevel(1)
    282 #endif
    283280
    284281          swap_level = 2
  • palm/trunk/SOURCE/time_integration.f90

    r2718 r2766  
    2525! -----------------
    2626! $Id$
     27! Removed preprocessor directive __chem
     28!
     29! 2718 2018-01-02 08:49:38Z maronga
    2730! Corrected "Former revisions" section
    2831!
     
    314317        ONLY:  calc_mean_profile
    315318
    316 #if defined( __chem )
    317319    USE chemistry_model_mod,                                                   &
    318320        ONLY:  chem_emissions, chem_species
     
    320322    USE chem_modules,                                                          &
    321323        ONLY:  nspec
    322 #endif
    323324
    324325    USE control_parameters,                                                    &
     
    693694          ENDIF
    694695          IF ( passive_scalar )  CALL exchange_horiz( s_p, nbgp )
    695 #if defined( __chem )
    696696          IF ( air_chemistry )  THEN
    697697             DO  n = 1, nspec     
     
    699699             ENDDO
    700700          ENDIF
    701 #endif
     701
    702702          CALL cpu_log( log_point(26), 'exchange-horiz-progn', 'stop' )
    703703
     
    10051005       ENDDO   ! Intermediate step loop
    10061006!
    1007 !--          If required, consider chemical emissions
    1008 !--          (todo (FK): Implement hourly call of emissions, using time_utc from
    1009 !--                      data_and_time_mod.f90;
    1010 !--                      move the CALL to appropriate location)
    1011 ! #if defined( __chem )
    1012 !              IF ( air_chemistry ) THEN
    1013 !                 CALL chem_emissions
    1014 !              ENDIF
    1015 ! #endif
     1007!--    If required, consider chemical emissions
     1008!--    (todo (FK): Implement hourly call of emissions, using time_utc from
     1009!--                data_and_time_mod.f90;
     1010!--                move the CALL to appropriate location)
     1011       IF ( air_chemistry ) THEN
     1012          CALL chem_emissions
     1013       ENDIF
    10161014!
    10171015!--    If required, do UV exposure calculations
  • palm/trunk/SOURCE/urban_surface_mod.f90

    r2765 r2766  
    2626! -----------------
    2727! $Id$
     28! Removed redundant commas, added some blanks
     29!
     30! 2765 2018-01-22 11:34:58Z maronga
    2831! Major bugfix in calculation of f_shf. Adjustment of roughness lengths in
    2932! building_pars
     
    58635866                                              start_index_on_file,             &
    58645867                                              end_index_on_file ,              &
    5865                                               nxlc, nysc,           ,          &
     5868                                              nxlc, nysc,                      &
    58665869                                              nys_on_file, nyn_on_file,        &
    5867                                               nxl_on_file,nxr_on_file )
     5870                                              nxl_on_file, nxr_on_file )
    58685871#else
    58695872                      IF ( kk == 1 )  THEN
     
    60226025                                              start_index_on_file,             &
    60236026                                              end_index_on_file ,              &
    6024                                               nxlc, nysc,           ,          &
     6027                                              nxlc, nysc,                      &
    60256028                                              nys_on_file, nyn_on_file,        &
    60266029                                              nxl_on_file,nxr_on_file )
     
    61356138                                              nxlf, nxrf, nysf, nynf,          &
    61366139                                              nys_on_file, nyn_on_file,        &
    6137                                               nxl_on_file,nxr_on_file )
     6140                                              nxl_on_file, nxr_on_file )
    61386141#endif
    61396142                   CASE ( 't_window_v(0)' )
     
    61526155                                              nxlf, nxrf, nysf, nynf,          &
    61536156                                              nys_on_file, nyn_on_file,        &
    6154                                               nxl_on_file,nxr_on_file )
     6157                                              nxl_on_file, nxr_on_file )
    61556158#else
    61566159                      IF ( kk == 1 )  THEN
     
    61816184                                              start_index_on_file,             &
    61826185                                              end_index_on_file ,              &
    6183                                               nxlc, nysc,           ,          &
     6186                                              nxlc, nysc,                      &
    61846187                                              nys_on_file, nyn_on_file,        &
    6185                                               nxl_on_file,nxr_on_file )
     6188                                              nxl_on_file, nxr_on_file )
    61866189#else
    61876190                      IF ( kk == 1 )  THEN
Note: See TracChangeset for help on using the changeset viewer.