This routine will provide the external pressure-gradient, either from analytical expressions or read-in from a file. The subroutine is called in get_all_obs() 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,fsecs use observations, only: init_saved_vars,read_obs use observations, only: pi,h_press,dpdx,dpdy use observations, only: AmpMu,AmpMv,PhaseMu,PhaseMv,PeriodM use observations, only: AmpSu,AmpSv,PhaseSu,PhaseSv,PeriodS use observations, only: PressConstU,PressConstV,PressHeight IMPLICIT NONEINPUT PARAMETERS:
integer, intent(in) :: method,unit,jul,secsREVISION HISTORY:
Original author(s): Karsten Bolding