Ignore:
Timestamp:
Jun 18, 2019 10:32:41 AM (5 years ago)
Author:
suehring
Message:

diagnostic output: Modularize diagnostic output, rename subroutines; formatting adjustments; allocate arrays only when required; add output of uu, vv, ww to enable variance calculation via temporal EC method; radiation: bugfix in masked data output; flow_statistics: Correct conversion to kinematic vertical scalar fluxes in case of pw-scheme and statistic regions

File:
1 edited

Legend:

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

    r4029 r4039  
    2525! -----------------
    2626! $Id$
     27! Rename subroutines in module for diagnostic quantities
     28!
     29! 4029 2019-06-14 14:04:35Z raasch
    2730! netcdf variable NF90_NOFILL is used as argument instead of "1" in calls to NF90_DEF_VAR_FILL
    2831!
     
    698701               skip_time_do3d, topography, num_leg, num_var_fl,                &
    699702               urban_surface
     703
     704    USE diagnostic_output_quantities_mod,                                      &
     705        ONLY:  doq_define_netcdf_grid           
    700706
    701707    USE grid_variables,                                                        &
     
    11241130                CASE ( 'e', 'nc', 'nr', 'p', 'pc', 'pr', 'prr',                &
    11251131                       'q', 'qc', 'ql', 'ql_c', 'ql_v', 'ql_vp', 'qr', 'qv',   &
    1126                        's', 'theta', 'thetal', 'thetav', 'ti' )
     1132                       's', 'theta', 'thetal', 'thetav' )
    11271133
    11281134                   grid_x = 'x'
     
    11691175                   ENDIF
    11701176
     1177                   IF ( .NOT. found )                                          &
     1178                      CALL doq_define_netcdf_grid( domask(mid,av,i), found,    &
     1179                                                   grid_x, grid_y, grid_z )
     1180                   
    11711181                   IF ( .NOT. found  .AND.  gust_module_enabled )  THEN
    11721182                      CALL gust_define_netcdf_grid( domask(mid,av,i), found,   &
     
    18701880                CASE ( 'e', 'nc', 'nr', 'p', 'pc', 'pr', 'prr',   &
    18711881                       'q', 'qc', 'ql', 'ql_c', 'ql_v', 'ql_vp', 'qr', 'qv',   &
    1872                        's', 'theta', 'thetal', 'thetav', 'ti' )
     1882                       's', 'theta', 'thetal', 'thetav' )
    18731883
    18741884                   grid_x = 'x'
     
    19741984                                                    grid_y, grid_z )
    19751985                   ENDIF
     1986                   
     1987                   IF ( .NOT. found )                                          &
     1988                      CALL doq_define_netcdf_grid( do3d(av,i), found, grid_x,  &
     1989                                                   grid_y, grid_z        )
    19761990                                                 
    19771991                   IF ( .NOT. found )  THEN
     
    28582872                             'ql_vp_xy', 'qr_xy', 'qv_xy',                     &
    28592873                             's_xy',                                           &
    2860                              'theta_xy', 'thetal_xy', 'thetav_xy', 'ti_xy' )
     2874                             'theta_xy', 'thetal_xy', 'thetav_xy' )
    28612875
    28622876                         grid_x = 'x'
     
    29442958                                                          grid_z )
    29452959                         ENDIF
    2946 
    2947 
     2960                         
     2961                         IF ( .NOT. found )                                    &
     2962                            CALL doq_define_netcdf_grid(                       &
     2963                                                    do2d(av,i), found, grid_x, &
     2964                                                    grid_y, grid_z              )
    29482965!
    29492966!--                      Check for user-defined quantities
     
    37623779                          'ql_xz', 'ql_c_xz', 'ql_v_xz', 'ql_vp_xz', 'qr_xz',  &
    37633780                          'qv_xz', 's_xz',                                     &
    3764                           'theta_xz', 'thetal_xz', 'thetav_xz', 'ti_xz' )
     3781                          'theta_xz', 'thetal_xz', 'thetav_xz' )
    37653782
    37663783                      grid_x = 'x'
     
    38373854                                                       grid_z )
    38383855                      ENDIF
     3856                     
     3857                      IF ( .NOT. found )                                       &
     3858                         CALL doq_define_netcdf_grid( do2d(av,i), found,       &
     3859                                                      grid_x, grid_y, grid_z )
    38393860
    38403861!
     
    46954716                                                       grid_z )
    46964717                      ENDIF
     4718                     
     4719                      IF ( .NOT. found )                                       &
     4720                         CALL doq_define_netcdf_grid(                          &
     4721                                                    do2d(av,i), found, grid_x, &
     4722                                                    grid_y, grid_z           )
    46974723!
    46984724!--                   Check for user-defined quantities
     
    60006026                       'q', 'qc', 'ql', 'ql_c', 'ql_v', 'ql_vp', 'qr', 'qv',   &
    60016027                       'rho_sea_water', 's', 'sa', &
    6002                        'theta', 'thetal', 'thetav', 'ti', 'u', 'v' )
     6028                       'theta', 'thetal', 'thetav', 'u', 'v' )
    60036029
    60046030                   grid_z = 'zu'
Note: See TracChangeset for help on using the changeset viewer.