<html><head></head><body bgcolor="#FFFFFF"><div>Arlindo,</div><div><br></div><div> I have a GPCP option in lats4d.sh; enter </div><div><br></div><div>% lats4d.sh -h</div><div><br></div><div>for the correct option.<br><br>
Sent from my iPhone</div><div><br>On Nov 8, 2011, at 4:59 PM, Arlindo Meque <<a href="mailto:mequitomz@yahoo.com.br">mequitomz@yahoo.com.br</a>> wrote:<br><br></div><div></div><blockquote type="cite"><div><div style="color:#000; background-color:#fff; font-family:times new roman, new york, times, serif;font-size:12pt">
<div class="MsoNormal">The lats4d.sh script worked properly. The only problem is that
the output from the lats4d.sh is a file with 144 points in x-direction and 73
points in y-direction while my observation file (gpcp) has 144 points in
x-direction and 72 points in y-direction, i.e, the number of points in y-direction
from the forecast and observation differ by 1.</div>
<div class="MsoNormal"> </div>
<div class="MsoNormal">I’ve tried again the re() <span style="mso-spacerun:yes"> </span>command but the but GrADS fail<span style="mso-spacerun:yes"> </span>to discern the coordinates:</div>
<div class="MsoNormal">ga-> sdfopen<span style="mso-spacerun:yes">
</span><a href="http://forecast.nc">forecast.nc</a><span style="mso-spacerun:yes"> </span></div>
<div class="MsoNormal">Scanning self-describing file:<span style="mso-spacerun:yes"> </span><a href="http://forecast.nc">forecast.nc</a></div>
<div class="MsoNormal">SDF file <a href="http://forecast.nc">forecast.nc</a> is open as file 1</div>
<div class="MsoNormal">LON set to 0 360 </div>
<div class="MsoNormal">LAT set to -90 90 </div>
<div class="MsoNormal">LEV set to 0 0 </div>
<div class="MsoNormal">Time values set: 2009:9:16:0 2009:9:16:0 </div>
<div class="MsoNormal">E set to 1 1 </div>
<div class="MsoNormal">ga-> define
gg=re(precip,144,linear,0,2.5,72,linear,-90,2.5,bl)</div>
<div class="MsoNormal">Define memory allocation size = 82944 bytes</div>
<div class="MsoNormal">ga-> set sdfwrite <a href="http://test.nc">test.nc</a></div>
<div class="MsoNormal">SDFWrite file name = <a href="http://test.nc">test.nc</a></div>
<div class="MsoNormal">SDFWrite will replace an existing file</div>
<div class="MsoNormal">SDFwrite file will have same number of dimensions as defined
variable</div>
<div class="MsoNormal">ga-> sdfwrite gg<span style="mso-spacerun:yes">
</span></div>
<div class="MsoNormal">Wrote variable gg to <a href="http://test.nc">test.nc</a></div>
<div class="MsoNormal">ga-> reinit</div>
<div class="MsoNormal">No hardcopy metafile open</div>
<div class="MsoNormal">All files closed; all defined objects released;</div>
<div class="MsoNormal">All GrADS attributes have been reinitialized</div>
<div class="MsoNormal">ga-> sdfopen <a href="http://test.nc">test.nc</a> </div>
<div class="MsoNormal">Scanning self-describing file:<span style="mso-spacerun:yes"> </span><a href="http://test.nc">test.nc</a></div>
<div class="MsoNormal">SDF file has no discernable time coordinate -- using default
values.</div>
<div class="MsoNormal">SDF file <a href="http://test.nc">test.nc</a> is open as file 1</div>
<div class="MsoNormal">LON set to 0 360 </div>
<div class="MsoNormal">LAT set to -90 87.5 </div>
<span style="font-size:11.0pt;line-height:115%;font-family:"Calibri","sans-serif";
mso-ascii-theme-font:minor-latin;mso-fareast-font-family:Calibri;mso-fareast-theme-font:
minor-latin;mso-hansi-theme-font:minor-latin;mso-bidi-font-family:"Times New Roman";
mso-bidi-theme-font:minor-bidi;mso-ansi-language:EN-US;mso-fareast-language:
EN-US;mso-bidi-language:AR-SA">LEV set to 0 0 </span>
<div class="MsoNormal"> </div>
<br><br><br><br><br><div style="font-family: times new roman, new york, times, serif; font-size: 12pt;"><div style="font-family: times new roman, new york, times, serif; font-size: 12pt;"><font size="2" face="Arial"><hr size="1">
<b><span style="font-weight:bold;">De:</span></b> Jennifer Adams <<a href="mailto:jma@cola.iges.org">jma@cola.iges.org</a>><br><b><span style="font-weight: bold;">Para:</span></b> GrADS Users Forum <<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>><br>
<b><span style="font-weight: bold;">Enviadas:</span></b> Terça-feira, 8 de Novembro de 2011 12:13<br><b><span style="font-weight: bold;">Assunto:</span></b> Re: [gradsusr] Enc: regriding problem<br></font><br><div id="yiv1685382208">
<div><br><div><div>On Nov 7, 2011, at 2:46 PM, Arlindo da Silva wrote:</div><br class="yiv1685382208Apple-interchange-newline"><blockquote type="cite">On Mon, Nov 7, 2011 at 10:11 AM, Arlindo Meque <span dir="ltr"><<a rel="nofollow" target="_blank" href="mailto:mequitomz@yahoo.com.br">mequitomz@yahoo.com.br</a>></span> wrote:<br>
<div class="yiv1685382208gmail_quote"><blockquote class="yiv1685382208gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div><div style="color:#000;background-color:#fff;font-family:times new roman, new york, times, serif;font-size:12pt;"><div style="font-family:times new roman, new york, times, serif;font-size:12pt;"><div style="font-family:times new roman, new york, times, serif;font-size:12pt;">
<div> Dear all,<br><br><br>1. I have a forecast file (1.875X1.25) with 192 points in x-direction and 145 points in y-direction. I am trying to regrid the file to match the observation (2.5X2.5) with 144 points in x-direction and 72 points on y-direction<br>
<br><br>2. I am able to create the regridded file, the problem is that the dimensions of the regridded file do not much those of the observation file, so I am not able to make the diffidence between the forecast and the observation.<br>
<br>3. Please find attached the forecast and observed file along with the regriding script.<br><br></div></div></div></div></div></blockquote><div><br></div><div>This more of a question for Jeniffer as I am not sure how sdfwrite figures out the coordinate variables of the output file (does it use the coordinates of the default file?).</div>
</div></blockquote><div>The sdfwrite copies the coordinate information from the default file at the time the 'define' command was invoked. </div><div>I have no idea what the re() function does when it returns the new grid. If you used lterp instead, you would have the right coordinate info because it's attached to the destination grid expression. </div>
<div><br></div><div>Since you are using opengrads, Arlindo's suggestion to use lats4d may also do what you need. </div><div>--Jennifer</div><div><br></div><br><blockquote type="cite"><div class="yiv1685382208gmail_quote">
<div> Now, a simple way to produce
a file interpolated to this particular grid is to use lats4d:</div>
<div><br></div><div>% lats4d.sh -i <a rel="nofollow" target="_blank" href="http://forecast.nc/">forecast.nc</a> -o fcstregridded -ncep2.5 -v</div><div><br></div><div>This is going to the standard NCEP 2.5x2.5 grid (which is not exactly what you had on your re() line). Are you sure your re() line is correct? It is kind of an unusual grid. In any event, if you would like to modify the output from ncep2.5 there is a simple way, but let me hear from you first.</div>
<div><br></div><div> (the other) Arlindo</div></div><div><br></div>-- <br>Arlindo da Silva<br><a rel="nofollow" target="_blank" href="mailto:dasilva@alum.mit.edu">dasilva@alum.mit.edu</a><br>
_______________________________________________<br>gradsusr mailing list<br><a rel="nofollow" target="_blank" href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br><a href="http://gradsusr.org/mailman/listinfo/gradsusr">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
</blockquote></div><br><div>
<span class="yiv1685382208Apple-style-span" style="border-collapse:separate;border-spacing:0px 0px;color:rgb(0, 0, 0);font-family:Helvetica;font-size:12px;font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:auto;text-indent:0px;text-transform:none;orphans:2;white-space:normal;widows:2;word-spacing:0px;"><span class="yiv1685382208Apple-style-span" style="border-collapse:separate;border-spacing:0px 0px;color:rgb(0, 0, 0);font-family:Helvetica;font-size:12px;font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:auto;text-indent:0px;text-transform:none;orphans:2;white-space:normal;widows:2;word-spacing:0px;"><span class="yiv1685382208Apple-style-span" style="border-collapse:separate;border-spacing:0px 0px;color:rgb(0, 0,
0);font-family:Helvetica;font-size:12px;font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-align:auto;text-indent:0px;text-transform:none;orphans:2;white-space:normal;widows:2;word-spacing:0px;"><div>
--</div><div>Jennifer M. Adams</div><div>IGES/COLA</div><div>4041 Powder Mill Road, Suite 302</div><div>Calverton, MD 20705</div><div><a rel="nofollow" target="_blank" href="mailto:jma@cola.iges.org">jma@cola.iges.org</a></div>
<div><br class="yiv1685382208khtml-block-placeholder"></div><br class="yiv1685382208Apple-interchange-newline"></span></span></span>
</div>
<br></div></div><br>_______________________________________________<br>gradsusr mailing list<br><a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br><a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br><br></div></div></div></div></blockquote><blockquote type="cite"><div><span>_______________________________________________</span><br><span>gradsusr mailing list</span><br><span><a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a></span><br>
<span><a href="http://gradsusr.org/mailman/listinfo/gradsusr">http://gradsusr.org/mailman/listinfo/gradsusr</a></span><br></div></blockquote></body></html>