correct bug in mesh check that did not plot zero differences in bed levels correctly, add more plots

.gitignore

@@ -266,3 +266,5 @@ ModelManifest.xml
 /RMAGEN/RCa18976
 /src/NEWRMGN-demo.F90
 /src/RMAGENV83e-demo.rc
+/RMAGEN/TRANSFER/RMAGENV84
+/RMAGEN/GINPT.exe

src/AREA.F90

@@ -205,16 +205,22 @@
           IF(EDIF(N) .GT. EMAX) EMAX=EDIF(N)
     400   CONTINUE
         ENDDO
-        NUMV=11
-        CONTUR(1)=0.
-        DO K=2,11
-          CONTUR(K)=EMAX/10.+CONTUR(K-1)
+        NUMV=13
+        CONTUR(1)=-0.5
+        DO K=2,13
+          CONTUR(K)=(EMAX+0.5)/12.+CONTUR(K-1)
         ENDDO
         
         DO N=1,NE
+!          IF(N .EQ. 46451) WRITE(155,*) N,EMAX,EDIF(N)
           IF(IMAT(N) .GT. 0) THEN
-            ICOL=EDIF(N)*10./EMAX+.999
+            IF(EDIF(N) .LT. 0.001) THEN
+              ICOL=1
+            ELSE
+              ICOL=EDIF(N)*12./EMAX+.999
+            ENDIF
             icll=itran(icol)
+!            IF(N .EQ. 46451) WRITE(155,*) N,ICOL,EMAX,EDIF(N)
             CALL FILLEMC(N,ICLL)
           ENDIF
         ENDDO

src/RMAGEN.F90

@@ -284,6 +284,24 @@
             YMAX = -YMIN
          ENDIF
 
+!         do j=1,ne
+!           if(nop(j,1) .eq. 0) cycle
+!           do k=1,8,2
+!             if(nop(j,k) .eq. 0) cycle
+!             if(wd(nop(j,k)) .gt. 27.) go to 7110
+!             if(wd(nop(j,k)) .lt. 22.) go to 7111
+             
+!           enddo
+!7110       continue    
+!           imat(j)=0
+!           IESKP(J)=1
+!           DO K=1,8
+!             NOP(J,K)=0
+!           ENDDO
+!7111       continue    
+!         enddo
+         
+         
          IF(NP .GT. 0) THEN
             DO 10  J=1,NP
                IF (CORD(J,1) .LT. VDX)  GOTO 10