Changes between Version 15 and Version 16 of Branches/Driver_Improvements
- Timestamp:
- 2016-05-27T08:49:17+02:00 (8 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
Branches/Driver_Improvements
v15 v16 36 36 37 37 It is important to understand that readim2.f90 is written in time steps and not physical time. There are 2 time stepping in the driver : 38 - n: is the list of indices of the time axis in the forcing file and readdim2.f90 will step from n-1 to n.39 - itau: is the time stepping of ORCHIDEE and takes 'split' values between n-1 and n.38 - ''n'' : is the list of indices of the time axis in the forcing file and readdim2.f90 will step from n-1 to n. 39 - ''itau'' : is the time stepping of ORCHIDEE and takes 'split' values between n-1 and n. 40 40 41 The readdim2.f90 will then interpolate using values at n-1 and n to obtain the forcing at itauin, which is the current time step of ORCHIDEE. Once itau has stepped through split values the forcing values of n will be copied to n-1 and the new values for nwill be read from the netCDF file.41 The readdim2.f90 will then interpolate using values at ''n-1'' and ''n'' to obtain the forcing at ''itauin'', which is the current time step of ORCHIDEE. Once ''itau'' has stepped through ''split'' values the forcing values of ''n'' will be copied to ''n-1'' and the new values for ''n'' will be read from the netCDF file. 42 42 43 All this is done independently from the physical time for which the values of the forcing at n-1 and n are valid over. This poses no problem for instantaneous variables in the forcing file as then the indexing space (n) and the physical time are identical.43 All this is done independently from the physical time for which the values of the forcing at ''n-1'' and ''n'' are valid over. This poses no problem for instantaneous variables in the forcing file as then the indexing space (''n'') and the physical time are identical. 44 44 45 The issue becomes more complicated when we consider fluxes which are averaged over an interval. We can only make the assumption that the time corresponding to index nis within the time interval of the average. But we in readim2.f90 we have no way of knowing if it is at the start, end or anywhere else.45 The issue becomes more complicated when we consider fluxes which are averaged over an interval. We can only make the assumption that the time corresponding to index ''n'' is within the time interval of the average. But we in readim2.f90 we have no way of knowing if it is at the start, end or anywhere else. 46 46 47 47 The graphic below explains probably better the above description. … … 50 50 51 51 52 The figure is especially meant to illustrate that the interpolation for fluxes assumes that the interval of validity of the fluxes at index n is assumed implicitly by readim2.f90 to be between n-1/2 and n+1/2.52 The figure is especially meant to illustrate that the interpolation for fluxes assumes that the interval of validity of the fluxes at index ''n'' is assumed implicitly by readim2.f90 to be between ''n-1/2'' and ''n+1/2''. 53 53 54 54 == Add an explicit time information within the forcing files ==