resetting initial conditions once a year

General scientific issues regarding ROMS

Moderators: arango, robertson

Post Reply
Message
Author
tasha
Posts: 15
Joined: Wed Nov 18, 2009 7:27 pm
Location: New York University

resetting initial conditions once a year

#1 Unread post by 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

User avatar
kate
Posts: 4091
Joined: Wed Jul 02, 2003 5:29 pm
Location: CFOS/UAF, USA

Re: resetting initial conditions once a year

#2 Unread post by 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
Posts: 15
Joined: Wed Nov 18, 2009 7:27 pm
Location: New York University

Re: resetting initial conditions once a year

#3 Unread post by 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

Post Reply