Hello,

I am using inverse distance weighting with gstat for temperature just to compare it with other kriging methods. For a temperature field, I have removed the elevation effect at the stations by normalizing all stations to the lowest with a lapse rate of -6.5 deg C /1000 m.

Now, I want to krig the normalized data and then add back the lapse rate multiplied by the difference in elevation to get the real temperatures at all real elevations. However, how can I do this before cross validation?

I only know cross validation with the following command where it does not allow manipulation of the kriged field before cross validation:

crossval_idw <- krige.cv(Temp93hr~1, DataCoordhr)

Thank you in advance for your help.

Cheers,

Cara

	[[alternative HTML version deleted]]

