I've compiled successfully the roms plot package and there are four programs, including ccnt, cnt, csec and sec in the PLT_BINDIR.
When I used this package to plot the ROMS output, it reported that:
Code: Select all
[cloud@igplogin ROMS_output]$ ./csec < csec.in
At line 576 of file csec.f90 (unit = 5, file = 'stdin')
Fortran runtime error: Bad integer for item 1 in list input
Here is my csec.in file:READ (stdinp,*) nfields
Code: Select all
2016 -1 : year and starting year-day (use yearday<0, for no time label)
ROMS 3.0
!Wind-Driven Upwelling/Downwelling over a Periodic Channel
ROMS output
Time-averaged fields
Potential Temperature (C), salinity, u, v, w, omega
1 NFIELDS: number of fields to plot. Line below, field(s) types:
122 field identification: FLDID(1:NFIELDS).
0 FRSTD : first day (>0) or time record (<0) to plot
0 LASTD : last day (>0) or time record (<0) to plot
0 DSKIP : plot every other DSKIP days (see below)
1 FINDX : forecast record to process, if any (atmospheric files only).
2 ISECPOS: section positions: [0] grid units, [1] Cartesian, [2] spherical
21 XSEC1 : starting section X-position (grid units, or km, or degrees).
110 YSEC1 : starting section Y-position (grid units, or km, or degrees).
21 XSEC2 : ending section X-position (grid units, or km, or degrees).
160 YSEC2 : ending section Y-position (grid units, or km, or degrees).
1.0 DELS : section horizontal grid spacing (km or degrees).
0 SZTOP : shallowest section depth to consider (meters).
150 SZBOT : deepest section depth to consider (meters).
16 KMS : number of refiment section levels between SZTOP and SZBOT.
0 VINTRP : vertical interpolation: [0] linear, [1] cubic splines.
250.0 PMIN : field minimum value for color palette (0.0 for default).
300.0 PMAX : field maximum value for color palette (0.0 for default).
1 ICNT : draw contours between color bands: 0=no, 1=yes.
0 IREF : secondary or reference field option (see below).
26 IDOVER : overlay field identification (for IREF=1,2 only).
0.0 RMIN : overlay field minimum value to consider (0.0 for default).
0.0 RMAX : overlay field maximum value to consider (0.0 for default).
1 NPAGE : number of plots per page (currently 1, 2, or 4).
1 IFILL : mask area below water column: [0] no, [1] yes.
F READGRD: logical switch to read in positions from grid NetCDF file.
F PLTLOGO: logical switch Logo.
T WRTHDR : logical switch to write out the plot header titles.
T WRTBLAB: logical switch to write out the plot bottom title.
T T WRTRANG: logical switch to write out data/overlay range values and CI.
F WRTFNAM: logical switch to write out input primary filename.
F WRTDATE: logical switch to write out current date.
/work/users/cloud/user/manh/anaconda3/plot/Data/varid.dat
/work/users/cloud/user/manh/anaconda3/plot/Palettes/jet120.pal
/work/users/cloud/user/manh/anaconda3/plot/Data/default.cnt
cean_avg1_00001.nc
ocean_hist1_00001.nc
Thank all,