Hi <span class="Apple-style-span" style="font-family: arial, sans-serif; font-size: 13px; border-collapse: collapse; ">Matthias,</span><div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;"><br>
</span></font></div><div><span class="Apple-style-span" style="font-size: 13px; "></span><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;">I am glad you liked it.</span></font></div>
<div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;"><br></span></font></div><div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;">I am using FORTRAN as it is convenient to me. You can use any other language which can handle binary files.</span></font></div>
<div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;"><br></span></font></div><div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;">Basically, what GIS (ArcView) gives us is the IDs of the grid points which are inside the country boundary. All my program does is assign 1s to those points and 0s to rest of the points in a domain considered (mine happens to be 86x46).</span></font></div>
<div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;"><br></span></font></div><div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;">In binary file, replace all values by 'undefined' values (-999.0) where masking file values are zero and save this as new file. Thats it.</span></font></div>
<div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;"><br></span></font></div><div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;">Attached is the RAR file which comprises the following.</span></font></div>
<div><span class="Apple-style-span" style="font-family: arial, sans-serif; border-collapse: collapse; ">ASCII</span></div><div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;"><br>
</span></font></div><div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;">masking file</span></font></div><div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;">unmasked-data file</span></font></div>
<div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;">masked-data file</span></font></div><div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;">(header file is also included to import it into ArcView)</span></font></div>
<div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;">control file - to display in grads</span></font></div><div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;"><br>
</span></font></div><div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;">BINARY</span></font></div><div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;"><br>
</span></font></div><div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;">unmasked data file</span></font></div><div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;">masked data file</span></font></div>
<div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;"> </span></font></div><div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;"><br>
</span></font></div><div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;">If you have any query please let me know.</span></font></div><div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;"><br>
</span></font></div><div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;">cheers!</span></font></div><div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;"><br>
</span></font></div><div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;">-rupak</span></font></div><div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;"><br>
</span></font></div><div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;"><br></span></font></div><div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;"><br>
</span></font><div class="gmail_quote">On Wed, Oct 28, 2009 at 4:47 PM, Matthias Fripp <span dir="ltr"><<a href="mailto:matthias.fripp@eci.ox.ac.uk">matthias.fripp@eci.ox.ac.uk</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div style="word-wrap:break-word">Wow, you used almost exactly the same approach I suggested (but with more Fortran and less Excel)! Glad to know I was on to something.<div><br></div><div>I'm not fluent in Fortran, so I had trouble deciphering the cr_mask_fl.f program that converts an ASCII list of points into an ASCII masking file that GRADS can read. Would you be able to send a sample of the final ASCII masking file that you produced and the command you use to load it into GRADS? Or if there's an obvious place to look in the instructions, that would be great too. So far I only know how to get GRADS to read local or remote binary files. It would be great to have a way to display data from ASCII/text files that I create myself.</div>
<div><br></div><div>Best,</div><div><br></div><div>Matthias</div><div><br><div><div><div></div><div class="h5"><div>On Oct 28, 2009, at 8:39 AM, Rupak Rajbhandari wrote:</div><br></div></div><blockquote type="cite"><div><div>
</div><div class="h5">Hi,<div><br></div><div>Creating a masking file is not a simple procedure. Here I am sending you a masking instruction I prepared (last year for small training group in Nepal). The masking file is for 0.1 resolution, so you need to edit a little. The attached RAR file consists of a instruction file and two fortran programs. In addition you will need ArcView GIS and India shape file (boundary polygon) and fortran compiler. You should first edit the program to suit your resolution. The rest should be clear from the instruction file.</div>
<div><br></div><div>Once the masking file has been created, simply replace the master file with undefined values where grid point point values are zero in the masking file.</div><div><br></div><div>Hope it helps!</div><div>
<br></div><div>-rupak</div><div><br></div><div><br></div><div><br></div><div> <br><br><div class="gmail_quote">On Tue, Oct 27, 2009 at 9:04 PM, Matthias Fripp <span dir="ltr"><<a href="mailto:matthias.fripp@eci.ox.ac.uk" target="_blank">matthias.fripp@eci.ox.ac.uk</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word"><div>It sounds like you will need to do something like this:</div><div><br></div><div>
1. Create a file listing all the grid points in or near India (e.g., within a simple latitude and longitude range that covers all of India and some surrounding areas).</div><div><br></div><div>2. Import this list into a GIS program and intersect it with India's land area.</div>
<div><br></div><div>3. Export a list of Indian grid points from the GIS program.</div><div><br></div><div>4. Import this list into GRADS.</div><div><br></div><div>Step 1 would be easy to do, e.g., in a text (.csv) file created in Excel, if your grid is regularly spaced.</div>
<div><br></div><div>Step 2 and 3 are easy if you have a GIS program and difficult otherwise. I could help you with that if you want.</div><div><br></div><div>I don't know how you would do Step 4. I can easily create a text file listing all the relevant points, but I don't know how you would reference that text file from GRADS to create a mask variable. Anyone else have ideas about this?</div>
<div><br></div><font color="#888888"><div>Matthias</div></font><div><br><div><div class="gmail_quote">On Tue, Oct 27, 2009 at 7:23 AM, Sachin Ghude <span dir="ltr"><<a href="mailto:sachinghude@tropmet.res.in" target="_blank">sachinghude@tropmet.res.in</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left-width:1px;border-left-style:solid;border-left-color:rgb(204, 204, 204);margin-top:0pt;margin-right:0pt;margin-bottom:0pt;margin-left:0.8ex;padding-left:1ex"> Hi,<br> <br>
I want to plot tropospheric NO2 only over the Indian grid points (not Nepal, Pakistan<br> and BanglaDesh). Is there any such utility so that it will enable me to plot data only<br> for India girds in GRADs. ( my grid size is 0.25 x 0.25 degrees).<br>
<br> please do the needful..<br> <br> bye..<br> <font color="#888888"><br> Sachin<br> </font></blockquote></div></div><br></div></div></blockquote></div><br></div> </div></div><span><masking.rar></span></blockquote>
</div><br></div></div></blockquote></div><br></div>