Ignore:
Timestamp:
Dec 14, 2017 5:12:51 PM (4 years ago)
Author:
kanani
Message:

Merge of branch palm4u into trunk

Location:
palm/trunk
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk

  • palm/trunk/SOURCE

  • palm/trunk/SOURCE/sum_up_3d_data.f90

    r2292 r2696  
    11!> @file sum_up_3d_data.f90
    22!------------------------------------------------------------------------------!
    3 ! This file is part of PALM.
     3! This file is part of the PALM model system.
    44!
    55! PALM is free software: you can redistribute it and/or modify it under the
     
    2525! -----------------
    2626! $Id$
     27! - Implementation of uv exposure model (FK)
     28! - output of diss_av, kh_av, km_av (turbulence_closure_mod) (TG)
     29! - Implementation of chemistry module (FK)
     30! - Workaround for sum-up usm arrays in case of restart runs, to avoid program
     31!   crash (MS)
     32!
     33! 2292 2017-06-20 09:51:42Z schwenkel
    2734! Implementation of new microphysic scheme: cloud_scheme = 'morrison'
    2835! includes two more prognostic equations for cloud drop concentration (nc) 
     
    144151
    145152    USE averaging,                                                             &
    146         ONLY:  e_av, lpt_av, lwp_av, nc_av, nr_av, ol_av, p_av, pc_av, pr_av,  &
    147                prr_av, precipitation_rate_av, pt_av, q_av, qc_av, ql_av,       &
    148                ql_c_av, ql_v_av, ql_vp_av, qr_av, qsws_av, qv_av, rho_ocean_av,&
    149                s_av, sa_av, shf_av, ssws_av, ts_av, u_av, us_av, v_av, vpt_av, &
    150                w_av, z0_av, z0h_av, z0q_av
     153        ONLY:  diss_av, e_av, kh_av, km_av, lpt_av, lwp_av, nc_av, nr_av,      &
     154               ol_av, p_av, pc_av, pr_av, prr_av, precipitation_rate_av, pt_av,&
     155               q_av, qc_av, ql_av, ql_c_av, ql_v_av, ql_vp_av, qr_av, qsws_av, &
     156               qv_av, rho_ocean_av, s_av, sa_av, shf_av, ssws_av, ts_av, u_av, &
     157               us_av, v_av, vpt_av, w_av, z0_av, z0h_av, z0q_av
     158#if defined( __chem )
     159    USE chemistry_model_mod,                                                   &
     160        ONLY:  chem_3d_data_averaging, chem_integrate, chem_species, nspec                                   
     161#endif
    151162
    152163    USE cloud_parameters,                                                      &
     
    154165
    155166    USE control_parameters,                                                    &
    156         ONLY:  average_count_3d, cloud_physics, doav, doav_n, land_surface,    &
    157                rho_surface, urban_surface, varnamelength
     167        ONLY:  air_chemistry, average_count_3d, cloud_physics, doav, doav_n,   &
     168               land_surface, rho_surface, urban_surface, uv_exposure,          &
     169               varnamelength
    158170
    159171    USE cpulog,                                                                &
     
    177189        ONLY:  surf_def_h, surf_lsm_h, surf_usm_h
    178190
     191    USE turbulence_closure_mod,                                                &
     192        ONLY:  tcm_3d_data_averaging
     193
    179194    USE urban_surface_mod,                                                     &
    180195        ONLY:  usm_average_3d_data
     196
     197    USE uv_exposure_model_mod,                                                &
     198        ONLY:  uvem_3d_data_averaging
    181199
    182200
     
    436454
    437455!
     456!--             Turbulence closure module
     457                CALL tcm_3d_data_averaging( 'allocate', doav(ii) )
     458
     459!
    438460!--             Land surface quantity
    439461                IF ( land_surface )  THEN
     
    445467                IF ( radiation )  THEN
    446468                   CALL radiation_3d_data_averaging( 'allocate', doav(ii) )
     469                ENDIF
     470
     471!
     472!--             Chemical quantity                                           
     473#if defined( __chem )               
     474                IF ( air_chemistry  .AND.  trimvar(1:3) == 'kc_')  THEN
     475                   CALL chem_3d_data_averaging( 'allocate', doav(ii) )
     476                ENDIF
     477#endif
     478
     479!
     480!--             UV exposure quantity
     481                IF ( uv_exposure  .AND.  trimvar(1:5) == 'uvem_')  THEN
     482                   CALL uvem_3d_data_averaging( 'allocate', doav(ii) )
    447483                ENDIF
    448484
     
    889925             ENDDO
    890926!             
    891 !--       Block of urban surface model outputs
     927!--       Block of urban surface model outputs.
     928!--       In case of urban surface variables it should be always checked
     929!--       if respective arrays are allocated, at least in case of a restart
     930!--       run, as usm arrays are not read from file at the moment.
    892931          CASE ( 'usm_output' )
     932             CALL usm_average_3d_data( 'allocate', doav(ii) )
    893933             CALL usm_average_3d_data( 'sum', doav(ii) )
    894934
    895935          CASE DEFAULT
     936!
     937!--          Turbulence closure module
     938             CALL tcm_3d_data_averaging( 'sum', doav(ii) )
     939
    896940!
    897941!--          Land surface quantity
     
    907951
    908952!
     953!--          Chemical quantity
     954#if defined( __chem )               
     955             IF ( air_chemistry  .AND.  trimvar(1:3) == 'kc_')  THEN
     956                CALL chem_3d_data_averaging( 'sum',doav(ii) )
     957             ENDIF
     958#endif
     959
     960!
     961!--          UV exposure quantity
     962             IF ( uv_exposure )  THEN
     963                CALL uvem_3d_data_averaging( 'sum', doav(ii) )
     964             ENDIF
     965
     966!
    909967!--          User-defined quantity
    910968             CALL user_3d_data_averaging( 'sum', doav(ii) )
Note: See TracChangeset for help on using the changeset viewer.