FESOM / fesom2

Multi-resolution ocean general circulation model.
http://fesom.de/
GNU General Public License v3.0
47 stars 48 forks source link

Refactoring improve4ogurses calendar #591

Closed patrickscholz closed 2 months ago

patrickscholz commented 4 months ago

Solve special problem here: Ozgurs wants to repeat only a part of the JRA55 forcing periodically (e.g repeat the period 1958-1967) through re-linking of the original forcing files.

JanStreffing commented 4 months ago

I'm not sure the visual inspection is enough for this. Would it be possible to run a 65 year cycle with this turned on and off? If results are bit-identical, I think it's good to go.

Or does the Draft status mean more work needs to be done?

patrickscholz commented 4 months ago

@JanStreffing Im still running some tests to check if everything fits together!

ogurses commented 4 months ago

I see that we turn leapyear off in FESOM config. Does this code patch consider the exclusion of 29th Feb. from the forcing data if the forcing file is on leapyear but FESOM clock not. In other words, do I have to drop 29th Feb. from the forcing files of leap year before I make any linking operation?

patrickscholz commented 3 months ago

@ogurses you were right in the previous solution the 29 Feb. was not fully skipped in the cache of the temporal interpolation. I changed that now! So when include_fleapyear=.false., the model checks if it finds a forcing file whos original time-axes represent a leapyear and than jumps over the 29. Feb when selecting the indices for the time interpolation cache. So now you should be able to truly only relink the original JRA55 forcing files to create your periodically repeating forcing. So no need to fully dublicate the forcing and delete leapyear 29. Feb from hand.

patrickscholz commented 3 months ago

There still seems to be some bias between this and the original refactoring branch, that do not seem to happen in the test case. Since the testcase seems to work with CORE2 data which do not have leap years.