Changeset 996 for palm/trunk/SOURCE


Ignore:
Timestamp:
Sep 7, 2012 10:41:47 AM (12 years ago)
Author:
raasch
Message:

parameter use_prior_plot1d_parameters removed; little reformatting

Location:
palm/trunk/SOURCE
Files:
9 edited

Legend:

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

    r979 r996  
    44! Current revisions:
    55! -----------------
    6 !
     6! little reformatting
    77
    88! Former revisions:
     
    259259
    260260!
    261 !--       Calculate the phase speeds for u,v, and w, first local and then
    262 !--       average parallel along the outflow boundary. 
    263           DO k = nzb+1, nzt+1
    264              DO i = nxl, nxr
     261!--       Calculate the phase speeds for u, v, and w, first local and then
     262!--       average along the outflow boundary.
     263          DO  k = nzb+1, nzt+1
     264             DO  i = nxl, nxr
    265265
    266266                denom = u_m_s(k,0,i) - u_m_s(k,1,i)
    267267
    268268                IF ( denom /= 0.0 )  THEN
    269                    c_u(k,i) = -c_max * ( u(k,0,i) - u_m_s(k,0,i) )             &
    270                               / ( denom * tsc(2) )
     269                   c_u(k,i) = -c_max * ( u(k,0,i) - u_m_s(k,0,i) ) / ( denom * tsc(2) )
    271270                   IF ( c_u(k,i) < 0.0 )  THEN
    272271                      c_u(k,i) = 0.0
     
    281280
    282281                IF ( denom /= 0.0 )  THEN
    283                    c_v(k,i) = -c_max * ( v(k,1,i) - v_m_s(k,1,i) )             &
    284                               / ( denom * tsc(2) )
     282                   c_v(k,i) = -c_max * ( v(k,1,i) - v_m_s(k,1,i) ) / ( denom * tsc(2) )
    285283                   IF ( c_v(k,i) < 0.0 )  THEN
    286284                      c_v(k,i) = 0.0
     
    295293
    296294                IF ( denom /= 0.0 )  THEN
    297                    c_w(k,i) = -c_max * ( w(k,0,i) - w_m_s(k,0,i) )             &
    298                               / ( denom * tsc(2) )
     295                   c_w(k,i) = -c_max * ( w(k,0,i) - w_m_s(k,0,i) ) / ( denom * tsc(2) )
    299296                   IF ( c_w(k,i) < 0.0 )  THEN
    300297                      c_w(k,i) = 0.0
     
    343340!
    344341!--       Calculate the new velocities
    345           DO k = nzb+1, nzt+1
    346              DO i = nxlg, nxrg
     342          DO  k = nzb+1, nzt+1
     343             DO  i = nxlg, nxrg
    347344                u_p(k,-1,i) = u(k,-1,i) - dt_3d * tsc(2) * c_u_m(k) *          &
    348345                                       ( u(k,-1,i) - u(k,0,i) ) * ddy
     
    401398
    402399!
    403 !--       Calculate the phase speeds for u,v, and w, first local and then
    404 !--       average parallel along the outflow boundary. 
    405           DO k = nzb+1, nzt+1
    406              DO i = nxl, nxr
     400!--       Calculate the phase speeds for u, v, and w, first local and then
     401!--       average along the outflow boundary.
     402          DO  k = nzb+1, nzt+1
     403             DO  i = nxl, nxr
    407404
    408405                denom = u_m_n(k,ny,i) - u_m_n(k,ny-1,i)
    409406
    410407                IF ( denom /= 0.0 )  THEN
    411                    c_u(k,i) = -c_max * ( u(k,ny,i) - u_m_n(k,ny,i) )           &
    412                               / ( denom * tsc(2) )
     408                   c_u(k,i) = -c_max * ( u(k,ny,i) - u_m_n(k,ny,i) ) / ( denom * tsc(2) )
    413409                   IF ( c_u(k,i) < 0.0 )  THEN
    414410                      c_u(k,i) = 0.0
     
    423419
    424420                IF ( denom /= 0.0 )  THEN
    425                    c_v(k,i) = -c_max * ( v(k,ny,i) - v_m_n(k,ny,i) )           &
    426                               / ( denom * tsc(2) )
     421                   c_v(k,i) = -c_max * ( v(k,ny,i) - v_m_n(k,ny,i) ) / ( denom * tsc(2) )
    427422                   IF ( c_v(k,i) < 0.0 )  THEN
    428423                      c_v(k,i) = 0.0
     
    437432
    438433                IF ( denom /= 0.0 )  THEN
    439                    c_w(k,i) = -c_max * ( w(k,ny,i) - w_m_n(k,ny,i) )           &
    440                               / ( denom * tsc(2) )
     434                   c_w(k,i) = -c_max * ( w(k,ny,i) - w_m_n(k,ny,i) ) / ( denom * tsc(2) )
    441435                   IF ( c_w(k,i) < 0.0 )  THEN
    442436                      c_w(k,i) = 0.0
     
    485479!
    486480!--       Calculate the new velocities
    487           DO k = nzb+1, nzt+1
    488              DO i = nxlg, nxrg
     481          DO  k = nzb+1, nzt+1
     482             DO  i = nxlg, nxrg
    489483                u_p(k,ny+1,i) = u(k,ny+1,i) - dt_3d * tsc(2) * c_u_m(k) *      &
    490484                                       ( u(k,ny+1,i) - u(k,ny,i) ) * ddy
     
    543537
    544538!
    545 !--       Calculate the phase speeds for u,v, and w, first local and then
    546 !--       average parallel along the outflow boundary. 
    547           DO k = nzb+1, nzt+1
    548              DO j = nys, nyn
     539!--       Calculate the phase speeds for u, v, and w, first local and then
     540!--       average along the outflow boundary.
     541          DO  k = nzb+1, nzt+1
     542             DO  j = nys, nyn
    549543
    550544                denom = u_m_l(k,j,1) - u_m_l(k,j,2)
    551545
    552546                IF ( denom /= 0.0 )  THEN
    553                    c_u(k,j) = -c_max * ( u(k,j,1) - u_m_l(k,j,1) )             &
    554                               / ( denom * tsc(2) )
     547                   c_u(k,j) = -c_max * ( u(k,j,1) - u_m_l(k,j,1) ) / ( denom * tsc(2) )
    555548                   IF ( c_u(k,j) < 0.0 )  THEN
    556549                      c_u(k,j) = 0.0
     
    565558
    566559                IF ( denom /= 0.0 )  THEN
    567                    c_v(k,j) = -c_max * ( v(k,j,0) - v_m_l(k,j,0) )             &
    568                               / ( denom * tsc(2) )
     560                   c_v(k,j) = -c_max * ( v(k,j,0) - v_m_l(k,j,0) ) / ( denom * tsc(2) )
    569561                   IF ( c_v(k,j) < 0.0 )  THEN
    570562                      c_v(k,j) = 0.0
     
    579571
    580572                IF ( denom /= 0.0 )  THEN
    581                    c_w(k,j) = -c_max * ( w(k,j,0) - w_m_l(k,j,0) )             &
    582                               / ( denom * tsc(2) )
     573                   c_w(k,j) = -c_max * ( w(k,j,0) - w_m_l(k,j,0) ) / ( denom * tsc(2) )
    583574                   IF ( c_w(k,j) < 0.0 )  THEN
    584575                      c_w(k,j) = 0.0
     
    627618!
    628619!--       Calculate the new velocities
    629           DO k = nzb+1, nzt+1
    630              DO i = nxlg, nxrg
     620          DO  k = nzb+1, nzt+1
     621             DO  i = nxlg, nxrg
    631622                u_p(k,j,0) = u(k,j,0) - dt_3d * tsc(2) * c_u_m(k) *            &
    632623                                       ( u(k,j,0) - u(k,j,1) ) * ddx
     
    685676
    686677!
    687 !--       Calculate the phase speeds for u,v, and w, first local and then
    688 !--       average parallel along the outflow boundary. 
    689           DO k = nzb+1, nzt+1
    690              DO j = nys, nyn
     678!--       Calculate the phase speeds for u, v, and w, first local and then
     679!--       average along the outflow boundary.
     680          DO  k = nzb+1, nzt+1
     681             DO  j = nys, nyn
    691682
    692683                denom = u_m_r(k,j,nx) - u_m_r(k,j,nx-1)
    693684
    694685                IF ( denom /= 0.0 )  THEN
    695                    c_u(k,j) = -c_max * ( u(k,j,nx) - u_m_r(k,j,nx) )           &
    696                               / ( denom * tsc(2) )
     686                   c_u(k,j) = -c_max * ( u(k,j,nx) - u_m_r(k,j,nx) ) / ( denom * tsc(2) )
    697687                   IF ( c_u(k,j) < 0.0 )  THEN
    698688                      c_u(k,j) = 0.0
     
    707697
    708698                IF ( denom /= 0.0 )  THEN
    709                    c_v(k,j) = -c_max * ( v(k,j,nx) - v_m_r(k,j,nx) )           &
    710                               / ( denom * tsc(2) )
     699                   c_v(k,j) = -c_max * ( v(k,j,nx) - v_m_r(k,j,nx) ) / ( denom * tsc(2) )
    711700                   IF ( c_v(k,j) < 0.0 )  THEN
    712701                      c_v(k,j) = 0.0
     
    721710
    722711                IF ( denom /= 0.0 )  THEN
    723                    c_w(k,j) = -c_max * ( w(k,j,nx) - w_m_r(k,j,nx) )           &
    724                               / ( denom * tsc(2) )
     712                   c_w(k,j) = -c_max * ( w(k,j,nx) - w_m_r(k,j,nx) ) / ( denom * tsc(2) )
    725713                   IF ( c_w(k,j) < 0.0 )  THEN
    726714                      c_w(k,j) = 0.0
     
    769757!
    770758!--       Calculate the new velocities
    771           DO k = nzb+1, nzt+1
    772              DO i = nxlg, nxrg
     759          DO  k = nzb+1, nzt+1
     760             DO  i = nxlg, nxrg
    773761                u_p(k,j,nx+1) = u(k,j,nx+1) - dt_3d * tsc(2) * c_u_m(k) *      &
    774762                                       ( u(k,j,nx+1) - u(k,j,nx) ) * ddx
  • palm/trunk/SOURCE/check_parameters.f90

    r979 r996  
    44! Current revisions:
    55! -----------------
    6 !
     6! little reformatting
    77!
    88! Former revisions:
     
    27122712             IF ( TRIM( var ) == 'u*'     )  unit = 'm/s'
    27132713             IF ( TRIM( var ) == 'z0*'    )  unit = 'm'
    2714              IF ( TRIM( var ) == 'z0h*'    )  unit = 'm'
     2714             IF ( TRIM( var ) == 'z0h*'   )  unit = 'm'
    27152715
    27162716
     
    29802980!-- potential temperature, check the width of the damping layer
    29812981    IF ( bc_lr /= 'cyclic' ) THEN
    2982        IF ( pt_damping_width < 0.0 .OR. pt_damping_width > REAL( nx * dx ) )  THEN
     2982       IF ( pt_damping_width < 0.0  .OR. pt_damping_width > REAL( nx * dx ) )  THEN
    29832983          message_string = 'pt_damping_width out of range'
    29842984          CALL message( 'check_parameters', 'PA0124', 1, 2, 0, 6, 0 )
     
    29872987
    29882988    IF ( bc_ns /= 'cyclic' )  THEN
    2989        IF ( pt_damping_width < 0.0 .OR. pt_damping_width > REAL( ny * dy ) )  THEN
     2989       IF ( pt_damping_width < 0.0  .OR. pt_damping_width > REAL( ny * dy ) )  THEN
    29902990          message_string = 'pt_damping_width out of range'
    29912991          CALL message( 'check_parameters', 'PA0124', 1, 2, 0, 6, 0 )
  • palm/trunk/SOURCE/init_1d_model.f90

    r979 r996  
    44! Current revisions:
    55! -----------------
    6 !
     6! little reformatting
    77!
    88! Former revisions:
     
    151151    usws1d = 0.0; usws1d_m = 0.0
    152152    vsws1d = 0.0; vsws1d_m = 0.0
    153     z01d = roughness_length
     153    z01d  = roughness_length
    154154    z0h1d = z0h_factor * z01d
    155155    IF ( humidity .OR. passive_scalar )  qs1d = 0.0
     
    448448!
    449449!--                Stable stratification
    450                    ts1d = kappa * ( pt_init(nzb+1) - pt_init(nzb) ) /      &
     450                   ts1d = kappa * ( pt_init(nzb+1) - pt_init(nzb) ) /       &
    451451                          ( LOG( zu(nzb+1) / z0h1d ) + 5.0 * rif1d(nzb+1) * &
    452452                                          ( zu(nzb+1) - z0h1d ) / zu(nzb+1) &
     
    462462!--                occur in the argument of the logarithm.
    463463                   IF ( a == 0.0  .OR.  b == 0.0 )  THEN
    464                       ts1d = kappa * ( pt_init(nzb+1) - pt_init(nzb) ) /      &
     464                      ts1d = kappa * ( pt_init(nzb+1) - pt_init(nzb) ) /       &
    465465                             ( LOG( zu(nzb+1) / z0h1d ) + 5.0 * rif1d(nzb+1) * &
    466466                                             ( zu(nzb+1) - z0h1d ) / zu(nzb+1) &
     
    581581!
    582582!--                Stable stratification
    583                    qs1d = kappa * ( q_init(nzb+1) - q_init(nzb) ) /        &
     583                   qs1d = kappa * ( q_init(nzb+1) - q_init(nzb) ) /         &
    584584                          ( LOG( zu(nzb+1) / z0h1d ) + 5.0 * rif1d(nzb+1) * &
    585585                                          ( zu(nzb+1) - z0h1d ) / zu(nzb+1) &
     
    595595!--                occur in the argument of the logarithm.
    596596                   IF ( a == 1.0  .OR.  b == 1.0 )  THEN
    597                       qs1d = kappa * ( q_init(nzb+1) - q_init(nzb) ) /        &
     597                      qs1d = kappa * ( q_init(nzb+1) - q_init(nzb) ) /         &
    598598                             ( LOG( zu(nzb+1) / z0h1d ) + 5.0 * rif1d(nzb+1) * &
    599599                                             ( zu(nzb+1) - z0h1d ) / zu(nzb+1) &
  • palm/trunk/SOURCE/init_3d_model.f90

    r979 r996  
    77! Current revisions:
    88! ------------------
    9 !
     9! little reformatting
    1010!
    1111! Former revisions:
     
    190190
    191191    INTEGER, DIMENSION(:,:), ALLOCATABLE ::  ngp_2dh_outer_l,  &
    192          ngp_2dh_s_inner_l
     192                                             ngp_2dh_s_inner_l
    193193
    194194    REAL ::  a, b
     
    213213              ngp_2dh_s_inner(nzb:nzt+1,0:statistic_regions),               &
    214214              ngp_2dh_s_inner_l(nzb:nzt+1,0:statistic_regions),             &
    215               rmask(nysg:nyng,nxlg:nxrg,0:statistic_regions),           &
     215              rmask(nysg:nyng,nxlg:nxrg,0:statistic_regions),               &
    216216              sums(nzb:nzt+1,pr_palm+max_pr_user),                          &
    217217              sums_l(nzb:nzt+1,pr_palm+max_pr_user,0:threads_per_task-1),   &
     
    225225              ts(nysg:nyng,nxlg:nxrg), tswst_1(nysg:nyng,nxlg:nxrg),  &
    226226              us(nysg:nyng,nxlg:nxrg), usws_1(nysg:nyng,nxlg:nxrg),   &
    227               uswst_1(nysg:nyng,nxlg:nxrg),                               &
    228               vsws_1(nysg:nyng,nxlg:nxrg),                                &
    229               vswst_1(nysg:nyng,nxlg:nxrg), z0(nysg:nyng,nxlg:nxrg),      &
     227              uswst_1(nysg:nyng,nxlg:nxrg),                           &
     228              vsws_1(nysg:nyng,nxlg:nxrg),                            &
     229              vswst_1(nysg:nyng,nxlg:nxrg), z0(nysg:nyng,nxlg:nxrg),  &
    230230              z0h(nysg:nyng,nxlg:nxrg) )
    231231
     
    242242    ENDIF
    243243
    244     ALLOCATE( d(nzb+1:nzta,nys:nyna,nxl:nxra),         &
     244    ALLOCATE( d(nzb+1:nzta,nys:nyna,nxl:nxra),     &
    245245              e_1(nzb:nzt+1,nysg:nyng,nxlg:nxrg),  &
    246246              e_2(nzb:nzt+1,nysg:nyng,nxlg:nxrg),  &
     
    256256              u_2(nzb:nzt+1,nysg:nyng,nxlg:nxrg),  &
    257257              u_3(nzb:nzt+1,nysg:nyng,nxlg:nxrg),  &
    258               v_1(nzb:nzt+1,nysg:nyng,nxlg:nxrg), &
    259               v_2(nzb:nzt+1,nysg:nyng,nxlg:nxrg), &
    260               v_3(nzb:nzt+1,nysg:nyng,nxlg:nxrg), &
     258              v_1(nzb:nzt+1,nysg:nyng,nxlg:nxrg),  &
     259              v_2(nzb:nzt+1,nysg:nyng,nxlg:nxrg),  &
     260              v_3(nzb:nzt+1,nysg:nyng,nxlg:nxrg),  &
    261261              w_1(nzb:nzt+1,nysg:nyng,nxlg:nxrg),  &
    262262              w_2(nzb:nzt+1,nysg:nyng,nxlg:nxrg),  &
     
    374374
    375375       IF ( passive_scalar ) THEN
    376           ALLOCATE ( sls(nzb:nzt+1,nysg:nyng,nxlg:nxrg),   &
     376          ALLOCATE ( sls(nzb:nzt+1,nysg:nyng,nxlg:nxrg),  &
    377377                     sec(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    378378       ENDIF
    379379
    380380       IF ( cthf /= 0.0 ) THEN
    381           ALLOCATE ( lai(nzb:nzt+1,nysg:nyng,nxlg:nxrg),   &
     381          ALLOCATE ( lai(nzb:nzt+1,nysg:nyng,nxlg:nxrg),  &
    382382                     canopy_heat_flux(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    383383       ENDIF
     
    423423                 c_w(nzb:nzt+1,nxlg:nxrg) )
    424424    ENDIF
    425     IF ( outflow_l  .OR.  outflow_r .OR. outflow_s .OR. outflow_n )  THEN
     425    IF ( outflow_l  .OR.  outflow_r  .OR.  outflow_s  .OR. outflow_n )  THEN
    426426       ALLOCATE( c_u_m_l(nzb:nzt+1), c_v_m_l(nzb:nzt+1), c_w_m_l(nzb:nzt+1) )                   
    427427       ALLOCATE( c_u_m(nzb:nzt+1), c_v_m(nzb:nzt+1), c_w_m(nzb:nzt+1) )
     
    601601!--          ---------  advection scheme: keep u and v zero one layer below
    602602!--                     the topography.
    603 !
    604 !--           Following was removed, because mirror boundary condition are
    605 !--           replaced by dirichlet boundary conditions
    606 !
    607 !             IF ( ibc_uv_b == 0 )  THEN
    608 !!
    609 !!--             Satisfying the Dirichlet condition with an extra layer below
    610 !!--             the surface where the u and v component change their sign.
    611 !                DO  i = nxl-1, nxr+1
    612 !                   DO  j = nys-1, nyn+1
    613 !                      IF ( nzb_u_inner(j,i) == 0 ) u(0,j,i) = -u(1,j,i)
    614 !                      IF ( nzb_v_inner(j,i) == 0 ) v(0,j,i) = -v(1,j,i)
    615 !                   ENDDO
    616 !                ENDDO
    617 !
    618 !             ELSE
    619603             IF ( ibc_uv_b == 1 )  THEN
    620604!
     
    10731057             IF ( zu(k) <= inflow_damping_height )  THEN
    10741058                inflow_damping_factor(k) = 1.0
    1075              ELSEIF ( zu(k) <= inflow_damping_height +  &
    1076                                inflow_damping_width )  THEN
    1077                 inflow_damping_factor(k) = 1.0 -                            &
    1078                                         ( zu(k) - inflow_damping_height ) / &
    1079                                         inflow_damping_width
     1059             ELSEIF ( zu(k) <= ( inflow_damping_height + inflow_damping_width ) )  THEN
     1060                inflow_damping_factor(k) = 1.0 -                               &
     1061                                           ( zu(k) - inflow_damping_height ) / &
     1062                                           inflow_damping_width
    10801063             ELSE
    10811064                inflow_damping_factor(k) = 0.0
     
    10971080             DO  i = nxlg, nxrg
    10981081                DO  j = nysg, nyng
    1099                    u  (nzb:nzb_u_inner(j,i),j,i) = 0.0
    1100                    v  (nzb:nzb_v_inner(j,i),j,i) = 0.0
    1101                    w  (nzb:nzb_w_inner(j,i),j,i) = 0.0
    1102                    e  (nzb:nzb_w_inner(j,i),j,i) = 0.0
    1103                    u_m(nzb:nzb_u_inner(j,i),j,i) = 0.0
    1104                    v_m(nzb:nzb_v_inner(j,i),j,i) = 0.0
    1105                    w_m(nzb:nzb_w_inner(j,i),j,i) = 0.0
    1106                    e_m(nzb:nzb_w_inner(j,i),j,i) = 0.0
    1107                    tu_m(nzb:nzb_u_inner(j,i),j,i) = 0.0
    1108                    tv_m(nzb:nzb_v_inner(j,i),j,i) = 0.0
    1109                    tw_m(nzb:nzb_w_inner(j,i),j,i) = 0.0
    1110                    te_m(nzb:nzb_w_inner(j,i),j,i) = 0.0
     1082                   u  (nzb:nzb_u_inner(j,i),j,i)   = 0.0
     1083                   v  (nzb:nzb_v_inner(j,i),j,i)   = 0.0
     1084                   w  (nzb:nzb_w_inner(j,i),j,i)   = 0.0
     1085                   e  (nzb:nzb_w_inner(j,i),j,i)   = 0.0
     1086                   u_m(nzb:nzb_u_inner(j,i),j,i)   = 0.0
     1087                   v_m(nzb:nzb_v_inner(j,i),j,i)   = 0.0
     1088                   w_m(nzb:nzb_w_inner(j,i),j,i)   = 0.0
     1089                   e_m(nzb:nzb_w_inner(j,i),j,i)   = 0.0
     1090                   tu_m(nzb:nzb_u_inner(j,i),j,i)  = 0.0
     1091                   tv_m(nzb:nzb_v_inner(j,i),j,i)  = 0.0
     1092                   tw_m(nzb:nzb_w_inner(j,i),j,i)  = 0.0
     1093                   te_m(nzb:nzb_w_inner(j,i),j,i)  = 0.0
    11111094                   tpt_m(nzb:nzb_w_inner(j,i),j,i) = 0.0
    11121095                ENDDO
     
    15541537    ptdf_x = 0.0
    15551538    ptdf_y = 0.0
    1556     IF ( bc_lr_dirrad .OR. bc_lr_dirneu )  THEN
    1557        DO i = nxl, nxr
     1539    IF ( bc_lr_dirrad  .OR. bc_lr_dirneu )  THEN
     1540       DO  i = nxl, nxr
    15581541          IF ( ( i * dx ) < pt_damping_width )  THEN
    15591542             ptdf_x(i) = pt_damping_factor * ( SIN( pi * 0.5 *        &
     
    15621545          ENDIF
    15631546       ENDDO
    1564     ELSEIF ( bc_lr_raddir .OR. bc_lr_neudir )  THEN
    1565        DO i = nxl, nxr
     1547    ELSEIF ( bc_lr_raddir  .OR. bc_lr_neudir )  THEN
     1548       DO  i = nxl, nxr
    15661549          IF ( ( i * dx ) > ( nx * dx - pt_damping_width ) )  THEN
    1567              ptdf_x(i) = pt_damping_factor * ( SIN( pi * 0.5 *                 &
    1568                          REAL( ( i - nx ) * dx + pt_damping_width ) / (        &
    1569                          REAL( pt_damping_width )                     ) ) )**2       
     1550             ptdf_x(i) = pt_damping_factor *                                      &
     1551                         SIN( pi * 0.5 * ( ( i - nx ) * dx + pt_damping_width ) / &
     1552                                         REAL( pt_damping_width ) )**2
    15701553          ENDIF
    15711554       ENDDO
    1572     ELSEIF ( bc_ns_dirrad .OR. bc_ns_dirneu )  THEN
    1573        DO j = nys, nyn
     1555    ELSEIF ( bc_ns_dirrad  .OR. bc_ns_dirneu )  THEN
     1556       DO  j = nys, nyn
    15741557          IF ( ( j * dy ) > ( ny * dy - pt_damping_width ) )  THEN
    1575              ptdf_y(j) = pt_damping_factor * ( SIN( pi * 0.5 *                 &
    1576                          REAL( ( j - ny ) * dy + pt_damping_width ) / (        &
    1577                          REAL( pt_damping_width )                     ) ) )**2       
     1558             ptdf_y(j) = pt_damping_factor *                                      &
     1559                         SIN( pi * 0.5 * ( ( j - ny ) * dy + pt_damping_width ) / &
     1560                                         REAL( pt_damping_width ) )**2
    15781561          ENDIF
    15791562       ENDDO
    1580     ELSEIF ( bc_ns_raddir .OR. bc_ns_neudir )  THEN
    1581        DO j = nys, nyn
     1563    ELSEIF ( bc_ns_raddir  .OR. bc_ns_neudir )  THEN
     1564       DO  j = nys, nyn
    15821565          IF ( ( j * dy ) < pt_damping_width )  THEN
    1583              ptdf_y(j) = pt_damping_factor * ( SIN( pi * 0.5 *        &
    1584                          REAL( pt_damping_width - j * dy ) / (        &
    1585                          REAL( pt_damping_width )            ) ) )**2       
     1566             ptdf_y(j) = pt_damping_factor *                             &
     1567                         SIN( pi * 0.5 * ( pt_damping_width - j * dy ) / &
     1568                                         REAL( pt_damping_width ) )**2
    15861569          ENDIF
    15871570       ENDDO
  • palm/trunk/SOURCE/init_grid.f90

    r979 r996  
    44! Current revisions:
    55! -----------------
    6 !
     6! little reformatting
    77!
    88! Former revisions:
     
    613613!-- steering the degradation of order of the applied advection scheme.
    614614!-- In case of non-cyclic lateral boundaries, the order of the advection
    615 !-- scheme is reduced at the lateral boundaries up to nzt.
     615!-- scheme have to be reduced up to nzt (required at the lateral boundaries).
    616616    nzb_max = MAXVAL( nzb_local )
    617617    IF ( inflow_l .OR. outflow_l .OR. inflow_r .OR. outflow_r .OR.    &
  • palm/trunk/SOURCE/init_masks.f90

    r979 r996  
    44! Current revisions:
    55! -----------------
    6 !
     6! little reformatting
    77!
    88! Former revisions:
     
    8080!-- Parallel mask output not yet tested
    8181    IF ( netcdf_data_format > 2 )  THEN
    82        message_string = 'NetCDF file formats '//                         &
    83                         '3 (NetCDF 4) and 4 (NetCDF 4 Classic model)'//  &
     82       message_string = 'netCDF file formats '//                         &
     83                        '3 (netCDF 4) and 4 (netCDF 4 Classic model)'//  &
    8484                        '&are currently not supported (not yet tested).'
    8585       CALL message( 'init_masks', 'PA0328', 1, 2, 0, 6, 0 )
     
    8888!
    8989!-- Store data output parameters for masked data output in few shared arrays
    90     DO mid = 1, masks
     90    DO  mid = 1, masks
    9191   
    9292       do_mask     (mid,:) = data_output_masks(mid,:)
     
    9696       mask      (mid,3,:) = mask_z(mid,:)
    9797       
    98        IF ( mask_x_loop(mid,1) == -1.0 .AND. mask_x_loop(mid,2) == -1.0  &
    99             .AND. mask_x_loop(mid,3) == -1.0 ) THEN
    100           mask_loop(mid,1,1:2) = -1.0
     98       IF ( mask_x_loop(mid,1) == -1.0  .AND. mask_x_loop(mid,2) == -1.0  &
     99            .AND.  mask_x_loop(mid,3) == -1.0 ) THEN
     100          mask_loop(mid,1,1:2) = -1.0
    101101          mask_loop(mid,1,3) = 0.0
    102102       ELSE
    103103          mask_loop(mid,1,:) = mask_x_loop(mid,:)
    104104       ENDIF
    105        IF ( mask_y_loop(mid,1) == -1.0 .AND. mask_y_loop(mid,2) == -1.0  &
    106             .AND. mask_y_loop(mid,3) == -1.0 ) THEN
    107           mask_loop(mid,2,1:2) = -1.0
     105       IF ( mask_y_loop(mid,1) == -1.0  .AND. mask_y_loop(mid,2) == -1.0  &
     106            .AND.  mask_y_loop(mid,3) == -1.0 ) THEN
     107          mask_loop(mid,2,1:2) = -1.0
    108108          mask_loop(mid,2,3) = 0.0
    109109       ELSE
    110110          mask_loop(mid,2,:) = mask_y_loop(mid,:)
    111111       ENDIF
    112        IF ( mask_z_loop(mid,1) == -1.0 .AND. mask_z_loop(mid,2) == -1.0  &
    113             .AND. mask_z_loop(mid,3) == -1.0 ) THEN
    114           mask_loop(mid,3,1:2) = -1.0
     112       IF ( mask_z_loop(mid,1) == -1.0  .AND. mask_z_loop(mid,2) == -1.0  &
     113            .AND.  mask_z_loop(mid,3) == -1.0 ) THEN
     114          mask_loop(mid,3,1:2) = -1.0
    115115          mask_loop(mid,3,3) = 0.0
    116116       ELSE
     
    158158!--    Check and set steering parameters for mask data output and averaging
    159159       i   = 1
    160        DO  WHILE ( do_mask(mid,i) /= ' '  .AND.  i <= 100 )
     160       DO WHILE ( do_mask(mid,i) /= ' '  .AND.  i <= 100 )
    161161!
    162162!--       Check for data averaging
     
    265265                     'for horizontal cross section'
    266266                CALL message( 'init_masks', 'PA0111', 1, 2, 0, 6, 0 )
    267 !                IF ( TRIM( var ) == 'lwp*'  .AND.  .NOT. cloud_physics )  THEN
    268 !                   WRITE ( message_string, * ) 'output of "', TRIM( var ), &
    269 !                         '" requires cloud_physics = .TRUE.'
    270 !                   CALL message( 'init_masks', 'PA0108', 1, 2, 0, 6, 0 )
    271 !                ENDIF
    272 !                IF ( TRIM( var ) == 'pra*'  .AND.  .NOT. precipitation )  THEN
    273 !                   WRITE ( message_string, * ) 'output of "', TRIM( var ), &
    274 !                         '" requires precipitation = .TRUE.'
    275 !                   CALL message( 'init_masks', 'PA0112', 1, 2, 0, 6, 0 )
    276 !                ENDIF
    277 !                IF ( TRIM( var ) == 'pra*'  .AND.  j == 1 )  THEN
    278 !                   WRITE ( message_string, * ) 'temporal averaging of ', &
    279 !                         ' precipitation amount "', TRIM( var ),         &
    280 !                         '" not possible'
    281 !                   CALL message( 'init_masks', 'PA0113', 1, 2, 0, 6, 0 )
    282 !                ENDIF
    283 !                IF ( TRIM( var ) == 'prr*'  .AND.  .NOT. precipitation )  THEN
    284 !                   WRITE ( message_string, * ) 'output of "', TRIM( var ), &
    285 !                         '" requires precipitation = .TRUE.'
    286 !                   CALL message( 'init_masks', 'PA0112', 1, 2, 0, 6, 0 )
    287 !                ENDIF
    288 !
    289 !                IF ( TRIM( var ) == 'u*'   )  unit = 'm/s'
    290 !                IF ( TRIM( var ) == 't*'   )  unit = 'K'
    291 !                IF ( TRIM( var ) == 'lwp*' )  unit = 'kg/kg*m'
    292 !                IF ( TRIM( var ) == 'pra*' )  unit = 'mm'
    293 !                IF ( TRIM( var ) == 'prr*' )  unit = 'mm/s'
    294 !                IF ( TRIM( var ) == 'z0*'  )  unit = 'm'
    295267
    296268             CASE ( 'p', 'pt', 'u', 'v', 'w' )
     
    404376!--       If at least part of the mask resides on the PE, send the index limits
    405377!--       for the target array, otherwise send -9999 to PE0.
    406           IF ( mask_size_l(mid,1) > 0 .AND.  mask_size_l(mid,2) > 0 .AND. &
    407                mask_size_l(mid,3) > 0  ) &
    408                THEN
     378          IF ( mask_size_l(mid,1) > 0  .AND.  mask_size_l(mid,2) > 0  .AND.  &
     379               mask_size_l(mid,3) > 0  )  THEN
    409380             ind(1) = mask_start_l(mid,1)
    410381             ind(2) = mask_start_l(mid,1) + mask_size_l(mid,1) - 1
     
    599570          DO  m = loop_begin, loop_end, loop_stride
    600571             count = count + 1
    601              IF ( m >= lb .AND. m <= ub )  THEN
     572             IF ( m >= lb  .AND. m <= ub )  THEN
    602573                IF ( count_l == 0 )  mask_start_l(mid,dim) = count
    603574                count_l = count_l + 1
  • palm/trunk/SOURCE/modules.f90

    r979 r996  
    44! Current revisions:
    55! -----------------
    6 !
     6! -use_prior_plot1d_parameters
    77!
    88! Former revisions:
     
    627627                stop_dt = .FALSE., synchronous_exchange = .FALSE., &
    628628                terminate_run = .FALSE., turbulent_inflow = .FALSE., &
    629                 use_prescribed_profile_data = .FALSE., &
    630                 use_prior_plot1d_parameters = .FALSE., use_reference = .FALSE.,&
     629                use_prescribed_profile_data = .FALSE., use_reference = .FALSE.,&
    631630                use_surface_fluxes = .FALSE., use_top_fluxes = .FALSE., &
    632631                use_ug_for_galilei_tr = .TRUE., use_upstream_for_tke = .FALSE.,&
  • palm/trunk/SOURCE/parin.f90

    r979 r996  
    44! Current revisions:
    55! -----------------
    6 !
     6! -use_prior_plot1d_parameters
    77!
    88! Former revisions:
     
    229229             skip_time_do2d_xz, skip_time_do2d_yz, skip_time_do3d, &
    230230             skip_time_domask, synchronous_exchange, termination_time_needed, &
    231              use_prior_plot1d_parameters, z_max_do2d
     231             z_max_do2d
    232232
    233233
     
    276276!
    277277!--       Read the control parameters for initialization.
    278 !--       The namelist "inipar" must be provided in the NAMELIST-file. If this
    279 !--       is not the case and the file contains - instead of "inipar" - any
    280 !--       other namelist, a read error is created on t3e and control is
    281 !--       transferred to the statement with label 10. Therefore, on t3e
    282 !--       machines one can not distinguish between errors produced by a wrong
    283 !--       "inipar" namelist or because this namelist is totally missing.
    284            READ ( 11, inipar, ERR=10, END=11 )
     278!--       The namelist "inipar" must be provided in the NAMELIST-file.
     279          READ ( 11, inipar, ERR=10, END=11 )
    285280
    286281#if defined ( __check )
     
    289284!--       used. The p3d file here must be closed and the p3df file for reading
    290285!--       3dpar is opened.
    291           IF ( check_restart == 1 ) THEN
     286          IF ( check_restart == 1 )  THEN
    292287             CALL close_file( 11 )
    293288             check_restart = 2
    294289             CALL check_open( 11 )             
    295290             initializing_actions = 'read_restart_data'
    296           END IF
     291          ENDIF
    297292#endif
    298            GOTO 12
     293          GOTO 12
    299294
    300295 10       message_string = 'errors in \$inipar &or no \$inipar-namelist ' // &
     
    309304!--       a prior run). All PEs are reading from file created by PE0 (see
    310305!--       check_open)
    311 
    312 
    313306 12       IF ( TRIM( initializing_actions ) == 'read_restart_data' )  THEN
    314307#if ! defined ( __check )
  • palm/trunk/SOURCE/poismg.f90

    r979 r996  
    77! Current revisions:
    88! -----------------
    9 !
     9! little reformatting
    1010!
    1111! Former revisions:
     
    12321232!--       outflow conditions have to be used on all PEs after the switch,
    12331233!--       because then they have the total domain.
    1234           IF ( bc_lr_dirrad .OR. bc_lr_dirneu )  THEN
     1234          IF ( bc_lr_dirrad  .OR. bc_lr_dirneu )  THEN
    12351235             inflow_l  = .TRUE.
    12361236             inflow_r  = .FALSE.
    12371237             outflow_l = .FALSE.
    12381238             outflow_r = .TRUE.
    1239           ELSEIF ( bc_lr_raddir .OR. bc_lr_neudir )  THEN
     1239          ELSEIF ( bc_lr_raddir  .OR. bc_lr_neudir )  THEN
    12401240             inflow_l  = .FALSE.
    12411241             inflow_r  = .TRUE.
     
    12441244          ENDIF
    12451245
    1246           IF ( bc_ns_dirrad .OR. bc_ns_dirneu )  THEN
     1246          IF ( bc_ns_dirrad  .OR. bc_ns_dirneu )  THEN
    12471247             inflow_n  = .TRUE.
    12481248             inflow_s  = .FALSE.
    12491249             outflow_n = .FALSE.
    12501250             outflow_s = .TRUE.
    1251           ELSEIF ( bc_ns_raddir .OR. bc_ns_neudir )  THEN
     1251          ELSEIF ( bc_ns_raddir  .OR. bc_ns_neudir )  THEN
    12521252             inflow_n  = .FALSE.
    12531253             inflow_s  = .TRUE.
     
    13191319
    13201320          IF ( pleft == MPI_PROC_NULL )  THEN
    1321              IF ( bc_lr_dirrad .OR. bc_lr_dirneu )  THEN
     1321             IF ( bc_lr_dirrad  .OR. bc_lr_dirneu )  THEN
    13221322                inflow_l  = .TRUE.
    1323              ELSEIF ( bc_lr_raddir .OR. bc_lr_neudir )  THEN
     1323             ELSEIF ( bc_lr_raddir  .OR. bc_lr_neudir )  THEN
    13241324                outflow_l = .TRUE.
    13251325             ENDIF
     
    13271327
    13281328          IF ( pright == MPI_PROC_NULL )  THEN
    1329              IF ( bc_lr_dirrad .OR. bc_lr_dirneu )  THEN
     1329             IF ( bc_lr_dirrad  .OR. bc_lr_dirneu )  THEN
    13301330                outflow_r = .TRUE.
    1331              ELSEIF ( bc_lr_raddir .OR. bc_lr_neudir )  THEN
     1331             ELSEIF ( bc_lr_raddir  .OR. bc_lr_neudir )  THEN
    13321332                inflow_r  = .TRUE.
    13331333             ENDIF
     
    13351335
    13361336          IF ( psouth == MPI_PROC_NULL )  THEN
    1337              IF ( bc_ns_dirrad .OR. bc_ns_dirneu )  THEN
     1337             IF ( bc_ns_dirrad  .OR. bc_ns_dirneu )  THEN
    13381338                outflow_s = .TRUE.
    1339              ELSEIF ( bc_ns_raddir .OR. bc_ns_neudir )  THEN
     1339             ELSEIF ( bc_ns_raddir  .OR. bc_ns_neudir )  THEN
    13401340                inflow_s  = .TRUE.
    13411341             ENDIF
     
    13431343
    13441344          IF ( pnorth == MPI_PROC_NULL )  THEN
    1345              IF ( bc_ns_dirrad .OR. bc_ns_dirneu )  THEN
     1345             IF ( bc_ns_dirrad  .OR. bc_ns_dirneu )  THEN
    13461346                inflow_n  = .TRUE.
    1347              ELSEIF ( bc_ns_raddir .OR. bc_ns_neudir )  THEN
     1347             ELSEIF ( bc_ns_raddir  .OR. bc_ns_neudir )  THEN
    13481348                outflow_n = .TRUE.
    13491349             ENDIF
Note: See TracChangeset for help on using the changeset viewer.