This routine will provide the internal pressure-gradients, either analytically prescribed or read from a file. The subroutine is called in the get_all_obs() subroutine as part of the main integration loop. The spatial interpolation is done via the reading routine and the temporal interpolation is done in this routine.
use time, only: time_diff,julian_day use observations, only: init_saved_vars,read_profiles use observations, only: dsdx,dsdy,dtdx,dtdy IMPLICIT NONEINPUT PARAMETERS:
integer, intent(in) :: method integer, intent(in) :: unit integer, intent(in) :: jul,secs integer, intent(in) :: nlev REALTYPE, intent(in) :: z(0:nlev)REVISION HISTORY:
Original author(s): Karsten Bolding