Changeset 2128


Ignore:
Timestamp:
Jan 23, 2017 3:00:03 PM (5 years ago)
Author:
suehring
Message:

Bugfix in topography in case of ocean simulations

File:
1 edited

Legend:

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

    r2101 r2128  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! Bugfix in setting topography from file in case of ocean simulations
    2323!
    2424! Former revisions:
     
    750750          DO  i = nxl, nxr
    751751             DO  j = nys, nyn
    752                 nzb_local(j,i) = MINLOC( ABS( zw - topo_height(j,i) ), 1 ) - 1
    753                 IF ( ABS( zw(nzb_local(j,i)  ) - topo_height(j,i) ) == &
    754                      ABS( zw(nzb_local(j,i)+1) - topo_height(j,i) )    )  &
    755                    nzb_local(j,i) = nzb_local(j,i) + 1
     752                IF ( .NOT. ocean )  THEN
     753                   nzb_local(j,i) = MINLOC( ABS( zw - topo_height(j,i) ), 1 ) - 1
     754                   IF ( ABS( zw(nzb_local(j,i)  ) - topo_height(j,i) ) ==      &
     755                        ABS( zw(nzb_local(j,i)+1) - topo_height(j,i) )    )    &
     756                      nzb_local(j,i) = nzb_local(j,i) + 1
     757                ELSE
     758                   nzb_local(j,i) = MINLOC( ABS( zw - zw(0)                    &
     759                                                    - topo_height(j,i) ), 1 ) - 1
     760                   IF ( ABS( zw(nzb_local(j,i)  ) - zw(0)                      &
     761                                                  - topo_height(j,i) )  ==     &
     762                        ABS( zw(nzb_local(j,i)+1) - zw(0)                      &
     763                                                  - topo_height(j,i) )    )    &
     764                      nzb_local(j,i) = nzb_local(j,i) + 1
     765                ENDIF
     766
    756767             ENDDO
    757768          ENDDO
Note: See TracChangeset for help on using the changeset viewer.