function postpr_flow (dirnow,dirfig,start_date);

stop_date = [datestr(now,'yyyymmdd') '000000'];

runids    = {'csm' , 'hol' , 'egm'};
%
% Matroos station names
%

loc_matroos = {'vlissingen'           , 'hoekvanholland'   , 'ijmuiden'       , 'den helder' , ...
               'terschelling noordzee', 'huibertgat'       };
%
% CSM station names
%

loc_csm     = {'VLISSGN '             , 'HOEKVHLD'         , 'IJMDBTHVN'      , 'DENHDR' , ...
               'TERSLNZE'             , 'HUIBGT'           };
%
% Get the observed water level data from Matroos
%
time_obs   (1:length(loc_matroos),1:14400) = NaN;
wl_obs     (1:length(loc_matroos),1:14400) = NaN;

for istat = 1: length(loc_matroos);
    t(1)  = NaN;
    wl(1) = NaN;
    try;
      [t, wl]=get_series('waterlevel'      , 'observed'          , ...
                                             loc_matroos{istat}  , [start_date '000000'], ...
                                             stop_date                               );
       pause(10);
       time_obs(istat,1:length(t))  = t;
       wl_obs  (istat,1:length(wl)) = wl;
    catch;
       clear t, wl;
       pause (600);
      [t, wl]=get_series('waterlevel'      , 'observed'          , ...
                                             loc_matroos{istat}  , [start_date '000000'], ...
                                             stop_date                               );
       pause(10);
       time_obs(istat,1:length(t))  = t;
       wl_obs  (istat,1:length(wl)) = wl;
    end;
end;

%
% Get the computed water level data from csm simulation
%

cd (dirnow);

for istat = 1: length(loc_csm);
   for irun = 1: length(runids);
      nfs     = vs_use(strcat('trih-',runids{irun}),'quiet');

      stations= vs_get('his-const'  ,'NAMST','quiet');

      nstat = strmatch (loc_csm(istat), stations);
      if (length(nstat) > 0);
         val   = qpread   (nfs,'water level','data',0,nstat);
         if (length(val.Time) > 0);
            time_cmp (istat,1:length(val.Time)) = val.Time;
            wl_cmp   (istat,1:length(val.Val))  = val.Val;
         end;
         clear val;
      end;
   end;
end;

figure ('Units','centimeters','Paperunits','centimeters','Position',[1  1   38  20], ...
        'PaperPosition',[1  1   38  20]);

sp(1) = subplot('Position',[ 0.075  0.55 0.25 0.40]);
sp(2) = subplot('Position',[ 0.375  0.55 0.25 0.40]);
sp(3) = subplot('Position',[ 0.675  0.55 0.25 0.40]);
sp(4) = subplot('Position',[ 0.075  0.05 0.25 0.40]);
sp(5) = subplot('Position',[ 0.375  0.05 0.25 0.40]);
sp(6) = subplot('Position',[ 0.675  0.05 0.25 0.40]);

for istat = 1: length(loc_matroos);
    axes (sp(istat));
    plot (time_cmp(istat,:),wl_cmp(istat,:),'r','LineWidth',2.0);
    hold on
    plot (time_obs(istat,:),wl_obs(istat,:),'k.','MarkerSize',10.);
    axis([datenum(start_date,'yyyymmdd') + 4   datenum(start_date,'yyyymmdd') + 9   -2.5    2.5]);
    set (gca,'LineWidth' ,3.0);
    set (gca,'XTick'     ,[datenum(start_date,'yyyymmdd') + 4 :2: datenum(start_date,'yyyymmdd') +  9]);
    set (gca,'XTickLabel',datestr(get(gca,'XTick')));
    set (gca,'FontSize',15);
    title  (loc_matroos(istat),'Fontsize',20);
    ylabel ('Waterlevel (m NAP)','Fontsize',15);
    legend ('Hind/Forecast','Observation');
end;
pngnam = [dirfig '\waterlevel'];
print (gcf,'-dpng',pngnam);

close;

for istat = 1: length(loc_matroos);
    figure ('Units','centimeters','Paperunits','centimeters','Position',[1  1   38  20], ...
            'PaperPosition',[1  1   38  20]);
    plot (time_cmp(istat,:),wl_cmp(istat,:),'r','LineWidth',2.0);
    hold on
    plot (time_obs(istat,:),wl_obs(istat,:),'k.','MarkerSize',10.);
    axis([datenum(start_date,'yyyymmdd') + 4   datenum(start_date,'yyyymmdd') + 9   -2.5    2.5]);
    set (gca,'LineWidth' ,3.0);
    set (gca,'XTick'     ,[datenum(start_date,'yyyymmdd') + 4 :1: datenum(start_date,'yyyymmdd') + 9]);
    set (gca,'XTickLabel',datestr(get(gca,'XTick')));
    set (gca,'FontSize',15);
    title  (loc_matroos(istat),'Fontsize',20);
    ylabel ('Waterlevel (m NAP)','Fontsize',15);
    legend ('Hind/Forecast','Observation');
    pngnam = [dirfig '\waterlevel_' loc_matroos{istat}];
    print (gcf,'-dpng',pngnam);
    close;
end;
