search

MODFLOW-USGS / mfusg

Public repository for MODFLOW-USG
https://water.usgs.gov/ogw/mfusg/
7 stars 4 forks source link

Bug with CLN CONST HEAD cell-by-cell flow output #32

Closed langevin-usgs closed 7 years ago

langevin-usgs commented 7 years ago

Hi Chris,

One more CLN bug (in output) is below. Need to write the CLN node number and not the global node number. I have changed this in the Beta version for writing to output listing file as well as to the binary CBB file.

Sorab

From: Bill Dolinar [mailto:bdolinar@aquaveo.com] Sent: Tuesday, March 14, 2017 12:02 PM To: Sorab Panday sp@gsi-net.com Subject: Bug with CLN CONST HEAD cell-by-cell flow output

Hi Sorab,

I ran across a bug in the CCF output for CLN CONST HEAD. The CLN node ID should be written but the value getting written is the global node ID (CLN ID + NODES). The section of code from GWF2CLNU1BDCHWR below shows the fix I've made in our code. The first change for the LST file and the second change for the CCF file. I'm not sure about the fix for the LST file, should it show the global node ID or the CLN ID?

Thanks, Bill

C6--------PRINT THE FLOW FOR THE CELL IF REQUESTED.
        IF(IBD.LT.0) THEN
          IF(IBDLBL.EQ.0) WRITE(IOUT,899) TEXT(1),KPER,KSTP
  899     FORMAT(1X,/1X,A,'   PERIOD',I3,'   STEP',I3)
C         WRITE(IOUT,910) N,RATE                                        !aq
          WRITE(IOUT,910) N-NODES,RATE                                  !aq
  910     FORMAT(1X,'NODE',I8,'   RATE',1PG15.6)
          IBDLBL=1
        END IF
C
C7------IF SAVING CELL-BY-CELL FLOW IN LIST, WRITE FLOW FOR CELL.
        IF(IBD.EQ.2)THEN
          SRATE = RATE
C         CALL UBDSVAU(ICLNCB,NODES,N,SRATE,IBOUND)                     !aq
          CALL UBDSVAU(ICLNCB,NODES,N-NODES,SRATE,IBOUND)               !aq
C--------------------------------------------------------------
        ENDIF
langevin-usgs commented 7 years ago

Fixed.