subroutine get_w_adv(method,unit,jul,secs)

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
    integer, intent(in)                 :: method,unit,jul,secs
    Original author(s): Karsten Bolding
    $Log: get_w_adv.F90,v $
    Revision 1.9  2010-09-17 12:53:50  jorn
    extensive code clean-up to ensure proper initialization and clean-up of all variables
    Revision 1.8  2006-11-27 09:25:18  kbk
    use logical var init_saved_vars to initialise saved variables
    Revision 1.7  2005/11/15 11:02:32  lars
    documentation finish for print
    Revision 1.6  2005/08/15 11:54:01  hb
    sequence of reading w_adv and w_height changed,
    w_adv_height0 introduced, documentation extended
    Revision 1.5  2005/06/27 13:44:07  kbk
    modified + removed traling blanks
    Revision 1.4  2003/03/28 09:20:35  kbk
    added new copyright to files
    Revision 1.3  2003/03/28 09:02:09  kbk
    removed tabs
    Revision 1.2  2003/03/10 08:51:57  gotm
    Improved documentation and cleaned up code
    Revision  2001/02/12 15:55:58  gotm
    initial import into CVS

Karsten Bolding 2012-01-24