Changeset 75 for palm/trunk


Ignore:
Timestamp:
Mar 22, 2007 9:54:05 AM (15 years ago)
Author:
raasch
Message:

preliminary update for changes concerning non-cyclic boundary conditions

Location:
palm/trunk
Files:
46 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/DOC/app/chapter_4.1.html

    r72 r75  
    8888<span style="font-weight: bold;">bc_lr</span> may
    8989also be
    90 assigned the values <span style="font-style: italic;">'dirichlet/neumann'</span>
    91 (inflow from left, outflow to the right) or <span style="font-style: italic;">'neumann/dirichlet'</span>
     90assigned the values <span style="font-style: italic;">'dirichlet/radiation'</span>
     91(inflow from left, outflow to the right) or <span style="font-style: italic;">'radiation/dirichlet'</span>
    9292(inflow from
    9393right, outflow to the left). This requires the multi-grid method to be
    9494used for solving the Poisson equation for perturbation pressure (see <a href="http://www.muk.uni-hannover.de/%7Eraasch/PALM_group/doc/app/chapter_4.2.html#psolver">psolver</a>)
    95 and it also requires cyclic boundary conditions along y (see<br> <a href="#bc_ns">bc_ns</a>).<br> <br>
     95and it also requires cyclic boundary conditions along y (see&nbsp;<a href="#bc_ns">bc_ns</a>).<br> <br>
    9696In case of these non-cyclic lateral boundaries, a Dirichlet condition
    9797is used at the inflow for all quantities (initial vertical profiles -
    9898see <a href="#initializing_actions">initializing_actions</a>
    9999- are fixed during the run) except u, to which a Neumann (zero
    100 gradient) condition is applied. At the outflow, a Neumann (zero
    101 gradient) condition is used for all quantities except v, which is set
    102 to its horizontal average along the outflow (e.g. v(k,:,nx+1) =
    103 average_along_y( v(k,:,nx)), and except w, which is set to zero
    104 (Dirichlet condition). These conditions ensure the velocity field to be
    105 free of divergence at the inflow and at the outflow. For perturbation
     100gradient) condition is applied. At the outflow, a radiation condition is used for all velocity components, while a Neumann (zero
     101gradient) condition is used for the scalars. For perturbation
    106102pressure Neumann (zero gradient) conditions are assumed both at the
    107103inflow and at the outflow.<br> <br>
     
    137133<span style="font-weight: bold;">bc_ns</span> may
    138134also be
    139 assigned the values <span style="font-style: italic;">'dirichlet/neumann'</span>
    140 (inflow from rear ("north"), outflow to the front ("south")) or <span style="font-style: italic;">'neumann/dirichlet'</span>
     135assigned the values <span style="font-style: italic;">'dirichlet/radiation'</span>
     136(inflow from rear ("north"), outflow to the front ("south")) or <span style="font-style: italic;">'radiation/dirichlet'</span>
    141137(inflow from front ("south"), outflow to the rear ("north")). This
    142138requires the multi-grid
     
    146142In case of these non-cyclic lateral boundaries, a Dirichlet condition
    147143is used at the inflow for all quantities (initial vertical profiles -
    148 see <a href="#initializing_actions">initializing_actions</a>
    149 - are fixed during the run) except v, to which a Neumann (zero
    150 gradient) condition is applied. At the outflow, a Neumann (zero
    151 gradient) condition is used for all quantities except u, which is set
    152 to its horizontal average along the outflow (e.g. u(k,ny+1,:) =
    153 average_along_x( u(k,ny,:)), and except w, which is set to zero
    154 (Dirichlet condition). These conditions ensure the velocity field to be
    155 free of divergence at the inflow and at the outflow. For perturbation
     144see <a href="chapter_4.1.html#initializing_actions">initializing_actions</a>
     145- are fixed during the run) except u, to which a Neumann (zero
     146gradient) condition is applied. At the outflow, a radiation condition is used for all velocity components, while a Neumann (zero
     147gradient) condition is used for the scalars. For perturbation
    156148pressure Neumann (zero gradient) conditions are assumed both at the
    157149inflow and at the outflow.<br> <br>
     
    428420Also cloud-top cooling by longwave radiation can be utilized (see <a href="#radiation">radiation</a>)<br> <b><br>
    429421cloud_physics =</b> <span style="font-style: italic;">.TRUE.
    430 </span>requires <a href="#moisture">moisture</a>
     422</span>requires&nbsp;<a href="#humidity">humidity</a>
    431423=<span style="font-style: italic;"> .TRUE.</span> .<br>
    432424Detailed information about the condensation scheme is given in the
     
    890882the so-called Blackadar mixing length is used (l = kappa * z / ( 1 +
    891883kappa * z / lambda ) with the limiting value lambda = 2.7E-4 * u_g / f).<br>
    892 </td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="moisture"></a><b>moisture</b></p>
     884</td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="humidity"></a><b>humidity</b></p>
    893885</td> <td style="vertical-align: top;">L</td>
    894886<td style="vertical-align: top;"><i>.F.</i></td>
     
    953945<br> <span style="font-weight: bold;">Important: </span>The&nbsp;
    954946upstream-spline scheme is not implemented for humidity and passive
    955 scalars (see <a href="#moisture">moisture</a>
     947scalars (see&nbsp;<a href="#humidity">humidity</a>
    956948and <a href="#passive_scalar">passive_scalar</a>)
    957949and requires the use of a 2d-domain-decomposition. The last conditions
     
    12061198With <span style="font-weight: bold;">passive_scalar</span>
    12071199switched
    1208 on, the simultaneous use of humidity (see <a href="#moisture">moisture</a>)
     1200on, the simultaneous use of humidity (see&nbsp;<a href="#humidity">humidity</a>)
    12091201is impossible.</p> </td> </tr> <tr> <td style="vertical-align: top;"> <p><a name="phi"></a><b>phi</b></p>
    12101202</td> <td style="vertical-align: top;">R</td>
     
    13671359<br> <b>_level</b></p> </td> <td style="vertical-align: top;">R (10)</td> <td style="vertical-align: top;"> <p><i>10 *</i>&nbsp;
    13681360<i>0.0</i></p> </td> <td style="vertical-align: top;"> <p>Height level from
    1369 which on the moisture gradient defined by <a href="#q_vertical_gradient">q_vertical_gradient</a>
     1361which on the humidity gradient defined by <a href="#q_vertical_gradient">q_vertical_gradient</a>
    13701362is effective (in m).&nbsp; </p> <p>The height levels
    13711363are to be assigned in ascending order. The
     
    15331525</p> <p style="margin-left: 40px;"><span style="font-weight: bold;">Important: </span>The&nbsp;
    15341526upstream-spline scheme is not implemented for humidity and passive
    1535 scalars (see <a href="#moisture">moisture</a>
     1527scalars (see&nbsp;<a href="#humidity">humidity</a>
    15361528and <a href="#passive_scalar">passive_scalar</a>)
    15371529and requires the use of a 2d-domain-decomposition. The last conditions
     
    18231815<i>&nbsp;</i><a href="#alpha_surface">alpha_surface</a>
    18241816= 0.0, <a href="#bc_lr">bc_lr</a> = <a href="#bc_ns">bc_ns</a> = <span style="font-style: italic;">'cyclic'</span>,&nbsp;<a style="" href="#galilei_transformation">galilei_transformation</a>
    1825 = <span style="font-style: italic;">.F.</span>,&nbsp;<a href="#cloud_physics">cloud_physics&nbsp;</a> = <span style="font-style: italic;">.F.</span>,&nbsp; <a href="#cloud_droplets">cloud_droplets</a> = <span style="font-style: italic;">.F.</span>,&nbsp; <a href="#moisture">moisture</a> = <span style="font-style: italic;">.F.</span>, and <a href="#prandtl_layer">prandtl_layer</a> = .T..<br>
     1817= <span style="font-style: italic;">.F.</span>,&nbsp;<a href="#cloud_physics">cloud_physics&nbsp;</a> = <span style="font-style: italic;">.F.</span>,&nbsp; <a href="#cloud_droplets">cloud_droplets</a> = <span style="font-style: italic;">.F.</span>,&nbsp;&nbsp;<a href="#humidity">humidity</a> = <span style="font-style: italic;">.F.</span>, and <a href="#prandtl_layer">prandtl_layer</a> = .T..<br>
    18261818<font color="#000000"><br>
    18271819Note that an inclined model domain requires the use of <span style="font-weight: bold;">topography</span> = <span style="font-style: italic;">'flat'</span> and a
  • palm/trunk/DOC/app/chapter_4.2.html

    r72 r75  
    394394= <span style="font-style: italic;">.TRUE.</span></td></tr><tr><td style="width: 106px; vertical-align: top;"><span style="font-style: italic;">pt</span></td><td style="width: 196px; vertical-align: top;">potential
    395395temperature<br></td><td style="vertical-align: top;">K</td><td style="vertical-align: top;"></td></tr><tr><td style="width: 106px; vertical-align: top;"><span style="font-style: italic;">q</span></td><td style="width: 196px; vertical-align: top;">specific humidity
    396 (or total water content, if cloud physics is switched on)</td><td style="vertical-align: top;">kg/kg</td><td style="vertical-align: top;">requires <a href="chapter_4.1.html#moisture">moisture</a> = <span style="font-style: italic;">.TRUE.</span></td></tr><tr><td style="width: 106px; vertical-align: top;"><span style="font-style: italic;">ql</span></td><td style="width: 196px; vertical-align: top;">liquid water
     396(or total water content, if cloud physics is switched on)</td><td style="vertical-align: top;">kg/kg</td><td style="vertical-align: top;">requires&nbsp;<a href="chapter_4.1.html#humidity">humidity</a> = <span style="font-style: italic;">.TRUE.</span></td></tr><tr><td style="width: 106px; vertical-align: top;"><span style="font-style: italic;">ql</span></td><td style="width: 196px; vertical-align: top;">liquid water
    397397content</td><td style="vertical-align: top;">kg/kg</td><td style="vertical-align: top;">requires <a href="chapter_4.1.html#cloud_physics">cloud_physics</a>
    398398= <span style="font-style: italic;">.TRUE.</span>
     
    414414is allowed</td></tr><tr><td style="width: 106px; vertical-align: top;"><span style="font-style: italic;">v</span></td><td style="width: 196px; vertical-align: top;">v-component of
    415415the velocity</td><td style="vertical-align: top;">m/s</td><td style="vertical-align: top;"></td></tr><tr><td style="width: 106px; vertical-align: top;"><span style="font-style: italic;">vpt</span></td><td style="width: 196px; vertical-align: top;">virtual potential
    416 temperature</td><td style="vertical-align: top;">K</td><td style="vertical-align: top;">requires <a href="chapter_4.1.html#moisture">moisture</a> = <span style="font-style: italic;">.TRUE.</span></td></tr><tr><td style="width: 106px; vertical-align: top;"><span style="font-style: italic;">w</span></td><td style="width: 196px; vertical-align: top;">w-component of
     416temperature</td><td style="vertical-align: top;">K</td><td style="vertical-align: top;">requires&nbsp;<a href="chapter_4.1.html#humidity">humidity</a> = <span style="font-style: italic;">.TRUE.</span></td></tr><tr><td style="width: 106px; vertical-align: top;"><span style="font-style: italic;">w</span></td><td style="width: 196px; vertical-align: top;">w-component of
    417417the velocity</td><td style="vertical-align: top;">m/s</td><td style="vertical-align: top;"></td></tr><tr><td style="vertical-align: top;"><span style="font-style: italic;">z0*</span></td><td style="vertical-align: top;">roughness length</td><td style="vertical-align: top;">m</td><td></td></tr></tbody></table><br>Multiple
    418418quantities can be assigned, e.g. <span style="font-weight: bold;">data_output</span>
  • palm/trunk/DOC/app/chapter_4.6.html

    r72 r75  
    652652<td style="vertical-align: middle;" width="57%"> <p>Variable
    653653to adjust the subdomain sizes in parallel runs.</p> </td> </tr>
    654 <tr> <td style="vertical-align: middle;"><b><a href="chapter_4.1.html#inflow_disturbance_begin"><b>inflow_disturbance_begin</b></a></b></td>
     654<tr><td><b><a href="chapter_4.1.html#humidity"><b>humidity</b></a></b></td><td>I</td><td>L</td><td><span style="font-style: italic;">.F.</span></td><td>Parameter
     655to switch on the prognostic equation for
     656specific
     657humidity q.</td></tr><tr> <td style="vertical-align: middle;"><b><a href="chapter_4.1.html#inflow_disturbance_begin"><b>inflow_disturbance_begin</b></a></b></td>
    655658<td style="vertical-align: middle;">I<br> </td>
    656659<td style="vertical-align: middle;">I<br> </td>
     
    749752which are
    750753to be created by the <span style="font-weight: bold;">dvrp</span>
    751 software. <br> </p> </td> </tr> <tr>
    752 <td style="vertical-align: middle;" width="15%"> <p><a href="chapter_4.1.html#moisture"><b>moisture</b></a></p>
    753 </td> <td style="vertical-align: middle;" width="5%">
    754 <p>I</p> </td> <td style="vertical-align: middle;" width="7%"> <p>L</p>
    755 </td> <td style="vertical-align: middle;" width="16%">
    756 <p><i>.F.</i></p> </td> <td style="vertical-align: middle;" width="57%"> <p>Parameter
    757 to switch on the prognostic equation for
    758 specific
    759 humidity q.</p> </td> </tr> <tr> <td style="vertical-align: middle;" width="15%"> <p><a href="chapter_4.1.html#momentum_advec"><b>momentum_advec</b></a></p>
     754software. <br> </p> </td> </tr>  <tr> <td style="vertical-align: middle;" width="15%"> <p><a href="chapter_4.1.html#momentum_advec"><b>momentum_advec</b></a></p>
    760755</td> <td style="vertical-align: middle;" width="5%">
    761756<p>I</p> </td> <td style="vertical-align: middle;" width="7%"> <p>C
     
    11051100(10)</p> </td> <td style="vertical-align: middle;" width="16%"> <p><i>10 *</i>&nbsp; <i>0.0</i></p>
    11061101</td> <td style="vertical-align: middle;" width="57%">
    1107 <p>Height level from which on the moisture gradient defined
     1102<p>Height level from which on the humidity gradient defined
    11081103by <a href="chapter_4.1.html#q_vertical_gradient">q_vertical_gradient</a>
    11091104is effective (in m).&nbsp;</p> </td> </tr> <tr>
  • palm/trunk/SCRIPTS/mbuild

    r69 r75  
    13951395          fi    # ENDE UEBERSETZUNG DER UTILITY-PROGRAMME
    13961396
     1397          rm -rf  ${remote_host}_last_make_protokoll
    13971398
    13981399             # NUR AUF EINEM HLRN-RECHNER UEBERSETZEN
  • palm/trunk/SCRIPTS/mrun

    r74 r75  
    26882688       # NAMELIST-DATEI MIT WERTEN VON ENVIRONMENT-VARIABLEN ERZEUGEN (ZU
    26892689       # LESEN VON PALM)
     2690       # DAZU VORHER DIE "GLOBAL REVISION" VON SUBVERSION ERMITTELN, FALLS
     2691       # DIE SOURCEN UNTER SUBVERSION-KONTROLLE STEHEN
     2692    global_revision=`svnversion $source_path`
     2693    global_revision="Rev: $global_revision"
     2694
    26902695    cat  >  ENVPAR  <<  %%END%%
    26912696 &envpar  run_identifier = '$fname', host = '$localhost',
    26922697          write_binary = '$write_binary', tasks_per_node = $tasks_per_node,
    2693           maximum_cpu_time_allowed = ${cpumax}. /
     2698          maximum_cpu_time_allowed = ${cpumax}.,
     2699          revision = '$global_revision' /
    26942700
    26952701%%END%%
  • palm/trunk/SOURCE/CURRENT_MODIFICATIONS

    r73 r75  
    3333Changed:
    3434-------
    35 General revision of non-cyclic horizontal boundary conditions: radiation boundary conditions are now used instead of Neumann conditions at the outflow (calculation needs velocity values for t-dt, which are stored on new arrays u_m_l, u_m_r, etc.), calculation of mean outflow is not needed any more, additional gridpoints along x and y (uxrp, vynp) are not needed any more, routine "boundary_conds" now operates on timelevel t+dt and is not split in two parts (main, uvw_outflow) any more, Neumann boundary conditions at inflow/outflow in case of non-cyclic boundary conditions for all 2d-arrays that are handled by exchange_horiz_2d
     35General revision of non-cyclic horizontal boundary conditions: radiation boundary conditions are now used instead of Neumann conditions at the outflow (calculation needs velocity values for t-dt, which are stored on new arrays u_m_l, u_m_r, etc.), calculation of mean outflow is not needed any more, volume flow control is added for the outflow boundary (currently only for the north boundary!!), additional gridpoints along x and y (uxrp, vynp) are not needed any more, routine "boundary_conds" now operates on timelevel t+dt and is not split in two parts (main, uvw_outflow) any more, Neumann boundary conditions at inflow/outflow in case of non-cyclic boundary conditions for all 2d-arrays that are handled by exchange_horiz_2d
    3636
    3737+age_m in particle_type
     
    4848q is not allowed to become negative (prognostic_equations).
    4949
     50poisfft_init is only called if fft-solver is switched on (init_pegrid).
     51
     52d3par-parameter moisture renamed to humidity.
     53
     54Subversion global revision number is read from mrun and added to the run description header and to the run control (_rc) file.
     55
    5056__vtk directives removed from main program.
    5157
    5258The uitility routine interpret_config reads PALM environment variables from NAMELIST instead using the system call GETENV.
    5359
    54 check_parameters, data_output_dvrp, data_output_ptseries, data_output_ts, exchange_horiz_2d, flow_statistics, header, init_particles, init_pegrid, init_3d_model, modules, palm, package_parin, parin, prognostic_equations, read_3d_binary, time_integration, write_3d_binary
     60advec_u_pw, advec_u_up, advec_v_pw, advec_v_up, asselin_filter, check_parameters, coriolis, data_output_dvrp, data_output_ptseries, data_output_ts, data_output_2d, data_output_3d, diffusion_u, diffusion_v, exchange_horiz, exchange_horiz_2d, flow_statistics, header, init_grid, init_particles, init_pegrid, init_rankine, init_pt_anomaly, init_1d_model, init_3d_model, modules, palm, package_parin, parin, poismg, prandtl_fluxes, pres, production_e, prognostic_equations, read_var_list, read_3d_binary, sor, swap_timelevel, time_integration, write_var_list, write_3d_binary
    5561
    5662
  • palm/trunk/SOURCE/advec_particles.f90

    r64 r75  
    99! + user_advec_particles, particles-package is now part of the defaut code,
    1010! array arguments in sendrecv calls have to refer to first element (1) due to
    11 ! mpich (mpiI) interface requirements
     11! mpich (mpiI) interface requirements,
     12! 2nd+3rd argument removed from exchange horiz
    1213! TEST: PRINT statements on unit 9 (commented out)
    1314!
     
    662663!
    663664!--    Lateral boundary conditions
    664        CALL exchange_horiz( de_dx, 0, 0 )
    665        CALL exchange_horiz( de_dy, 0, 0 )
    666        CALL exchange_horiz( de_dz, 0, 0 )
    667        CALL exchange_horiz( diss,  0, 0 )
     665       CALL exchange_horiz( de_dx )
     666       CALL exchange_horiz( de_dy )
     667       CALL exchange_horiz( de_dz )
     668       CALL exchange_horiz( diss  )
    668669
    669670!
  • palm/trunk/SOURCE/advec_u_pw.f90

    r4 r75  
    44! Actual revisions:
    55! -----------------
    6 !
     6! uxrp eliminated
    77!
    88! Former revisions:
     
    5454       gu = 2.0 * u_gtrans
    5555       gv = 2.0 * v_gtrans
    56        DO  i = nxl, nxr+uxrp
     56       DO  i = nxl, nxr
    5757          DO  j = nys, nyn
    5858             DO  k = nzb_u_inner(j,i)+1, nzt
  • palm/trunk/SOURCE/advec_u_up.f90

    r4 r75  
    44! Actual revisions:
    55! -----------------
    6 !
     6! uxrp eliminated
    77!
    88! Former revisions:
     
    5353
    5454
    55        DO  i = nxl, nxr+uxrp
     55       DO  i = nxl, nxr
    5656          DO  j = nys, nyn
    5757             DO  k = nzb_u_inner(j,i)+1, nzt
  • palm/trunk/SOURCE/advec_v_pw.f90

    r4 r75  
    44! Actual revisions:
    55! -----------------
    6 !
     6! vynp eliminated
    77!
    88! Former revisions:
     
    5656       gv = 2.0 * v_gtrans
    5757       DO  i = nxl, nxr
    58           DO  j = nys, nyn+vynp
     58          DO  j = nys, nyn
    5959             DO  k = nzb_v_inner(j,i)+1, nzt
    6060                tend(k,j,i) = tend(k,j,i) - 0.25 * (                           &
  • palm/trunk/SOURCE/advec_v_up.f90

    r4 r75  
    44! Actual revisions:
    55! -----------------
    6 !
     6! vynp eliminated
    77!
    88! Former revisions:
     
    5353
    5454       DO  i = nxl, nxr
    55           DO  j = nys, nyn+vynp
     55          DO  j = nys, nyn
    5656             DO  k = nzb_v_inner(j,i)+1, nzt
    5757!
  • palm/trunk/SOURCE/asselin_filter.f90

    r4 r75  
    44! Actual revisions:
    55! -----------------
    6 !
     6! moisture renamed humidity
    77!
    88! Former revisions:
     
    7272          ENDIF
    7373
    74           IF ( ( moisture .OR. passive_scalar )  .AND. &
     74          IF ( ( humidity .OR. passive_scalar )  .AND. &
    7575             scalar_advec /= 'bc-scheme' )  THEN
    7676             DO  k = nzb_2d(j,i), nzt+1
     
    9696    ENDIF
    9797
    98     IF ( ( moisture .OR. passive_scalar )  .AND. &
     98    IF ( ( humidity .OR. passive_scalar )  .AND. &
    9999         scalar_advec /= 'bc-scheme' )  THEN
    100100       q = q + asselin_filter_factor * ( q_p - 2.0 * q + q_m )
  • palm/trunk/SOURCE/boundary_conds.f90

    r73 r75  
    44! Actual revisions:
    55! -----------------
    6 ! The "main" part sets conditions for time level t+dt insteat of level t,
    7 ! outflow boundary conditions changed from Neumann to radiation condition
     6! The "main" part sets conditions for time level t+dt instead of level t,
     7! outflow boundary conditions changed from Neumann to radiation condition,
     8! uxrp, vynp eliminated, moisture renamed humidity
    89!
    910! Former revisions:
     
    126127!--    Boundary conditions for total water content or scalar,
    127128!--    bottom and surface boundary (see also temperature)
    128        IF ( moisture  .OR.  passive_scalar )  THEN
    129 !
    130 !--       Surface conditions for constant_moisture_flux
     129       IF ( humidity  .OR.  passive_scalar )  THEN
     130!
     131!--       Surface conditions for constant_humidity_flux
    131132          IF ( ibc_q_b == 0 ) THEN
    132133             DO  i = nxl-1, nxr+1
     
    154155          v_p(:,nys,:) = v_p(:,nys-1,:)
    155156       ELSEIF ( inflow_n ) THEN
    156           v_p(:,nyn+vynp,:) = v_p(:,nyn+vynp+1,:)
     157          v_p(:,nyn,:) = v_p(:,nyn+1,:)
    157158       ELSEIF ( inflow_l ) THEN
    158159          u_p(:,:,nxl) = u_p(:,:,nxl-1)
    159160       ELSEIF ( inflow_r ) THEN
    160           u_p(:,:,nxr+uxrp) = u_p(:,:,nxr+uxrp+1)
     161          u_p(:,:,nxr) = u_p(:,:,nxr+1)
    161162       ENDIF
    162163
     
    166167          pt_p(:,nys-1,:)     = pt_p(:,nys,:)
    167168          IF ( .NOT. constant_diffusion     )  e_p(:,nys-1,:) = e_p(:,nys,:)
    168           IF ( moisture .OR. passive_scalar )  q_p(:,nys-1,:) = q_p(:,nys,:)
     169          IF ( humidity .OR. passive_scalar )  q_p(:,nys-1,:) = q_p(:,nys,:)
    169170       ELSEIF ( outflow_n )  THEN
    170171          pt_p(:,nyn+1,:)     = pt_p(:,nyn,:)
    171172          IF ( .NOT. constant_diffusion     )  e_p(:,nyn+1,:) = e_p(:,nyn,:)
    172           IF ( moisture .OR. passive_scalar )  q_p(:,nyn+1,:) = q_p(:,nyn,:)
     173          IF ( humidity .OR. passive_scalar )  q_p(:,nyn+1,:) = q_p(:,nyn,:)
    173174       ELSEIF ( outflow_l )  THEN
    174175          pt_p(:,:,nxl-1)     = pt_p(:,:,nxl)
    175176          IF ( .NOT. constant_diffusion     )  e_p(:,:,nxl-1) = e_p(:,:,nxl)
    176           IF ( moisture .OR. passive_scalar )  q_p(:,:,nxl-1) = q_p(:,:,nxl)     
     177          IF ( humidity .OR. passive_scalar )  q_p(:,:,nxl-1) = q_p(:,:,nxl)
    177178       ELSEIF ( outflow_r )  THEN
    178179          pt_p(:,:,nxr+1)     = pt_p(:,:,nxr)
    179180          IF ( .NOT. constant_diffusion     )  e_p(:,:,nxr+1) = e_p(:,:,nxr)
    180           IF ( moisture .OR. passive_scalar )  q_p(:,:,nxr+1) = q_p(:,:,nxr)
     181          IF ( humidity .OR. passive_scalar )  q_p(:,:,nxr+1) = q_p(:,:,nxr)
    181182       ENDIF
    182183
    183184    ENDIF
    184185
    185     IF ( range == 'outflow_uvw' ) THEN
    186 !
    187 !--    Radiation boundary condition for the velocities at the respective outflow
    188        IF ( outflow_s ) THEN
    189 !          v(:,nys-1,:) = v(:,nys,:)
    190 !          w(:,nys-1,:) = 0.0
    191 !!
    192 !!--       Compute the mean horizontal wind parallel to and within the outflow
    193 !!--       wall and use this as boundary condition for u
    194 !#if defined( __parallel )
    195 !          CALL MPI_ALLREDUCE( uvmean_outflow_l, uvmean_outflow, nzt-nzb+2, &
    196 !                              MPI_REAL, MPI_SUM, comm1dx, ierr )
    197 !          uvmean_outflow = uvmean_outflow / ( nx + 1.0 )
    198 !#else
    199 !          uvmean_outflow = uvmean_outflow_l / ( nx + 1.0 )
    200 !#endif
    201 !          DO  k = nzb, nzt+1
    202 !             u(k,nys-1,:) = uvmean_outflow(k)
    203 !          ENDDO
    204        ENDIF
    205 
    206        IF ( outflow_n  .AND.                                                 &
    207             intermediate_timestep_count == intermediate_timestep_count_max ) &
    208        THEN
    209 
    210           c_max = dy / dt_3d
    211 
    212           DO i = nxl-1, nxr+1
    213              DO k = nzb+1, nzt+1
    214 
    215 !
    216 !--             First calculate the phase speeds for u,v, and w
    217                 denom = u_m_n(k,ny,i,-2) - u_m_n(k,ny-1,i,-2)
    218 
    219                 IF ( denom /= 0.0 )  THEN
    220                    c_u = -c_max * ( u_m_n(k,ny,i,-1)-u_m_n(k,ny,i,-2) ) / denom
    221                    IF ( c_u < 0.0 )  THEN
    222                       c_u = 0.0
    223                    ELSEIF ( c_u > c_max )  THEN
    224                       c_u = c_max
    225                    ENDIF
    226                 ELSE
     186!
     187!-- Radiation boundary condition for the velocities at the respective outflow
     188    IF ( outflow_s  .AND.                                                 &
     189         intermediate_timestep_count == intermediate_timestep_count_max ) &
     190    THEN
     191
     192       c_max = dy / dt_3d
     193
     194       DO i = nxl-1, nxr+1
     195          DO k = nzb+1, nzt+1
     196
     197!
     198!--          First calculate the phase speeds for u,v, and w
     199             denom = u_m_s(k,0,i) - u_m_s(k,1,i)
     200
     201             IF ( denom /= 0.0 )  THEN
     202                c_u = -c_max * ( u(k,0,i) - u_m_s(k,0,i) ) / denom
     203                IF ( c_u < 0.0 )  THEN
     204                   c_u = 0.0
     205                ELSEIF ( c_u > c_max )  THEN
    227206                   c_u = c_max
    228207                ENDIF
    229 
    230                 denom = v_m_n(k,ny,i,-2) - v_m_n(k,ny-1,i,-2)
    231 
    232                 IF ( denom /= 0.0 )  THEN
    233                    c_v = -c_max * ( v_m_n(k,ny,i,-1)-v_m_n(k,ny,i,-2) ) / denom
    234                    IF ( c_v < 0.0 )  THEN
    235                       c_v = 0.0
    236                    ELSEIF ( c_v > c_max )  THEN
    237                       c_v = c_max
    238                    ENDIF
    239                 ELSE
     208             ELSE
     209                c_u = c_max
     210             ENDIF
     211             denom = v_m_s(k,0,i) - v_m_s(k,1,i)
     212
     213             IF ( denom /= 0.0 )  THEN
     214                c_v = -c_max * ( v(k,0,i) - v_m_s(k,0,i) ) / denom
     215                IF ( c_v < 0.0 )  THEN
     216                   c_v = 0.0
     217                ELSEIF ( c_v > c_max )  THEN
    240218                   c_v = c_max
    241219                ENDIF
    242 
    243                 denom = w_m_n(k,ny,i,-2) - w_m_n(k,ny-1,i,-2)
    244 
    245                 IF ( denom /= 0.0 )  THEN
    246                    c_w = -c_max * ( w_m_n(k,ny,i,-1)-w_m_n(k,ny,i,-2) ) / denom
    247                    IF ( c_w < 0.0 )  THEN
    248                       c_w = 0.0
    249                    ELSEIF ( c_w > c_max )  THEN
    250                       c_w = c_max
    251                    ENDIF
    252                 ELSE
     220             ELSE
     221                c_v = c_max
     222             ENDIF
     223
     224             denom = w_m_s(k,0,i) - w_m_s(k,1,i)
     225
     226             IF ( denom /= 0.0 )  THEN
     227                c_w = -c_max * ( w(k,0,i) - w_m_s(k,0,i) ) / denom
     228                IF ( c_w < 0.0 )  THEN
     229                   c_w = 0.0
     230                ELSEIF ( c_w > c_max )  THEN
    253231                   c_w = c_max
    254232                ENDIF
    255 
    256 !
    257 !--             Calculate the new velocities
    258                 u(k,ny+1,i) = u_m_n(k,ny+1,i,-1) - dt_3d * c_u * &
    259                               ( u_m_n(k,ny+1,i,-1) - u_m_n(k,ny,i,-1) ) * ddy
    260 
    261                 v(k,ny+1,i) = v_m_n(k,ny+1,i,-1) - dt_3d * c_v * &
    262                               ( v_m_n(k,ny+1,i,-1) - v_m_n(k,ny,i,-1) ) * ddy
    263 
    264                 w(k,ny+1,i) = w_m_n(k,ny+1,i,-1) - dt_3d * c_w * &
    265                               ( w_m_n(k,ny+1,i,-1) - w_m_n(k,ny,i,-1) ) * ddy
    266 
    267 !
    268 !--             Swap timelevels for the next timestep
    269                 u_m_n(k,:,i,-2) = u_m_n(k,:,i,-1)
    270                 u_m_n(k,:,i,-1) = u(k,ny-1:ny+1,i)
    271                 v_m_n(k,:,i,-2) = v_m_n(k,:,i,-1)
    272                 v_m_n(k,:,i,-1) = v(k,ny-1:ny+1,i)
    273                 w_m_n(k,:,i,-2) = w_m_n(k,:,i,-1)
    274                 w_m_n(k,:,i,-1) = w(k,ny-1:ny+1,i)
    275 
    276              ENDDO
    277           ENDDO
    278 
    279 !
    280 !--       Bottom boundary at the outflow
    281           IF ( ibc_uv_b == 0 )  THEN
    282              u(nzb,ny+1,:) = -u(nzb+1,ny+1,:)
    283              v(nzb,ny+1,:) = -v(nzb+1,ny+1,:) 
    284           ELSE                   
    285              u(nzb,ny+1,:) =  u(nzb+1,ny+1,:)
    286              v(nzb,ny+1,:) =  v(nzb+1,ny+1,:)
    287           ENDIF
    288           w(nzb,ny+1,:) = 0.0
    289 
    290 !
    291 !--       Top boundary at the outflow
    292           IF ( ibc_uv_t == 0 )  THEN
    293              u(nzt+1,ny+1,:) = ug(nzt+1)
    294              v(nzt+1,ny+1,:) = vg(nzt+1)
    295           ELSE
    296              u(nzt+1,ny+1,:) = u(nzt,nyn+1,:)
    297              v(nzt+1,ny+1,:) = v(nzt,nyn+1,:)
    298           ENDIF
    299           w(nzt:nzt+1,ny+1,:) = 0.0
    300 
    301        ENDIF
    302 
    303        IF ( outflow_l ) THEN
    304 !          u(:,:,nxl-1) = u(:,:,nxl)
    305 !          w(:,:,nxl-1) = 0.0
    306 !
    307 !--       Compute the mean horizontal wind parallel to and within the outflow
    308 !--       wall and use this as boundary condition for v
    309 !#if defined( __parallel )
    310 !          CALL MPI_ALLREDUCE( uvmean_outflow_l, uvmean_outflow, nzt-nzb+2, &
    311 !                              MPI_REAL, MPI_SUM, comm1dy, ierr )
    312 !          uvmean_outflow = uvmean_outflow / ( ny + 1.0 )
    313 !#else
    314 !          uvmean_outflow = uvmean_outflow_l / ( ny + 1.0 )
    315 !#endif
    316 !          DO  k = nzb, nzt+1
    317 !             v(k,:,nxl-1) = uvmean_outflow(k)
    318 !          ENDDO   
    319 !
    320        ENDIF
    321 
    322        IF ( outflow_r  .AND.                                                 &
    323             intermediate_timestep_count == intermediate_timestep_count_max ) &
    324        THEN
    325 
    326           c_max = dx / dt_3d
    327 
    328           DO j = nys-1, nyn+1
    329              DO k = nzb+1, nzt+1
    330 
    331 !
    332 !--             First calculate the phase speeds for u,v, and w
    333                 denom = u_m_r(k,j,nx,-2) - u_m_r(k,j,nx-1,-2)
    334 
    335                 IF ( denom /= 0.0 )  THEN
    336                    c_u = -c_max * ( u_m_r(k,j,nx,-1)-u_m_r(k,j,nx,-2) ) / denom
    337                    IF ( c_u < 0.0 )  THEN
    338                       c_u = 0.0
    339                    ELSEIF ( c_u > c_max )  THEN
    340                       c_u = c_max
    341                    ENDIF
    342                 ELSE
     233             ELSE
     234                c_w = c_max
     235             ENDIF
     236
     237!
     238!--          Calculate the new velocities
     239             u_p(k,-1,i) = u(k,-1,i) + dt_3d * c_u * &
     240                                       ( u(k,-1,i) - u(k,0,i) ) * ddy
     241
     242             v_p(k,-1,i) = v(k,-1,i) + dt_3d * c_v * &
     243                                       ( v(k,-1,i) - v_m_s(k,0,i) ) * ddy
     244
     245             w_p(k,-1,i) = w(k,-1,i) + dt_3d * c_w * &
     246                                       ( w(k,-1,i) - w(k,0,i) ) * ddy
     247
     248!
     249!--          Save old timelevels for the next timestep
     250             u_m_s(k,:,i) = u(k,-1:1,i)
     251             v_m_s(k,:,i) = v(k,-1:1,i)
     252             w_m_s(k,:,i) = w(k,-1:1,i)
     253
     254          ENDDO
     255       ENDDO
     256
     257!
     258!--    Bottom boundary at the outflow
     259       IF ( ibc_uv_b == 0 )  THEN
     260          u_p(nzb,-1,:) = -u_p(nzb+1,-1,:)
     261          v_p(nzb,-1,:) = -v_p(nzb+1,-1,:) 
     262       ELSE                   
     263          u_p(nzb,-1,:) =  u_p(nzb+1,-1,:)
     264          v_p(nzb,-1,:) =  v_p(nzb+1,-1,:)
     265       ENDIF
     266       w_p(nzb,ny+1,:) = 0.0
     267
     268!
     269!--    Top boundary at the outflow
     270       IF ( ibc_uv_t == 0 )  THEN
     271          u_p(nzt+1,-1,:) = ug(nzt+1)
     272          v_p(nzt+1,-1,:) = vg(nzt+1)
     273       ELSE
     274          u_p(nzt+1,-1,:) = u(nzt,-1,:)
     275          v_p(nzt+1,-1,:) = v(nzt,-1,:)
     276       ENDIF
     277       w_p(nzt:nzt+1,-1,:) = 0.0
     278
     279    ENDIF
     280
     281    IF ( outflow_n  .AND.                                                 &
     282         intermediate_timestep_count == intermediate_timestep_count_max ) &
     283    THEN
     284
     285       c_max = dy / dt_3d
     286
     287       DO i = nxl-1, nxr+1
     288          DO k = nzb+1, nzt+1
     289
     290!
     291!--          First calculate the phase speeds for u,v, and w
     292             denom = u_m_n(k,ny,i) - u_m_n(k,ny-1,i)
     293
     294             IF ( denom /= 0.0 )  THEN
     295                c_u = -c_max * ( u(k,ny,i) - u_m_n(k,ny,i) ) / denom
     296                IF ( c_u < 0.0 )  THEN
     297                   c_u = 0.0
     298                ELSEIF ( c_u > c_max )  THEN
    343299                   c_u = c_max
    344300                ENDIF
    345 
    346                 denom = v_m_r(k,j,nx,-2) - v_m_r(k,j,nx-1,-2)
    347 
    348                 IF ( denom /= 0.0 )  THEN
    349                    c_v = -c_max * ( v_m_r(k,j,nx,-1)-v_m_r(k,j,nx,-2) ) / denom
    350                    IF ( c_v < 0.0 )  THEN
    351                       c_v = 0.0
    352                    ELSEIF ( c_v > c_max )  THEN
    353                       c_v = c_max
    354                    ENDIF
    355                 ELSE
     301             ELSE
     302                c_u = c_max
     303             ENDIF
     304
     305             denom = v_m_n(k,ny,i) - v_m_n(k,ny-1,i)
     306
     307             IF ( denom /= 0.0 )  THEN
     308                c_v = -c_max * ( v(k,ny,i) - v_m_n(k,ny,i) ) / denom
     309                IF ( c_v < 0.0 )  THEN
     310                   c_v = 0.0
     311                ELSEIF ( c_v > c_max )  THEN
    356312                   c_v = c_max
    357313                ENDIF
    358 
    359                 denom = w_m_r(k,j,nx,-2) - w_m_r(k,j,nx-1,-2)
    360 
    361                 IF ( denom /= 0.0 )  THEN
    362                    c_w = -c_max * ( w_m_r(k,j,nx,-1)-w_m_n(k,j,nx,-2) ) / denom
    363                    IF ( c_w < 0.0 )  THEN
    364                       c_w = 0.0
    365                    ELSEIF ( c_w > c_max )  THEN
    366                       c_w = c_max
    367                    ENDIF
    368                 ELSE
     314             ELSE
     315                c_v = c_max
     316             ENDIF
     317
     318             denom = w_m_n(k,ny,i) - w_m_n(k,ny-1,i)
     319
     320             IF ( denom /= 0.0 )  THEN
     321                c_w = -c_max * ( w(k,ny,i) - w_m_n(k,ny,i) ) / denom
     322                IF ( c_w < 0.0 )  THEN
     323                   c_w = 0.0
     324                ELSEIF ( c_w > c_max )  THEN
    369325                   c_w = c_max
    370326                ENDIF
    371 
    372 !
    373 !--             Calculate the new velocities
    374                 u(k,j,nx+1) = u_m_r(k,j,nx+1,-1) - dt_3d * c_u * &
    375                               ( u_m_r(k,j,nx+1,-1) - u_m_r(k,j,nx,-1) ) * ddx
    376 
    377                 v(k,j,nx+1) = v_m_r(k,j,nx+1,-1) - dt_3d * c_v * &
    378                               ( v_m_r(k,j,nx+1,-1) - v_m_r(k,j,nx,-1) ) * ddx
    379 
    380                 w(k,j,nx+1) = w_m_r(k,j,nx+1,-1) - dt_3d * c_w * &
    381                               ( w_m_r(k,j,nx+1,-1) - w_m_r(k,j,nx,-1) ) * ddx
    382 
    383 !
    384 !--             Swap timelevels for the next timestep
    385                 u_m_r(k,j,:,-2) = u_m_r(k,j,:,-1)
    386                 u_m_r(k,j,:,-1) = u(k,j,nx-1:nx+1)
    387                 v_m_r(k,j,:,-2) = v_m_r(k,j,:,-1)
    388                 v_m_r(k,j,:,-1) = v(k,j,nx-1:nx+1)
    389                 w_m_r(k,j,:,-2) = w_m_r(k,j,:,-1)
    390                 w_m_r(k,j,:,-1) = w(k,j,nx-1:nx+1)
    391 
    392              ENDDO
    393           ENDDO
    394 
    395 !
    396 !--       Bottom boundary at the outflow
    397           IF ( ibc_uv_b == 0 )  THEN
    398              u(nzb,ny+1,:) = -u(nzb+1,ny+1,:)
    399              v(nzb,ny+1,:) = -v(nzb+1,ny+1,:) 
    400           ELSE                   
    401              u(nzb,ny+1,:) =  u(nzb+1,ny+1,:)
    402              v(nzb,ny+1,:) =  v(nzb+1,ny+1,:)
    403           ENDIF
    404           w(nzb,ny+1,:) = 0.0
    405 
    406 !
    407 !--       Top boundary at the outflow
    408           IF ( ibc_uv_t == 0 )  THEN
    409              u(nzt+1,ny+1,:) = ug(nzt+1)
    410              v(nzt+1,ny+1,:) = vg(nzt+1)
    411           ELSE
    412              u(nzt+1,ny+1,:) = u(nzt,nyn+1,:)
    413              v(nzt+1,ny+1,:) = v(nzt,nyn+1,:)
    414           ENDIF
    415           w(nzt:nzt+1,ny+1,:) = 0.0
    416 
    417        ENDIF
     327             ELSE
     328                c_w = c_max
     329             ENDIF
     330
     331!
     332!--          Calculate the new velocities
     333             u_p(k,ny+1,i) = u(k,ny+1,i) - dt_3d * c_u * &
     334                                           ( u(k,ny+1,i) - u(k,ny,i) ) * ddy
     335
     336             v_p(k,ny+1,i) = v(k,ny+1,i) - dt_3d * c_v * &
     337                                           ( v(k,ny+1,i) - v(k,ny,i) ) * ddy
     338
     339             w_p(k,ny+1,i) = w(k,ny+1,i) - dt_3d * c_w * &
     340                                           ( w(k,ny+1,i) - w(k,ny,i) ) * ddy
     341
     342!
     343!--          Swap timelevels for the next timestep
     344             u_m_n(k,:,i) = u(k,ny-1:ny+1,i)
     345             v_m_n(k,:,i) = v(k,ny-1:ny+1,i)
     346             w_m_n(k,:,i) = w(k,ny-1:ny+1,i)
     347
     348          ENDDO
     349       ENDDO
     350
     351!
     352!--    Bottom boundary at the outflow
     353       IF ( ibc_uv_b == 0 )  THEN
     354          u_p(nzb,ny+1,:) = -u_p(nzb+1,ny+1,:)
     355          v_p(nzb,ny+1,:) = -v_p(nzb+1,ny+1,:) 
     356       ELSE                   
     357          u_p(nzb,ny+1,:) =  u_p(nzb+1,ny+1,:)
     358          v_p(nzb,ny+1,:) =  v_p(nzb+1,ny+1,:)
     359       ENDIF
     360       w_p(nzb,ny+1,:) = 0.0
     361
     362!
     363!--    Top boundary at the outflow
     364       IF ( ibc_uv_t == 0 )  THEN
     365          u_p(nzt+1,ny+1,:) = ug(nzt+1)
     366          v_p(nzt+1,ny+1,:) = vg(nzt+1)
     367       ELSE
     368          u_p(nzt+1,ny+1,:) = u_p(nzt,nyn+1,:)
     369          v_p(nzt+1,ny+1,:) = v_p(nzt,nyn+1,:)
     370       ENDIF
     371       w_p(nzt:nzt+1,ny+1,:) = 0.0
     372
     373    ENDIF
     374
     375    IF ( outflow_l  .AND.                                                 &
     376         intermediate_timestep_count == intermediate_timestep_count_max ) &
     377    THEN
     378
     379       c_max = dx / dt_3d
     380
     381       DO j = nys-1, nyn+1
     382          DO k = nzb+1, nzt+1
     383
     384!
     385!--          First calculate the phase speeds for u,v, and w
     386             denom = u_m_l(k,j,0) - u_m_l(k,j,1)
     387
     388             IF ( denom /= 0.0 )  THEN
     389                c_u = -c_max * ( u(k,j,0) - u_m_r(k,j,0) ) / denom
     390                IF ( c_u > 0.0 )  THEN
     391                   c_u = 0.0
     392                ELSEIF ( c_u < -c_max )  THEN
     393                   c_u = -c_max
     394                ENDIF
     395             ELSE
     396                c_u = -c_max
     397             ENDIF
     398
     399             denom = v_m_l(k,j,0) - v_m_l(k,j,1)
     400
     401             IF ( denom /= 0.0 )  THEN
     402                c_v = -c_max * ( v(k,j,0) - v_m_l(k,j,0) ) / denom
     403                IF ( c_v < 0.0 )  THEN
     404                   c_v = 0.0
     405                ELSEIF ( c_v > c_max )  THEN
     406                   c_v = c_max
     407                ENDIF
     408             ELSE
     409                c_v = c_max
     410             ENDIF
     411
     412             denom = w_m_l(k,j,0) - w_m_l(k,j,1)
     413
     414             IF ( denom /= 0.0 )  THEN
     415                c_w = -c_max * ( w(k,j,0) - w_m_l(k,j,0) ) / denom
     416                IF ( c_w < 0.0 )  THEN
     417                   c_w = 0.0
     418                ELSEIF ( c_w > c_max )  THEN
     419                   c_w = c_max
     420                ENDIF
     421             ELSE
     422                c_w = c_max
     423             ENDIF
     424
     425!
     426!--          Calculate the new velocities
     427             u_p(k,j,-1) = u(k,j,-1) + dt_3d * c_u * &
     428                                       ( u(k,j,-1) - u(k,j,0) ) * ddx
     429
     430             v_p(k,j,-1) = v(k,j,-1) + dt_3d * c_v * &
     431                                       ( v(k,j,-1) - v(k,j,0) ) * ddx
     432
     433             w_p(k,j,-1) = w(k,j,-1) + dt_3d * c_w * &
     434                                       ( w(k,j,-1) - w(k,j,0) ) * ddx
     435
     436!
     437!--          Swap timelevels for the next timestep
     438             u_m_l(k,j,:) = u(k,j,-1:1)
     439             v_m_l(k,j,:) = v(k,j,-1:1)
     440             w_m_l(k,j,:) = w(k,j,-1:1)
     441
     442          ENDDO
     443       ENDDO
     444
     445!
     446!--    Bottom boundary at the outflow
     447       IF ( ibc_uv_b == 0 )  THEN
     448          u_p(nzb,:,-1) = -u_p(nzb+1,:,-1)
     449          v_p(nzb,:,-1) = -v_p(nzb+1,:,-1) 
     450       ELSE                   
     451          u_p(nzb,:,-1) =  u_p(nzb+1,:,-1)
     452          v_p(nzb,:,-1) =  v_p(nzb+1,:,-1)
     453       ENDIF
     454       w_p(nzb,:,-1) = 0.0
     455
     456!
     457!--    Top boundary at the outflow
     458       IF ( ibc_uv_t == 0 )  THEN
     459          u_p(nzt+1,:,-1) = ug(nzt+1)
     460          v_p(nzt+1,:,-1) = vg(nzt+1)
     461       ELSE
     462          u_p(nzt+1,:,-1) = u_p(nzt,:,-1)
     463          v_p(nzt+1,:,-1) = v_p(nzt,:,-1)
     464       ENDIF
     465       w_p(nzt:nzt+1,:,-1) = 0.0
     466
     467    ENDIF
     468
     469    IF ( outflow_r  .AND.                                                 &
     470         intermediate_timestep_count == intermediate_timestep_count_max ) &
     471    THEN
     472
     473       c_max = dx / dt_3d
     474
     475       DO j = nys-1, nyn+1
     476          DO k = nzb+1, nzt+1
     477
     478!
     479!--          First calculate the phase speeds for u,v, and w
     480             denom = u_m_r(k,j,nx) - u_m_r(k,j,nx-1)
     481
     482             IF ( denom /= 0.0 )  THEN
     483                c_u = -c_max * ( u(k,j,nx) - u_m_r(k,j,nx) ) / denom
     484                IF ( c_u < 0.0 )  THEN
     485                   c_u = 0.0
     486                ELSEIF ( c_u > c_max )  THEN
     487                   c_u = c_max
     488                ENDIF
     489             ELSE
     490                c_u = c_max
     491             ENDIF
     492
     493             denom = v_m_r(k,j,nx) - v_m_r(k,j,nx-1)
     494
     495             IF ( denom /= 0.0 )  THEN
     496                c_v = -c_max * ( v(k,j,nx) - v_m_r(k,j,nx) ) / denom
     497                IF ( c_v < 0.0 )  THEN
     498                   c_v = 0.0
     499                ELSEIF ( c_v > c_max )  THEN
     500                   c_v = c_max
     501                ENDIF
     502             ELSE
     503                c_v = c_max
     504             ENDIF
     505
     506             denom = w_m_r(k,j,nx) - w_m_r(k,j,nx-1)
     507
     508             IF ( denom /= 0.0 )  THEN
     509                c_w = -c_max * ( w(k,j,nx) - w_m_r(k,j,nx) ) / denom
     510                IF ( c_w < 0.0 )  THEN
     511                   c_w = 0.0
     512                ELSEIF ( c_w > c_max )  THEN
     513                   c_w = c_max
     514                ENDIF
     515             ELSE
     516                c_w = c_max
     517             ENDIF
     518
     519!
     520!--          Calculate the new velocities
     521             u_p(k,j,nx+1) = u(k,j,nx+1) - dt_3d * c_u * &
     522                                           ( u(k,j,nx+1) - u(k,j,nx) ) * ddx
     523
     524             v_p(k,j,nx+1) = v(k,j,nx+1) - dt_3d * c_v * &
     525                                           ( v(k,j,nx+1) - v(k,j,nx) ) * ddx
     526
     527             w_p(k,j,nx+1) = w(k,j,nx+1) - dt_3d * c_w * &
     528                                           ( w(k,j,nx+1) - w(k,j,nx) ) * ddx
     529
     530!
     531!--          Swap timelevels for the next timestep
     532             u_m_r(k,j,:) = u(k,j,nx-1:nx+1)
     533             v_m_r(k,j,:) = v(k,j,nx-1:nx+1)
     534             w_m_r(k,j,:) = w(k,j,nx-1:nx+1)
     535
     536          ENDDO
     537       ENDDO
     538
     539!
     540!--    Bottom boundary at the outflow
     541       IF ( ibc_uv_b == 0 )  THEN
     542          u_p(nzb,:,nx+1) = -u_p(nzb+1,:,nx+1)
     543          v_p(nzb,:,nx+1) = -v_p(nzb+1,:,nx+1) 
     544       ELSE                   
     545          u_p(nzb,:,nx+1) =  u_p(nzb+1,:,nx+1)
     546          v_p(nzb,:,nx+1) =  v_p(nzb+1,:,nx+1)
     547       ENDIF
     548       w_p(nzb,:,nx+1) = 0.0
     549
     550!
     551!--    Top boundary at the outflow
     552       IF ( ibc_uv_t == 0 )  THEN
     553          u_p(nzt+1,:,nx+1) = ug(nzt+1)
     554          v_p(nzt+1,:,nx+1) = vg(nzt+1)
     555       ELSE
     556          u_p(nzt+1,:,nx+1) = u_p(nzt,:,nx+1)
     557          v_p(nzt+1,:,nx+1) = v_p(nzt,:,nx+1)
     558       ENDIF
     559       w(nzt:nzt+1,:,nx+1) = 0.0
    418560
    419561    ENDIF
  • palm/trunk/SOURCE/check_parameters.f90

    r73 r75  
    66! "by_user" allowed as initializing action, -data_output_ts,
    77! leapfrog with non-flat topography not allowed any more, loop_optimization
    8 ! and pt_reference are checked,
     8! and pt_reference are checked, moisture renamed humidity,
    99! output of precipitation amount/rate and roughnes length + check
    1010! possible negative humidities are avoided in initial profile,
    11 ! dirichlet/neumann changed to dirichlet/radiation, etc.
     11! dirichlet/neumann changed to dirichlet/radiation, etc.,
     12! revision added to run_description_header
    1213!
    1314! Former revisions:
     
    8081    run_time = time(1:2)//':'//time(3:4)//':'//time(5:6)
    8182
    82     WRITE ( run_description_header, '(A,2X,A,A,A,I2.2,2X,A,A,2X,A,1X,A)' )  &
    83               TRIM( version ),'run: ', TRIM( run_identifier ), '.', &
     83    WRITE ( run_description_header, '(A,2X,A,2X,A,A,A,I2.2,2X,A,A,2X,A,1X,A)' )&
     84              TRIM( version ), TRIM( revision ), &
     85              'run: ', TRIM( run_identifier ), '.', &
    8486              runnr, 'host: ', TRIM( host ), run_date, run_time
    8587
     
    131133          WRITE( action, '(A)' )  'cloud_droplets = .TRUE.'
    132134       ENDIF
    133        IF ( moisture )  THEN
    134           WRITE( action, '(A)' )  'moisture = .TRUE.'
     135       IF ( humidity )  THEN
     136          WRITE( action, '(A)' )  'humidity = .TRUE.'
    135137       ENDIF
    136138       IF ( .NOT. prandtl_layer )  THEN
     
    353355    ENDIF
    354356
    355     IF ( cloud_physics  .AND.  .NOT. moisture )  THEN
     357    IF ( cloud_physics  .AND.  .NOT. humidity )  THEN
    356358       IF ( myid == 0 )  PRINT*, '+++ check_parameters: cloud_physics =', &
    357359                                 cloud_physics, ' is not allowed with ',  &
    358                                  'moisture =', moisture
     360                                 'humidity =', humidity
    359361       CALL local_stop
    360362    ENDIF
     
    367369    ENDIF
    368370
    369     IF ( moisture  .AND.  sloping_surface )  THEN
    370        IF ( myid == 0 )  PRINT*, '+++ check_parameters: moisture = TRUE', &
     371    IF ( humidity  .AND.  sloping_surface )  THEN
     372       IF ( myid == 0 )  PRINT*, '+++ check_parameters: humidity = TRUE', &
    371373                                 'and hang = TRUE are not',               &
    372374                                 ' allowed simultaneously'
     
    374376    ENDIF
    375377
    376     IF ( moisture  .AND.  scalar_advec == 'ups-scheme' )  THEN
     378    IF ( humidity  .AND.  scalar_advec == 'ups-scheme' )  THEN
    377379       IF ( myid == 0 )  PRINT*, '+++ check_parameters: UPS-scheme', &
    378                                  'is not implemented for moisture'
     380                                 'is not implemented for humidity'
    379381       CALL local_stop       
    380382    ENDIF
    381383
    382     IF ( passive_scalar  .AND.  moisture )  THEN
    383        IF ( myid == 0 )  PRINT*, '+++ check_parameters: moisture = TRUE and', &
     384    IF ( passive_scalar  .AND.  humidity )  THEN
     385       IF ( myid == 0 )  PRINT*, '+++ check_parameters: humidity = TRUE and', &
    384386                                 'passive_scalar = TRUE is not allowed ',     &
    385387                                 'simultaneously'
     
    410412       v_init  = vg_surface
    411413       pt_init = pt_surface
    412        IF ( moisture )        q_init = q_surface
     414       IF ( humidity )        q_init = q_surface
    413415       IF ( passive_scalar )  q_init = s_surface
    414416
     
    552554       ENDIF
    553555
    554        IF ( moisture  .OR.  passive_scalar )  THEN
     556       IF ( humidity  .OR.  passive_scalar )  THEN
    555557
    556558          i = 1
     
    760762          CALL local_stop
    761763       ENDIF
    762        IF ( conserve_volume_flow )  THEN
    763           IF ( myid == 0 )  THEN
    764              PRINT*, '+++ check_parameters:'
    765              PRINT*, '    non-cyclic lateral boundaries do not allow', &
    766                           ' conserve_volume_flow = .T.'
    767           ENDIF
    768           CALL local_stop
    769        ENDIF
     764!       IF ( conserve_volume_flow )  THEN
     765!          IF ( myid == 0 )  THEN
     766!             PRINT*, '+++ check_parameters:'
     767!             PRINT*, '    non-cyclic lateral boundaries do not allow', &
     768!                          ' conserve_volume_flow = .T.'
     769!          ENDIF
     770!          CALL local_stop
     771!       ENDIF
    770772    ENDIF
    771773
     
    907909
    908910!
    909 !-- In case of moisture or passive scalar, set boundary conditions for total
     911!-- In case of humidity or passive scalar, set boundary conditions for total
    910912!-- water content / scalar
    911     IF ( moisture  .OR.  passive_scalar ) THEN
    912        IF ( moisture )  THEN
     913    IF ( humidity  .OR.  passive_scalar ) THEN
     914       IF ( humidity )  THEN
    913915          sq = 'q'
    914916       ELSE
     
    942944!
    943945!--    A given surface humidity implies Dirichlet boundary condition for
    944 !--    moisture. In this case specification of a constant water flux is
     946!--    humidity. In this case specification of a constant water flux is
    945947!--    forbidden.
    946948       IF ( ibc_q_b == 0  .AND.  constant_waterflux )  THEN
     
    15521554
    15531555          CASE ( 'w"qv"' )
    1554              IF ( moisture  .AND.  .NOT. cloud_physics ) &
     1556             IF ( humidity  .AND.  .NOT. cloud_physics ) &
    15551557             THEN
    15561558                dopr_index(i) = 48
    15571559                hom(:,2,48,:) = SPREAD( zw, 2, statistic_regions+1 )
    1558              ELSEIF( moisture .AND. cloud_physics ) THEN
     1560             ELSEIF( humidity .AND. cloud_physics ) THEN
    15591561                dopr_index(i) = 51
    15601562                hom(:,2,51,:) = SPREAD( zw, 2, statistic_regions+1 )
     
    15641566                           data_output_pr(i),                          &
    15651567                           '    is not implemented for cloud_physics = FALSE', &
    1566                            '    and                    moisture      = FALSE'
     1568                           '    and                    humidity      = FALSE'
    15671569                ENDIF
    15681570                CALL local_stop                   
     
    15701572
    15711573          CASE ( 'w*qv*' )
    1572              IF ( moisture  .AND.  .NOT. cloud_physics ) &
     1574             IF ( humidity  .AND.  .NOT. cloud_physics ) &
    15731575             THEN
    15741576                dopr_index(i) = 49
    15751577                hom(:,2,49,:) = SPREAD( zw, 2, statistic_regions+1 )
    1576              ELSEIF( moisture .AND. cloud_physics ) THEN
     1578             ELSEIF( humidity .AND. cloud_physics ) THEN
    15771579                dopr_index(i) = 52
    15781580                hom(:,2,52,:) = SPREAD( zw, 2, statistic_regions+1 )
     
    15821584                           data_output_pr(i),                                  &
    15831585                           '    is not implemented for cloud_physics = FALSE', &
    1584                            '                       and moisture      = FALSE'
     1586                           '                       and humidity      = FALSE'
    15851587                ENDIF
    15861588                CALL local_stop                   
     
    15881590
    15891591          CASE ( 'wqv' )
    1590              IF ( moisture  .AND.  .NOT. cloud_physics ) &
     1592             IF ( humidity  .AND.  .NOT. cloud_physics ) &
    15911593             THEN
    15921594                dopr_index(i) = 50
    15931595                hom(:,2,50,:) = SPREAD( zw, 2, statistic_regions+1 )
    1594              ELSEIF( moisture .AND. cloud_physics ) THEN
     1596             ELSEIF( humidity .AND. cloud_physics ) THEN
    15951597                dopr_index(i) = 53
    15961598                hom(:,2,53,:) = SPREAD( zw, 2, statistic_regions+1 )
     
    16001602                           data_output_pr(i),                                  &
    16011603                           '    is not implemented for cloud_physics = FALSE', &
    1602                            '                       and moisture      = FALSE'
     1604                           '                       and humidity      = FALSE'
    16031605                ENDIF
    16041606                CALL local_stop                   
     
    18191821
    18201822          CASE ( 'q', 'vpt' )
    1821              IF ( .NOT. moisture )  THEN
     1823             IF ( .NOT. humidity )  THEN
    18221824                IF ( myid == 0 )  THEN
    18231825                   PRINT*, '+++ check_parameters: output of "', TRIM( var ), &
    1824                                 '" requires moisture = .TRUE.'
     1826                                '" requires humidity = .TRUE.'
    18251827                ENDIF
    18261828                CALL local_stop
  • palm/trunk/SOURCE/coriolis.f90

    r4 r75  
    44! Actual revisions:
    55! -----------------
    6 !
     6! uxrp, vynp eliminated
    77!
    88! Former revisions:
     
    5656!--       u-component
    5757          CASE ( 1 )
    58              DO  i = nxl, nxr+uxrp
     58             DO  i = nxl, nxr
    5959                DO  j = nys, nyn
    6060                   DO  k = nzb_u_inner(j,i)+1, nzt
     
    7474          CASE ( 2 )
    7575             DO  i = nxl, nxr
    76                 DO  j = nys, nyn+vynp
     76                DO  j = nys, nyn
    7777                   DO  k = nzb_v_inner(j,i)+1, nzt
    7878                      tend(k,j,i) = tend(k,j,i) - f *     ( 0.25 *          &
  • palm/trunk/SOURCE/data_output_2d.f90

    r72 r75  
    44! Actual revisions:
    55! -----------------
    6 ! Output of precipitation amount/rate and roughness length
     6! Output of precipitation amount/rate and roughness length,
     7! 2nd+3rd argument removed from exchange horiz
    78!
    89! Former revisions:
     
    208209                IF ( av == 0 )  THEN
    209210                   tend = prt_count
    210                    CALL exchange_horiz( tend, 0, 0 )
     211                   CALL exchange_horiz( tend )
    211212                   DO  i = nxl-1, nxr+1
    212213                      DO  j = nys-1, nyn+1
     
    218219                   resorted = .TRUE.
    219220                ELSE
    220                    CALL exchange_horiz( pc_av, 0, 0 )
     221                   CALL exchange_horiz( pc_av )
    221222                   to_be_resorted => pc_av
    222223                ENDIF
     
    243244                      ENDDO
    244245                   ENDDO
    245                    CALL exchange_horiz( tend, 0, 0 )
     246                   CALL exchange_horiz( tend )
    246247                   DO  i = nxl-1, nxr+1
    247248                      DO  j = nys-1, nyn+1
     
    253254                   resorted = .TRUE.
    254255                ELSE
    255                    CALL exchange_horiz( pr_av, 0, 0 )
     256                   CALL exchange_horiz( pr_av )
    256257                   to_be_resorted => pr_av
    257258                ENDIF
  • palm/trunk/SOURCE/data_output_3d.f90

    r4 r75  
    44! Actual revisions:
    55! -----------------
    6 !
     6! 2nd+3rd argument removed from exchange horiz
    77!
    88! Former revisions:
     
    129129             IF ( av == 0 )  THEN
    130130                tend = prt_count
    131                 CALL exchange_horiz( tend, 0, 0 )
     131                CALL exchange_horiz( tend )
    132132                DO  i = nxl-1, nxr+1
    133133                   DO  j = nys-1, nyn+1
     
    139139                resorted = .TRUE.
    140140             ELSE
    141                 CALL exchange_horiz( pc_av, 0, 0 )
     141                CALL exchange_horiz( pc_av )
    142142                to_be_resorted => pc_av
    143143             ENDIF
     
    164164                   ENDDO
    165165                ENDDO
    166                 CALL exchange_horiz( tend, 0, 0 )
     166                CALL exchange_horiz( tend )
    167167                DO  i = nxl-1, nxr+1
    168168                   DO  j = nys-1, nyn+1
     
    174174                resorted = .TRUE.
    175175             ELSE
    176                 CALL exchange_horiz( pr_av, 0, 0 )
     176                CALL exchange_horiz( pr_av )
    177177                to_be_resorted => pr_av
    178178             ENDIF
  • palm/trunk/SOURCE/data_output_dvrp.f90

    r60 r75  
    3232! Actual revisions:
    3333! -----------------
    34 ! Particles-package is now part of the default code.
     34! Particles-package is now part of the default code,
     35! moisture renamed humidity
    3536! TEST: write statements
    3637!
     
    391392
    392393             CASE ( 'q', 'q_xy', 'q_xz', 'q_yz' )
    393                 IF ( moisture  .OR.  passive_scalar )  THEN
     394                IF ( humidity  .OR.  passive_scalar )  THEN
    394395                   DO  i = nxl, nxr+1
    395396                      DO  j = nys, nyn+1
     
    401402                ELSE
    402403                   IF ( myid == 0 )  THEN
    403                       PRINT*, '+++ data_output_dvrp: if moisture/passive_scalar = ', &
     404                      PRINT*, '+++ data_output_dvrp: if humidity/passive_scalar = ', &
    404405                              'FALSE output of ', output_variable,            &
    405406                              'is not provided'
  • palm/trunk/SOURCE/diffusion_u.f90

    r57 r75  
    55! -----------------
    66! Wall functions now include diabatic conditions, call of routine wall_fluxes,
    7 ! z0 removed from argument list
     7! z0 removed from argument list, uxrp eliminated
    88!
    99! Former revisions:
     
    6565       REAL, DIMENSION(:,:),   POINTER ::  usws
    6666       REAL, DIMENSION(:,:,:), POINTER ::  km, u, v, w
    67        REAL, DIMENSION(nzb:nzt+1,nys:nyn,nxl:nxr+uxrp) ::  usvs
     67       REAL, DIMENSION(nzb:nzt+1,nys:nyn,nxl:nxr) ::  usvs
    6868
    6969!
     
    7171!--    if neccessary
    7272       IF ( topography /= 'flat' )  THEN
    73           CALL wall_fluxes( usvs, 1.0, 0.0, 0.0, 0.0, uxrp, 0, nzb_u_inner, &
     73          CALL wall_fluxes( usvs, 1.0, 0.0, 0.0, 0.0, nzb_u_inner, &
    7474                            nzb_u_outer, wall_u )
    7575       ENDIF
    7676
    77        DO  i = nxl, nxr+uxrp
     77       DO  i = nxl, nxr
    7878          DO  j = nys,nyn
    7979!
  • palm/trunk/SOURCE/diffusion_v.f90

    r57 r75  
    55! -----------------
    66! Wall functions now include diabatic conditions, call of routine wall_fluxes,
    7 ! z0 removed from argument list
     7! z0 removed from argument list, vynp eliminated
    88!
    99! Former revisions:
     
    6363       REAL, DIMENSION(:,:),   POINTER ::  vsws
    6464       REAL, DIMENSION(:,:,:), POINTER ::  km, u, v, w
    65        REAL, DIMENSION(nzb:nzt+1,nys:nyn+vynp,nxl:nxr) ::  vsus
     65       REAL, DIMENSION(nzb:nzt+1,nys:nyn,nxl:nxr) ::  vsus
    6666
    6767!
     
    6969!--    if neccessary
    7070       IF ( topography /= 'flat' )  THEN
    71           CALL wall_fluxes( vsus, 0.0, 1.0, 0.0, 0.0, 0, vynp, nzb_v_inner, &
     71          CALL wall_fluxes( vsus, 0.0, 1.0, 0.0, 0.0, nzb_v_inner, &
    7272                            nzb_v_outer, wall_v )
    7373       ENDIF
    7474
    7575       DO  i = nxl, nxr
    76           DO  j = nys, nyn+vynp
     76          DO  j = nys, nyn
    7777!
    7878!--          Compute horizontal diffusion
  • palm/trunk/SOURCE/diffusion_w.f90

    r57 r75  
    6969!--    walls, if neccessary
    7070       IF ( topography /= 'flat' )  THEN
    71           CALL wall_fluxes( wsus, 0.0, 0.0, 0.0, 1.0, 0, 0, nzb_w_inner, &
     71          CALL wall_fluxes( wsus, 0.0, 0.0, 0.0, 1.0, nzb_w_inner, &
    7272                            nzb_w_outer, wall_w_x )
    73           CALL wall_fluxes( wsvs, 0.0, 0.0, 1.0, 0.0, 0, 0, nzb_w_inner, &
     73          CALL wall_fluxes( wsvs, 0.0, 0.0, 1.0, 0.0, nzb_w_inner, &
    7474                            nzb_w_outer, wall_w_y )
    7575       ENDIF
  • palm/trunk/SOURCE/disturb_field.f90

    r4 r75  
    1  SUBROUTINE disturb_field( nzb_uv_inner, dist1, field, xrp, ynp )
     1 SUBROUTINE disturb_field( nzb_uv_inner, dist1, field )
    22
    33!------------------------------------------------------------------------------!
    44! Actual revisions:
    55! -----------------
    6 !
     6! xrp, ynp eliminated, 2nd+3rd argument removed from exchange horiz
    77!
    88! Former revisions:
     
    3636    IMPLICIT NONE
    3737
    38     INTEGER ::  i, j, k, xrp, ynp
     38    INTEGER ::  i, j, k
    3939    INTEGER ::  nzb_uv_inner(nys-1:nyn+1,nxl-1:nxr+1)
    4040
    4141    REAL    ::  randomnumber,                             &
    4242                dist1(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1), &
    43                 field(nzb:nzt+1,nys-1:nyn+ynp+1,nxl-1:nxr+xrp+1)
     43                field(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1)
    4444    REAL, DIMENSION(:,:,:), ALLOCATABLE ::  dist2
    4545
     
    9494!
    9595!-- Exchange of ghost points for the random perturbation
    96     CALL exchange_horiz( dist1, 0, 0 )
     96    CALL exchange_horiz( dist1 )
    9797
    9898!
     
    114114!-- Exchange of ghost points for the filtered perturbation.
    115115!-- Afterwards, filter operation and exchange of ghost points are repeated.
    116     CALL exchange_horiz( dist2, 0, 0 )
     116    CALL exchange_horiz( dist2 )
    117117    DO  i = nxl, nxr
    118118       DO  j = nys, nyn
     
    125125       ENDDO
    126126    ENDDO
    127     CALL exchange_horiz( dist1, 0, 0 )
     127    CALL exchange_horiz( dist1 )
    128128
    129129!
  • palm/trunk/SOURCE/exchange_horiz.f90

    r4 r75  
    1  SUBROUTINE exchange_horiz( ar, xrp, ynp )
     1 SUBROUTINE exchange_horiz( ar )
    22
    33!------------------------------------------------------------------------------!
    44! Actual revisions:
    55! -----------------
    6 !
     6! Special cases for additional gridpoints along x or y in case of non-cyclic
     7! boundary conditions are not regarded any more
    78!
    89! Former revisions:
     
    3233    IMPLICIT NONE
    3334
    34     INTEGER ::  xrp, ynp
    35 
    3635#if defined( __parallel )
    37     INTEGER                               ::  typexz
    3836    INTEGER, DIMENSION(4)                 ::  req
    3937    INTEGER, DIMENSION(MPI_STATUS_SIZE,4) ::  wait_stat
    4038#endif
    4139
    42     REAL ::  ar(nzb:nzt+1,nys-1:nyn+ynp+1,nxl-1:nxr+xrp+1)
     40    REAL ::  ar(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1)
    4341
    4442
     
    5957
    6058    ELSE
     59
    6160       req = 0
    6261!
     
    7675               ar(nzb,nys-1,nxl-1), ngp_yz(grid_level), MPI_REAL, pleft,  1, &
    7776                          comm2d, req(4), ierr )
    78        call MPI_Waitall (4,req,wait_stat,ierr)
     77       CALL MPI_WAITALL( 4, req, wait_stat, ierr )
     78
    7979    ENDIF
    8080
     
    9090
    9191    ELSE
    92 !
    93 !--    Set the MPI data type, which depends on the size of the array
    94 !--    (the v array has an additional gridpoint along y in case of non-cyclic
    95 !--    boundary conditions)
    96        IF ( ynp == 0 )  THEN
    97           typexz = type_xz(grid_level)
    98        ELSE
    99           typexz = type_xz_p
    100        ENDIF
    10192
    10293       req = 0
    10394!
    10495!--    Send front boundary, receive rear one
    105        CALL MPI_ISEND( ar(nzb,nys,nxl-1),   1, typexz, psouth, 0, comm2d, &
    106                        req(1), ierr )
    107        CALL MPI_IRECV( ar(nzb,nyn+1,nxl-1), 1, typexz, pnorth, 0, comm2d, &
    108                        req(2), ierr )
     96       CALL MPI_ISEND( ar(nzb,nys,nxl-1),   1, type_xz(grid_level), psouth, 0, &
     97                       comm2d, req(1), ierr )
     98       CALL MPI_IRECV( ar(nzb,nyn+1,nxl-1), 1, type_xz(grid_level), pnorth, 0, &
     99                       comm2d, req(2), ierr )
    109100!
    110101!--    Send rear boundary, receive front one
    111        CALL MPI_ISEND( ar(nzb,nyn,nxl-1),   1, typexz, pnorth, 1, comm2d, &
    112                        req(3), ierr )
    113        CALL MPI_IRECV( ar(nzb,nys-1,nxl-1), 1, typexz, psouth, 1, comm2d, &
    114                        req(4), ierr )
    115        call MPI_Waitall (4,req,wait_stat,ierr)
     102       CALL MPI_ISEND( ar(nzb,nyn,nxl-1),   1, type_xz(grid_level), pnorth, 1, &
     103                       comm2d, req(3), ierr )
     104       CALL MPI_IRECV( ar(nzb,nys-1,nxl-1), 1, type_xz(grid_level), psouth, 1, &
     105                       comm2d, req(4), ierr )
     106       call MPI_WAITALL( 4, req, wait_stat, ierr )
     107
    116108    ENDIF
    117109
  • palm/trunk/SOURCE/flow_statistics.f90

    r51 r75  
    55! -----------------
    66! Collection of time series quantities moved from routine flow_statistics to
    7 ! here, routine user_statistics is called for each statistic region
     7! here, routine user_statistics is called for each statistic region,
     8! moisture renamed humidity
    89!
    910! Former revisions:
     
    115116!--    total water content, specific humidity and liquid water potential
    116117!--    temperature
    117        IF ( moisture )  THEN
     118       IF ( humidity )  THEN
    118119          !$OMP DO
    119120          DO  i = nxl, nxr
     
    164165             sums_l(:,2,0) = sums_l(:,2,0) + sums_l(:,2,i)
    165166             sums_l(:,4,0) = sums_l(:,4,0) + sums_l(:,4,i)
    166              IF ( moisture )  THEN
     167             IF ( humidity )  THEN
    167168                sums_l(:,41,0) = sums_l(:,41,0) + sums_l(:,41,i)
    168169                sums_l(:,44,0) = sums_l(:,44,0) + sums_l(:,44,i)
     
    187188       CALL MPI_ALLREDUCE( sums_l(nzb,4,0), sums(nzb,4), nzt+2-nzb, MPI_REAL, &
    188189                           MPI_SUM, comm2d, ierr )
    189        IF ( moisture ) THEN
     190       IF ( humidity ) THEN
    190191          CALL MPI_ALLREDUCE( sums_l(nzb,44,0), sums(nzb,44), nzt+2-nzb, &
    191192                              MPI_REAL, MPI_SUM, comm2d, ierr )
     
    208209       sums(:,2) = sums_l(:,2,0)
    209210       sums(:,4) = sums_l(:,4,0)
    210        IF ( moisture ) THEN
     211       IF ( humidity ) THEN
    211212          sums(:,44) = sums_l(:,44,0)
    212213          sums(:,41) = sums_l(:,41,0)
     
    231232!
    232233!--    Humidity and cloud parameters
    233        IF ( moisture ) THEN
     234       IF ( humidity ) THEN
    234235          sums(:,44) = sums(:,44) / ngp_2dh_outer(:,sr)
    235236          sums(:,41) = sums(:,41) / ngp_2dh_outer(:,sr)
     
    365366!
    366367!--             Buoyancy flux, water flux (humidity flux) w"q"
    367                 IF ( moisture ) THEN
     368                IF ( humidity ) THEN
    368369                   sums_l(k,45,tn) = sums_l(k,45,tn)                           &
    369370                                         - 0.5 * ( kh(k,j,i) + kh(k+1,j,i) )   &
     
    407408                sums_l(nzb,61,tn) = sums_l(nzb,61,tn) + &
    408409                                    0.0 * rmask(j,i,sr)           ! v"pt"
    409                 IF ( moisture )  THEN
     410                IF ( humidity )  THEN
    410411                   sums_l(nzb,48,tn) = sums_l(nzb,48,tn) + &
    411412                                       qsws(j,i) * rmask(j,i,sr)  ! w"q" (w"qv")
     
    437438                sums_l(nzt,61,tn) = sums_l(nzt,61,tn) + &
    438439                                    0.0 * rmask(j,i,sr)           ! v"pt"
    439                 IF ( moisture )  THEN
     440                IF ( humidity )  THEN
    440441                   sums_l(nzt,48,tn) = sums_l(nzt,48,tn) + &
    441442                                       qswst(j,i) * rmask(j,i,sr)  ! w"q" (w"qv")
     
    499500!--             Buoyancy flux, water flux, humidity flux and liquid water
    500501!--             content
    501                 IF ( moisture )  THEN
     502                IF ( humidity )  THEN
    502503                   pts = 0.5 * ( vpt(k,j,i)   - hom(k,1,44,sr) + &
    503504                                 vpt(k+1,j,i) - hom(k+1,1,44,sr) )
  • palm/trunk/SOURCE/header.f90

    r63 r75  
    55! -----------------
    66! Output of netcdf_64bit_3d, particles-package is now part of the default code,
    7 ! output of the loop optimization method.
     7! output of the loop optimization method, moisture renamed humidity,
     8! output of subversion revision number
    89!
    910! Former revisions:
     
    5253    CHARACTER (LEN=10) ::  coor_chr, host_chr
    5354    CHARACTER (LEN=16) ::  begin_chr
     55    CHARACTER (LEN=21) ::  ver_rev
    5456    CHARACTER (LEN=40) ::  output_format
    5557    CHARACTER (LEN=70) ::  char1, char2, coordinates, gradients, dopr_chr, &
     
    9597!-- Run-identification, date, time, host
    9698    host_chr = host(1:10)
    97     WRITE ( io, 100 )  version, TRIM( run_classification ), run_date, &
     99    ver_rev = TRIM( version ) // '  ' // TRIM( revision )
     100    WRITE ( io, 100 )  ver_rev, TRIM( run_classification ), run_date, &
    98101                       run_identifier, run_time, runnr, ADJUSTR( host_chr )
    99102#if defined( __parallel )
     
    217220       WRITE ( io, 123 )  rayleigh_damping_height, rayleigh_damping_factor
    218221    ENDIF
    219     IF ( moisture )  THEN
     222    IF ( humidity )  THEN
    220223       IF ( .NOT. cloud_physics )  THEN
    221224          WRITE ( io, 129 )
     
    370373    ENDIF
    371374
    372     IF ( moisture  .OR.  passive_scalar )  THEN
    373        IF ( moisture )  THEN
     375    IF ( humidity  .OR.  passive_scalar )  THEN
     376       IF ( humidity )  THEN
    374377          IF ( ibc_q_b == 0 )  THEN
    375378             runten = 'q(0)     = q_surface'
     
    405408          IF ( random_heatflux )  WRITE ( io, 307 )
    406409       ENDIF
    407        IF ( moisture  .AND.  constant_waterflux )  THEN
     410       IF ( humidity  .AND.  constant_waterflux )  THEN
    408411          WRITE ( io, 311 ) surface_waterflux
    409412       ENDIF
     
    418421          WRITE ( io, 306 )  top_heatflux
    419422       ENDIF
    420        IF ( moisture  .OR.  passive_scalar )  THEN
     423       IF ( humidity  .OR.  passive_scalar )  THEN
    421424          WRITE ( io, 315 )
    422425       ENDIF
     
    426429       WRITE ( io, 305 )  zu(1), roughness_length, kappa, rif_min, rif_max
    427430       IF ( .NOT. constant_heatflux )  WRITE ( io, 308 )
    428        IF ( moisture  .AND.  .NOT. constant_waterflux )  THEN
     431       IF ( humidity  .AND.  .NOT. constant_waterflux )  THEN
    429432          WRITE ( io, 312 )
    430433       ENDIF
     
    955958!-- Initial humidity profile
    956959!-- Building output strings, starting with surface humidity
    957     IF ( moisture  .OR.  passive_scalar )  THEN
     960    IF ( humidity  .OR.  passive_scalar )  THEN
    958961       WRITE ( temperatures, '(E8.1)' )  q_surface
    959962       gradients = '--------'
     
    978981       ENDDO
    979982
    980        IF ( moisture )  THEN
     983       IF ( humidity )  THEN
    981984          WRITE ( io, 421 )  TRIM( coordinates ), TRIM( temperatures ), &
    982985                             TRIM( gradients ), TRIM( slices )
     
    10211024       WRITE ( io, 475 )  pt_surface_initial_change
    10221025    ENDIF
    1023     IF ( moisture  .AND.  q_surface_initial_change /= 0.0 )  THEN
     1026    IF ( humidity  .AND.  q_surface_initial_change /= 0.0 )  THEN
    10241027       WRITE ( io, 476 )  q_surface_initial_change       
    10251028    ENDIF
     
    11131116
    11141117 99 FORMAT (1X,78('-'))
    1115 100 FORMAT (/10X,'****************',11X,28('-')/                &
    1116             10X,'*  ',A12,'*',11X,A/                            &
    1117             10X,'****************',11X,28('-')//                &
     1118100 FORMAT (/1X,'*************************',11X,28('-')/        &
     1119            1X,'* ',A,' *',11X,A/                               &
     1120            1X,'*************************',11X,28('-')//        &
    11181121            ' Date:            ',A8,11X,'Run:       ',A20/      &
    11191122            ' Time:            ',A8,11X,'Run-No.:   ',I2.2/     &
  • palm/trunk/SOURCE/init_1d_model.f90

    r46 r75  
    44! Actual revisions:
    55! -----------------
    6 ! Bugfix: preset of tendencies te_em, te_um, te_vm
     6! Bugfix: preset of tendencies te_em, te_um, te_vm,
     7! moisture renamed humidity
    78!
    89! Former revisions:
     
    125126    vsws1d = 0.0; vsws1d_m = 0.0
    126127    z01d = roughness_length
    127     IF ( moisture .OR. passive_scalar )  qs1d = 0.0
     128    IF ( humidity .OR. passive_scalar )  qs1d = 0.0
    128129
    129130!
     
    216217                kmzm = 0.5 * ( km1d_m(k-1) + km1d_m(k) )
    217218                kmzp = 0.5 * ( km1d_m(k) + km1d_m(k+1) )
    218                 IF ( .NOT. moisture )  THEN
     219                IF ( .NOT. humidity )  THEN
    219220                   pt_0 = pt_init(k)
    220221                   flux =  ( pt_init(k+1)-pt_init(k-1) ) * dd2zu(k)
     
    256257             kmzm = 0.5 * ( km1d_m(k-1) + km1d_m(k) )
    257258             kmzp = 0.5 * ( km1d_m(k) + km1d_m(k+1) )
    258              IF ( .NOT. moisture )  THEN
     259             IF ( .NOT. humidity )  THEN
    259260                pt_0 = pt_init(k)
    260261                flux =  ( pt_init(k+1)-pt_init(k-1) ) * dd2zu(k)
     
    450451
    451452             IF ( prandtl_layer )  THEN
    452                 IF ( .NOT. moisture )  THEN
     453                IF ( .NOT. humidity )  THEN
    453454                   pt_0 = pt_init(nzb+1)
    454455                   flux = ts1d
     
    462463
    463464             DO  k = nzb_diff, nzt
    464                 IF ( .NOT. moisture )  THEN
     465                IF ( .NOT. humidity )  THEN
    465466                   pt_0 = pt_init(k)
    466467                   flux = ( pt_init(k+1) - pt_init(k-1) ) * dd2zu(k)
     
    544545                e1d(nzb) = e1d(nzb+1)
    545546
    546                 IF ( moisture .OR. passive_scalar ) THEN
     547                IF ( humidity .OR. passive_scalar ) THEN
    547548!
    548549!--                Compute q*
  • palm/trunk/SOURCE/init_3d_model.f90

    r73 r75  
    99! Arrays for radiation boundary conditions are allocated (u_m_l, u_m_r, etc.),
    1010! bugfix for cases with the outflow damping layer extending over more than one
    11 ! subdomain,
    12 ! New initializing action "by_user" calls user_init_3d_model,
     11! subdomain, moisture renamed humidity,
     12! new initializing action "by_user" calls user_init_3d_model,
    1313! precipitation_amount/rate, ts_value are allocated, +module netcdf_control,
    1414! initial velocities at nzb+1 are regarded for volume
    1515! flow control in case they have been set zero before (to avoid small timesteps)
     16! -uvmean_outflow, uxrp, vynp eliminated
    1617!
    1718! Former revisions:
     
    7677              sums_divnew_l(0:statistic_regions),                           &
    7778              sums_divold_l(0:statistic_regions) )
    78     ALLOCATE( rdf(nzb+1:nzt), uvmean_outflow(nzb:nzt+1),                    &
    79               uvmean_outflow_l(nzb:nzt+1) )
     79    ALLOCATE( rdf(nzb+1:nzt) )
    8080    ALLOCATE( hom_sum(nzb:nzt+1,var_hom,0:statistic_regions),               &
    8181              ngp_2dh_outer(nzb:nzt+1,0:statistic_regions),                 &
     
    105105    ENDIF
    106106
    107     ALLOCATE( d(nzb+1:nzta,nys:nyna,nxl:nxra),             &
    108               e_1(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1),      &
    109               e_2(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1),      &
    110               e_3(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1),      &
    111               kh_1(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1),     &
    112               km_1(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1),     &
    113               p(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1),        &
    114               pt_1(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1),     &
    115               pt_2(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1),     &
    116               pt_3(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1),     &
    117               tend(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1),     &
    118               u_1(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1+uxrp), &
    119               u_2(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1+uxrp), &
    120               u_3(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1+uxrp), &
    121               v_1(nzb:nzt+1,nys-1:nyn+1+vynp,nxl-1:nxr+1), &
    122               v_2(nzb:nzt+1,nys-1:nyn+1+vynp,nxl-1:nxr+1), &
    123               v_3(nzb:nzt+1,nys-1:nyn+1+vynp,nxl-1:nxr+1), &
    124               w_1(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1),      &
    125               w_2(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1),      &
     107    ALLOCATE( d(nzb+1:nzta,nys:nyna,nxl:nxra),         &
     108              e_1(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1),  &
     109              e_2(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1),  &
     110              e_3(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1),  &
     111              kh_1(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1), &
     112              km_1(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1), &
     113              p(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1),    &
     114              pt_1(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1), &
     115              pt_2(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1), &
     116              pt_3(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1), &
     117              tend(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1), &
     118              u_1(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1), &
     119              u_2(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1), &
     120              u_3(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1), &
     121              v_1(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1), &
     122              v_2(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1), &
     123              v_3(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1), &
     124              w_1(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1),  &
     125              w_2(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1),  &
    126126              w_3(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
    127127
     
    131131    ENDIF
    132132
    133     IF ( moisture  .OR.  passive_scalar ) THEN
    134 !
    135 !--    2D-moisture/scalar arrays
     133    IF ( humidity  .OR.  passive_scalar ) THEN
     134!
     135!--    2D-humidity/scalar arrays
    136136       ALLOCATE ( qs(nys-1:nyn+1,nxl-1:nxr+1),     &
    137137                  qsws_1(nys-1:nyn+1,nxl-1:nxr+1), &
     
    143143       ENDIF
    144144!
    145 !--    3D-moisture/scalar arrays
     145!--    3D-humidity/scalar arrays
    146146       ALLOCATE( q_1(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1), &
    147147                 q_2(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1), &
     
    149149
    150150!
    151 !--    3D-arrays needed for moisture only
    152        IF ( moisture )  THEN
     151!--    3D-arrays needed for humidity only
     152       IF ( humidity )  THEN
    153153          ALLOCATE( vpt_1(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
    154154
     
    212212!-- conditions
    213213    IF ( outflow_l )  THEN
    214        ALLOCATE( u_m_l(nzb:nzt+1,nys-1:nyn+1,-1:1,-2:-1), &
    215                  v_m_l(nzb:nzt+1,nys-1:nyn+1,-1:1,-2:-1), &
    216                  w_m_l(nzb:nzt+1,nys-1:nyn+1,-1:1,-2:-1) )
     214       ALLOCATE( u_m_l(nzb:nzt+1,nys-1:nyn+1,-1:1), &
     215                 v_m_l(nzb:nzt+1,nys-1:nyn+1,-1:1), &
     216                 w_m_l(nzb:nzt+1,nys-1:nyn+1,-1:1) )
    217217    ENDIF
    218218    IF ( outflow_r )  THEN
    219        ALLOCATE( u_m_r(nzb:nzt+1,nys-1:nyn+1,nx-1:nx+1,-2:-1), &
    220                  v_m_r(nzb:nzt+1,nys-1:nyn+1,nx-1:nx+1,-2:-1), &
    221                  w_m_r(nzb:nzt+1,nys-1:nyn+1,nx-1:nx+1,-2:-1) )
     219       ALLOCATE( u_m_r(nzb:nzt+1,nys-1:nyn+1,nx-1:nx+1), &
     220                 v_m_r(nzb:nzt+1,nys-1:nyn+1,nx-1:nx+1), &
     221                 w_m_r(nzb:nzt+1,nys-1:nyn+1,nx-1:nx+1) )
    222222    ENDIF
    223223    IF ( outflow_s )  THEN
    224        ALLOCATE( u_m_s(nzb:nzt+1,-1:1,nxl-1:nxr+1,-2:-1), &
    225                  v_m_s(nzb:nzt+1,-1:1,nxl-1:nxr+1,-2:-1), &
    226                  w_m_s(nzb:nzt+1,-1:1,nxl-1:nxr+1,-2:-1) )
     224       ALLOCATE( u_m_s(nzb:nzt+1,-1:1,nxl-1:nxr+1), &
     225                 v_m_s(nzb:nzt+1,-1:1,nxl-1:nxr+1), &
     226                 w_m_s(nzb:nzt+1,-1:1,nxl-1:nxr+1) )
    227227    ENDIF
    228228    IF ( outflow_n )  THEN
    229        ALLOCATE( u_m_n(nzb:nzt+1,ny-1:ny+1,nxl-1:nxr+1,-2:-1), &
    230                  v_m_n(nzb:nzt+1,ny-1:ny+1,nxl-1:nxr+1,-2:-1), &
    231                  w_m_n(nzb:nzt+1,ny-1:ny+1,nxl-1:nxr+1,-2:-1) )
     229       ALLOCATE( u_m_n(nzb:nzt+1,ny-1:ny+1,nxl-1:nxr+1), &
     230                 v_m_n(nzb:nzt+1,ny-1:ny+1,nxl-1:nxr+1), &
     231                 w_m_n(nzb:nzt+1,ny-1:ny+1,nxl-1:nxr+1) )
    232232    ENDIF
    233233
     
    249249       w_m  => w_1;   w  => w_2;   w_p  => w_3;   tw_m  => w_3
    250250
    251        IF ( moisture  .OR.  passive_scalar )  THEN
     251       IF ( humidity  .OR.  passive_scalar )  THEN
    252252          qsws_m  => qsws_1;   qsws  => qsws_2
    253253          qswst_m => qswst_1;  qswst => qswst_2
    254254          q_m    => q_1;     q    => q_2;     q_p => q_3;     tq_m => q_3
    255           IF ( moisture )        vpt_m  => vpt_1;   vpt  => vpt_2
     255          IF ( humidity )        vpt_m  => vpt_1;   vpt  => vpt_2
    256256          IF ( cloud_physics )   ql   => ql_1
    257257          IF ( cloud_droplets )  THEN
     
    276276       w     => w_1;   w_p  => w_2;   tw_m  => w_3;   w_m  => w_3
    277277
    278        IF ( moisture  .OR.  passive_scalar )  THEN
     278       IF ( humidity  .OR.  passive_scalar )  THEN
    279279          qsws   => qsws_1
    280280          qswst  => qswst_1
    281281          q      => q_1;     q_p  => q_2;     tq_m => q_3;    q_m => q_3
    282           IF ( moisture )        vpt  => vpt_1
     282          IF ( humidity )        vpt  => vpt_1
    283283          IF ( cloud_physics )   ql   => ql_1
    284284          IF ( cloud_droplets )  THEN
     
    309309                km(:,j,i) = km1d
    310310                pt(:,j,i) = pt_init
    311              ENDDO
    312           ENDDO
    313           DO  i = nxl-1, nxr+uxrp+1
    314              DO  j = nys-1, nyn+1
    315311                u(:,j,i)  = u1d
    316              ENDDO
    317           ENDDO
    318           DO  i = nxl-1, nxr+1
    319              DO  j = nys-1, nyn+vynp+1
    320312                v(:,j,i)  = v1d
    321313             ENDDO
    322314          ENDDO
    323315
    324           IF ( moisture  .OR.  passive_scalar )  THEN
     316          IF ( humidity  .OR.  passive_scalar )  THEN
    325317             DO  i = nxl-1, nxr+1
    326318                DO  j = nys-1, nyn+1
     
    368360!--       This could actually be computed more accurately in the 1D model.
    369361!--       Update when opportunity arises!
    370           IF ( moisture  .OR.  passive_scalar )  qs = 0.0
     362          IF ( humidity  .OR.  passive_scalar )  qs = 0.0
    371363
    372364!
     
    419411             DO  j = nys-1, nyn+1
    420412                pt(:,j,i) = pt_init
    421              ENDDO
    422           ENDDO
    423           DO  i = nxl-1, nxr+uxrp+1
    424              DO  j = nys-1, nyn+1
    425413                u(:,j,i)  = u_init
    426              ENDDO
    427           ENDDO
    428           DO  i = nxl-1, nxr+1
    429              DO  j = nys-1, nyn+vynp+1
    430414                v(:,j,i)  = v_init
    431415             ENDDO
    432416          ENDDO
     417
    433418!
    434419!--       Set initial horizontal velocities at the lowest computational grid levels
     
    458443          ENDIF
    459444
    460           IF ( moisture  .OR.  passive_scalar )  THEN
     445          IF ( humidity  .OR.  passive_scalar )  THEN
    461446             DO  i = nxl-1, nxr+1
    462447                DO  j = nys-1, nyn+1
     
    482467          usws = 0.0
    483468          vsws = 0.0
    484           IF ( moisture  .OR.  passive_scalar ) qs = 0.0
     469          IF ( humidity  .OR.  passive_scalar ) qs = 0.0
    485470
    486471!
     
    499484!
    500485!--    Calculate virtual potential temperature
    501        IF ( moisture ) vpt = pt * ( 1.0 + 0.61 * q )
     486       IF ( humidity ) vpt = pt * ( 1.0 + 0.61 * q )
    502487
    503488!
     
    515500
    516501
    517        IF ( moisture )  THEN
     502       IF ( humidity )  THEN
    518503!
    519504!--       Store initial profile of total water content, virtual potential
     
    564549!
    565550!--       Determine the near-surface water flux
    566           IF ( moisture  .OR.  passive_scalar )  THEN
     551          IF ( humidity  .OR.  passive_scalar )  THEN
    567552             IF ( constant_waterflux )  THEN
    568553                qsws   = surface_waterflux
     
    585570             IF ( ASSOCIATED( tswst_m ) )  tswst_m = tswst
    586571
    587              IF ( moisture  .OR.  passive_scalar )  THEN
     572             IF ( humidity  .OR.  passive_scalar )  THEN
    588573                qswst = 0.0
    589574                IF ( ASSOCIATED( qswst_m ) )  qswst_m = qswst
     
    610595          ENDIF
    611596
    612           IF ( moisture  .OR.  passive_scalar )  THEN
     597          IF ( humidity  .OR.  passive_scalar )  THEN
    613598             IF ( .NOT. constant_waterflux )  THEN
    614599                qsws   = 0.0
     
    709694!--    If required, change the surface humidity/scalar at the start of the 3D
    710695!--    run
    711        IF ( ( moisture .OR. passive_scalar ) .AND. &
     696       IF ( ( humidity .OR. passive_scalar ) .AND. &
    712697            q_surface_initial_change /= 0.0 )  THEN
    713698          q(nzb,:,:) = q(nzb,:,:) + q_surface_initial_change
     
    722707!--    remove the divergences from the velocity field
    723708       IF ( create_disturbances )  THEN
    724           CALL disturb_field( nzb_u_inner, tend, u, uxrp,    0 )
    725           CALL disturb_field( nzb_v_inner, tend, v,    0, vynp )
     709          CALL disturb_field( nzb_u_inner, tend, u )
     710          CALL disturb_field( nzb_v_inner, tend, v )
    726711          n_sor = nsor_ini
    727712          CALL pres
     
    746731       e_p = e; pt_p = pt; u_p = u; v_p = v; w_p = w
    747732
    748        IF ( moisture  .OR.  passive_scalar )  THEN
     733       IF ( humidity  .OR.  passive_scalar )  THEN
    749734          IF ( ASSOCIATED( q_m ) )  q_m = q
    750735          IF ( timestep_scheme(1:5) == 'runge' )  tq_m = 0.0
    751736          q_p = q
    752           IF ( moisture  .AND.  ASSOCIATED( vpt_m ) )  vpt_m = vpt
     737          IF ( humidity  .AND.  ASSOCIATED( vpt_m ) )  vpt_m = vpt
    753738       ENDIF
    754739
     
    756741!--    Initialize old timelevels needed for radiation boundary conditions
    757742       IF ( outflow_l )  THEN
    758           u_m_l(:,:,:,-2) = u(:,:,-1:1)
    759           v_m_l(:,:,:,-2) = v(:,:,-1:1)
    760           w_m_l(:,:,:,-2) = w(:,:,-1:1)
    761           u_m_l(:,:,:,-1) = u(:,:,-1:1)
    762           v_m_l(:,:,:,-1) = v(:,:,-1:1)
    763           w_m_l(:,:,:,-1) = w(:,:,-1:1)
     743          u_m_l(:,:,:) = u(:,:,-1:1)
     744          v_m_l(:,:,:) = v(:,:,-1:1)
     745          w_m_l(:,:,:) = w(:,:,-1:1)
    764746       ENDIF
    765747       IF ( outflow_r )  THEN
    766           u_m_r(:,:,:,-2) = u(:,:,nx-1:nx+1)
    767           v_m_r(:,:,:,-2) = v(:,:,nx-1:nx+1)
    768           w_m_r(:,:,:,-2) = w(:,:,nx-1:nx+1)
    769           u_m_r(:,:,:,-1) = u(:,:,nx-1:nx+1)
    770           v_m_r(:,:,:,-1) = v(:,:,nx-1:nx+1)
    771           w_m_r(:,:,:,-1) = w(:,:,nx-1:nx+1)
     748          u_m_r(:,:,:) = u(:,:,nx-1:nx+1)
     749          v_m_r(:,:,:) = v(:,:,nx-1:nx+1)
     750          w_m_r(:,:,:) = w(:,:,nx-1:nx+1)
    772751       ENDIF
    773752       IF ( outflow_s )  THEN
    774           u_m_s(:,:,:,-2) = u(:,-1:1,:)
    775           v_m_s(:,:,:,-2) = v(:,-1:1,:)
    776           w_m_s(:,:,:,-2) = w(:,-1:1,:)
    777           u_m_s(:,:,:,-1) = u(:,-1:1,:)
    778           v_m_s(:,:,:,-1) = v(:,-1:1,:)
    779           w_m_s(:,:,:,-1) = w(:,-1:1,:)
     753          u_m_s(:,:,:) = u(:,-1:1,:)
     754          v_m_s(:,:,:) = v(:,-1:1,:)
     755          w_m_s(:,:,:) = w(:,-1:1,:)
    780756       ENDIF
    781757       IF ( outflow_n )  THEN
    782           u_m_n(:,:,:,-2) = u(:,ny-1:ny+1,:)
    783           v_m_n(:,:,:,-2) = v(:,ny-1:ny+1,:)
    784           w_m_n(:,:,:,-2) = w(:,ny-1:ny+1,:)
    785           u_m_n(:,:,:,-1) = u(:,ny-1:ny+1,:)
    786           v_m_n(:,:,:,-1) = v(:,ny-1:ny+1,:)
    787           w_m_n(:,:,:,-1) = w(:,ny-1:ny+1,:)
     758          u_m_n(:,:,:) = u(:,ny-1:ny+1,:)
     759          v_m_n(:,:,:) = v(:,ny-1:ny+1,:)
     760          w_m_n(:,:,:) = w(:,ny-1:ny+1,:)
    788761       ENDIF
    789762
     
    803776!--    including ghost points)
    804777       e_p = e; pt_p = pt; u_p = u; v_p = v; w_p = w
    805        IF ( moisture  .OR.  passive_scalar )  q_p = q
     778       IF ( humidity  .OR.  passive_scalar )  q_p = q
    806779
    807780    ELSE
  • palm/trunk/SOURCE/init_grid.f90

    r49 r75  
    44! Actual revisions:
    55! -----------------
    6 ! storage of topography height arrays zu_s_inner and zw_s_inner
     6! storage of topography height arrays zu_s_inner and zw_s_inner,
     7! 2nd+3rd argument removed from exchange horiz
    78!
    89! Former revisions:
     
    823824!
    824825!-- Need to set lateral boundary conditions for l_wall
    825     CALL exchange_horiz( l_wall, 0, 0 )
     826    CALL exchange_horiz( l_wall )
    826827
    827828    DEALLOCATE( corner_nl, corner_nr, corner_sl, corner_sr, nzb_local, &
  • palm/trunk/SOURCE/init_pegrid.f90

    r73 r75  
    55! -----------------
    66! uxrp, vynp eliminated,
    7 ! dirichlet/neumann changed to dirichlet/radiation, etc.
     7! dirichlet/neumann changed to dirichlet/radiation, etc.,
     8! poisfft_init is only called if fft-solver is switched on
    89!
    910! Former revisions:
     
    493494    nzta = nz
    494495    nnz  = nz
    495 
    496 !
    497 !-- For non-cyclic boundaries extend array u (v) by one gridpoint
    498 !    IF ( bc_lr /= 'cyclic' )  uxrp = 1
    499 !    IF ( bc_ns /= 'cyclic' )  vynp = 1
    500496
    501497!
     
    757753          inflow_r  = .TRUE.
    758754       ENDIF
    759 !       uxrp      = 1
    760755    ENDIF
    761756
     
    774769          outflow_n = .TRUE.
    775770       ENDIF
    776 !       vynp      = 1
    777     ENDIF
    778 
    779 !
    780 !-- Additional MPI derived data type for the exchange of ghost points along x
    781 !-- needed in case of non-cyclic boundary conditions along y on the northmost
    782 !-- processors (for the exchange of the enlarged v array)
    783     IF ( bc_ns /= 'cyclic'  .AND.  pnorth == MPI_PROC_NULL )  THEN
    784        ngp_yz_p = ( nzt - nzb + 2 ) * ( nyn + vynp - nys + 3 )
    785        CALL MPI_TYPE_VECTOR( nxr-nxl+3, nzt-nzb+2, ngp_yz_p, &
    786                              MPI_REAL, type_xz_p, ierr )
    787        CALL MPI_TYPE_COMMIT( type_xz_p, ierr )
    788     ENDIF
     771    ENDIF
     772
    789773#else
    790774    IF ( bc_lr == 'dirichlet/radiation' )  THEN
    791775       inflow_l  = .TRUE.
    792776       outflow_r = .TRUE.
    793 !       uxrp      = 1
    794777    ELSEIF ( bc_lr == 'radiation/dirichlet' )  THEN
    795778       outflow_l = .TRUE.
     
    803786       outflow_n = .TRUE.
    804787       inflow_s  = .TRUE.
    805 !       vynp      = 1
    806788    ENDIF
    807789#endif
     
    809791    IF ( psolver == 'poisfft_hybrid' )  THEN
    810792       CALL poisfft_hybrid_ini
    811     ELSE
     793    ELSEIF ( psolver == 'poisfft' )  THEN
    812794       CALL poisfft_init
    813795    ENDIF
  • palm/trunk/SOURCE/init_pt_anomaly.f90

    r39 r75  
    44! Actual revisions:
    55! -----------------
    6 !
     6! 2nd+3rd argument removed from exchange horiz
    77!
    88! Former revisions:
     
    7070!
    7171!-- Exchange of boundary values for temperature
    72     CALL exchange_horiz( pt, 0, 0 )
     72    CALL exchange_horiz( pt )
    7373
    7474
  • palm/trunk/SOURCE/init_rankine.f90

    r4 r75  
    44! Actual revisions:
    55! -----------------
    6 !
     6! uxrp, vynp eliminated, 2nd+3rd argument removed from exchange horiz
    77!
    88! Former revisions:
     
    129129!
    130130!-- Exchange of boundary values for the velocities.
    131     CALL exchange_horiz( u, uxrp,    0 )
    132     CALL exchange_horiz( v,    0, vynp )
     131    CALL exchange_horiz( u )
     132    CALL exchange_horiz( v )
    133133!
    134134!-- Make velocity field nondivergent.
  • palm/trunk/SOURCE/modules.f90

    r73 r75  
    99! +loop_optimization, netcdf_64bit_3d, zu_s_inner, zw_w_inner, id_var_zusi_*,
    1010! id_var_zwwi_*, ts_value, u_nzb_p1_for_vfc, v_nzb_p1_for_vfc, pt_reference,
    11 ! use_pt_reference, precipitation_amount_interval
    12 ! +age_m in particle_type
    13 ! -data_output_ts, dots_n
    14 ! arrays dots_label and dots_unit now dimensioned with dots_max
     11! use_pt_reference, precipitation_amount_interval, revision
     12! +age_m in particle_type, moisture renamed humidity,
     13! -data_output_ts, dots_n, uvmean_outflow, uxrp, vynp,
     14! arrays dots_label and dots_unit now dimensioned with dots_max,
     15! setting of palm version moved to main program
    1516!
    1617! Former revisions:
     
    8182    REAL, DIMENSION(:), ALLOCATABLE ::                                         &
    8283          ddzu, dd2zu, dzu, ddzw, dzw, km_damp_x, km_damp_y, l_grid, pt_init,  &
    83           q_init, rdf, ug, uvmean_outflow, uvmean_outflow_l, u_init,           &
    84           u_nzb_p1_for_vfc, vg, v_init, v_nzb_p1_for_vfc, zu, zw
     84          q_init, rdf, ug, u_init, u_nzb_p1_for_vfc, vg, v_init,               &
     85          v_nzb_p1_for_vfc, zu, zw
    8586
    8687    REAL, DIMENSION(:,:), ALLOCATABLE ::                                       &
     
    9697
    9798    REAL, DIMENSION(:,:,:), ALLOCATABLE ::                                     &
    98           d, diss, l_wall, tend
     99          d, diss, l_wall, tend, u_m_l, u_m_n, u_m_r, u_m_s, v_m_l, v_m_n,     &
     100          v_m_r, v_m_s, w_m_l, w_m_n, w_m_r, w_m_s
    99101
    100102    REAL, DIMENSION(:,:,:), ALLOCATABLE, TARGET ::                             &
     
    111113          vpt, vpt_m, w, w_m, w_p
    112114
    113     REAL, DIMENSION(:,:,:,:), ALLOCATABLE ::  rif_wall, u_m_l, u_m_n, u_m_r,   &
    114           u_m_s, v_m_l, v_m_n, v_m_r, v_m_s, w_m_l, w_m_n, w_m_r, w_m_s
    115                                                    
     115    REAL, DIMENSION(:,:,:,:), ALLOCATABLE ::  rif_wall
    116116
    117117    SAVE
     
    217217    CHARACTER (LEN=8)   ::  run_date, run_time
    218218    CHARACTER (LEN=9)   ::  simulated_time_chr
    219     CHARACTER (LEN=12)  ::  version = 'PALM   3.2'
     219    CHARACTER (LEN=12)  ::  version = ' ', revision = ' '
    220220    CHARACTER (LEN=16)  ::  loop_optimization = 'default', &
    221221                            momentum_advec = 'pw-scheme', &
     
    236236    CHARACTER (LEN=40)  ::  avs_data_file, topography = 'flat'
    237237    CHARACTER (LEN=64)  ::  host
    238     CHARACTER (LEN=80)  ::  log_message, run_description_header, run_identifier
    239     CHARACTER (LEN=100) ::  initializing_actions = ' '
     238    CHARACTER (LEN=80)  ::  log_message, run_identifier
     239    CHARACTER (LEN=100) ::  initializing_actions = ' ', run_description_header
    240240
    241241    CHARACTER (LEN=7),  DIMENSION(100) ::  do3d_comp_prec = ' '
     
    293293                first_call_advec_particles = .TRUE., &
    294294                force_print_header = .FALSE., galilei_transformation = .FALSE.,&
    295                 inflow_l = .FALSE., inflow_n = .FALSE., inflow_r = .FALSE., &
    296                 inflow_s = .FALSE., iso2d_output = .FALSE., &
    297                 mg_switch_to_pe0 = .FALSE., moisture = .FALSE., &
     295                humidity = .FALSE., inflow_l = .FALSE., inflow_n = .FALSE., &
     296                inflow_r = .FALSE., inflow_s = .FALSE., iso2d_output = .FALSE.,&
     297                mg_switch_to_pe0 = .FALSE., &
    298298                netcdf_output = .FALSE., netcdf_64bit = .FALSE., &
    299299                netcdf_64bit_3d = .TRUE., &
     
    528528    INTEGER ::  ngp_sums, nnx, nx = 0, nxa, nxl, nxr, nxra, nny, ny = 0, nya,  &
    529529                nyn, nyna, nys, nnz, nz = 0, nza, nzb, nzb_diff, nzt, nzta,    &
    530                 nzt_diff, uxrp = 0, vynp = 0
     530                nzt_diff
    531531
    532532    INTEGER, DIMENSION(:), ALLOCATABLE ::                                      &
     
    889889#if defined( __parallel )
    890890    INTEGER ::  comm1dx, comm1dy, comm2d, comm_palm, ierr, myidx, myidy,       &
    891                 ndim = 2, ngp_y, ngp_yz_p, pleft, pnorth, pright, psouth,      &
     891                ndim = 2, ngp_y, pleft, pnorth, pright, psouth,                &
    892892                sendrecvcount_xy, sendrecvcount_yz, sendrecvcount_zx,          &
    893893                sendrecvcount_zyd, sendrecvcount_yxd,                          &
    894                 type_x, type_x_int, type_xz_p, ibuf(12), pcoord(2), pdims(2),  &
     894                type_x, type_x_int, ibuf(12), pcoord(2), pdims(2),             &
    895895                status(MPI_STATUS_SIZE)
    896896
  • palm/trunk/SOURCE/palm.f90

    r70 r75  
    55! -----------------
    66! __vtk directives removed, write_particles is called only in case of particle
    7 ! advection switched on, open unit 9 for debug output
     7! advection switched on, open unit 9 for debug output,
     8! setting of palm version moved from modules to here
    89!
    910! Former revisions:
     
    5758    CHARACTER (LEN=1) ::  cdum
    5859    INTEGER           ::  i, run_description_header_i(80)
     60
     61    version = 'PALM 3.2'
    5962
    6063#if defined( __parallel )
  • palm/trunk/SOURCE/parin.f90

    r72 r75  
    55! -----------------
    66! +dt_max, netcdf_64bit_3d, precipitation_amount_interval in d3par,
    7 ! +loop_optimization, pt_reference in inipar, -data_output_ts
     7! +loop_optimization, pt_reference in inipar, -data_output_ts,
     8! moisture renamed humidity
    89!
    910! Former revisions:
     
    5859                       dz_stretch_factor, dz_stretch_level, e_min, &
    5960                       end_time_1d, fft_method, galilei_transformation, &
    60                        grid_matching, inflow_disturbance_begin, &
     61                       grid_matching, humidity, inflow_disturbance_begin, &
    6162                       inflow_disturbance_end, initializing_actions, &
    6263                       km_constant, km_damp_max, long_filter_factor, &
    63                        loop_optimization, mixing_length_1d, moisture, &
     64                       loop_optimization, mixing_length_1d, &
    6465                       momentum_advec, netcdf_precision, npex, npey, nsor_ini, &
    6566                       nx, ny, nz, omega, outflow_damping_width, &
     
    117118
    118119
    119     NAMELIST /envpar/  host, maximum_cpu_time_allowed, run_identifier, &
    120                        tasks_per_node, write_binary
     120    NAMELIST /envpar/  host, maximum_cpu_time_allowed, revision, &
     121                       run_identifier, tasks_per_node, write_binary
    121122
    122123
  • palm/trunk/SOURCE/poismg.f90

    r4 r75  
    88! Actual revisions:
    99! -----------------
    10 !
     10! 2nd+3rd argument removed from exchange horiz
    1111!
    1212! Former revisions:
     
    5959!
    6060!-- Some boundaries have to be added to divergence array
    61     CALL exchange_horiz( d, 0, 0 )
     61    CALL exchange_horiz( d )
    6262    d(nzb,:,:) = d(nzb+1,:,:)
    6363
     
    164164!
    165165!-- Horizontal boundary conditions
    166     CALL exchange_horiz( r, 0, 0 )
     166    CALL exchange_horiz( r )
    167167
    168168    IF ( bc_lr /= 'cyclic' )  THEN
     
    257257!
    258258!-- Horizontal boundary conditions
    259     CALL exchange_horiz( f_mg, 0, 0 )
     259    CALL exchange_horiz( f_mg )
    260260
    261261    IF ( bc_lr /= 'cyclic' )  THEN
     
    355355!
    356356!-- Horizontal boundary conditions
    357     CALL exchange_horiz( temp, 0, 0 )
     357    CALL exchange_horiz( temp )
    358358
    359359    IF ( bc_lr /= 'cyclic' )  THEN
     
    627627!
    628628!--       Horizontal boundary conditions
    629           CALL exchange_horiz( p_mg, 0, 0 )
     629          CALL exchange_horiz( p_mg )
    630630
    631631          IF ( bc_lr /= 'cyclic' )  THEN
  • palm/trunk/SOURCE/prandtl_fluxes.f90

    r4 r75  
    44! Actual revisions:
    55! -----------------
    6 !
     6! moisture renamed humidity
    77!
    88! Former revisions:
     
    9797!
    9898!-- Compute z_p/L (corresponds to the Richardson-flux number)
    99     IF ( .NOT. moisture ) THEN
     99    IF ( .NOT. humidity ) THEN
    100100       !$OMP PARALLEL DO PRIVATE( k, z_p )
    101101       DO  i = nxl-1, nxr+1
     
    274274!
    275275!-- If required compute q*
    276     IF ( moisture  .OR.  passive_scalar )  THEN
     276    IF ( humidity  .OR.  passive_scalar )  THEN
    277277       IF ( constant_waterflux )  THEN
    278278!
     
    332332    CALL exchange_horiz_2d( usws )
    333333    CALL exchange_horiz_2d( vsws )
    334     IF ( moisture  .OR.  passive_scalar )  CALL exchange_horiz_2d( qsws )
     334    IF ( humidity  .OR.  passive_scalar )  CALL exchange_horiz_2d( qsws )
    335335
    336336!
     
    347347!
    348348!-- Compute the vertical water/scalar flux
    349     IF ( .NOT. constant_heatflux .AND. ( moisture .OR. passive_scalar ) ) THEN
     349    IF ( .NOT. constant_heatflux .AND. ( humidity .OR. passive_scalar ) ) THEN
    350350       !$OMP PARALLEL DO
    351351       DO  i = nxl-1, nxr+1
  • palm/trunk/SOURCE/pres.f90

    r73 r75  
    44! Actual revisions:
    55! -----------------
    6 !
     6! Volume flow control for non-cyclic boundary conditions added (currently only
     7! for the north boundary!!), 2nd+3rd argument removed from exchange horiz
    78!
    89! Former revisions:
     
    5455       ALLOCATE( d(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1) )
    5556    ENDIF
     57
     58!
     59!-- Conserve the volume flow at the outflow in case of non-cyclic lateral
     60!-- boundary conditions
     61    IF ( conserve_volume_flow  .AND.  bc_ns == 'radiation/dirichlet')  THEN
     62
     63       volume_flow(2)   = 0.0
     64       volume_flow_l(2) = 0.0
     65
     66       IF ( nyn == ny )  THEN
     67          j = ny+1
     68          DO  i = nxl, nxr
     69!
     70!--          Sum up the volume flow through the north boundary
     71             DO  k = nzb_2d(j,i) + 1, nzt
     72                volume_flow_l(2) = volume_flow_l(2) + v(k,j,i) * dzu(k)
     73             ENDDO
     74          ENDDO
     75       ENDIF
     76#if defined( __parallel )   
     77       CALL MPI_ALLREDUCE( volume_flow_l(2), volume_flow(2), 1, MPI_REAL, &
     78                           MPI_SUM, comm1dx, ierr )   
     79#else
     80       volume_flow = volume_flow_l 
     81#endif
     82       volume_flow_offset(2) = ( volume_flow_initial(2) - volume_flow(2) )    &
     83                               / volume_flow_area(2)                         
     84
     85       IF ( outflow_n )  THEN
     86          j = nyn+1
     87          DO  i = nxl, nxr
     88             DO  k = nzb_v_inner(j,i) + 1, nzt
     89                v(k,j,i) = v(k,j,i) + volume_flow_offset(2)
     90             ENDDO
     91          ENDDO
     92       ENDIF
     93
     94       CALL exchange_horiz( v )
     95
     96    ENDIF
     97
    5698
    5799!
     
    322364!
    323365!--    Exchange boundaries for p
    324        CALL exchange_horiz( tend, 0, 0 )
     366       CALL exchange_horiz( tend )
    325367     
    326368    ELSEIF ( psolver == 'sor' )  THEN
     
    365407!-- Correction of the provisional velocities with the current perturbation
    366408!-- pressure just computed
    367     IF ( bc_lr /= 'cyclic'  .OR.  bc_ns /= 'cyclic' )  uvmean_outflow_l = 0.0
    368     IF ( conserve_volume_flow )  THEN
     409    IF ( conserve_volume_flow  .AND. &
     410         ( bc_lr == 'cyclic'  .OR.  bc_ns == 'cyclic' ) )  THEN
    369411       volume_flow_l(1) = 0.0
    370412       volume_flow_l(2) = 0.0
     
    393435
    394436!
    395 !--       Sum up the horizontal velocity along the outflow plane (in case
    396 !--       of non-cyclic boundary conditions). The respective mean velocity
    397 !--       is calculated from this in routine boundary_conds.
    398 !          IF ( outflow_l  .AND.  i == nxl )  THEN
    399 !             !$OMP CRITICAL
    400 !             DO  k = nzb, nzt+1
    401 !                uvmean_outflow_l(k) = uvmean_outflow_l(k) + v(k,j,nxl)
    402 !             ENDDO
    403 !             !$OMP END CRITICAL
    404 !          ELSEIF ( outflow_r  .AND.  i == nxr )  THEN
    405 !             !$OMP CRITICAL
    406 !             DO  k = nzb, nzt+1
    407 !                uvmean_outflow_l(k) = uvmean_outflow_l(k) + v(k,j,nxr)
    408 !             ENDDO
    409 !             !$OMP END CRITICAL
    410 !          ELSEIF ( outflow_s  .AND.  j == nys )  THEN
    411 !             !$OMP CRITICAL
    412 !             DO  k = nzb, nzt+1
    413 !                uvmean_outflow_l(k) = uvmean_outflow_l(k) + u(k,nys,i)
    414 !             ENDDO
    415 !             !$OMP END CRITICAL
    416 !          ELSEIF ( outflow_n  .AND.  j == nyn )  THEN
    417 !             !$OMP CRITICAL
    418 !             DO  k = nzb, nzt+1
    419 !                uvmean_outflow_l(k) = uvmean_outflow_l(k) + u(k,nyn,i)
    420 !             ENDDO
    421 !             !$OMP END CRITICAL
    422 !          ENDIF
    423 
    424 !
    425437!--       Sum up the volume flow through the right and north boundary
    426           IF ( conserve_volume_flow  .AND.  i == nx )  THEN
     438          IF ( conserve_volume_flow  .AND.  bc_lr == 'cyclic'  .AND. &
     439               i == nx )  THEN
    427440             !$OMP CRITICAL
    428441             DO  k = nzb_2d(j,i) + 1, nzt
     
    431444             !$OMP END CRITICAL
    432445          ENDIF
    433           IF ( conserve_volume_flow  .AND.  j == ny )  THEN
     446          IF ( conserve_volume_flow  .AND.  bc_ns == 'cyclic'  .AND. &
     447               j == ny )  THEN
    434448             !$OMP CRITICAL
    435449             DO  k = nzb_2d(j,i) + 1, nzt
     
    445459!
    446460!-- Conserve the volume flow
    447     IF ( conserve_volume_flow )  THEN
     461    IF ( conserve_volume_flow  .AND. &
     462         ( bc_lr == 'cyclic'  .OR.  bc_ns == 'cyclic' ) )  THEN
    448463
    449464#if defined( __parallel )   
     
    461476       DO  i = nxl, nxr
    462477          DO  j = nys, nyn
    463              DO  k = nzb_u_inner(j,i) + 1, nzt
    464                 u(k,j,i) = u(k,j,i) + volume_flow_offset(1)
    465              ENDDO
    466              DO  k = nzb_v_inner(j,i) + 1, nzt
    467                 v(k,j,i) = v(k,j,i) + volume_flow_offset(2)
    468              ENDDO
     478             IF ( bc_lr == 'cyclic' )  THEN
     479                DO  k = nzb_u_inner(j,i) + 1, nzt
     480                   u(k,j,i) = u(k,j,i) + volume_flow_offset(1)
     481                ENDDO
     482             ENDIF
     483             IF ( bc_ns == 'cyclic' )  THEN
     484                DO  k = nzb_v_inner(j,i) + 1, nzt
     485                   v(k,j,i) = v(k,j,i) + volume_flow_offset(2)
     486                ENDDO
     487             ENDIF
    469488          ENDDO
    470489       ENDDO
     
    475494!
    476495!-- Exchange of boundaries for the velocities
    477     CALL exchange_horiz( u, uxrp,    0 )
    478     CALL exchange_horiz( v,    0, vynp )
    479     CALL exchange_horiz( w,    0,    0 )
     496    CALL exchange_horiz( u )
     497    CALL exchange_horiz( v )
     498    CALL exchange_horiz( w )
    480499
    481500!
  • palm/trunk/SOURCE/production_e.f90

    r73 r75  
    55! -----------------
    66! Wall functions now include diabatic conditions, call of routine wall_fluxes_e,
    7 ! reference temperature pt_reference can be used in buoyancy term
     7! reference temperature pt_reference can be used in buoyancy term,
     8! moisture renamed humidity
    89!
    910! Former revisions:
     
    358359!
    359360!--       Calculate TKE production by buoyancy
    360           IF ( .NOT. moisture )  THEN
     361          IF ( .NOT. humidity )  THEN
    361362
    362363             IF ( use_pt_reference )  THEN
     
    756757!
    757758!--    Calculate TKE production by buoyancy
    758        IF ( .NOT. moisture )  THEN
     759       IF ( .NOT. humidity )  THEN
    759760
    760761          IF ( use_pt_reference )  THEN
  • palm/trunk/SOURCE/prognostic_equations.f90

    r73 r75  
    55! -----------------
    66! checking for negative q and limiting for positive values,
    7 ! z0 removed from arguments in calls of diffusion_u/v/w,
    8 ! subroutine names changed to .._noopt, .._cache, and .._vector
     7! z0 removed from arguments in calls of diffusion_u/v/w, uxrp, vynp eliminated,
     8! subroutine names changed to .._noopt, .._cache, and .._vector,
     9! moisture renamed humidity
    910!
    1011! Former revisions:
     
    102103!-- global communication
    103104    CALL calc_mean_pt_profile( pt, 4 )
    104     IF ( moisture )  CALL calc_mean_pt_profile( vpt, 44 )
     105    IF ( humidity )  CALL calc_mean_pt_profile( vpt, 44 )
    105106
    106107!
     
    117118!
    118119!-- u-tendency terms with no communication
    119     DO  i = nxl, nxr+uxrp
     120    DO  i = nxl, nxr
    120121       DO  j = nys, nyn
    121122!
     
    186187!-- v-tendency terms with no communication
    187188    DO  i = nxl, nxr
    188        DO  j = nys, nyn+vynp
     189       DO  j = nys, nyn
    189190!
    190191!--       Tendency terms
     
    273274          ENDIF
    274275          CALL coriolis( i, j, 3 )
    275           IF ( .NOT. moisture )  THEN
     276          IF ( .NOT. humidity )  THEN
    276277             CALL buoyancy( i, j, pt, 3, 4 )
    277278          ELSE
     
    410411!
    411412!-- If required, compute prognostic equation for total water content / scalar
    412     IF ( moisture  .OR.  passive_scalar )  THEN
     413    IF ( humidity  .OR.  passive_scalar )  THEN
    413414
    414415       CALL cpu_log( log_point(29), 'q/s-equation', 'start' )
     
    550551             IF ( scalar_advec == 'bc-scheme'  .AND.  &
    551552                  .NOT. use_upstream_for_tke )  THEN
    552                 IF ( .NOT. moisture )  THEN
     553                IF ( .NOT. humidity )  THEN
    553554                   CALL diffusion_e( i, j, ddzu, dd2zu, ddzw, diss, e, km, &
    554555                                     l_grid, pt, rif, tend, zu )
     
    575576                IF ( tsc(2) == 2.0  .AND.  timestep_scheme(1:8) == 'leapfrog' )&
    576577                THEN
    577                    IF ( .NOT. moisture )  THEN
     578                   IF ( .NOT. humidity )  THEN
    578579                      CALL diffusion_e( i, j, ddzu, dd2zu, ddzw, diss, e_m, &
    579580                                        km_m, l_grid, pt_m, rif_m, tend, zu )
     
    583584                   ENDIF
    584585                ELSE
    585                    IF ( .NOT. moisture )  THEN
     586                   IF ( .NOT. humidity )  THEN
    586587                      CALL diffusion_e( i, j, ddzu, dd2zu, ddzw, diss, e, km, &
    587588                                        l_grid, pt, rif, tend, zu )
     
    661662!-- global communication
    662663    CALL calc_mean_pt_profile( pt, 4 )
    663     IF ( moisture )  CALL calc_mean_pt_profile( vpt, 44 )
     664    IF ( humidity )  CALL calc_mean_pt_profile( vpt, 44 )
    664665    IF ( .NOT. constant_diffusion )  CALL production_e_init
    665666
     
    669670!$OMP PARALLEL private (i,j,k)
    670671!$OMP DO
    671     DO  i = nxl, nxr+uxrp       ! Additional levels for non cyclic boundary
    672        DO  j = nys, nyn+vynp    ! conditions are included
     672    DO  i = nxl, nxr
     673       DO  j = nys, nyn
    673674!
    674675!--       Tendency terms for u-velocity component
     
    789790             ENDIF
    790791             CALL coriolis( i, j, 3 )
    791              IF ( .NOT. moisture )  THEN
     792             IF ( .NOT. humidity )  THEN
    792793                CALL buoyancy( i, j, pt, 3, 4 )
    793794             ELSE
     
    881882!--          If required, compute prognostic equation for total water content /
    882883!--          scalar
    883              IF ( moisture  .OR.  passive_scalar )  THEN
     884             IF ( humidity  .OR.  passive_scalar )  THEN
    884885
    885886!
     
    954955                IF ( tsc(2) == 2.0  .AND.  timestep_scheme(1:8) == 'leapfrog' )&
    955956                THEN
    956                    IF ( .NOT. moisture )  THEN
     957                   IF ( .NOT. humidity )  THEN
    957958                      CALL diffusion_e( i, j, ddzu, dd2zu, ddzw, diss, e_m, &
    958959                                        km_m, l_grid, pt_m, rif_m, tend, zu )
     
    962963                   ENDIF
    963964                ELSE
    964                    IF ( .NOT. moisture )  THEN
     965                   IF ( .NOT. humidity )  THEN
    965966                      CALL diffusion_e( i, j, ddzu, dd2zu, ddzw, diss, e, km, &
    966967                                        l_grid, pt, rif, tend, zu )
     
    10321033!-- global communication
    10331034    CALL calc_mean_pt_profile( pt, 4 )
    1034     IF ( moisture )  CALL calc_mean_pt_profile( vpt, 44 )
     1035    IF ( humidity )  CALL calc_mean_pt_profile( vpt, 44 )
    10351036
    10361037!
     
    10681069!
    10691070!-- Prognostic equation for u-velocity component
    1070     DO  i = nxl, nxr+uxrp
     1071    DO  i = nxl, nxr
    10711072       DO  j = nys, nyn
    10721073          DO  k = nzb_u_inner(j,i)+1, nzt
     
    10851086    IF ( timestep_scheme(1:5) == 'runge' )  THEN
    10861087       IF ( intermediate_timestep_count == 1 )  THEN
    1087           DO  i = nxl, nxr+uxrp
     1088          DO  i = nxl, nxr
    10881089             DO  j = nys, nyn
    10891090                DO  k = nzb_u_inner(j,i)+1, nzt
     
    10941095       ELSEIF ( intermediate_timestep_count < &
    10951096                intermediate_timestep_count_max )  THEN
    1096           DO  i = nxl, nxr+uxrp
     1097          DO  i = nxl, nxr
    10971098             DO  j = nys, nyn
    10981099                DO  k = nzb_u_inner(j,i)+1, nzt
     
    11401141!-- Prognostic equation for v-velocity component
    11411142    DO  i = nxl, nxr
    1142        DO  j = nys, nyn+vynp
     1143       DO  j = nys, nyn
    11431144          DO  k = nzb_v_inner(j,i)+1, nzt
    11441145             v_p(k,j,i) = ( 1.0-tsc(1) ) * v_m(k,j,i) + tsc(1) * v(k,j,i) +    &
     
    11571158       IF ( intermediate_timestep_count == 1 )  THEN
    11581159          DO  i = nxl, nxr
    1159              DO  j = nys, nyn+vynp
     1160             DO  j = nys, nyn
    11601161                DO  k = nzb_v_inner(j,i)+1, nzt
    11611162                   tv_m(k,j,i) = tend(k,j,i)
     
    11661167                intermediate_timestep_count_max )  THEN
    11671168          DO  i = nxl, nxr
    1168              DO  j = nys, nyn+vynp
     1169             DO  j = nys, nyn
    11691170                DO  k = nzb_v_inner(j,i)+1, nzt
    11701171                   tv_m(k,j,i) = -9.5625 * tend(k,j,i) + 5.3125 * tv_m(k,j,i)
     
    12061207    ENDIF
    12071208    CALL coriolis( 3 )
    1208     IF ( .NOT. moisture )  THEN
     1209    IF ( .NOT. humidity )  THEN
    12091210       CALL buoyancy( pt, 3, 4 )
    12101211    ELSE
     
    13551356!
    13561357!-- If required, compute prognostic equation for total water content / scalar
    1357     IF ( moisture  .OR.  passive_scalar )  THEN
     1358    IF ( humidity  .OR.  passive_scalar )  THEN
    13581359
    13591360       CALL cpu_log( log_point(29), 'q/s-equation', 'start' )
     
    14931494       IF ( scalar_advec == 'bc-scheme'  .AND.  .NOT. use_upstream_for_tke )  &
    14941495       THEN
    1495           IF ( .NOT. moisture )  THEN
     1496          IF ( .NOT. humidity )  THEN
    14961497             CALL diffusion_e( ddzu, dd2zu, ddzw, diss, e, km, l_grid, pt, &
    14971498                               rif, tend, zu )
     
    15161517          ENDIF
    15171518          IF ( tsc(2) == 2.0  .AND.  timestep_scheme(1:8) == 'leapfrog' )  THEN
    1518              IF ( .NOT. moisture )  THEN
     1519             IF ( .NOT. humidity )  THEN
    15191520                CALL diffusion_e( ddzu, dd2zu, ddzw, diss, e_m, km_m, l_grid, &
    15201521                                  pt_m, rif_m, tend, zu )
     
    15241525             ENDIF
    15251526          ELSE
    1526              IF ( .NOT. moisture )  THEN
     1527             IF ( .NOT. humidity )  THEN
    15271528                CALL diffusion_e( ddzu, dd2zu, ddzw, diss, e, km, l_grid, pt, &
    15281529                                  rif, tend, zu )
  • palm/trunk/SOURCE/read_var_list.f90

    r63 r75  
    44! Actual revisions:
    55! -----------------
    6 ! +loop_optimization, pt_reference
     6! +loop_optimization, pt_reference, moisture renamed humidity
    77!
    88! Former revisions:
     
    212212          CASE ( 'mixing_length_1d' )
    213213             READ ( 13 )  mixing_length_1d
    214           CASE ( 'moisture' )
    215              READ ( 13 )  moisture
     214          CASE ( 'humidity' )
     215             READ ( 13 )  humidity
    216216          CASE ( 'momentum_advec' )
    217217             READ ( 13 )  momentum_advec
  • palm/trunk/SOURCE/sor.f90

    r4 r75  
    44! Actual revisions:
    55! -----------------
    6 !
     6! 2nd+3rd argument removed from exchange horiz
    77!
    88! Former revisions:
     
    9898!
    9999!--    Exchange of boundary values for p.
    100        CALL exchange_horiz( p, 0, 0 )
     100       CALL exchange_horiz( p )
    101101
    102102!
     
    143143!
    144144!--    Exchange of boundary values for p.
    145        CALL exchange_horiz( p, 0, 0 )
     145       CALL exchange_horiz( p )
    146146
    147147!
  • palm/trunk/SOURCE/swap_timelevel.f90

    r39 r75  
    44! Actual revisions:
    55! -----------------
    6 !
     6! moisture renamed humidity
    77!
    88! Former revisions:
     
    5656                e_m => e_1;  e => e_2;  e_p => e_3
    5757             ENDIF
    58              IF ( moisture  .OR.  passive_scalar )  THEN
     58             IF ( humidity  .OR.  passive_scalar )  THEN
    5959                q_m => q_1;  q => q_2;  q_p => q_3
    6060             ENDIF
     
    6969                e_m => e_2;  e => e_3;  e_p => e_1
    7070             ENDIF
    71              IF ( moisture  .OR.  passive_scalar )  THEN
     71             IF ( humidity  .OR.  passive_scalar )  THEN
    7272                q_m => q_2;  q => q_3;  q_p => q_1
    7373             ENDIF
     
    8282                e_m => e_3;  e => e_1;  e_p => e_2
    8383             ENDIF
    84              IF ( moisture  .OR.  passive_scalar )  THEN
     84             IF ( humidity  .OR.  passive_scalar )  THEN
    8585                q_m => q_3;  q => q_1;  q_p => q_2
    8686             ENDIF
     
    105105                e => e_1;  e_p => e_2
    106106             ENDIF
    107              IF ( moisture  .OR.  passive_scalar )  THEN
     107             IF ( humidity  .OR.  passive_scalar )  THEN
    108108                q => q_1;  q_p => q_2
    109109             ENDIF
     
    119119                   vsws_m => vsws_1;  vsws => vsws_2
    120120                   shf_m  => shf_1;   shf  => shf_2
    121                    IF ( moisture  .OR.  passive_scalar )  THEN
     121                   IF ( humidity  .OR.  passive_scalar )  THEN
    122122                      qsws_m => qsws_1;  qsws => qsws_2
    123123                   ENDIF
     
    128128                IF ( use_top_fluxes )  THEN
    129129                   tswst_m => tswst_1;  tswst => tswst_2
    130                    IF ( moisture  .OR.  passive_scalar )  THEN
     130                   IF ( humidity  .OR.  passive_scalar )  THEN
    131131                      qswst_m => qswst_1;  qswst => qswst_2
    132132                   ENDIF
     
    134134             ENDIF
    135135
    136              IF ( moisture )  THEN
     136             IF ( humidity )  THEN
    137137                vpt_m => vpt_1;  vpt => vpt_2
    138138             ENDIF
     
    151151                e => e_2;  e_p => e_1
    152152             ENDIF
    153              IF ( moisture  .OR.  passive_scalar )  THEN
     153             IF ( humidity  .OR.  passive_scalar )  THEN
    154154                q => q_2;  q_p => q_1
    155155             ENDIF
     
    164164                   vsws_m => vsws_2;  vsws => vsws_1
    165165                   shf_m  => shf_2;   shf  => shf_1
    166                    IF ( moisture  .OR.  passive_scalar )  THEN
     166                   IF ( humidity  .OR.  passive_scalar )  THEN
    167167                      qsws_m => qsws_2;  qsws => qsws_1
    168168                   ENDIF
     
    173173                IF ( use_top_fluxes )  THEN
    174174                   tswst_m  => tswst_2;  tswst => tswst_1
    175                    IF ( moisture  .OR.  passive_scalar )  THEN
     175                   IF ( humidity  .OR.  passive_scalar )  THEN
    176176                      qswst_m => qswst_2;  qswst => qswst_1
    177177                   ENDIF
     
    179179             ENDIF
    180180
    181              IF ( moisture )  THEN
     181             IF ( humidity )  THEN
    182182                vpt_m => vpt_2;  vpt => vpt_1
    183183             ENDIF
  • palm/trunk/SOURCE/time_integration.f90

    r73 r75  
    77! and counters,
    88! calls of prognostic_equations_.. changed to .._noopt, .._cache, and
    9 ! .._vector, these calls are now controlled by switch loop_optimization
     9! .._vector, these calls are now controlled by switch loop_optimization,
     10! uxrp, vynp eliminated, 2nd+3rd argument removed from exchange horiz,
     11! moisture renamed humidity
    1012!
    1113! Former revisions:
     
    136138!--       Exchange of ghost points (lateral boundary conditions)
    137139          CALL cpu_log( log_point(26), 'exchange-horiz-progn', 'start' )
    138           CALL exchange_horiz( u_p, uxrp,    0 )
    139           CALL exchange_horiz( v_p,    0, vynp )
    140           CALL exchange_horiz( w_p,    0,    0 )
    141           CALL exchange_horiz( pt_p,   0,    0 )
    142           IF ( .NOT. constant_diffusion       )  CALL exchange_horiz( e_p, 0, 0)
    143           IF ( moisture  .OR.  passive_scalar )  CALL exchange_horiz( q_p, 0, 0)
     140          CALL exchange_horiz( u_p )
     141          CALL exchange_horiz( v_p )
     142          CALL exchange_horiz( w_p )
     143          CALL exchange_horiz( pt_p )
     144          IF ( .NOT. constant_diffusion       )  CALL exchange_horiz( e_p )
     145          IF ( humidity  .OR.  passive_scalar )  CALL exchange_horiz( q_p )
    144146          IF ( cloud_droplets )  THEN
    145              CALL exchange_horiz( ql,    0, 0 )
    146              CALL exchange_horiz( ql_c,  0, 0 )
    147              CALL exchange_horiz( ql_v,  0, 0 )
    148              CALL exchange_horiz( ql_vp, 0, 0 )
     147             CALL exchange_horiz( ql )
     148             CALL exchange_horiz( ql_c )
     149             CALL exchange_horiz( ql_v )
     150             CALL exchange_horiz( ql_vp )
    149151          ENDIF
    150152
     
    182184             IF ( time_disturb >= dt_disturb )  THEN
    183185                IF ( hom(nzb+5,1,var_hom,0) < disturbance_energy_limit )  THEN
    184                    CALL disturb_field( nzb_u_inner, tend, u, uxrp, 0    )
    185                    CALL disturb_field( nzb_v_inner, tend, v ,   0, vynp )
     186                   CALL disturb_field( nzb_u_inner, tend, u )
     187                   CALL disturb_field( nzb_v_inner, tend, v )
    186188                ELSEIF ( bc_lr /= 'cyclic'  .OR.  bc_ns /= 'cyclic' )  THEN
    187189!
     
    189191!--                near the inflow throughout the whole simulation
    190192                   dist_range = 1
    191                    CALL disturb_field( nzb_u_inner, tend, u, uxrp, 0    )
    192                    CALL disturb_field( nzb_v_inner, tend, v ,   0, vynp )
     193                   CALL disturb_field( nzb_u_inner, tend, u )
     194                   CALL disturb_field( nzb_v_inner, tend, v )
    193195                   dist_range = 0
    194196                ENDIF
     
    208210!--       In case of a non-cyclic lateral wall, set the boundary conditions for
    209211!--       the velocities at the outflow
    210           IF ( bc_lr /= 'cyclic'  .OR.  bc_ns /= 'cyclic' )  THEN
    211              CALL boundary_conds( 'outflow_uvw' )
    212           ENDIF
     212!          IF ( bc_lr /= 'cyclic'  .OR.  bc_ns /= 'cyclic' )  THEN
     213!             CALL boundary_conds( 'outflow_uvw' )
     214!          ENDIF
    213215
    214216!
    215217!--       If required, compute virtuell potential temperature
    216           IF ( moisture ) CALL compute_vpt
     218          IF ( humidity ) CALL compute_vpt
    217219
    218220!
     
    235237!--          Compute the diffusion coefficients
    236238             CALL cpu_log( log_point(17), 'diffusivities', 'start' )
    237              IF ( .NOT. moisture ) THEN
     239             IF ( .NOT. humidity ) THEN
    238240                CALL diffusivities( pt )
    239241             ELSE
  • palm/trunk/SOURCE/wall_fluxes.f90

    r60 r75  
    3737! Call for all grid points
    3838!------------------------------------------------------------------------------!
    39     SUBROUTINE wall_fluxes( wall_flux, a, b, c1, c2, ixp, jyp, nzb_uvw_inner, &
     39    SUBROUTINE wall_fluxes( wall_flux, a, b, c1, c2, nzb_uvw_inner, &
    4040                            nzb_uvw_outer, wall )
    4141
     
    4848       IMPLICIT NONE
    4949
    50        INTEGER ::  i, ixp, j, jyp, k, wall_index
     50       INTEGER ::  i, j, k, wall_index
    5151
    5252       INTEGER, DIMENSION(nys-1:nyn+1,nxl-1:nxr+1) ::  nzb_uvw_inner, &
     
    5555       REAL ::  pts, pt_i, rifs, u_i, v_i, us_wall, vel_total, ws, wspts
    5656
    57        REAL, DIMENSION(nys-1:nyn+1,nxl-1:nxr+1)           ::  wall
    58        REAL, DIMENSION(nzb:nzt+1,nys:nyn+jyp,nxl:nxr+ixp) ::  wall_flux
     57       REAL, DIMENSION(nys-1:nyn+1,nxl-1:nxr+1)   ::  wall
     58       REAL, DIMENSION(nzb:nzt+1,nys:nyn,nxl:nxr) ::  wall_flux
    5959
    6060
     
    6363       wall_index = NINT( a+ 2*b + 3*c1 + 4*c2 )
    6464
    65        DO  i = nxl, nxr+ixp
    66           DO  j = nys, nyn+jyp
     65       DO  i = nxl, nxr
     66          DO  j = nys, nyn
    6767
    6868             IF ( wall(j,i) /= 0.0 )  THEN
  • palm/trunk/SOURCE/write_3d_binary.f90

    r73 r75  
    55! -----------------
    66! +precipitation_amount, precipitation_rate_av, rif_wall, u_m_l, u_m_r, etc.,
    7 ! z0_av
     7! z0_av, moisture renamed humidity
    88!
    99! Former revisions:
     
    107107    ENDIF
    108108    WRITE ( 14 )  'pt_m                ';  WRITE ( 14 )  pt_m
    109     IF ( moisture  .OR. passive_scalar )  THEN
     109    IF ( humidity  .OR. passive_scalar )  THEN
    110110       WRITE ( 14 )  'q                   ';  WRITE ( 14 )  q
    111111       IF ( ALLOCATED( q_av ) )  THEN
     
    209209       WRITE ( 14 )  'v_m_s               ';  WRITE ( 14 )  v_m_s
    210210    ENDIF
    211     IF ( moisture )  THEN
     211    IF ( humidity )  THEN
    212212       WRITE ( 14 )  'vpt                 ';  WRITE ( 14 )  vpt
    213213       IF ( ALLOCATED( vpt_av ) )  THEN
  • palm/trunk/SOURCE/write_var_list.f90

    r63 r75  
    44! Actual revisions:
    55! -----------------
    6 ! +loop_optimization, pt_refrence
     6! +loop_optimization, pt_refrence, moisture renamed humidity
    77!
    88! Former revisions:
     
    178178    WRITE ( 14 )  'mixing_length_1d              '
    179179    WRITE ( 14 )  mixing_length_1d
    180     WRITE ( 14 )  'moisture                      '
    181     WRITE ( 14 )  moisture
     180    WRITE ( 14 )  'humidity                      '
     181    WRITE ( 14 )  humidity
    182182    WRITE ( 14 )  'momentum_advec                '
    183183    WRITE ( 14 )  momentum_advec
Note: See TracChangeset for help on using the changeset viewer.