Changeset 3833 for palm/trunk/SOURCE


Ignore:
Timestamp:
Mar 28, 2019 3:04:04 PM (3 years ago)
Author:
forkel
Message:

removed USE chem_gasphase_mod from chem_modules, apply USE chem_gasphase for nvar, nspec, cs_mech and spc_names instead

Location:
palm/trunk/SOURCE
Files:
12 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/Makefile

    r3745 r3833  
    874874        modules.o
    875875chem_modules.o: \
    876         chem_gasphase_mod.o \
    877876        mod_kinds.o \
    878877        modules.o
     
    14401439        bulk_cloud_model_mod.o \
    14411440        chemistry_model_mod.o \
    1442         chem_modules.o \
     1441        chem_gasphase_mod.o \
    14431442        mod_kinds.o \
    14441443        mod_particle_attributes.o \
     
    15111510        buoyancy.o \
    15121511        chemistry_model_mod.o \
     1512        chem_gasphase_mod.o \
    15131513        chem_modules.o \
    15141514        chem_photolysis_mod.o \
     
    16261626        basic_constants_and_equations_mod.o \
    16271627        bulk_cloud_model_mod.o \
     1628        chem_gasphase_mod.o \
    16281629        chem_modules.o \
    16291630        exchange_horiz_2d.o \
     
    16341635        urban_surface_mod.o
    16351636surface_mod.o: \
     1637        chem_gasphase_mod.o \
    16361638        chem_modules.o \
    16371639        model_1d_mod.o \
     
    16701672        chem_emissions_mod.o \
    16711673        chemistry_model_mod.o \
     1674        chem_gasphase_mod.o \
    16721675        chem_modules.o \
    16731676        cpulog_mod.o \
     
    18471850        cpulog_mod.o \
    18481851        chemistry_model_mod.o \
    1849         chem_modules.o \
     1852        chem_gasphase_mod.o \
    18501853        mod_kinds.o \
    18511854        modules.o \
  • palm/trunk/SOURCE/chem_gasphase_mod.f90

    r3820 r3833  
    6060  !SAVE  ! note: occurs again in automatically generated code ...
    6161
    62 !  PUBLIC :: IERR_NAMES
    63  
    64 ! PUBLIC :: SPC_NAMES,EQN_NAMES,EQN_TAGS,REQ_HET,REQ_AEROSOL,REQ_PHOTRAT &
    65 !         ,REQ_MCFCT,IP_MAX,jname
    66 
     62! Public variables
     63  PUBLIC :: atol
    6764  PUBLIC :: cs_mech
    68   PUBLIC :: eqn_names, phot_names, spc_names
     65  PUBLIC :: eqn_names
     66  PUBLIC :: fakt
    6967  PUBLIC :: nmaxfixsteps
    70   PUBLIC :: atol, rtol
    71   PUBLIC :: nspec, nreact
     68  PUBLIC :: nphot
     69  PUBLIC :: nreact
     70  PUBLIC :: nspec
     71  PUBLIC :: nvar
     72  PUBLIC :: qvap
     73  PUBLIC :: phot
     74  PUBLIC :: phot_names
     75  PUBLIC :: rconst
     76  PUBLIC :: rtol
     77  PUBLIC :: spc_names
    7278  PUBLIC :: temp
    73   PUBLIC :: qvap
    74   PUBLIC :: fakt
    75   PUBLIC :: phot
    76   PUBLIC :: rconst
    77   PUBLIC :: nvar
    78   PUBLIC :: nphot
    79   PUBLIC :: vl_dim                     ! PUBLIC to ebable other MODULEs to distiguish between scalar and vec
    80  
    81   PUBLIC :: initialize, integrate, update_rconst
     79  PUBLIC :: vl_dim                     !< PUBLIC to enable other MODULEs to distiguish between scalar and vec
     80 
     81! Public routines
    8282  PUBLIC :: chem_gasphase_integrate
     83  PUBLIC :: get_mechanism_name
     84  PUBLIC :: initialize
    8385  PUBLIC :: initialize_kpp_ctrl
    84   PUBLIC :: get_mechanism_name
     86  PUBLIC :: integrate
     87  PUBLIC :: update_rconst
    8588
    8689! END OF MODULE HEADER TEMPLATE
     
    114117!
    115118! File                 : chem_gasphase_mod_Parameters.f90
    116 ! Time                 : Wed Mar 27 09:51:40 2019
     119! Time                 : Thu Mar 28 15:59:30 2019
    117120! Working directory    : /home/forkel-r/palmstuff/work/trunk20190327/UTIL/chemistry/gasphase_preproc/tmp_kpp4palm
    118121! Equation file        : chem_gasphase_mod.kpp
     
    191194!
    192195! File                 : chem_gasphase_mod_Global.f90
    193 ! Time                 : Wed Mar 27 09:51:40 2019
     196! Time                 : Thu Mar 28 15:59:30 2019
    194197! Working directory    : /home/forkel-r/palmstuff/work/trunk20190327/UTIL/chemistry/gasphase_preproc/tmp_kpp4palm
    195198! Equation file        : chem_gasphase_mod.kpp
     
    256259!
    257260! File                 : chem_gasphase_mod_JacobianSP.f90
    258 ! Time                 : Wed Mar 27 09:51:40 2019
     261! Time                 : Thu Mar 28 15:59:30 2019
    259262! Working directory    : /home/forkel-r/palmstuff/work/trunk20190327/UTIL/chemistry/gasphase_preproc/tmp_kpp4palm
    260263! Equation file        : chem_gasphase_mod.kpp
     
    300303!
    301304! File                 : chem_gasphase_mod_Monitor.f90
    302 ! Time                 : Wed Mar 27 09:51:40 2019
     305! Time                 : Thu Mar 28 15:59:30 2019
    303306! Working directory    : /home/forkel-r/palmstuff/work/trunk20190327/UTIL/chemistry/gasphase_preproc/tmp_kpp4palm
    304307! Equation file        : chem_gasphase_mod.kpp
     
    362365!
    363366! File                 : chem_gasphase_mod_Initialize.f90
    364 ! Time                 : Wed Mar 27 09:51:40 2019
     367! Time                 : Thu Mar 28 15:59:30 2019
    365368! Working directory    : /home/forkel-r/palmstuff/work/trunk20190327/UTIL/chemistry/gasphase_preproc/tmp_kpp4palm
    366369! Equation file        : chem_gasphase_mod.kpp
     
    388391!
    389392! File                 : chem_gasphase_mod_Integrator.f90
    390 ! Time                 : Wed Mar 27 09:51:40 2019
     393! Time                 : Thu Mar 28 15:59:30 2019
    391394! Working directory    : /home/forkel-r/palmstuff/work/trunk20190327/UTIL/chemistry/gasphase_preproc/tmp_kpp4palm
    392395! Equation file        : chem_gasphase_mod.kpp
     
    446449!
    447450! File                 : chem_gasphase_mod_LinearAlgebra.f90
    448 ! Time                 : Wed Mar 27 09:51:40 2019
     451! Time                 : Thu Mar 28 15:59:30 2019
    449452! Working directory    : /home/forkel-r/palmstuff/work/trunk20190327/UTIL/chemistry/gasphase_preproc/tmp_kpp4palm
    450453! Equation file        : chem_gasphase_mod.kpp
     
    473476!
    474477! File                 : chem_gasphase_mod_Jacobian.f90
    475 ! Time                 : Wed Mar 27 09:51:40 2019
     478! Time                 : Thu Mar 28 15:59:30 2019
    476479! Working directory    : /home/forkel-r/palmstuff/work/trunk20190327/UTIL/chemistry/gasphase_preproc/tmp_kpp4palm
    477480! Equation file        : chem_gasphase_mod.kpp
     
    500503!
    501504! File                 : chem_gasphase_mod_Function.f90
    502 ! Time                 : Wed Mar 27 09:51:40 2019
     505! Time                 : Thu Mar 28 15:59:30 2019
    503506! Working directory    : /home/forkel-r/palmstuff/work/trunk20190327/UTIL/chemistry/gasphase_preproc/tmp_kpp4palm
    504507! Equation file        : chem_gasphase_mod.kpp
     
    529532!
    530533! File                 : chem_gasphase_mod_Rates.f90
    531 ! Time                 : Wed Mar 27 09:51:40 2019
     534! Time                 : Thu Mar 28 15:59:30 2019
    532535! Working directory    : /home/forkel-r/palmstuff/work/trunk20190327/UTIL/chemistry/gasphase_preproc/tmp_kpp4palm
    533536! Equation file        : chem_gasphase_mod.kpp
     
    555558!
    556559! File                 : chem_gasphase_mod_Util.f90
    557 ! Time                 : Wed Mar 27 09:51:40 2019
     560! Time                 : Thu Mar 28 15:59:30 2019
    558561! Working directory    : /home/forkel-r/palmstuff/work/trunk20190327/UTIL/chemistry/gasphase_preproc/tmp_kpp4palm
    559562! Equation file        : chem_gasphase_mod.kpp
     
    589592  INTEGER, PARAMETER, PUBLIC :: nkppctrl = 20
    590593  !
    591   INTEGER, DIMENSION(nkppctrl), PUBLIC     :: icntrl = 0
     594  ! steering PARAMETERs for chemistry solver (see kpp domumentation)
     595  INTEGER, DIMENSION(nkppctrl), PUBLIC      :: icntrl = 0
    592596  REAL(dp), DIMENSION(nkppctrl), PUBLIC     :: rcntrl = 0.0_dp
     597  ! t_steps: fixed time steps in vector mode
    593598  REAL(dp), DIMENSION(nmaxfixsteps), PUBLIC :: t_steps = 0.0_dp
    594599
  • palm/trunk/SOURCE/chem_modules.f90

    r3827 r3833  
    2727! -----------------
    2828! $Id$
     29! removed USE chem_gasphase_mod
     30!
     31! 3827 2019-03-27 17:20:32Z forkel
    2932! some formatting  and reordering (ecc)
    3033!
     
    7578 MODULE chem_modules
    7679
    77     USE chem_gasphase_mod,                                                     &   
    78         ONLY: cs_mech, nspec, nvar, spc_names
    79 
    8080    USE control_parameters,                                                    &
    8181        ONLY: varnamelength
     
    8484
    8585    IMPLICIT NONE
    86 
    87     PUBLIC cs_mech
    88     PUBLIC nspec
    89     PUBLIC nvar
    90     PUBLIC spc_names
    9186
    9287    CHARACTER (LEN=20)                ::  bc_cs_b             = 'dirichlet'         !< namelist parameter
  • palm/trunk/SOURCE/chemistry_model_mod.f90

    r3824 r3833  
    2727! -----------------
    2828! $Id: chemistry_model_mod.f90 3784 2019-03-05 14:16:20Z banzhafs
     29! added cs_mech to USE chem_gasphase_mod
     30!
     31! 3784 2019-03-05 14:16:20Z banzhafs
    2932! renamed get_mechanismname to get_mechanism_name
    3033! renamed do_emiss to emissions_anthropogenic and do_depo to deposition_dry (ecc)
     
    257260! ------------
    258261!> Chemistry model for PALM-4U
     262!> @todo Extend chem_species type by nspec and nvar as addititional elements
    259263!> @todo Adjust chem_rrd_local to CASE structure of others modules. It is not
    260264!>       allowed to use the chemistry model in a precursor run and additionally
    261265!>       not using it in a main run
    262266!> @todo Update/clean-up todo list! (FK)
    263 !> @todo Set proper fill values (/= 0) for chem output arrays! (FK)
    264267!> @todo Add routine chem_check_parameters, add checks for inconsistent or
    265268!>       unallowed parameter settings.
    266269!>       CALL of chem_check_parameters from check_parameters. (FK)
    267 !> @todo Make routine chem_header available, CALL from header.f90
    268 !>       (see e.g. how it is done in routine lsm_header in
    269 !>        land_surface_model_mod.f90). chem_header should include all setup
    270 !>        info about chemistry parameter settings. (FK)
    271270!> @todo Implement turbulent inflow of chem spcs in inflow_turbulence.
    272271!> @todo Separate boundary conditions for each chem spcs to be implemented
     
    307306
    308307    USE chem_gasphase_mod,                                                                         &
    309          ONLY:  atol, chem_gasphase_integrate, get_mechanism_name, nkppctrl, nmaxfixsteps,         &
    310          nphot, nreact, nspec, nvar, phot_names, rtol, spc_names, t_steps, vl_dim
     308         ONLY:  atol, chem_gasphase_integrate, cs_mech, get_mechanism_name, nkppctrl,              &
     309         nmaxfixsteps, nphot, nreact, nspec, nvar, phot_names, rtol, spc_names, t_steps, vl_dim
    311310
    312311    USE chem_modules
  • palm/trunk/SOURCE/multi_agent_system_mod.f90

    r3766 r3833  
    2525! -----------------
    2626! $Id$
     27! replaced nspec by nvar: only variable species should bconsidered, fixed species are not relevant
     28!
     29! 3766 2019-02-26 16:23:41Z raasch
    2730! save attribute added to local targets to avoid outlive pointer target warning
    2831!
     
    30153018
    30163019!       USE chem_gasphase_mod,                                                  &
    3017 !           ONLY:  nspec
     3020!           ONLY:  nvar
    30183021
    30193022!       USE chemistry_model_mod,                                                &
     
    31083111! !--    Get indices of PM10 and PM2.5 species, if active
    31093112!        IF ( air_chemistry ) THEN
    3110 !           DO il = 1, nspec
     3113!           DO il = 1, nvar
    31113114! print*,chem_species(il)%name
    31123115! !              IF ( spec_name(1:4) == 'PM10' ) THEN 
  • palm/trunk/SOURCE/pmc_interface_mod.f90

    r3822 r3833  
    2525! -----------------
    2626! $Id$
     27! replaced USE chem_modules by USE chem_gasphase_mod
     28!
     29! 3822 2019-03-27 13:10:23Z hellstea
    2730! Temporary increase of the vertical dimension of the parent-grid arrays and
    2831! workarrc_t is cancelled as unnecessary.
     
    397400               roughness_length, topography, volume_flow
    398401
    399     USE chem_modules,                                                          &
     402    USE chem_gasphase_mod,                                                          &
    400403        ONLY:  nspec
    401404
  • palm/trunk/SOURCE/prognostic_equations.f90

    r3820 r3833  
    2525! -----------------
    2626! $Id$
     27! added USE chem_gasphase_mod for nvar, nspec and spc_names
     28!
     29! 3820 2019-03-27 11:53:41Z forkel
    2730! renamed do_depo to deposition_dry (ecc)
    2831!
     
    378381               deposition_dry
    379382
     383    USE chem_gasphase_mod,                                                     &
     384        ONLY:  nspec, nvar, spc_names
     385
    380386    USE chem_photolysis_mod,                                                   &
    381387        ONLY:  photolysis_control
     
    383389    USE chemistry_model_mod,                                                   &
    384390        ONLY:  chem_boundary_conds, chem_depo, chem_integrate,                 &
    385                chem_prognostic_equations, chem_species,                        &
    386                nspec, nvar, spc_names
     391               chem_prognostic_equations, chem_species
    387392
    388393    USE control_parameters,                                                    &
  • palm/trunk/SOURCE/salsa_mod.f90

    r3787 r3833  
    2626! -----------------
    2727! $Id$
     28! added USE chem_gasphase_mod for nvar, nspec and spc_names
     29!
     30! 3787 2019-03-07 08:43:54Z raasch
    2831! unused variables removed
    2932!
     
    9699        ONLY:  c_p, g, p_0, pi, r_d
    97100 
     101    USE chem_gasphase_mod,                                                     &
     102        ONLY:  nspec, nvar, spc_names
     103
    98104    USE chemistry_model_mod,                                                   &
    99         ONLY:  chem_species, nspec, nvar, spc_names
     105        ONLY:  chem_species
    100106
    101107    USE chem_modules,                                                          &
  • palm/trunk/SOURCE/surface_layer_fluxes_mod.f90

    r3787 r3833  
    2626! -----------------
    2727! $Id$
     28! added USE chem_gasphase_mod
     29!
     30! 3787 2019-03-07 08:43:54Z raasch
    2831! unused variables removed
    2932!
     
    255258    USE basic_constants_and_equations_mod,                                     &
    256259        ONLY:  g, kappa, lv_d_cp, pi, rd_d_rv
     260
     261    USE chem_gasphase_mod,                                                     &
     262        ONLY:  nvar
    257263
    258264    USE chem_modules,                                                          &
  • palm/trunk/SOURCE/surface_mod.f90

    r3772 r3833  
    2626! -----------------
    2727! $Id$
     28! added USE chem_gasphase_mod (chem_modules will not transport nvar and nspec anymore)
     29!
     30! 3772 2019-02-28 15:51:57Z suehring
    2831! small change in the todo's
    2932!
     
    270273        ONLY:  heatflux_input_conversion, momentumflux_input_conversion,       &
    271274               rho_air, rho_air_zw, zu, zw, waterflux_input_conversion
     275
     276    USE chem_gasphase_mod,                                                     &
     277        ONLY:  nvar, spc_names
    272278
    273279    USE chem_modules
  • palm/trunk/SOURCE/time_integration.f90

    r3820 r3833  
    2525! -----------------
    2626! $Id$
     27! added USE chem_gasphase_mod, replaced nspec by nvar since fixed compounds are not integrated
     28!
     29! 3820 2019-03-27 11:53:41Z forkel
    2730! renamed do_emiss to emissions_anthropogenic (ecc)
    2831!
     
    489492        ONLY:  chem_emissions_setup
    490493
     494    USE chem_gasphase_mod,                                                                         &
     495        ONLY:  nvar
     496
    491497    USE chem_modules,                                                                              &
    492         ONLY:  bc_cs_t_val, cs_name, emissions_anthropogenic, nspec, nspec_out
     498        ONLY:  bc_cs_t_val, cs_name, emissions_anthropogenic, nspec_out
    493499
    494500    USE chemistry_model_mod,                                                                       &
     
    844850           bc_q_t_val  = ( q_init(nzt+1) - q_init(nzt) ) / dzu(nzt+1)
    845851           IF ( air_chemistry )  THEN
    846               DO  lsp = 1, nspec
     852              DO  lsp = 1, nvar
    847853                 bc_cs_t_val = (  chem_species(lsp)%conc_pr_init(nzt+1)                            &
    848854                                - chem_species(lsp)%conc_pr_init(nzt) )                            &
     
    10071013          IF ( passive_scalar )  CALL exchange_horiz( s_p, nbgp )
    10081014          IF ( air_chemistry )  THEN
    1009              DO  lsp = 1, nspec
     1015             DO  lsp = 1, nvar
    10101016                CALL exchange_horiz( chem_species(lsp)%conc_p, nbgp )
    10111017!
     
    11101116
    11111117                IF ( air_chemistry )  THEN
    1112                    DO  n = 1, nspec     
     1118                   DO  n = 1, nvar     
    11131119                      CALL exchange_horiz( chem_species(n)%conc, nbgp )
    11141120                   ENDDO
  • palm/trunk/SOURCE/virtual_measurement_mod.f90

    r3766 r3833  
    2525! -----------------
    2626! $Id$
     27! renamed nvar to nmeas, replaced USE chem_modules by USE chem_gasphase_mod and nspec by nvar
     28!
     29! 3766 2019-02-26 16:23:41Z raasch
    2730! unused variables removed
    2831!
     
    7679        ONLY:  q, pt, u, v, w, zu, zw
    7780
    78     USE chem_modules,                                                          &
    79         ONLY:  nspec
     81    USE chem_gasphase_mod,                                                     &
     82        ONLY:  nvar
    8083
    8184    USE chemistry_model_mod,                                                   &
     
    132135       INTEGER(iwp) ::  ns_tot = 0      !< total number of observation coordinates, for atmospheric measurements
    133136       INTEGER(iwp) ::  ntraj           !< number of trajectories of a measurement
    134        INTEGER(iwp) ::  nvar            !< number of measured variables (atmosphere + soil)
     137       INTEGER(iwp) ::  nmeas           !< number of measured variables (atmosphere + soil)
    135138       
    136139       INTEGER(iwp) ::  ns_soil = 0     !< number of observation coordinates on subdomain, for soil measurements
     
    551554!--    for a NULL to get the correct character length in order to compare
    552555!--    them with the list of allowed variables.
    553        vmea(l)%nvar   = 0
     556       vmea(l)%nmeas  = 0
    554557       DO ll = 1, SIZE( measured_variables_file )
    555558          IF ( measured_variables_file(ll)(1:1) /= CHAR(0)  .AND.              &
     
    569572                IF ( measured_variables_file(ll)(1:len_char) ==                &
    570573                     TRIM( list_allowed_variables(lll) ) )  THEN
    571                    vmea(l)%nvar = vmea(l)%nvar + 1
    572                    measured_variables(vmea(l)%nvar) =                          &
     574                   vmea(l)%nmeas = vmea(l)%nmeas + 1
     575                   measured_variables(vmea(l)%nmeas) =                         &
    573576                                       measured_variables_file(ll)(1:len_char)
    574577                ENDIF
     
    578581!
    579582!--    Allocate array for the measured variables names for the respective site.
    580        ALLOCATE( vmea(l)%measured_vars_name(1:vmea(l)%nvar) )
    581 
    582        DO  ll = 1, vmea(l)%nvar
     583       ALLOCATE( vmea(l)%measured_vars_name(1:vmea(l)%nmeas) )
     584
     585       DO  ll = 1, vmea(l)%nmeas
    583586          vmea(l)%measured_vars_name(ll) = TRIM( measured_variables(ll) )
    584587       ENDDO
     
    587590!--    chemistry mechanism.
    588591!        IF ( air_chemistry )  THEN
    589 !           DO  ll = 1, vmea(l)%nvar
     592!           DO  ll = 1, vmea(l)%nmeas
    590593!              chem_include = .FALSE.
    591 !              DO  n = 1, nspec
     594!              DO  n = 1, nvar
    592595!                 IF ( TRIM( vmea(l)%measured_vars_name(ll) ) ==                 &
    593596!                      TRIM( chem_species(n)%name ) )  chem_include = .TRUE.
     
    609612!--    as this would exceed memory requirements, particularly for trajectory
    610613!--    measurements.
    611        IF ( vmea(l)%nvar > 0 )  THEN
     614       IF ( vmea(l)%nmeas > 0 )  THEN
    612615!
    613616!--       For stationary measurements UTM coordinates are just one value and
     
    855858!
    856859!--       Allocate array to save the sampled values.
    857           ALLOCATE( vmea(l)%measured_vars(1:vmea(l)%ns,1:vmea(l)%nvar) )
     860          ALLOCATE( vmea(l)%measured_vars(1:vmea(l)%ns,1:vmea(l)%nmeas) )
    858861         
    859862          IF ( vmea(l)%soil_sampling )                                         &
    860863             ALLOCATE( vmea(l)%measured_vars_soil(1:vmea(l)%ns_soil,           &
    861                                                   1:vmea(l)%nvar) )
     864                                                  1:vmea(l)%nmeas) )
    862865!
    863866!--       Initialize with _FillValues
    864           vmea(l)%measured_vars(1:vmea(l)%ns,1:vmea(l)%nvar) = vmea(l)%fillout
     867          vmea(l)%measured_vars(1:vmea(l)%ns,1:vmea(l)%nmeas) = vmea(l)%fillout
    865868          IF ( vmea(l)%soil_sampling )                                         &
    866              vmea(l)%measured_vars_soil(1:vmea(l)%ns_soil,1:vmea(l)%nvar) =    &
     869             vmea(l)%measured_vars_soil(1:vmea(l)%ns_soil,1:vmea(l)%nmeas) =   &
    867870                                                                vmea(l)%fillout
    868871!
     
    955958                 WRITE ( 27 )  vmea(l)%ns_tot
    956959                 WRITE ( 27 )  'number of measured variables      '
    957                  WRITE ( 27 )  vmea(l)%nvar
     960                 WRITE ( 27 )  vmea(l)%nmeas
    958961                 WRITE ( 27 )  'variables                         '
    959962                 WRITE ( 27 )  vmea(l)%measured_vars_name(:)
     
    10131016!--              Skip binary writing if no observation points are defined on PE
    10141017                 IF ( vmea(l)%ns < 1  .AND.  vmea(l)%ns_soil < 1)  CYCLE                 
    1015                  DO  n = 1, vmea(l)%nvar
     1018                 DO  n = 1, vmea(l)%nmeas
    10161019                    WRITE( 27 )  vmea(l)%measured_vars_name(n)
    10171020                    IF ( vmea(l)%soil_sampling  .AND.                           &
     
    11061109!
    11071110!--     Loop over all variables measured at this site. 
    1108         DO  n = 1, vmea(l)%nvar
     1111        DO  n = 1, vmea(l)%nmeas
    11091112       
    11101113           SELECT CASE ( TRIM( vmea(l)%measured_vars_name(n) ) )
     
    12631266!--                 Run loop over all chemical species, if the measured
    12641267!--                 variable matches the interal name, sample the variable.
    1265                     DO  nn = 1, nspec                   
     1268                    DO  nn = 1, nvar                   
    12661269                       IF ( TRIM( chem_vars(1,ind_chem) ) ==                   &
    12671270                            TRIM( chem_species(nn)%name ) )  THEN                           
Note: See TracChangeset for help on using the changeset viewer.