Forum

  1. Home
  2. 3. Archive
  3. DELWAQ
  4. RE: DynDepth and z-layers

RE: DynDepth and z-layers

Rusty Holleman, modified 9 Years ago.

DynDepth and z-layers

Youngling Post: 1 Join Date: 8/3/15 Recent Posts
I am coupling unstructured z-layer hydrodynamics to DelWaq. In an effort to save disk space, I have DelWaq calculate the depth of segments based on surface area and volume via the DynDepth process. When a surface segment dries out, though, the depth is left as the last depth before it dried out. This depth is still included in the totdepth process, such that totaldepth and localdepth are overestimated. I believe that a fix is to add lines 7-8 below to ddepth.f:

1

2
...
3
DEPTH = VOLUME / SURF
4
!
5
      PMSA (IP3 ) = DEPTH
6
!
7
      ELSE ! This might fix for dry z-levels
8
         PMSA (IP3 ) = 0 ! dry cells get 0 depth
9
      ENDIF
10
!
11
      IFLUX = IFLUX + NOFLUX
12
      IP1   = IP1   + INCREM (  1 )
13
...


Am I missing some bigger problem, or is there any problem with this "fix"? I'm using r5740 - any chance this has been addressed in a more recent version?

thank you--
Michelle Jeuken, modified 6 Years ago.

RE: DynDepth and z-layers

Padawan Posts: 28 Join Date: 1/21/13 Recent Posts
<p>Hi Rusty, You are right about this! It seems that totaldepth doesn&#39;t skip inactive calls( any more?), but dyndepth does. This way old values remain in memory. Your solution works fine, but I think it is better to remove the check on the cell being active, so depth is always calculated, even if the volume/depth is below the threshold: IF (BTEST(IKNMRK(ISEG),0)) THEN I have committed this to the trunk as revision 6071 Thank you for posting this problem!</p>
tyron woodley, modified 5 Years ago.

RE: DynDepth and z-layers