resetting initial conditions once a year

[tasha]

Hi all,
What is the most straightforward way to reset certain tracers back to their initial state (from the *ini*.nc file) once per year?

I tried copying and modifying get_state and calling it from nl_ocean.h. It calls get_state okay, but the reset values don't show up in ocean_his.nc. Also the program crashes now after running a little while (Blowing up, usually @ Tindex = 31 in debug mode and segmentation fault in parallel mode)

Thanks for your help!
Tasha

[kate]

Rather than modifying the ROMS code, I would create a new "initial" file based on the current state of the restart file with the tracers copied over.

[tasha]

Thank you for the suggestion, Kate. I could edit the new "initial" file and run the model again, but I was hoping to run the model for 15 yrs in a row, which would not be ideal to edit the file again and again.

I ended up getting my code modification to work, I just had to fix one little bug.

Thanks, Tasha