[MITgcm-support] Exf package and interpolation in time
Martin Losch
Martin.Losch at awi.de
Fri Jan 8 02:26:54 EST 2010
Hi George,
Hi Georges,
your hfluxperiod has to be 100hours (hfluxperiod = 360000.,), if you
have 100hourly data. The code reads new data every 100h and
interpolates in-between. You configuration makes the code believe that
you have data every 150s, so it reads it every 150s (your time step)
and then runs out of data.
Martin
On Jan 7, 2010, at 6:46 PM, Georges Djoumna wrote:
> Dear all,
> I'm trying to use the EXF package with interpolation in time, I
> have 100 hourly observed data at three point in the computational
> domain,
> I've done the interpolation in space ( I've created a binary file
> containing hourly observed data) and I want the code to do the
> linear interpolation in time. The time step for the simulation is
> 150s.
>
> In data.exf, I've set the startdate1, startdate2 and the period for
> each forcing . For example
> hfluxstartdate1 = 20050401,
> hfluxstartdate2 = 00000,
> hfluxperiod = 150.0,
>
> In data.cal I have
> #
> # *******************
> # Calendar Parameters
> # *******************
> &CAL_NML
> # TheCalendar='model',
> TheCalendar='gregorian',
> startDate_1=20050401,
> startDate_2=000000,
> &
>
>
> when the time step is close to the number of observation value, the
> code stop because the binary file contains 100 observed values.
>
> I set #defined verbose in EXF_OPTION and gets
> (PID.TID 0000.0001) exf_GetFFieldsRec:
> (PID.TID 0000.0001) MDS_READ_FIELD: opening global file:
> ZWindMet_Erie2002_IDWSpace_2Km.bin
> (PID.TID 0000.0001) exf_GetFFieldsRec:
> (PID.TID 0000.0001) exf_GetFFieldsRec: first, changed, fac: F
> T 0.10000000D+01
> (PID.TID 0000.0001) exf_GetFFieldsRec: myiter, count0, count1:
> 1 2 3
> (PID.TID 0000.0001) exf_GetFFieldsRec:
> (PID.TID 0000.0001) MDS_READ_FIELD: opening global file:
> MWindMet_Erie2002_IDWSpace_2Km.bin
> (PID.TID 0000.0001) exf_GetFFieldsRec:
> (PID.TID 0000.0001) exf_GetFFieldsRec: first, changed, fac: F
> T 0.10000000D+01
> (PID.TID 0000.0001) exf_GetFFieldsRec: myiter, count0, count1:
> 1 2 3
> (PID.TID 0000.0001) exf_GetFFieldsRec:
>
> I doing something wrong somewhere, but I don't know what exactly.
> Any help is welcome.
> Thanks so much.
>
> Georges Djoumna
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
More information about the MITgcm-support
mailing list