Updated rhosigchi.f with a warning if the x and y dims are too big

Author: cecgresonant

prog/rhosigchi.f

@@ -187,6 +187,14 @@ PROGRAM RhoSigChi
           Fi(i)=0.
           Ff(i)=0.
         ENDDO
+        IF(X_dim.GT.512)THEN
+           X_dim=512
+           WRITE(6,*)'WARNING, more than 512 channels on x-axis'
+        ENDIF
+        IF(Y_dim.GT.512)THEN
+           Y_dim=512
+           WRITE(6,*)'WARNING, more than 512 channels on y-axis'
+        ENDIF
         DO j=0,Y_dim-1
           DO i=0,X_dim-1
             IF(rmatrix(i,j).LT.eps) rmatrix(i,j) = eps
@@ -234,6 +242,14 @@ PROGRAM RhoSigChi
           Fi(i)=0.
           Ff(i)=0.
         ENDDO
+        IF(X_dim.GT.512)THEN
+           X_dim=512
+           WRITE(6,*)'WARNING, more than 512 channels on x-axis'
+        ENDIF
+        IF(Y_dim.GT.512)THEN
+           Y_dim=512
+           WRITE(6,*)'WARNING, more than 512 channels on y-axis'
+        ENDIF        
         DO j=0,Y_dim-1
           DO i=0,X_dim-1
             IF(rmatrix(i,j).LT.eps) rmatrix(i,j) = eps
@@ -281,6 +297,14 @@ PROGRAM RhoSigChi
           Fi(i)=0.
           Ff(i)=0.
         ENDDO
+        IF(X_dim.GT.512)THEN
+           X_dim=512
+           WRITE(6,*)'WARNING, more than 512 channels on x-axis'
+        ENDIF
+        IF(Y_dim.GT.512)THEN
+           Y_dim=512
+           WRITE(6,*)'WARNING, more than 512 channels on y-axis'
+        ENDIF        
         DO j=0,Y_dim-1
           DO i=0,X_dim-1
             IF(rmatrix(i,j).LT.eps) rmatrix(i,j) = eps
@@ -1086,6 +1110,7 @@ SUBROUTINE R2Dim_matrix
       device=20
       OPEN(device,FILE=filex,ACCESS='SEQUENTIAL',status='old',ERR=901)
 
+
 C Read all 10 lines of the file header
       DO i=1,10
         READ(device,100,ERR=800)text