This routine is responsible for providing sane values to `observed' vertical velocity which will then be applied for vertical advection of mean flow properties. A height and a vertical velocity value are either set to constant values or read from a file. The height will be assigned to be the position of maximum vertical velocity, and the vertical profiles of vertical velocity will be then constructed in such a way that the velocity is linearly decreasing away from this height, with zero values at the surface and the bottom. The subroutine is called in the get_all_obs() subroutine as part of the main integration loop. In case of observations from file the temporal interpolation is done in this routine.
use time, only: time_diff,julian_day use observations, only: init_saved_vars,read_obs use observations, only: w_adv,w_adv0,w_adv_height0,w_height IMPLICIT NONEINPUT PARAMETERS:
integer, intent(in) :: method,unit,jul,secsREVISION HISTORY:
Original author(s): Karsten Bolding