[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