Changeset 2369 for palm/trunk


Ignore:
Timestamp:
Aug 22, 2017 3:20:37 PM (7 years ago)
Author:
suehring
Message:

Bugfix, no masking of topography in eqn_seawater

File:
1 edited

Legend:

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

    r2233 r2369  
    2525! -----------------
    2626! $Id$
     27! Bugfix, do not mask topography here, since density becomes zero, leading to
     28! division by zero in production_e
     29!
     30! 2233 2017-05-30 18:08:54Z suehring
    2731!
    2832! 2232 2017-05-30 17:47:52Z suehring
     
    123127           ONLY:  hyp, prho, pt_p, rho_ocean, sa_p
    124128       USE indices,                                                            &
    125            ONLY:  nxl, nxr, nyn, nys, nzb, nzt, wall_flags_0
     129           ONLY:  nxl, nxr, nyn, nys, nzb, nzt
    126130
    127131       USE surface_mod,                                                        &
     
    138142       INTEGER(iwp) ::  surf_s  !< Start index of surface elements at (j,i)-gridpoint
    139143
    140        REAL(wp) ::  flag   !< flag to mask topography grid points
    141144       REAL(wp) ::  pden   !<
    142145       REAL(wp) ::  pnom   !<
     
    183186                       den(10)*sa15*pt2
    184187!
    185 !--             Predetermine flag to mask topography
    186                 flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_0(k,j,i), 0 ) )
    187 !
    188188!--             Potential density (without pressure terms)
    189                 prho(k,j,i) = pnom / pden * flag
     189                prho(k,j,i) = pnom / pden
    190190
    191191                pnom = pnom +             nom(8)*p1      + nom(9)*p1*pt2  +    &
     
    197197!
    198198!--             In-situ density
    199                 rho_ocean(k,j,i) = pnom / pden * flag
     199                rho_ocean(k,j,i) = pnom / pden
    200200
    201201             ENDDO
     
    242242           
    243243       USE indices,                                                            &
    244            ONLY:  nzb, nzt, wall_flags_0
     244           ONLY:  nzb, nzt
    245245
    246246       USE surface_mod,                                                        &
     
    257257       INTEGER(iwp) ::  surf_s  !< Start index of surface elements at (j,i)-gridpoint
    258258
    259        REAL(wp) ::  flag   !< flag to mask topography grid points
    260259       REAL(wp) ::  pden   !<
    261260       REAL(wp) ::  pnom   !<
     
    298297                 den(10)*sa15*pt2
    299298!
    300 !--       Predetermine flag to mask topography
    301           flag = MERGE( 1.0_wp, 0.0_wp, BTEST( wall_flags_0(k,j,i), 0 ) )
    302 !
    303299!--       Potential density (without pressure terms)
    304           prho(k,j,i) = pnom / pden * flag
     300          prho(k,j,i) = pnom / pden
    305301
    306302          pnom = pnom +             nom(8)*p1      + nom(9)*p1*pt2  +          &
     
    311307!
    312308!--       In-situ density
    313           rho_ocean(k,j,i) = pnom / pden * flag
     309          rho_ocean(k,j,i) = pnom / pden
    314310
    315311
Note: See TracChangeset for help on using the changeset viewer.