[MITgcm-support] External forcing for deltaT < 1 sec
Jody Klymak
jklymak at uvic.ca
Fri Apr 11 23:24:01 EDT 2014
On Apr 11, 2014, at 13:09 PM, Dustin Carroll <dcarroll at uoregon.edu> wrote:
> I tried using the basic periodicExternalForcing but that didn’t properly loop through my forcing files properly when deltaT < 1. Can anyone confirm that periodicExternalForcing works with deltaT < 1 or recommend a method for time dependent bc's with these small time steps?
Hi Dustin,
Not sure if I understand "didn't properly loop". If deltaT>1. does it "properly loop"? The source code you want to look at is get_periodic_interval.F, maybe that has a clue.
One thing to be careful about, which has got me in the past, is that your first record in your OBET_hr.bin etc is actually centered at 0.5*externForcingPeriod (=1800s), and that at t=0, the boundary condition is actually a linear interpolation between the records at -1800 s and at +1800 s i.e. between the *last* record in your file and the first. If you are doing something with a linear initial ramp, and the last value you have in your hourly forcing is large, your t=0 startup will be a lot bumpier than you were hoping.
The justification for that convention is that a lot of folks force with monthly climatologies etc, which are most appropriately centered on the period.
Cheers, Jody
More information about the MITgcm-support
mailing list