<div dir="ltr"><div>Hi Geeta,</div><div>It would be difficult to answer this question without a data file and the script used to plot these figures. Is it possible to share data and the NCL script?</div><div>Cheers,</div><div>Rashed<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Sep 17, 2020 at 8:05 PM 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>
This issue was sorted out a few days ago. Now I am finding a small problem which is the shift of the position of the clouds.<br>
I mean when I plot the data for regional domain, clouds seem to be shifted by about 2 degree compared to the original domain plot (attachments) .<br>
<br>
<br>
Requesting once again to sort out this issue too. <br>
Geeta.<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
On Friday, 11 September, 2020, 10:10:16 pm IST, Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>> wrote:<br>
<br>
<br>
<br>
<br>
<br>
Are the latitude and longitudes on the file?<br>
1) If "yes", just read as indicated below.2) If "no", the projection information must be used to generate the latitude and longitude information<br>
------------------------<br>
The following functions may be useful when the latitudes and longitudes are available<br>
<br>
region_ind<br>
<br>
getind_latlon2d<br>
<br>
<br>
-------------------------<br>
f = addfile("3DIMG_23JUN2018_0300_L1C_ASIA_MER.h5" , "r")<br>
<br>
glon2d = f->...longitude_variable ...<br>
glat2d = f->...latitude_variable...<br>
<br>
printMinMax(glat2d, 0)<br>
printMinMax(glon2d, 0)<br>
<br>
latS = 10 ; California [*rough*]<br>
latN = 20<br>
lonW = 70<br>
lonE = 75<br>
<br>
<br>
ji = region_ind (glat2d,glon2d, latS, latN, lonW, lonE)<br>
<br>
jStrt = ji(0) ; lat start<br>
jLast = ji(1) ; lat last<br>
iStrt = ji(2) ; lon start<br>
iLast = ji(3) ; lon last<br>
<br>
LAT2D = glat2d(jStrt:jLast,iStrt:iLast)<br>
LON2D = glon2d(jStrt:jLast,iStrt:iLast)<br>
<br>
printMinMax(LAT2D, 0)<br>
printMinMax(LON2D, 0)<br>
<br>
; read data just for the region of interest<br>
x = f->IMG_TIR1(0,jStrt:jLast,iStrt:iLast)<br>
printVarSummary(x)<br>
<br>
<br>
printMinMax(x)<br>
<br>
<br>
<br>
<br>
<br>
On Fri, Sep 11, 2020 at 8:41 AM Geeta Geeta via ncl-talk <<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a>> wrote:<br>
the NxM grid is the larger grid, of the size 55N and 61E.<br>
<br>
will it be minI = 70* 61/(105.5-44.5) ??<br>
<br>
<br>
<br>
does not seem to be correct, I am not achieving anything<br>
<br>
<br>
<br>
<br>
<br>
Geeta.<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
On Friday, 11 September, 2020, 07:48:54 pm IST, Rick Brownrigg <<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>> wrote:<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
Hi,<br>
<br>
No, not in meters -- you're trying to find the indices of the grid cells that correspond to the corners of your smaller domain. For instance, what are the values i,j that correspond to the corner (-10N,44.5E)<br>
<br>
<br>
<br>
RB<br>
<br>
<br>
<br>
On Fri, Sep 11, 2020 at 8:11 AM Geeta Geeta <<a href="mailto:geetag54@yahoo.com" target="_blank">geetag54@yahoo.com</a>> wrote:<br>
<br>
Soory Mr Rick. It is not clear to me.<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
4 corners of large domain are (-10N,44.5E) , (-10,105.5), (45.5,105.5), (45.5,44.5).<br>
<br>
<br>
<br>
I want small domain to be 10-20N/70-75E.<br>
<br>
<br>
<br>
So I have to calculate the X Y (distance in meters)?.<br>
<br>
<br>
<br>
kindly clarify .<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
Geeta.<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
On Friday, 11 September, 2020, 07:20:57 pm IST, Rick Brownrigg <<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>> wrote:<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
Hi,<br>
<br>
<br>
<br>
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>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
minI = tointeger( minLonSmallerDomain * (M -0) / (maxLon - minLon) ) ; truncating<br>
<br>
<br>
<br>
maxI = tointeger( maxLonSmallerDomain * (M -0) / (maxLon - minLon) )...similarly for Y...<br>
<br>
<br>
<br>
and then:<br>
<br>
<br>
<br>
Plot(0) = gsn_csm_contour_map(wks,TIR1_INVERTED(0,minJ:maxJ,minI:maxI),rescount)<br>
<br>
<br>
<br>
and don't forget to substitute appropriate values into the mpLeft/RightCornerLat/LonF resources.<br>
<br>
<br>
<br>
Does that make sense?<br>
<br>
<br>
<br>
Rick<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
On Fri, Sep 11, 2020 at 4:22 AM Geeta Geeta via ncl-talk <<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a>> wrote:<br>
<br>
<br>
<br>
Hi.<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
I have satellite data which has a resultion of 4 km (TIR1) and there is projection information provided in the datafile.<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
integer Projection_Information ( DIM_000 )<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
DIMENSION_LIST : <ARRAY of 0 elements><br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
false_easting : 0<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
false_northing : 0<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
grid_mapping_name : mercatorpor Brightness Temperaturerature<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
longitude_of_projection_origin : 75<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
lower_left_lat_lon(degrees) : ( -10, 44.5 )<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
lower_right_lat_lon(degrees) : ( -10, 105.5 )<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
semi_major_axis : 6378137<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
semi_minor_axis : 6356752.3142<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
standard_parallel : 17.75<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
upper_left_lat_lon(degrees) : ( 45.5, 44.5 )<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
upper_left_xy(meters) : ( -3234623.003937, 5401854.420193 )<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
upper_right_lat_lon(degrees) : ( 45.5, 105.5 )<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
The lat/lon are not provided.<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<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>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
My questions are 1. I am Not able to plot the data for a smaller domain of my ineterst<br>
<br>
<br>
<br>
<br>
<br>
<br>
<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>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
the Images are attached.<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
Request someone to help me.<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
Geeta._______________________________________________<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
ncl-talk mailing list<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a><br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
List instructions, subscriber options, unsubscribe:<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<a href="https://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">https://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br>
_______________________________________________<br>
<br>
ncl-talk mailing list<br>
<br>
<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a><br>
<br>
List instructions, subscriber options, unsubscribe:<br>
<br>
<a href="https://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">https://mailman.ucar.edu/mailman/listinfo/ncl-talk</a>_______________________________________________<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>