Changeset 4519 for palm


Ignore:
Timestamp:
May 5, 2020 5:33:30 PM (7 months ago)
Author:
suehring
Message:

Add missing computation of passive scalar scaling parameter

File:
1 edited

Legend:

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

    r4370 r4519  
    2626! -----------------
    2727! $Id$
     28! Add missing computation of passive scalar scaling parameter
     29!
     30! 4370 2020-01-10 14:00:44Z raasch
    2831! bugfix: openacc porting for vector version of OL calculation added
    2932!
     
    15861589                surf%ss(m) = -surf%ssws(m) / ( surf%us(m) + 1E-30_wp )
    15871590             ENDDO
     1591          ELSEIF ( .NOT. surf_vertical )  THEN
     1592
     1593             !$OMP PARALLEL DO PRIVATE( i, j, k, z_mo )
     1594             DO  m = 1, surf%ns   
     1595                i    = surf%i(m)
     1596                j    = surf%j(m)
     1597                k    = surf%k(m)
     1598                z_mo = surf%z_mo(m)
     1599
     1600                surf%ss(m) = kappa * ( s(k,j,i) - s(k-1,j,i) )              &
     1601                                  / ( LOG( z_mo / surf%z0h(m) )             &
     1602                                      - psi_h( z_mo / surf%ol(m) )          &
     1603                                      + psi_h( surf%z0h(m) / surf%ol(m) ) )
     1604
     1605             ENDDO
    15881606          ENDIF
    15891607!
Note: See TracChangeset for help on using the changeset viewer.