Code: Select all

```
+902 # ifdef Q_PSOURCE
+903 !
+904 ! Apply mass point sources - Volume influx.
+905 !
+906 DO is=1,Nsrc
+907 i=Isrc(is)
+908 IF (LtracerSrc(itrc,ng)) THEN
+909 IF (((IstrR.le.i).and.(i.le.IendR)).and. &
+910 & ((JstrR.le.j).and.(j.le.JendR))) THEN
+911 IF (j.eq.Jsrc(is)) THEN
+912 DO k=1,N(ng)
+913 FC(i,k)=FC(i,k)+0.5_r8* &
+914 & (Qsrc(is,k )*Tsrc(is,k ,itrc)+ &
+915 & Qsrc(is,k+1)*Tsrc(is,k+1,itrc))
+916 END DO
+917 END IF
+918 END IF
+919 END IF
+920 END DO
+921 # endif
```

Could someone please let me know how to remedy this bug - can we simply not carry-out the averaging on lines 914-915 in step3d_t.F?