Changeset 2088


Ignore:
Timestamp:
Dec 19, 2016 4:30:25 PM (7 years ago)
Author:
suehring
Message:

Bugfixes in initial salinity profile and generic topography definition in case of ocean simulations

Location:
palm/trunk/SOURCE
Files:
2 edited

Legend:

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

    r2085 r2088  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! Bugfix in initial salinity profile
    2323!
    2424! Former revisions:
     
    503503    LOGICAL     ::  found                            !<
    504504   
     505    REAL(wp)    ::  dum                              !<
    505506    REAL(wp)    ::  gradient                         !<
    506507    REAL(wp)    ::  remote = 0.0_wp                  !<
     
    14651466          CALL init_vertical_profiles( sa_vertical_gradient_level_ind,          &
    14661467                                       sa_vertical_gradient_level,              &
    1467                                        sa_vertical_gradient, s_init,            &
    1468                                        sa_surface, -999.0_wp )
     1468                                       sa_vertical_gradient, sa_init,           &
     1469                                       sa_surface, dum )
    14691470       ENDIF
    14701471
  • palm/trunk/SOURCE/init_grid.f90

    r2038 r2088  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! Bugfix in generic topography in case of ocean simulations
    2323!
    2424! Former revisions:
     
    574574          blx = NINT( building_length_x / dx )
    575575          bly = NINT( building_length_y / dy )
    576           bh  = MINLOC( ABS( zw - building_height ), 1 ) - 1
     576          IF ( .NOT. ocean )  THEN
     577             bh  = MINLOC( ABS( zw - building_height ), 1 ) - 1
     578          ELSE
     579             bh  = MINLOC( ABS( zw - zw(0) - building_height ), 1 ) - 1
     580          ENDIF
     581
    577582          IF ( ABS( zw(bh  ) - building_height ) == &
    578583               ABS( zw(bh+1) - building_height )    )  bh = bh + 1
     
    640645          ENDIF
    641646
    642           ch  = MINLOC( ABS( zw - canyon_height ), 1 ) - 1
     647          IF ( .NOT. ocean )  THEN
     648             ch  = MINLOC( ABS( zw - canyon_height ), 1 ) - 1
     649          ELSE
     650             ch  = MINLOC( ABS( zw - zw(0) - canyon_height ), 1 ) - 1
     651          ENDIF
     652
    643653          IF ( ABS( zw(ch  ) - canyon_height ) == &
    644654               ABS( zw(ch+1) - canyon_height )    )  ch = ch + 1
Note: See TracChangeset for help on using the changeset viewer.