Hi,
I'm trying to calculate the area-weighted mean inside OpenIFS using the following:
IFS_array(klon)
IFS_lat(klon)
sumx = 0.0
sumwgt = 0.0
Using the cosine weights, and the latitude is defined in radians:
do jl = kidia,kfdia
sumx = sumx + IFS_array(jl) * cos(IFS_lat(jl))
sumwgt = sumwgt + cos(IFS_lat(jl))
end do
area_weighted_mean = sumx/sumwgt
I'm wondering is this method is ok to use in OpenIFS given the reduced gaussian grid?
Appreciate your thoughts on this,
/M