To whom it may concern:
I am using the plot package provided by roms in fortran and I would like to know if there is something missing in the csec file, since one of the options is
1 ISECPOS: section positions: [0] grid units, [1] Cartesian, [2] spherical
but there is no option in the code to read isecpos=1, only 0 by default which you should enter an integer, and isecpos=2, spherical coordinates. Please see the section of the code below. Is it that I am using an old version of the plot codes or there is something missing?
i would appreciate your help.
thanks
Rafael
READ (stdinp,*) isecpos
!
! Read in the coordinates of the two points defining the section
! (grid units, or kilometers, or degrees). Convert longitudes to
! EAST-longitudes, if appropriate.
!
READ (stdinp,*) xsec1
READ (stdinp,*) ysec1
READ (stdinp,*) xsec2
READ (stdinp,*) ysec2
IF (isecpos.eq.2) THEN
IF (((xsec1.le.c0).and.(xsec2.ge.c0)).or. &
& ((xsec1.le.c0).and.(xsec2.le.c0)).or. &
& ((xsec1.gt.c0).and.(xsec2.gt.c0))) THEN
eastlon=.false.
ELSE
eastlon=.true.
IF (xsec1.lt.c0) xsec1=c360+xsec1
IF (xsec2.lt.c0) xsec2=c360+xsec2
END IF
END IF