Ignore:
Timestamp:
Oct 26, 2015 4:17:44 PM (9 years ago)
Author:
maronga
Message:

various bugfixes and modifications of the atmosphere-land-surface-radiation interaction. Completely re-written routine to calculate surface fluxes (surface_layer_fluxes.f90) that replaces prandtl_fluxes. Minor formatting corrections and renamings

File:
1 edited

Legend:

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

    r1683 r1691  
    1414! PALM. If not, see <http://www.gnu.org/licenses/>.
    1515!
    16 ! Copyright 1997-2014 Leibniz Universitaet Hannover
     16! Copyright 1997-2015 Leibniz Universitaet Hannover
    1717!--------------------------------------------------------------------------------!
    1818!
    1919! Current revisions:
    2020! ------------------
    21 !
     21! Added output of radiative heating rates for RRTMG. Corrected output of
     22! radiative fluxes
    2223!
    2324! Former revisions:
     
    116117#if defined( __netcdf )
    117118
    118     USE arrays_3d,                                                              &
     119    USE arrays_3d,                                                             &
    119120        ONLY:  zu, zw
    120121
    121     USE constants,                                                              &
     122    USE constants,                                                             &
    122123        ONLY:  pi
    123124
    124     USE control_parameters,                                                     &
    125         ONLY:  averaging_interval, averaging_interval_pr, averaging_interval_sp,&
    126         data_output_pr,  domask,  dopr_n,dopr_time_count, dopts_time_count,     &
    127         dots_time_count, dosp_time_count, do2d, do2d_xz_time_count, do3d,       &
    128         do2d_yz_time_count, mask_size, do2d_xy_time_count, do3d_time_count,     &
    129         domask_time_count, mask_i_global, mask_j_global,mask_k_global,          &
    130         message_string, mid, netcdf_data_format, netcdf_precision, ntdim_2d_xy, &
    131         ntdim_2d_xz, ntdim_2d_yz, ntdim_3d, nz_do3d, prt_time_count,            &
    132         run_description_header, section, simulated_time, topography
    133 
    134     USE grid_variables,                                                         &
     125    USE control_parameters,                                                    &
     126        ONLY:  averaging_interval, averaging_interval_pr,                      &
     127               averaging_interval_sp, data_output_pr,  domask,  dopr_n,        &
     128               dopr_time_count, dopts_time_count, dots_time_count,             &
     129               dosp_time_count, do2d, do2d_xz_time_count, do3d,                &
     130               do2d_yz_time_count, mask_size, do2d_xy_time_count,              &
     131               do3d_time_count, domask_time_count, mask_i_global,              &
     132               mask_j_global, mask_k_global, message_string, mid,              &
     133               netcdf_data_format, netcdf_precision, ntdim_2d_xy,              &
     134               ntdim_2d_xz, ntdim_2d_yz, ntdim_3d, nz_do3d, prt_time_count,    &
     135              run_description_header, section, simulated_time, topography
     136
     137    USE grid_variables,                                                        &
    135138        ONLY:  dx, dy, zu_s_inner, zw_w_inner
    136139
    137     USE indices,                                                                &
     140    USE indices,                                                               &
    138141        ONLY:  nx, ny, nz ,nzb, nzt
    139142
     
    150153    USE pegrid
    151154
    152     USE particle_attributes,                                                    &
     155    USE particle_attributes,                                                   &
    153156        ONLY:  maximum_number_of_particles, number_of_particle_groups
    154157
    155     USE profil_parameter,                                                       &
    156         ONLY:  crmax, cross_profiles, dopr_index,profile_columns, profile_rows
    157 
    158     USE radiation_model_mod,                                                    &
    159         ONLY:  rad_lw_in, rad_lw_out, rad_sw_in, rad_sw_out
    160 
    161     USE spectrum,                                                               &
     158    USE profil_parameter,                                                      &
     159        ONLY:  crmax, cross_profiles, dopr_index, profile_columns, profile_rows
     160
     161    USE radiation_model_mod,                                                   &
     162        ONLY:  rad_lw_in, rad_lw_out, rad_lw_cs_hr, rad_lw_hr,                 &
     163               rad_sw_in, rad_sw_out, rad_sw_cs_hr, rad_sw_hr
     164
     165
     166    USE spectrum,                                                              &
    162167        ONLY:  comp_spectra_level, data_output_sp, spectra_direction
    163168
    164     USE statistics,                                                             &
     169    USE statistics,                                                            &
    165170        ONLY:  hom, statistic_regions
    166171
     
    531536!
    532537!--             Most variables are defined on the scalar grid
    533                 CASE ( 'e', 'lpt', 'nr', 'p', 'pc', 'pr', 'prr', 'pt', 'q',&
    534                        'qc', 'ql', 'ql_c', 'ql_v', 'ql_vp', 'qr', 'qv',    &
    535                        'rho', 's', 'sa', 'vpt', 'rad_lw_in', 'rad_lw_out', &
    536                        'rad_sw_in', 'rad_sw_out' )
     538                CASE ( 'e', 'lpt', 'nr', 'p', 'pc', 'pr', 'prr', 'pt', 'q',    &
     539                       'qc', 'ql', 'ql_c', 'ql_v', 'ql_vp', 'qr', 'qv',        &
     540                       'rad_lw_cs_hr', 'rad_lw_hr', 'rad_sw_cs_hr',            &
     541                       'rad_sw_hr', 'rho', 's', 'sa', 'vpt' )
    537542
    538543                   grid_x = 'x'
     
    555560!
    556561!--             w grid
    557                 CASE ( 'w' )
     562                CASE ( 'rad_lw_in', 'rad_lw_out', 'rad_sw_in', 'rad_sw_out',   &
     563                       'w' )
    558564
    559565                   grid_x = 'x'
     
    562568!
    563569!--             soil grid
    564                 CASE ( 't_soil', 'm_soil' )
     570                CASE ( 'm_soil', 't_soil' )
    565571
    566572                   grid_x = 'x'
     
    11121118!
    11131119!--             Most variables are defined on the scalar grid
    1114                 CASE ( 'e', 'lpt', 'nr', 'p', 'pc', 'pr', 'prr', 'pt', 'q',   &
    1115                        'qc', 'ql', 'ql_c', 'ql_v', 'ql_vp', 'qr', 'qv', 'rho',&
    1116                        's', 'sa', 'vpt' , 'rad_lw_in', 'rad_lw_out',          &
    1117                        'rad_sw_in', 'rad_sw_out' )
     1120                CASE ( 'e', 'lpt', 'nr', 'p', 'pc', 'pr', 'prr', 'pt', 'q',    &
     1121                       'qc', 'ql', 'ql_c', 'ql_v', 'ql_vp', 'qr', 'qv', 'rho', &
     1122                       's', 'sa', 'vpt' , 'rad_lw_cs_hr', 'rad_lw_hr',         &
     1123                       'rad_sw_cs_hr', 'rad_sw_hr' )
    11181124
    11191125                   grid_x = 'x'
     
    11361142!
    11371143!--             w grid
    1138                 CASE ( 'w' )
     1144                CASE ( 'rad_lw_in', 'rad_lw_out', 'rad_sw_in', 'rad_sw_out',   &
     1145                       'w' )
    11391146
    11401147                   grid_x = 'x'
     
    11431150!
    11441151!--             soil grid
    1145                 CASE ( 't_soil', 'm_soil' )
     1152                CASE ( 'm_soil', 't_soil' )
    11461153
    11471154                   grid_x = 'x'
     
    17651772!
    17661773!--                   Most variables are defined on the zu grid
    1767                       CASE ( 'e_xy', 'lpt_xy', 'nr_xy', 'p_xy', 'pc_xy', 'pr_xy',&
    1768                              'prr_xy', 'pt_xy', 'q_xy', 'qc_xy', 'ql_xy',        &
    1769                              'ql_c_xy', 'ql_v_xy', 'ql_vp_xy', 'qr_xy', 'qv_xy', &
    1770                              'rho_xy', 's_xy', 'sa_xy', 'vpt_xy', 'rad_lw_in_xy',&
    1771                              'rad_lw_out_xy', 'rad_sw_in_xy', 'rad_sw_out_xy' )
     1774                      CASE ( 'e_xy', 'lpt_xy', 'nr_xy', 'p_xy', 'pc_xy',       &
     1775                             'pr_xy', 'prr_xy', 'pt_xy', 'q_xy', 'qc_xy',      &
     1776                             'ql_xy', 'ql_c_xy', 'ql_v_xy', 'ql_vp_xy',        &
     1777                             'qr_xy', 'qv_xy', 'rad_lw_cs_hr_xy',              &
     1778                             'rad_lw_hr_xy', 'rad_sw_cs_hr_xy', 'rad_sw_hr_xy',&
     1779                             'rho_xy', 's_xy', 'sa_xy', 'vpt_xy' )
    17721780
    17731781                         grid_x = 'x'
     
    17901798!
    17911799!--                   w grid
    1792                       CASE ( 'w_xy' )
     1800                      CASE ( 'rad_lw_in_xy', 'rad_lw_out_xy', 'rad_sw_in_xy',  &
     1801                             'rad_sw_out_xy' , 'w_xy' )
    17931802
    17941803                         grid_x = 'x'
     
    17971806!
    17981807!--                   soil grid
    1799                       CASE ( 't_soil_xy', 'm_soil_xy' )
     1808                      CASE ( 'm_soil_xy', 't_soil_xy' )
    18001809                         grid_x = 'x'
    18011810                         grid_y = 'y'
     
    24522461                          'prr_xz', 'pt_xz', 'q_xz', 'qc_xz', 'ql_xz',         &
    24532462                          'ql_c_xz', 'ql_v_xz', 'ql_vp_xz', 'qr_xz', 'qv_xz',  &
    2454                           'rho_xz', 's_xz', 'sa_xz', 'vpt_xz' , 'rad_lw_in_xz',&
    2455                              'rad_lw_out_xz', 'rad_sw_in_xz', 'rad_sw_out_xz' )
     2463                          'rad_lw_cs_hr_xz', 'rad_lw_hr_xz',                   &
     2464                          'rad_sw_cs_hr_xz', 'rad_sw_hr_xz''rho_xz', 's_xz',   &
     2465                          'sa_xz', 'vpt_xz' )
    24562466
    24572467                      grid_x = 'x'
     
    24742484!
    24752485!--                w grid
    2476                    CASE ( 'w_xz' )
     2486                   CASE ( 'rad_lw_in_xz', 'rad_lw_out_xz', 'rad_sw_in_xz',     &
     2487                          'rad_sw_out_xz', 'w_xz' )
    24772488
    24782489                      grid_x = 'x'
     
    24822493!
    24832494!--                soil grid
    2484                    CASE ( 't_soil_xz', 'm_soil_xz' )
     2495                   CASE ( 'm_soil_xz', 't_soil_xz' )
    24852496
    24862497                      grid_x = 'x'
     
    31263137!
    31273138!--                Most variables are defined on the zu grid
    3128                    CASE ( 'e_yz', 'lpt_yz', 'nr_yz', 'p_yz', 'pc_yz', 'pr_yz',&
    3129                           'prr_yz', 'pt_yz', 'q_yz', 'qc_yz', 'ql_yz',        &
    3130                           'ql_c_yz', 'ql_v_yz', 'ql_vp_yz', 'qr_yz', 'qv_yz', &
    3131                           'rho_yz', 's_yz', 'sa_yz', 'vpt_yz', 'rad_lw_in_yz',&
    3132                              'rad_lw_out_yz', 'rad_sw_in_yz', 'rad_sw_out_yz' )
     3139                   CASE ( 'e_yz', 'lpt_yz', 'nr_yz', 'p_yz', 'pc_yz', 'pr_yz', &
     3140                          'prr_yz', 'pt_yz', 'q_yz', 'qc_yz', 'ql_yz',         &
     3141                          'ql_c_yz', 'ql_v_yz', 'ql_vp_yz', 'qr_yz', 'qv_yz',  &
     3142                          'rad_lw_cs_hr_yz', 'rad_lw_hr_yz',                   &
     3143                          'rad_sw_cs_hr_yz', 'rad_sw_hr_yz''rho_yz', 's_yz',   &
     3144                          'sa_yz', 'vpt_yz' )
    31333145
    31343146                      grid_x = 'x'
     
    31513163!
    31523164!--                w grid
    3153                    CASE ( 'w_yz' )
     3165                   CASE ( 'rad_lw_in_yz', 'rad_lw_out_yz', 'rad_sw_in_yz',     &
     3166                          'rad_sw_out_yz', 'w_yz' )
    31543167
    31553168                      grid_x = 'x'
     
    31583171!
    31593172!--                soil grid
    3160                    CASE ( 't_soil_yz', 'm_soil_yz' )
     3173                   CASE ( 'm_soil_yz', 't_soil_yz' )
    31613174
    31623175                      grid_x = 'x'
Note: See TracChangeset for help on using the changeset viewer.