function [Xgrid,Ygrid]=obs_sample(InpObs,OutObs,Tstr,Tend);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Copyright (c) 2006 Rutgers University.                                    %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%                                                                           %
% function [Xgrid,Ygrid]=obs_ijpos(gname,obslon,obslat)                     %
%                                                                           %
% This function computes the fractional (I,J) model RHO-grid locations      %
% at the observations (obslon,obslat) locations.                            %
%                                                                           %
% On Input:                                                                 %
%                                                                           %
%    gname       NetCDF grid file name (character string).                  %
%    obslon      Observation longitude (positive, degrees_east).            %
%    obslat      Observation latitude (positive, degrees_north).            %
%    Ccorrection Small correction for curvilinear grids (0: no, 1: yes).    %
%                                                                           %
% On Ouput:                                                                 %
%                                                                           %
%    Xgrid       Model fractional X-grid locations at RHO-points.           %
%    Ygrid       Model fractional Y-grid locations at RHO-points.           %
%                                                                           %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
