Changeset 20 for palm/trunk


Ignore:
Timestamp:
Feb 26, 2007 12:12:32 AM (15 years ago)
Author:
raasch
Message:

new parameter use_top_fluxes, Bugfix: ddzw dimensioned 1:nzt+1

Location:
palm/trunk/SOURCE
Files:
9 edited

Legend:

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

    r19 r20  
    806806    ENDIF
    807807
    808     IF ( surface_heatflux == 9999999.9 )  constant_heatflux = .FALSE.
    809     IF ( top_heatflux     == 9999999.9 )  THEN
    810        constant_top_heatflux = .FALSE.
    811     ELSE
    812        use_top_fluxes = .TRUE.    ! because this is currently the only choice
    813     ENDIF
     808    IF ( surface_heatflux == 9999999.9 )  constant_heatflux     = .FALSE.
     809    IF ( top_heatflux     == 9999999.9 )  constant_top_heatflux = .FALSE.
    814810
    815811!
  • palm/trunk/SOURCE/diffusion_e.f90

    r19 r20  
    44! Actual revisions:
    55! -----------------
     6! Bugfix: ddzw dimensioned 1:nzt"+1"
    67! Calculation extended for gridpoint nzt
    78!
     
    5152       INTEGER ::  i, j, k
    5253       REAL            ::  dpt_dz, l_stable, phi_m
    53        REAL            ::  ddzu(1:nzt+1), dd2zu(1:nzt), ddzw(1:nzt), &
     54       REAL            ::  ddzu(1:nzt+1), dd2zu(1:nzt), ddzw(1:nzt+1), &
    5455                           l_grid(1:nzt), zu(0:nzt+1)
    5556       REAL, DIMENSION(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) :: diss, tend
     
    166167       INTEGER         ::  i, j, k
    167168       REAL            ::  dpt_dz, l_stable, phi_m
    168        REAL            ::  ddzu(1:nzt+1), dd2zu(1:nzt), ddzw(1:nzt), &
     169       REAL            ::  ddzu(1:nzt+1), dd2zu(1:nzt), ddzw(1:nzt+1), &
    169170                           l_grid(1:nzt), zu(0:nzt+1)
    170171       REAL, DIMENSION(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) ::  diss, tend
  • palm/trunk/SOURCE/diffusion_s.f90

    r19 r20  
    44! Actual revisions:
    55! -----------------
     6! Bugfix: ddzw dimensioned 1:nzt"+1"
    67! Calculation extended for gridpoint nzt, fluxes can be given at top,
    7 ! +s_flux_t in parameter list, s_flux renamed s_flux_b
     8! +s_flux_t in parameter list, s_flux renamed s_flux_b,
    89!
    910! Former revisions:
     
    5051       INTEGER ::  i, j, k
    5152       REAL    ::  vertical_gridspace
    52        REAL    ::  ddzu(1:nzt+1), ddzw(1:nzt)
     53       REAL    ::  ddzu(1:nzt+1), ddzw(1:nzt+1)
    5354       REAL    ::  tend(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1)
    5455       REAL, DIMENSION(:,:),   POINTER ::  s_flux_b, s_flux_t
     
    134135                tend(k,j,i) = tend(k,j,i)                                     &
    135136                                       + ( - s_flux_t(j,i)                    &
    136                                          - 0.5 * ( kh(k-1,j,i)+kh(k,j,i) )    &
    137                                                * ( s(k,j,i)-s(k-1,j,i) )      &
    138                                                * ddzu(k)                      &
     137                                           - 0.5 * ( kh(k-1,j,i)+kh(k,j,i) )  &
     138                                                 * ( s(k,j,i)-s(k-1,j,i) )    &
     139                                                 * ddzu(k)                    &
    139140                                         ) * ddzw(k)
    140141
     
    161162       INTEGER ::  i, j, k
    162163       REAL    ::  vertical_gridspace
    163        REAL    ::  ddzu(1:nzt+1), ddzw(1:nzt)
     164       REAL    ::  ddzu(1:nzt+1), ddzw(1:nzt+1)
    164165       REAL    ::  tend(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1)
    165166       REAL, DIMENSION(:,:),   POINTER ::  s_flux_b, s_flux_t
  • palm/trunk/SOURCE/diffusion_u.f90

    r4 r20  
    44! Actual revisions:
    55! -----------------
    6 !
     6! Bugfix: ddzw dimensioned 1:nzt"+1"
    77!
    88! Former revisions:
     
    5252       INTEGER ::  i, j, k
    5353       REAL    ::  kmym_x, kmym_y, kmyp_x, kmyp_y, kmzm, kmzp, usvs
    54        REAL    ::  ddzu(1:nzt+1), ddzw(1:nzt), km_damp_y(nys-1:nyn+1)
     54       REAL    ::  ddzu(1:nzt+1), ddzw(1:nzt+1), km_damp_y(nys-1:nyn+1)
    5555       REAL    ::  z0(nys-1:nyn+1,nxl-1:nxr+1)
    5656       REAL    ::  tend(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1)
     
    201201       INTEGER ::  i, j, k
    202202       REAL    ::  kmym_x, kmym_y, kmyp_x, kmyp_y, kmzm, kmzp, usvs
    203        REAL    ::  ddzu(1:nzt+1), ddzw(1:nzt), km_damp_y(nys-1:nyn+1)
     203       REAL    ::  ddzu(1:nzt+1), ddzw(1:nzt+1), km_damp_y(nys-1:nyn+1)
    204204       REAL    ::  z0(nys-1:nyn+1,nxl-1:nxr+1)
    205205       REAL    ::  tend(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1)
  • palm/trunk/SOURCE/diffusion_v.f90

    r4 r20  
    44! Actual revisions:
    55! -----------------
    6 !
     6! Bugfix: ddzw dimensioned 1:nzt"+1"
    77!
    88! Former revisions:
     
    5252       INTEGER ::  i, j, k
    5353       REAL    ::  kmxm_x, kmxm_y, kmxp_x, kmxp_y, kmzm, kmzp, vsus
    54        REAL    ::  ddzu(1:nzt+1), ddzw(1:nzt), km_damp_x(nxl-1:nxr+1)
     54       REAL    ::  ddzu(1:nzt+1), ddzw(1:nzt+1), km_damp_x(nxl-1:nxr+1)
    5555       REAL    ::  z0(nys-1:nyn+1,nxl-1:nxr+1)
    5656       REAL    ::  tend(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1)
     
    201201       INTEGER ::  i, j, k
    202202       REAL    ::  kmxm_x, kmxm_y, kmxp_x, kmxp_y, kmzm, kmzp, vsus
    203        REAL    ::  ddzu(1:nzt+1), ddzw(1:nzt), km_damp_x(nxl-1:nxr+1)
     203       REAL    ::  ddzu(1:nzt+1), ddzw(1:nzt+1), km_damp_x(nxl-1:nxr+1)
    204204       REAL    ::  z0(nys-1:nyn+1,nxl-1:nxr+1)
    205205       REAL    ::  tend(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1)
  • palm/trunk/SOURCE/diffusion_w.f90

    r4 r20  
    44! Actual revisions:
    55! -----------------
    6 !
     6! Bugfix: ddzw dimensioned 1:nzt"+1"
    77!
    88! Former revisions:
     
    5151       REAL    ::  kmxm_x, kmxm_z, kmxp_x, kmxp_z, kmym_y, kmym_z, kmyp_y, &
    5252                   kmyp_z, wsus, wsvs
    53        REAL    ::  ddzu(1:nzt+1), ddzw(1:nzt), km_damp_x(nxl-1:nxr+1), &
     53       REAL    ::  ddzu(1:nzt+1), ddzw(1:nzt+1), km_damp_x(nxl-1:nxr+1), &
    5454                   km_damp_y(nys-1:nyn+1)
    5555       REAL    ::  z0(nys-1:nyn+1,nxl-1:nxr+1)
     
    199199       REAL    ::  kmxm_x, kmxm_z, kmxp_x, kmxp_z, kmym_y, kmym_z, kmyp_y, &
    200200                   kmyp_z, wsus, wsvs
    201        REAL    ::  ddzu(1:nzt+1), ddzw(1:nzt), km_damp_x(nxl-1:nxr+1), &
     201       REAL    ::  ddzu(1:nzt+1), ddzw(1:nzt+1), km_damp_x(nxl-1:nxr+1), &
    202202                   km_damp_y(nys-1:nyn+1)
    203203       REAL    ::  z0(nys-1:nyn+1,nxl-1:nxr+1)
  • palm/trunk/SOURCE/parin.f90

    r19 r20  
    44! Actual revisions:
    55! -----------------
    6 ! +top_heatflux in inipar
     6! +top_heatflux, use_top_fluxes in inipar
    77!
    88! Former revisions:
     
    7474                       ug_vertical_gradient, ug_vertical_gradient_level, &
    7575                       ups_limit_e, ups_limit_pt, ups_limit_u, ups_limit_v, &
    76                        ups_limit_w, use_surface_fluxes, use_ug_for_galilei_tr, &
    77                        use_upstream_for_tke, vg_surface, vg_vertical_gradient, &
     76                       ups_limit_w, use_surface_fluxes, use_top_fluxes, &
     77                       use_ug_for_galilei_tr, use_upstream_for_tke, &
     78                       vg_surface, vg_vertical_gradient, &
    7879                       vg_vertical_gradient_level, wall_adjustment, &
    7980                       wall_heatflux
  • palm/trunk/SOURCE/read_var_list.f90

    r19 r20  
    44! Actual revisions:
    55! -----------------
    6 ! +top_heatflux
     6! +top_heatflux, use_top_fluxes
    77!
    88! Former revisions:
     
    372372          CASE ( 'use_surface_fluxes' )
    373373             READ ( 13 )  use_surface_fluxes
     374          CASE ( 'use_top_fluxes' )
     375             READ ( 13 )  use_top_fluxes
    374376          CASE ( 'use_ug_for_galilei_tr' )
    375377             READ ( 13 )  use_ug_for_galilei_tr
  • palm/trunk/SOURCE/write_var_list.f90

    r19 r20  
    44! Actual revisions:
    55! -----------------
    6 ! +top_heatflux
     6! +top_heatflux, use_top_fluxes
    77!
    88! Former revisions:
     
    338338    WRITE ( 14 )  'use_surface_fluxes            '
    339339    WRITE ( 14 )  use_surface_fluxes
     340    WRITE ( 14 )  'use_top_fluxes                '
     341    WRITE ( 14 )  use_top_fluxes
    340342    WRITE ( 14 )  'use_ug_for_galilei_tr         '
    341343    WRITE ( 14 )  use_ug_for_galilei_tr
Note: See TracChangeset for help on using the changeset viewer.