Ignore:
Timestamp:
Dec 18, 2019 11:55:56 AM (23 months ago)
Author:
motisi
Message:

Introduction of wall_flags_total_0, which currently sets bits based on static topography information used in wall_flags_static_0

File:
1 edited

Legend:

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

    r4329 r4346  
    2525! -----------------
    2626! $Id$
     27! Introduction of wall_flags_total_0, which currently sets bits based on static
     28! topography information used in wall_flags_static_0
     29!
     30! 4329 2019-12-10 15:46:36Z motisi
    2731! Renamed wall_flags_0 to wall_flags_static_0
    2832!
     
    8993
    9094    USE indices,                                                               &
    91         ONLY:  advc_flags_s, nxl, nxr, nyn, nys, nzb, nzt, wall_flags_static_0
     95        ONLY:  advc_flags_s, nxl, nxr, nyn, nys, nzb, nzt, wall_flags_total_0
    9296
    9397    USE surface_mod,                                                           &
     
    971975    USE indices,                                                               &
    972976        ONLY: nxl, nxlg, nxr, nxrg, nyn, nyng, nys, nysg, nzb, nzt,            &
    973               wall_flags_static_0
     977              wall_flags_total_0
    974978
    975979    IMPLICIT NONE
     
    10391043             DO  k = nzb_do, nzt_do
    10401044                local_pf(i,j,k) = MERGE( to_be_resorted(k,j,i),                &
    1041                                          REAL( fill_value, KIND = wp ),        &
    1042                                          BTEST( wall_flags_static_0(k,j,i), flag_nr ) )
     1045                                   REAL( fill_value, KIND = wp ),              &
     1046                                   BTEST( wall_flags_total_0(k,j,i), flag_nr ) )
    10431047             ENDDO
    10441048          ENDDO
     
    10661070    USE indices,                                                               &
    10671071        ONLY:  nxl, nxlg, nxr, nxrg, nyn, nyng, nys, nysg, nzb, nzt,           &
    1068                wall_flags_static_0
     1072               wall_flags_total_0
    10691073
    10701074    IMPLICIT NONE
     
    11311135             DO  k = nzb_do, nzt_do
    11321136                local_pf(i,j,k) = MERGE( to_be_resorted(k,j,i),                &
    1133                                          REAL( fill_value, KIND = wp ),        &
    1134                                          BTEST( wall_flags_static_0(k,j,i), flag_nr ) )
     1137                                   REAL( fill_value, KIND = wp ),              &
     1138                                   BTEST( wall_flags_total_0(k,j,i), flag_nr ) )
    11351139             ENDDO
    11361140          ENDDO
     
    16761680          DO  j = nys, nyn
    16771681             DO  k = nzb+1, nzt
    1678                 sa_p(k,j,i) = sa(k,j,i) + ( dt_3d * ( sbt * tend(k,j,i) +      &
    1679                                                       tsc(3) * tsa_m(k,j,i) )  &
    1680                                                   - tsc(5) * rdf_sc(k) *       &
    1681                                                     ( sa(k,j,i) - sa_init(k) ) &
    1682                                           )                                    &
    1683                                             * MERGE( 1.0_wp, 0.0_wp,           &
    1684                                                BTEST( wall_flags_static_0(k,j,i), 0 ) &
     1682                sa_p(k,j,i) = sa(k,j,i) + ( dt_3d * ( sbt * tend(k,j,i) +            &
     1683                                                      tsc(3) * tsa_m(k,j,i) )        &
     1684                                                  - tsc(5) * rdf_sc(k) *             &
     1685                                                    ( sa(k,j,i) - sa_init(k) )       &
     1686                                          )                                          &
     1687                                            * MERGE( 1.0_wp, 0.0_wp,                 &
     1688                                               BTEST( wall_flags_total_0(k,j,i), 0 ) &
    16851689                                                   )
    16861690                IF ( sa_p(k,j,i) < 0.0_wp )  sa_p(k,j,i) = 0.1_wp * sa(k,j,i)
     
    18291833                                             * ( sa(k,j,i) - sa_init(k) )      &
    18301834                                    ) * MERGE( 1.0_wp, 0.0_wp,                 &
    1831                                                BTEST( wall_flags_static_0(k,j,i), 0 ) )
     1835                                         BTEST( wall_flags_total_0(k,j,i), 0 ) )
    18321836
    18331837          IF ( sa_p(k,j,i) < 0.0_wp )  sa_p(k,j,i) = 0.1_wp * sa(k,j,i)
Note: See TracChangeset for help on using the changeset viewer.