Changeset 124 for palm/trunk/SOURCE


Ignore:
Timestamp:
Oct 19, 2007 3:47:46 PM (17 years ago)
Author:
raasch
Message:

further bugfixes for the ocean part

Location:
palm/trunk/SOURCE
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/CURRENT_MODIFICATIONS

    r121 r124  
    2525------
    2626
     27Several bugfixes in the ocean part: Initial density rho is calculated
     28(init_ocean). Error in initializing u_init and v_init removed
     29(check_parameters). Calculation of density flux now starts from
     30nzb+1 (production_e).
     31
    2732Bugfix: pleft/pright changed to pnorth/psouth in sendrecv of particle tail
    2833numbers along y, small bugfixes in the SGS part (advec_particles)
     
    3035Bugfix: model_string needed a default value (combine_plot_fields)
    3136
    32 Bugfix: initializing u_init and v_init in case of ocean runs (check_parameters)
    3337
    34 advec_particles, check_parameters, combine_plot_fields
     38advec_particles, check_parameters, combine_plot_fields, init_ocean, production_e
  • palm/trunk/SOURCE/init_ocean.f90

    r98 r124  
    44! Actual revisions:
    55! -----------------
    6 !
     6! Bugfix: Initial density rho is calculated
    77!
    88! Former revisions:
     
    7373    prho_reference = prho_reference / ( zu(nzt) - zu(nzb) )
    7474
     75!
     76!-- Calculate the initial potential density, based on the initial
     77!-- temperature and salinity profile
     78    CALL eqn_state_seawater
     79
    7580
    7681 END SUBROUTINE init_ocean
  • palm/trunk/SOURCE/production_e.f90

    r110 r124  
    44! Actual revisions:
    55! -----------------
    6 !
     6! Bugfix: calculation of density flux in the ocean now starts from nzb+1
    77!
    88! Former revisions:
     
    382382!--                the bottom and top surface layer
    383383                   DO  j = nys, nyn
    384                       DO  k = nzb_s_inner(j,i), nzt
     384                      DO  k = nzb_s_inner(j,i)+1, nzt
    385385                         tend(k,j,i) = tend(k,j,i) +                    &
    386386                                       kh(k,j,i) * g / prho_reference * &
     
    419419!--                the bottom and top surface layer
    420420                   DO  j = nys, nyn
    421                       DO  k = nzb_s_inner(j,i), nzt
     421                      DO  k = nzb_s_inner(j,i)+1, nzt
    422422                         tend(k,j,i) = tend(k,j,i) -                &
    423423                                       kh(k,j,i) * g / rho(k,j,i) * &
     
    814814!--             So far in the ocean no special treatment of density flux in the
    815815!--             bottom and top surface layer
    816                 k= nzb_s_inner(j,i)
    817                 tend(k,j,i) = tend(k,j,i) + kh(k,j,i) * g / prho_reference * &
    818                      ( rho(k+1,j,i) - rho(k,j,i) ) * dzu(k+1)
    819816                DO  k = nzb_s_inner(j,i)+1, nzt
    820817                   tend(k,j,i) = tend(k,j,i) + kh(k,j,i) * g / prho_reference * &
     
    847844!--             So far in the ocean no special treatment of density flux in the
    848845!--             bottom and top surface layer
    849                 DO  k = nzb_s_inner(j,i), nzt
     846                DO  k = nzb_s_inner(j,i)+1, nzt
    850847                   tend(k,j,i) = tend(k,j,i) + kh(k,j,i) * g / rho(k,j,i) * &
    851848                                      ( rho(k+1,j,i) - rho(k-1,j,i) ) * dd2zu(k)
Note: See TracChangeset for help on using the changeset viewer.