Dear Sravan,<br><br>It appears that there is a mismatch between your .ctl file and fortran declarations. <br><br><span lang="EN">You have this as a fortran declaration: parameter(nx=129,ny=144,nz=171,ns=30,nt=41,nd=5)</span><br>
<br>but your .ctl file contains 331 x-points and 306 y-points (check your z-levels as well). In using pdef, ARWPost is regridding to a new resolution, which will be reflected in your .ctl file.<br><br>matei<br><br><b><font style="color: rgb(255, 102, 0);" color="#888888">==============================</font><font style="color: rgb(255, 102, 0);" color="#888888">====</font><font style="color: rgb(255, 102, 0);" color="#888888">====</font></b><br>
<b>Matei Georgescu, Ph.D.</b><br><i style="color: rgb(255, 102, 0);">School of Mathematical and Statistical Sciences and<br>Center for Environmental Fluid Dynamics</i><br><i><span style="color: rgb(255, 102, 0);">Arizona State University</span></i><br>
<b>Phone</b>:
(480) 965-3770<br><b>Fax</b>: (480) 965-8119<br><b>E-mail</b>: <a href="mailto:Matei.Georgescu@asu.edu" target="_blank">Matei.Georgescu@asu.edu</a><br><b><font style="color: rgb(255, 102, 0);" color="#888888">==============================</font><font style="color: rgb(255, 102, 0);" color="#888888">====</font><font style="color: rgb(255, 102, 0);" color="#888888">====</font></b><br>
<br><div class="gmail_quote">On Tue, Oct 26, 2010 at 10:08 PM, Sravan <span dir="ltr"><<a href="mailto:sravan_sase@yahoo.in">sravan_sase@yahoo.in</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<div><div style="font-family: times new roman,new york,times,serif; font-size: 12pt; color: rgb(0, 0, 0);"><div></div>
<div>Dear Users</div>
<div>I am using ARW post to reate .ctl and .dat files for the GrADS.</div>
<div>Now I want to extract some perameters from the .dat file for some particular locations.</div>
<div>I have used a fortran script to extract the perameters. But it is not picking the one which I want.</div>
<div>I am giving the .ctl file for your reference.</div><span lang="EN">
<p><font color="#ff0000">dset ^dom3.dat</font></p>
<p><font color="#0000ff">options byteswapped</font></p>
<p><font color="#ff0000">undef 1.e30</font></p>
<p><font color="#ff0000">title OUTPUT FROM WRF V3.2.1 MODEL</font></p>
<p><font color="#ff0000">pdef 129 144 lcc 31.672 78.135 65.000 72.500 30.00000 10.00000 77.00000 9000.000 9000.000</font></p>
<p><font color="#ff0000">xdef 331 linear 71.47070 0.04054054</font></p>
<p><font color="#ff0000">ydef 306 linear 25.33539 0.04054054</font></p>
<p><font color="#ff0000">zdef 19 levels </font></p>
<p><font color="#ff0000" size="1">1000.00000</font></p>
<p><font color="#ff0000" size="1">950.00000</font></p>
<p><font color="#ff0000" size="1">900.00000</font></p>
<p><font color="#ff0000" size="1">850.00000</font></p>
<p><font color="#ff0000" size="1">800.00000</font></p>
<p><font color="#ff0000" size="1">750.00000</font></p>
<p><font color="#ff0000" size="1">700.00000</font></p>
<p><font color="#ff0000" size="1">650.00000</font></p>
<p><font color="#ff0000" size="1">600.00000</font></p>
<p><font color="#ff0000" size="1">550.00000</font></p>
<p><font color="#ff0000" size="1">500.00000</font></p>
<p><font color="#ff0000" size="1">450.00000</font></p>
<p><font color="#ff0000" size="1">400.00000</font></p>
<p><font color="#ff0000" size="1">350.00000</font></p>
<p><font color="#ff0000" size="1">300.00000</font></p>
<p><font color="#ff0000" size="1">250.00000</font></p>
<p><font color="#ff0000" size="1">200.00000</font></p>
<p><font color="#ff0000" size="1">150.00000</font></p>
<p><font color="#ff0000" size="1">100.00000</font></p>
<p><font color="#ff0000">tdef 41 linear 00Z26OCT2010 180MN </font></p>
<p><font color="#ff0000">VARS 27</font></p>
<p><font color="#ff0000" size="1">LU_INDEX 1 0 LAND USE CATEGORY (-)</font></p>
<p><font color="#ff0000" size="1">U 19 0 x-wind component (m s-1)</font></p>
<p><font color="#ff0000" size="1">V 19 0 y-wind component (m s-1)</font></p>
<p><font color="#ff0000" size="1">W 19 0 z-wind component (m s-1)</font></p>
<p><font color="#ff0000" size="1">PH 19 0 perturbation geopotential (m2 s-2)</font></p>
<p><font color="#ff0000" size="1">PHB 19 0 base-state geopotential (m2 s-2)</font></p>
<p><font color="#ff0000" size="1">T 19 0 perturbation potential temperature (theta-t0) (K)</font></p>
<p><font color="#ff0000" size="1">MU 1 0 perturbation dry air mass in column (Pa)</font></p>
<p><font color="#ff0000" size="1">MUB 1 0 base state dry air mass in column (Pa)</font></p>
<p><font color="#ff0000" size="1">P 19 0 perturbation pressure (Pa)</font></p>
<p><font color="#ff0000" size="1">PB 19 0 BASE STATE PRESSURE (Pa)</font></p>
<p><font color="#ff0000" size="1">T2 1 0 TEMP at 2 M (K)</font></p>
<p><font color="#ff0000" size="1">TH2 1 0 POT TEMP at 2 M (K)</font></p>
<p><font color="#ff0000" size="1">PSFC 1 0 SFC PRESSURE (Pa)</font></p>
<p><font color="#ff0000" size="1">U10 1 0 U at 10 M (m s-1)</font></p>
<p><font color="#ff0000" size="1">V10 1 0 V at 10 M (m s-1)</font></p>
<p><font color="#ff0000" size="1">LANDMASK 1 0 LAND MASK (1 FOR LAND, 0 FOR WATER) (-)</font></p>
<p><font color="#ff0000" size="1">VEGFRA 1 0 VEGETATION FRACTION (-)</font></p>
<p><font color="#ff0000" size="1">SST 1 0 SEA SURFACE TEMPERATURE (K)</font></p>
<p><font color="#ff0000" size="1">HGT 1 0 Terrain Height (m)</font></p>
<p><font color="#ff0000" size="1">TSK 1 0 SURFACE SKIN TEMPERATURE (K)</font></p>
<p><font color="#ff0000" size="1">RAINC 1 0 ACCUMULATED TOTAL CUMULUS PRECIPITATION (mm)</font></p>
<p><font color="#ff0000" size="1">RAINNC 1 0 ACCUMULATED TOTAL GRID SCALE PRECIPITATION (mm)</font></p>
<p><font color="#ff0000" size="1">XLAT 1 0 LATITUDE, SOUTH IS NEGATIVE (degree_north)</font></p>
<p><font color="#ff0000" size="1">XLONG 1 0 LONGITUDE, WEST IS NEGATIVE (degree_east)</font></p>
<p><font color="#ff0000" size="1">TMN 1 0 SOIL TEMPERATURE AT LOWER BOUNDARY (K)</font></p>
<p><font color="#ff0000" size="1">PBLH 1 0 PBL HEIGHT (m)</font></p>
<p><font color="#ff0000" size="1">ENDVARS</font></p>
<div></div></span>And the fortran code which I have used to extract the variables is given below.</div><span lang="EN">
<p>parameter(nx=129,ny=144,nz=171,ns=30,nt=41,nd=5)</p>
<p>dimension x(nx,ny,nz),dum(nx,ny)</p></span><span lang="EN">
<p>open(12,file='./dom3.dat',form='unformatted',access='direct',recl=nx*ny*4)</p><span lang="EN">
<p>do 5 it = 1,nt</p>
<p>do 10 iz = 1,nz</p>
<p>read (12) dum</p>
<p>do 15 ix = 1,nx</p>
<p>do 16 iy = 1,ny</p>
<p>x(ix,iy,iz) = dum(ix,iy)</p>
<p>16 continue</p>
<p>15 continue</p>
<p>10 continue</p><span lang="EN">
<p>do 99 k=1,10</p>
<p>p1(k,it) = x(ix,iy,166)+x(ix,iy,167)</p>
<div></div></span></span></span>write (*,*) p1(k,it)</div>
<div>99 continue</div>
<div>stop</div>
<div>end</div>
<div> </div>
<div>I am using a 64 bit Linux Machine with pgf77 compiler.</div>
<div>Can any one clarify what <font color="#0000ff">options byteswapped means?</font></div>
<div>Kindly Help Me</div>
<div> </div>
<div>Regards</div>
<div>Sravan<br> </div>############################<br>M.SRAVAN KUMAR<br>SCIENTIST<br>RDC SASE<br>CHANDIGARH - 160036<br>INDIA<br>############################
<div></div><br><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></blockquote></div>