Changeset 3780 for palm/trunk/SOURCE


Ignore:
Timestamp:
Mar 5, 2019 11:19:45 AM (3 years ago)
Author:
forkel
Message:

removed read from unit 10 in chemistry_model_mod.f90, added get_mechanismname

Location:
palm/trunk/SOURCE
Files:
4 edited

Legend:

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

    r3655 r3780  
    6565!         ,REQ_MCFCT,IP_MAX,jname
    6666
     67  PUBLIC :: cs_mech
    6768  PUBLIC :: eqn_names, phot_names, spc_names
    6869  PUBLIC :: nmaxfixsteps
     
    8182  PUBLIC :: chem_gasphase_integrate
    8283  PUBLIC :: initialize_kpp_ctrl
     84  PUBLIC :: get_mechanismname
    8385
    8486! END OF MODULE HEADER TEMPLATE
     
    113115!
    114116! File                 : chem_gasphase_mod_Parameters.f90
    115 ! Time                 : Fri Nov 30 13:52:19 2018
    116 ! Working directory    : /home/forkel-r/palmstuff/work/trunk20181130/UTIL/chemistry/gasphase_preproc/tmp_kpp4palm
     117! Time                 : Tue Mar  5 11:50:53 2019
     118! Working directory    : /home/forkel-r/palmstuff/work/trunk20190305/UTIL/chemistry/gasphase_preproc/tmp_kpp4palm
    117119! Equation file        : chem_gasphase_mod.kpp
    118120! Output root filename : chem_gasphase_mod
     
    190192!
    191193! File                 : chem_gasphase_mod_Global.f90
    192 ! Time                 : Fri Nov 30 13:52:19 2018
    193 ! Working directory    : /home/forkel-r/palmstuff/work/trunk20181130/UTIL/chemistry/gasphase_preproc/tmp_kpp4palm
     194! Time                 : Tue Mar  5 11:50:53 2019
     195! Working directory    : /home/forkel-r/palmstuff/work/trunk20190305/UTIL/chemistry/gasphase_preproc/tmp_kpp4palm
    194196! Equation file        : chem_gasphase_mod.kpp
    195197! Output root filename : chem_gasphase_mod
     
    236238  REAL(kind=dp):: fakt
    237239
     240! CS_MECH for check of mechanism name with namelist
     241  CHARACTER(len=30):: cs_mech
    238242
    239243! INLINED global variable declarations
     
    256260!
    257261! File                 : chem_gasphase_mod_JacobianSP.f90
    258 ! Time                 : Fri Nov 30 13:52:19 2018
    259 ! Working directory    : /home/forkel-r/palmstuff/work/trunk20181130/UTIL/chemistry/gasphase_preproc/tmp_kpp4palm
     262! Time                 : Tue Mar  5 11:50:53 2019
     263! Working directory    : /home/forkel-r/palmstuff/work/trunk20190305/UTIL/chemistry/gasphase_preproc/tmp_kpp4palm
    260264! Equation file        : chem_gasphase_mod.kpp
    261265! Output root filename : chem_gasphase_mod
     
    300304!
    301305! File                 : chem_gasphase_mod_Monitor.f90
    302 ! Time                 : Fri Nov 30 13:52:19 2018
    303 ! Working directory    : /home/forkel-r/palmstuff/work/trunk20181130/UTIL/chemistry/gasphase_preproc/tmp_kpp4palm
     306! Time                 : Tue Mar  5 11:50:53 2019
     307! Working directory    : /home/forkel-r/palmstuff/work/trunk20190305/UTIL/chemistry/gasphase_preproc/tmp_kpp4palm
    304308! Equation file        : chem_gasphase_mod.kpp
    305309! Output root filename : chem_gasphase_mod
     
    362366!
    363367! File                 : chem_gasphase_mod_Initialize.f90
    364 ! Time                 : Fri Nov 30 13:52:19 2018
    365 ! Working directory    : /home/forkel-r/palmstuff/work/trunk20181130/UTIL/chemistry/gasphase_preproc/tmp_kpp4palm
     368! Time                 : Tue Mar  5 11:50:53 2019
     369! Working directory    : /home/forkel-r/palmstuff/work/trunk20190305/UTIL/chemistry/gasphase_preproc/tmp_kpp4palm
    366370! Equation file        : chem_gasphase_mod.kpp
    367371! Output root filename : chem_gasphase_mod
     
    388392!
    389393! File                 : chem_gasphase_mod_Integrator.f90
    390 ! Time                 : Fri Nov 30 13:52:19 2018
    391 ! Working directory    : /home/forkel-r/palmstuff/work/trunk20181130/UTIL/chemistry/gasphase_preproc/tmp_kpp4palm
     394! Time                 : Tue Mar  5 11:50:53 2019
     395! Working directory    : /home/forkel-r/palmstuff/work/trunk20190305/UTIL/chemistry/gasphase_preproc/tmp_kpp4palm
    392396! Equation file        : chem_gasphase_mod.kpp
    393397! Output root filename : chem_gasphase_mod
     
    446450!
    447451! File                 : chem_gasphase_mod_LinearAlgebra.f90
    448 ! Time                 : Fri Nov 30 13:52:19 2018
    449 ! Working directory    : /home/forkel-r/palmstuff/work/trunk20181130/UTIL/chemistry/gasphase_preproc/tmp_kpp4palm
     452! Time                 : Tue Mar  5 11:50:53 2019
     453! Working directory    : /home/forkel-r/palmstuff/work/trunk20190305/UTIL/chemistry/gasphase_preproc/tmp_kpp4palm
    450454! Equation file        : chem_gasphase_mod.kpp
    451455! Output root filename : chem_gasphase_mod
     
    473477!
    474478! File                 : chem_gasphase_mod_Jacobian.f90
    475 ! Time                 : Fri Nov 30 13:52:19 2018
    476 ! Working directory    : /home/forkel-r/palmstuff/work/trunk20181130/UTIL/chemistry/gasphase_preproc/tmp_kpp4palm
     479! Time                 : Tue Mar  5 11:50:53 2019
     480! Working directory    : /home/forkel-r/palmstuff/work/trunk20190305/UTIL/chemistry/gasphase_preproc/tmp_kpp4palm
    477481! Equation file        : chem_gasphase_mod.kpp
    478482! Output root filename : chem_gasphase_mod
     
    500504!
    501505! File                 : chem_gasphase_mod_Function.f90
    502 ! Time                 : Fri Nov 30 13:52:19 2018
    503 ! Working directory    : /home/forkel-r/palmstuff/work/trunk20181130/UTIL/chemistry/gasphase_preproc/tmp_kpp4palm
     506! Time                 : Tue Mar  5 11:50:53 2019
     507! Working directory    : /home/forkel-r/palmstuff/work/trunk20190305/UTIL/chemistry/gasphase_preproc/tmp_kpp4palm
    504508! Equation file        : chem_gasphase_mod.kpp
    505509! Output root filename : chem_gasphase_mod
     
    529533!
    530534! File                 : chem_gasphase_mod_Rates.f90
    531 ! Time                 : Fri Nov 30 13:52:19 2018
    532 ! Working directory    : /home/forkel-r/palmstuff/work/trunk20181130/UTIL/chemistry/gasphase_preproc/tmp_kpp4palm
     535! Time                 : Tue Mar  5 11:50:53 2019
     536! Working directory    : /home/forkel-r/palmstuff/work/trunk20190305/UTIL/chemistry/gasphase_preproc/tmp_kpp4palm
    533537! Equation file        : chem_gasphase_mod.kpp
    534538! Output root filename : chem_gasphase_mod
     
    555559!
    556560! File                 : chem_gasphase_mod_Util.f90
    557 ! Time                 : Fri Nov 30 13:52:19 2018
    558 ! Working directory    : /home/forkel-r/palmstuff/work/trunk20181130/UTIL/chemistry/gasphase_preproc/tmp_kpp4palm
     561! Time                 : Tue Mar  5 11:50:53 2019
     562! Working directory    : /home/forkel-r/palmstuff/work/trunk20190305/UTIL/chemistry/gasphase_preproc/tmp_kpp4palm
    559563! Equation file        : chem_gasphase_mod.kpp
    560564! Output root filename : chem_gasphase_mod
     
    661665    MODULE PROCEDURE   kppdecomp
    662666  END INTERFACE        kppdecomp
     667 
     668  INTERFACE            get_mechanismname
     669    MODULE PROCEDURE   get_mechanismname
     670  END INTERFACE        get_mechanismname
    663671 
    664672  INTERFACE            chem_gasphase_integrate
     
    22642272  END SUBROUTINE kppdecomp                                           
    22652273 
     2274SUBROUTINE get_mechanismname                                       
     2275                                                                   
     2276  IMPLICIT NONE                                                     
     2277
     2278! Set cs_mech for check with mechanism name from namelist
     2279    cs_mech = 'phstatp'
     2280                                                                   
     2281  RETURN                                                           
     2282END SUBROUTINE get_mechanismname                                   
     2283                                                                   
     2284 
    22662285SUBROUTINE chem_gasphase_integrate (time_step_len, conc, tempi, qvapi, fakti, photo, ierrf, xnacc, xnrej, istatus, l_debug, pe, &
    22672286                     icntrl_i, rcntrl_i)
     
    23492368                                                                   
    23502369  RETURN                                                           
    2351 END SUBROUTINE chem_gasphase_integrate                                       
     2370END SUBROUTINE chem_gasphase_integrate                             
    23522371
    23532372END MODULE chem_gasphase_mod
  • palm/trunk/SOURCE/chem_modules.f90

    r3652 r3780  
    2727! -----------------
    2828! $Id$
     29! added cs_mech
     30!
     31!
     32! 3652 2019-01-07 15:29:59Z forkel
    2933! parameter chem_mechanism added (basit)
    3034!
     
    6670
    6771    USE chem_gasphase_mod,                                                     &   
    68         ONLY: nspec, nvar, spc_names
     72        ONLY: cs_mech, nspec, nvar, spc_names
    6973
    7074    USE control_parameters,                                                    &
     
    7983    IMPLICIT NONE
    8084
     85    PUBLIC cs_mech
    8186    PUBLIC nspec
    8287    PUBLIC nvar
  • palm/trunk/SOURCE/chemistry_model_mod.f90

    r3767 r3780  
    2727! -----------------
    2828! $Id$
     29! Removed READ from unit 10, added CALL get_mechanismname
     30!
     31!
     32! 3767 2019-02-27 08:18:02Z raasch
    2933! unused variable for file index removed from rrd-subroutines parameter list
    3034!
     
    278282
    279283    USE control_parameters,                                                                        &
    280          ONLY:  bc_lr_cyc, bc_ns_cyc, dt_3d, humidity, initializing_actions, message_string,        &
    281          omega, tsc, intermediate_timestep_count, intermediate_timestep_count_max,           &
     284         ONLY:  bc_lr_cyc, bc_ns_cyc, dt_3d, humidity, initializing_actions, message_string,       &
     285         omega, tsc, intermediate_timestep_count, intermediate_timestep_count_max,                 &
    282286         max_pr_user, timestep_scheme, use_prescribed_profile_data, ws_scheme_sca         
    283287
     
    286290
    287291    USE chem_gasphase_mod,                                                                         &
    288          ONLY:  nspec, spc_names, nkppctrl, nmaxfixsteps, t_steps, chem_gasphase_integrate,         &
    289          vl_dim, nvar, nreact,  atol, rtol, nphot, phot_names
     292         ONLY:  atol, chem_gasphase_integrate, get_mechanismname, nkppctrl, nmaxfixsteps,          &
     293         nphot, nreact, nspec, nvar, phot_names, rtol, spc_names, t_steps, vl_dim
    290294
    291295    USE chem_modules
     
    11491153    INTEGER (iwp) ::  lsp_usr      !< running index for user defined chem spcs
    11501154    INTEGER (iwp) ::  lsp          !< running index for chem spcs.
    1151     CHARACTER (LEN=30)       ::  cs_mech,a1,b1,string
    1152 
    1153 
    1154     OPEN (10,FILE="chem_gasphase_mod.f90")   !get the chem_mechanism name from the file.
    1155     READ (10, 100) a1,b1,string
    1156     cs_mech = trim(string(16:))
    1157  100    FORMAT(a)
    1158         CLOSE(10)
    1159 
    11601155!
    11611156!-- check for chemical reactions status
     
    11901185    ENDIF
    11911186!-- check for chemical mechanism used
     1187    CALL get_mechanismname
     1188    WRITE(06,*) 'cs_mech ',cs_mech
     1189    WRITE(06,*) 'chem_mechanism ', chem_mechanism
    11921190    IF (chem_mechanism /= trim(cs_mech) )  THEN
    1193        message_string = 'Incorrect chemical mechanism selected, please check spelling and/or chem_gasphase_mod'
     1191       message_string = 'Incorrect chemistry mechanism selected, check spelling in namelist and/or chem_gasphase_mod'
    11941192       CALL message( 'chem_check_parameters', 'CM0462', 1, 2, 0, 6, 0 )
    11951193    ENDIF
  • palm/trunk/SOURCE/salsa_mod.f90

    r3767 r3780  
    8988988   FORMAT  (/'    Aerosol bin subrange limits (in metres): ',  3(ES10.2E3), / &
    899899              '    Number of size bins for each aerosol subrange: ', 2I3,/     &
    900               '    Aerosol bin limits (in metres): ', *(ES10.2E3))
     900              '    Aerosol bin limits (in metres): ', 9(ES10.2E3))
    90190115  FORMAT   ('    Initial number concentration in bins at the lowest level',  &
    902               ' (#/m**3):', *(ES10.2E3))       
     902              ' (#/m**3):', 9(ES10.2E3))       
    90390313  FORMAT  (/'    Number of chemical components used: ', I1,/                 &
    904904              '       Species: ',7(A6),/                                       &
Note: See TracChangeset for help on using the changeset viewer.