Ignore:
Timestamp:
Aug 24, 2007 3:10:38 PM (14 years ago)
Author:
letzel
Message:
  • Improved coupler: evaporation - salinity-flux coupling for humidity = .T.,

avoid MPI hangs when coupled runs terminate, add DOC/app/chapter_3.8;

  • Optional calculation of km and kh from initial TKE e_init;
  • Default initialization of km,kh = 0.00001 for ocean = .T.;
  • Allow data_output_pr= q, wq, w"q", w*q* for humidity = .T.;
  • Bugfix: Rayleigh damping for ocean fixed.
File:
1 edited

Legend:

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

    r106 r108  
    88! u_0 and v_0 are calculated for nxr+1, nyn+1 also (otherwise these values are
    99! not available in case of non-cyclic boundary conditions)
     10! Bugfix for ocean density flux at bottom
    1011!
    1112! Former revisions:
     
    810811!--             So far in the ocean no special treatment of density flux in the
    811812!--             bottom and top surface layer
    812                 DO  k = nzb_s_inner(j,i), nzt
     813                k= nzb_s_inner(j,i)
     814                tend(k,j,i) = tend(k,j,i) + kh(k,j,i) * g / prho_reference * &
     815                     ( rho(k+1,j,i) - rho(k,j,i) ) * dzu(k+1)
     816                DO  k = nzb_s_inner(j,i)+1, nzt
    813817                   tend(k,j,i) = tend(k,j,i) + kh(k,j,i) * g / prho_reference * &
    814818                                      ( rho(k+1,j,i) - rho(k-1,j,i) ) * dd2zu(k)
Note: See TracChangeset for help on using the changeset viewer.