Hello NCAR-users,
I think there is a bug in the contour tools when the special value is 0.
In this case, contour lines are generated in the special region.
Is there an error in the interpolation routine ?
I let you below a small program which shows it.
Run it with spval = 0. and then with any other value and see the graphic output.
Do you know any solution to avoid this problem ?
Thanks for your help.
Jean-Francois Foccroulle
Institut d'Astronomie et de Geophysique
2, chemin du cyclotron
1348 Louvain-la-Neuve
PROGRAM EXMPL9
IMPLICIT NONE
REAL spval,ZDAT(40,40)
C
WRITE(6,*)'spval ?'
READ(5,*)spval
C Generate dummy data.
CALL GENDAT (ZDAT,40,40,spval)
C Open GKS.
CALL OPNGKS
CALL CPSETR ('SPV - OUT-OF-RANGE VALUE',spval)
CALL CPCNRC (ZDAT,40,40,40,0.,0.,-11.,3,0,-682)
CALL FRAME
CALL CLSGKS
STOP
END
SUBROUTINE GENDAT (DATA,M,N,spval)
IMPLICIT NONE
INTEGER I,J,M,N
REAL DATA(M,N),spval
C This is a routine to generate test data for two-dimensional graphics
C routines.
DO I=1,M
DO J=1,N
DATA(I,J) = (i-1)*j
ENDDO
ENDDO
DO I=1,M
DATA(I,N) = spval
ENDDO
RETURN
END
This archive was generated by hypermail 2b29 : Wed Jun 28 2000 - 09:45:31 MDT