null

Message Boards

DynDepth and z-layers

RH
Rusty Holleman, modified 3 Years ago.

DynDepth and z-layers

Youngling Posts: 4 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:


...
DEPTH = VOLUME / SURF
!
      PMSA (IP3 ) = DEPTH
!
      ELSE ! This might fix for dry z-levels
         PMSA (IP3 ) = 0 ! dry cells get 0 depth
      ENDIF
!
      IFLUX = IFLUX + NOFLUX
      IP1   = IP1   + INCREM (  1 )
...


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--
Michel Jeuken, modified 1 Year ago.

RE: DynDepth and z-layers

Jedi Knight Posts: 156 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>
RH
Rusty Holleman, modified 3 Years ago.

RE: DynDepth and z-layers

Youngling Posts: 4 Join Date: 8/3/15 Recent Posts
Thanks, Michel - I'll switch to your fix.

best,
Rusty