[AMBER] Ptraj and new netcdf Restart files

From: Joshua Adelman <jla65.pitt.edu>
Date: Thu, 31 May 2012 13:55:05 -0400

I am starting to use the new netcdf restart files available in Amber 12, and am running into a slight problem when using ptraj. In my ptraj file, I have:

trajin parent.restrt
trajin seg.mdcrd

To read in the coordinates from the restart file and then the trajectory data. This worked fine when I was using the ASCII format. However now using the netcdf format (ntxo=2), I get the following errors in my ptraj log:
--------
PTRAJ: trajin parent.restrt
  Checking coordinates: parent.restrt
NETCDF file:
Yikes! netcdfGetDimensionInfo on frame, error: NetCDF: Invalid dimension ID or name

. . .

PTRAJ: trajout pcoord.nc netcdf nobox
  parent.restrt: 379119632 frames.
  seg.mdcrd: 1 frames.

. . .

Error detected upon reading NetCDF coordinates on set 1: NetCDF: Start+count exceeds dimension bound 100%
--------

I would assume that the restart file should have only a single frame (The seg.mdcrd is also expected to have only a single frame).

When I open the parent.restrt and segment.mdcrd files using the netcdf4-python module in python, I see that the 'coordinates' array has different shapes for each. For the seg.mdcrd:

In [4]: f.variables['coordinates'][:].shape
Out[4]: (1, 70366, 3)

For parent.restrt:

In [8]: f.variables['coordinates'][:].shape
Out[8]: (70366, 3)

I'm not sure if this makes a difference or is a possible problem. I'm a fairly new Amber user and am not sure if this is potentially a bug or if I'm doing something incorrectly. Suggestions would be most appreciated.

Thanks,
Josh




_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Thu May 31 2012 - 11:00:02 PDT
Custom Search