[MITgcm-support] Exf package and interpolation in time

[Georges Djoumna]

Hi all,

I am still facing problems with time interpolation in EXF package.

I've done a test with 3-hourly observed data for the wind. I've  done  
a run for 3 days with
ustressperiod     = 259200.0,
Deltat=150,
( I do interpolation in space)

In exf_set_gen.F  it is written
Interpolate linearly onto the  time.
              genfld(i,j,bi,bj) = exf_inscal_gen * (
                                fac * gen0(i,j,bi,bj) +
                    (exf_one - fac) * gen1(i,j,bi,bj) )

so I print out gen0, gen1 and genfld.
gen0 is always the observed data at t=0, and gen1 is always zero  
which are not correct.

My question is how to fill in  gen1 with correct values?

I've checked the binary file for uwindfile and it looks correct.

Suggestions are all welcome.

Thanks

Georges Djoumna 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20100118/593f3db5/attachment.htm>