Ignore:
Timestamp:
Mar 26, 2012 2:21:38 PM (13 years ago)
Author:
suehring
Message:

Formatting adjustments.

File:
1 edited

Legend:

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

    r861 r862  
    2121! boundary_flags and logicals steering the degradation are removed.
    2222! Empty SUBROUTINE local_diss removed.
     23!
     24! Further formatting adjustments.
    2325!
    2426! Former revisions:
     
    274276       IMPLICIT NONE
    275277
    276        INTEGER ::  i, i_omp, j, k, tn, k_ppp, k_pp, k_mm
    277        REAL    ::  flux_d, diss_d, u_comp, v_comp, div
     278       INTEGER ::  i, i_omp, j, k, k_mm, k_pp, k_ppp,  tn
     279       REAL    ::  diss_d, div, flux_d, u_comp, v_comp
    278280       REAL, DIMENSION(:,:,:), POINTER    :: sk
    279        REAL, DIMENSION(nzb:nzt+1)         :: flux_t, diss_t, flux_r, diss_r,  &
    280                                              flux_n, diss_n
    281        REAL, DIMENSION(nzb+1:nzt,0:threads_per_task-1) :: swap_flux_y_local,  &
    282                                                           swap_diss_y_local
     281       REAL, DIMENSION(nzb:nzt+1)         :: diss_n, diss_r, diss_t, flux_n,  &
     282                                             flux_r, flux_t
     283       REAL, DIMENSION(nzb+1:nzt,0:threads_per_task-1) :: swap_diss_y_local,  &
     284                                                          swap_flux_y_local
    283285       REAL, DIMENSION(nzb+1:nzt,nys:nyn,0:threads_per_task-1) ::             &
    284                                                           swap_flux_x_local,  &
    285                                                           swap_diss_x_local
     286                                                          swap_diss_x_local,  &
     287                                                          swap_flux_x_local
    286288       CHARACTER (LEN = *), INTENT(IN)    :: sk_char
    287289
     
    671673       IMPLICIT NONE
    672674
    673        INTEGER ::  i, i_omp, j, k, tn, k_ppp, k_pp, k_mm
    674        REAL    ::  gu, gv, flux_d, diss_d, u_comp_l, v_comp, w_comp, div
    675        REAL, DIMENSION(nzb:nzt+1) ::  flux_t, diss_t, flux_r, diss_r,        &
    676                                       flux_n, diss_n, u_comp
     675       INTEGER ::  i, i_omp, j, k, k_mm, k_pp, k_ppp, tn
     676       REAL    ::  diss_d, div, flux_d, gu, gv, u_comp_l, v_comp, w_comp
     677       REAL, DIMENSION(nzb:nzt+1) :: diss_n, diss_r, diss_t, flux_n, flux_r,  &
     678                                     flux_t, u_comp
    677679
    678680       gu = 2.0 * u_gtrans
     
    10581060       IMPLICIT NONE
    10591061
    1060        INTEGER  ::  i, i_omp, j, k, tn, k_ppp, k_pp, k_mm
    1061        REAL     ::  gu, gv, flux_d, diss_d, u_comp, v_comp_l, w_comp, div
    1062        REAL, DIMENSION(nzb:nzt+1)  ::  flux_t, diss_t, flux_n,                &
    1063                                        diss_n, flux_r, diss_r, v_comp
     1062       INTEGER  ::  i, i_omp, j, k, k_mm, k_pp, k_ppp, tn
     1063       REAL     ::  diss_d, div, flux_d, gu, gv, u_comp, v_comp_l, w_comp
     1064       REAL, DIMENSION(nzb:nzt+1)  :: diss_n, diss_r, diss_t, flux_n, flux_r,  &
     1065                                      flux_t, v_comp
    10641066
    10651067       gu = 2.0 * u_gtrans
     
    14511453       IMPLICIT NONE
    14521454
    1453        INTEGER ::  i, i_omp, j, k, tn, k_ppp, k_pp, k_mm
    1454        REAL    ::  gu, gv, flux_d, diss_d, u_comp, v_comp, w_comp, div
    1455        REAL, DIMENSION(nzb:nzt+1)  ::  flux_t, diss_t, flux_r, diss_r, flux_n, &
    1456                                        diss_n
     1455       INTEGER ::  i, i_omp, j, k, k_mm, k_pp, k_ppp, tn
     1456       REAL    ::  diss_d, div, flux_d, gu, gv, u_comp, v_comp, w_comp
     1457       REAL, DIMENSION(nzb:nzt+1)  :: diss_n, diss_r, diss_t, flux_n, flux_r, &
     1458                                      flux_t
    14571459
    14581460       gu = 2.0 * u_gtrans
     
    14651467          DO  k = nzb+1, nzb_max
    14661468
    1467              v_comp      = v(k+1,j,i) + v(k,j,i) - gv
     1469             v_comp         = v(k+1,j,i) + v(k,j,i) - gv
    14681470             flux_s_w(k,tn) = v_comp * (                                      &
    14691471                         ( 37.0 * IBITS(wall_flags_0(k,j,i),32,1) * adv_mom_5 &
     
    15191521          DO  k = nzb+1, nzb_max
    15201522
    1521              u_comp         = u(k+1,j,i) + u(k,j,i) - gu
     1523             u_comp           = u(k+1,j,i) + u(k,j,i) - gu
    15221524             flux_l_w(k,j,tn) = u_comp * (                                     &
    15231525                          ( 37.0 * IBITS(wall_flags_0(k,j,i),29,1) * adv_mom_5 &
     
    15671569
    15681570       ENDIF
    1569 
    1570        flux_t(0) = 0.0
    1571        diss_t(0) = 0.0
    1572        flux_d    = 0.0
    1573        diss_d    = 0.0
     1571!
     1572!--    The lower flux has to be calculated explicetely for the tendency at
     1573!--    the first w-level. For topography wall this is done implicitely by
     1574!--    wall_flags_0.
     1575       k         = nzb + 1
     1576       w_comp    = w(k,j,i) + w(k-1,j,i)
     1577       flux_t(0) = w_comp       * ( w(k,j,i) + w(k-1,j,i) ) * adv_mom_1
     1578       diss_t(0) = -ABS(w_comp) * ( w(k,j,i) - w(k-1,j,i) ) * adv_mom_1
     1579       flux_d    = flux_t(0)
     1580       diss_d    = diss_t(0)
    15741581!
    15751582!--    Now compute the fluxes and tendency terms for the horizontal
     
    16751682                             ( w(k_ppp,j,i) - w(k_mm,j,i) )                  &
    16761683                                         )
     1684
    16771685!
    16781686!--       Calculate the divergence of the velocity field. A respective
     
    18171825       IMPLICIT NONE
    18181826
    1819        INTEGER ::  i, j, k, tn = 0, k_ppp, k_pp, k_mm
     1827       INTEGER ::  i, j, k, k_mm, k_pp, k_ppp, tn = 0
    18201828       REAL, DIMENSION(:,:,:), POINTER ::  sk
    1821        REAL    :: flux_d, diss_d, u_comp, v_comp, div
    1822        REAL, DIMENSION(nzb:nzt)  ::  flux_r, diss_r, flux_n, diss_n, flux_t,  &
    1823                                      diss_t
    1824        REAL, DIMENSION(nzb+1:nzt) :: swap_flux_y_local, swap_diss_y_local
    1825        REAL, DIMENSION(nzb+1:nzt,nys:nyn) :: swap_flux_x_local,               &
    1826                                              swap_diss_x_local
    1827        CHARACTER (LEN = *), INTENT(IN) :: sk_char
     1829       REAL    :: diss_d, div, flux_d, u_comp, v_comp
     1830       REAL, DIMENSION(nzb:nzt)   :: diss_n, diss_r, diss_t, flux_n, flux_r,  &
     1831                                     flux_t
     1832       REAL, DIMENSION(nzb+1:nzt) :: swap_diss_y_local, swap_flux_y_local
     1833       REAL, DIMENSION(nzb+1:nzt,nys:nyn) :: swap_diss_x_local,               &
     1834                                             swap_flux_x_local
     1835       CHARACTER (LEN = *), INTENT(IN)    :: sk_char
    18281836
    18291837!
     
    21972205       IMPLICIT NONE
    21982206
    2199        INTEGER ::  i, j, k, tn = 0,  k_ppp, k_pp, k_mm
    2200        REAL    ::  gu, gv, flux_d, diss_d, v_comp, w_comp, div
    2201        REAL, DIMENSION(nzb+1:nzt) :: swap_flux_y_local_u, swap_diss_y_local_u
    2202        REAL, DIMENSION(nzb+1:nzt,nys:nyn) :: swap_flux_x_local_u,           &
    2203                                              swap_diss_x_local_u
    2204        REAL, DIMENSION(nzb:nzt) :: flux_t, diss_t, flux_r, diss_r, flux_n,  &
    2205                                      diss_n, u_comp
     2207       INTEGER ::  i, j, k, k_mm, k_pp, k_ppp, tn = 0
     2208       REAL    ::  diss_d, div, flux_d, gu, gv, v_comp, w_comp
     2209       REAL, DIMENSION(nzb+1:nzt) :: swap_diss_y_local_u, swap_flux_y_local_u
     2210       REAL, DIMENSION(nzb+1:nzt,nys:nyn) :: swap_diss_x_local_u,              &
     2211                                             swap_flux_x_local_u
     2212       REAL, DIMENSION(nzb:nzt) :: diss_n, diss_r, diss_t, flux_n, flux_r,     &
     2213                                   flux_t, u_comp
    22062214 
    22072215       gu = 2.0 * u_gtrans
     
    25912599
    25922600
    2593        INTEGER ::  i, j, k, tn = 0, k_ppp, k_pp, k_mm
    2594        REAL    ::  gu, gv, flux_d, diss_d, u_comp, w_comp, div
    2595        REAL, DIMENSION(nzb+1:nzt) :: swap_flux_y_local_v, swap_diss_y_local_v
    2596        REAL, DIMENSION(nzb+1:nzt,nys:nyn) :: swap_flux_x_local_v,             &
    2597                                              swap_diss_x_local_v
    2598        REAL, DIMENSION(nzb:nzt) :: flux_t, diss_t, flux_n, diss_n, flux_r,    &
    2599                                    diss_r, v_comp
     2601       INTEGER ::  i, j, k, k_mm, k_pp, k_ppp, tn = 0
     2602       REAL    ::  diss_d, div, flux_d, gu, gv, u_comp, w_comp
     2603       REAL, DIMENSION(nzb+1:nzt) :: swap_diss_y_local_v, swap_flux_y_local_v
     2604       REAL, DIMENSION(nzb+1:nzt,nys:nyn) :: swap_diss_x_local_v,             &
     2605                                             swap_flux_x_local_v
     2606       REAL, DIMENSION(nzb:nzt) :: diss_n, diss_r, diss_t, flux_n, flux_r,    &
     2607                                   flux_t, v_comp
    26002608
    26012609       gu = 2.0 * u_gtrans
     
    29933001       IMPLICIT NONE
    29943002
    2995        INTEGER ::  i, j, k, tn = 0, k_ppp, k_pp, k_mm
    2996        REAL    ::  gu, gv, flux_d, diss_d, u_comp, v_comp, w_comp, div
    2997        REAL, DIMENSION(nzb:nzt)    ::  flux_t, diss_t
    2998        REAL, DIMENSION(nzb+1:nzt)  ::  flux_r, diss_r, flux_n, diss_n,        &
    2999                                        swap_flux_y_local_w,                   &
    3000                                        swap_diss_y_local_w
    3001        REAL, DIMENSION(nzb+1:nzt,nys:nyn) :: swap_flux_x_local_w,             &
    3002                                              swap_diss_x_local_w
     3003       INTEGER ::  i, j, k, k_mm, k_pp, k_ppp, tn = 0
     3004       REAL    ::  diss_d, div, flux_d, gu, gv, u_comp, v_comp, w_comp
     3005       REAL, DIMENSION(nzb:nzt)    ::  diss_t, flux_t
     3006       REAL, DIMENSION(nzb+1:nzt)  ::  diss_n, diss_r, flux_n, flux_r,        &
     3007                                       swap_diss_y_local_w,                   &
     3008                                       swap_flux_y_local_w
     3009       REAL, DIMENSION(nzb+1:nzt,nys:nyn) :: swap_diss_x_local_w,             &
     3010                                             swap_flux_x_local_w
    30033011 
    30043012       gu = 2.0 * u_gtrans
     
    31133121          DO  j = nys, nyn
    31143122
    3115              flux_t(0) = 0.0
    3116              diss_t(0) = 0.0
    3117              flux_d    = 0.0
    3118              diss_d    = 0.0
     3123!
     3124!--          The lower flux has to be calculated explicetely for the tendency
     3125!--          at the first w-level. For topography wall this is done implicitely
     3126!--          by wall_flags_0.
     3127             k         = nzb + 1
     3128             w_comp    = w(k,j,i) + w(k-1,j,i)
     3129             flux_t(0) = w_comp       * ( w(k,j,i) + w(k-1,j,i) ) * adv_mom_1
     3130             diss_t(0) = -ABS(w_comp) * ( w(k,j,i) - w(k-1,j,i) ) * adv_mom_1
     3131             flux_d    = flux_t(0)
     3132             diss_d    = diss_t(0)
    31193133
    31203134             DO  k = nzb+1, nzb_max
Note: See TracChangeset for help on using the changeset viewer.