<div dir="ltr"><div>Hi,</div><div><br></div><div>I *think* the trick here is to compute the i,j indices of the grid that correspond to your smaller domain. Since it's a mercator projection, and the original domain is relatively small (60x60 degs), I would think you could get decent results by using a simple linear scaling to find those values. So if the grid is NxM (LatxLon), and you extract the min/max lat/lons from the attributes into variables minLon, maxLon, etc.,  then something like:<br></div><div><br></div><div>  minI = tointeger( minLonSmallerDomain * (M -0) / (maxLon - minLon) )  ; truncating<br></div><div>  maxI = tointeger( maxLonSmallerDomain * (M -0) / (maxLon - minLon) )</div><div>...similarly for Y...</div><div><br></div><div>and then:</div><div><br></div><div> Plot(0)                           =       gsn_csm_contour_map(wks,TIR1_INVERTED(0,<span style="color:rgb(255,0,0)">minJ:maxJ,minI:maxI</span>),rescount)</div><div><br></div><div>and don't forget to substitute appropriate values into the mpLeft/RightCornerLat/LonF resources.</div><div><br></div><div>Does that make sense?</div><div><br></div><div>Rick</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Sep 11, 2020 at 4:22 AM Geeta Geeta via ncl-talk <<a href="mailto:ncl-talk@mailman.ucar.edu">ncl-talk@mailman.ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi.<br>
I have satellite data which has a resultion of 4 km (TIR1) and there is projection information provided in the datafile.<br>
integer Projection_Information ( DIM_000 )<br>
        DIMENSION_LIST :    <ARRAY of 0 elements><br>
        false_easting :      0<br>
        false_northing :      0<br>
        grid_mapping_name :    mercatorpor Brightness Temperaturerature<br>
        longitude_of_projection_origin :      75<br>
        lower_left_lat_lon(degrees) :    (  -10, 44.5 )<br>
        lower_right_lat_lon(degrees) :    (  -10, 105.5 )<br>
        semi_major_axis :    6378137<br>
        semi_minor_axis :    6356752.3142<br>
        standard_parallel :    17.75<br>
        upper_left_lat_lon(degrees) :    ( 45.5, 44.5 )<br>
        upper_left_xy(meters) :    ( -3234623.003937, 5401854.420193 )<br>
        upper_right_lat_lon(degrees) :    ( 45.5, 105.5 )<br>
The lat/lon are not provided.<br>
<br>
Plotting the TIR1 (count), radiance and Brightness temperature for the entire domain is happening by using information from the Projection info. This INFO is used from the Line numbers 109-114 of the code. <br>
<br>
My questions are 1. I am Not able to plot the data for a smaller domain of my ineterst<br>
                          2. since I am unable to attch the lat/lon arrays to the variables, I am unable to plot the data for a smaller domain and unable to write the data for a smaller domain of my intererst.<br>
the Images are attached.<br>
<br>
<br>
Request someone to help me.<br>
<br>
Geeta._______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="https://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">https://mailman.ucar.edu/mailman/listinfo/ncl-talk</a></blockquote></div>