Hi, 

I am trying to get a plot of Star formation rate in yt for RAMSES Data set and I am using the sample code given for Enzo instead 

import yt
import numpy as np
from yt.data_objects.particle_filters import add_particle_filter
from matplotlib import pyplot as plt
def formed_star(pfilter, data):
    filter = data["io", "creation_time"] > 0
    return filter
add_particle_filter("formed_star", function=formed_star, filtered_type='io',
                    requires=["creation_time"])
filename = "/Users/xxxxxx/Desktop/out/info_00014.txt"
ds = yt.load(filename)
ds.add_particle_filter('formed_star')
ad = ds.all_data()
masses = ad['formed_star', 'particle_mass'].in_units('Msun')
formation_time = ad['formed_star', 'creation_time'].in_units('yr')
time_range = [0, 5e9] # years
n_bins = 1000
hist, bins = np.histogram(formation_time, bins=n_bins, range=time_range,)
inds = np.digitize(formation_time, bins=bins)

time = (bins[:-1] + bins[1:])/2
sfr = np.array([masses[inds == j+1].sum()/(bins[j+1]-bins[j])
                for j in range(len(time))])
sfr[sfr == 0] = np.nan
plt.plot(time/1e6, sfr)
plt.xlabel('Time  [Myr]')

However, I am ending up with this error:

The fields

('io', 'creation_time'),

required by the "formed_star" particle filter, are not defined for this dataset. 

 

Could anyone please suggest how I could access this "formed_star"/ an equivalent field for RAMSES. Thanks!