[gradsusr] ARW Post for Grads

Thu Oct 28 05:29:10 EDT 2010

Dear Matei,
I want to inform you that the no. of points in x and y directions are as per the 
pdef in the ctl file.
Even the record length and all are perfectly matching. I have computed the no of 
levels also as per the control file.
One more thing is if I compute the same parameters from the WRF2GrADS output it 
is picking up the values very fine.
The problem is with the ARW post only. Even I have fixed the lat and lon values 
fix and picked up the 

24hr accumulated (rainc+rainnc) values. It is giving some exponentional values 
for some time steps. 

What can be done in this regard.

Thanks
Sravan

 ############################
M.SRAVAN KUMAR
SCIENTIST
RDC SASE
CHANDIGARH - 160036
INDIA
############################ 

________________________________
From: Matei Georgescu <Matei.Georgescu at asu.edu>
To: GrADS Users Forum <gradsusr at gradsusr.org>
Sent: Wed, 27 October, 2010 8:25:18 PM
Subject: Re: [gradsusr] ARW Post for Grads

Dear Sravan,

It appears that there is a mismatch between your .ctl file and fortran 
declarations. 

You have this as a fortran declaration: 
parameter(nx=129,ny=144,nz=171,ns=30,nt=41,nd=5)

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.

matei

======================================
Matei Georgescu, Ph.D.
School of Mathematical and Statistical Sciences and
Center for Environmental Fluid Dynamics
Arizona State University
Phone: (480) 965-3770
Fax: (480) 965-8119
E-mail: Matei.Georgescu at asu.edu
======================================

On Tue, Oct 26, 2010 at 10:08 PM, Sravan <sravan_sase at yahoo.in> wrote:

Dear Users
>I am using ARW post to reate .ctl and .dat files for the GrADS.
>Now I want to extract some perameters from the .dat file for some particular 
>locations.
>I have used a fortran script to extract the perameters. But it is not picking 
>the one which I want.
>I am giving the .ctl file for your reference.dset ^dom3.dat
>options byteswapped
>undef 1.e30
>title OUTPUT FROM WRF V3.2.1 MODEL
>pdef 129 144 lcc 31.672 78.135 65.000 72.500 30.00000 10.00000 77.00000 9000.000 
>9000.000
>xdef 331 linear 71.47070 0.04054054
>ydef 306 linear 25.33539 0.04054054
>zdef 19 levels 
>1000.00000
>950.00000
>900.00000
>850.00000
>800.00000
>750.00000
>700.00000
>650.00000
>600.00000
>550.00000
>500.00000
>450.00000
>400.00000
>350.00000
>300.00000
>250.00000
>200.00000
>150.00000
>100.00000
>tdef 41 linear 00Z26OCT2010 180MN 
>VARS 27
>LU_INDEX 1 0 LAND USE CATEGORY (-)
>U 19 0 x-wind component (m s-1)
>V 19 0 y-wind component (m s-1)
>W 19 0 z-wind component (m s-1)
>PH 19 0 perturbation geopotential (m2 s-2)
>PHB 19 0 base-state geopotential (m2 s-2)
>T 19 0 perturbation potential temperature (theta-t0) (K)
>MU 1 0 perturbation dry air mass in column (Pa)
>MUB 1 0 base state dry air mass in column (Pa)
>P 19 0 perturbation pressure (Pa)
>PB 19 0 BASE STATE PRESSURE (Pa)
>T2 1 0 TEMP at 2 M (K)
>TH2 1 0 POT TEMP at 2 M (K)
>PSFC 1 0 SFC PRESSURE (Pa)
>U10 1 0 U at 10 M (m s-1)
>V10 1 0 V at 10 M (m s-1)
>LANDMASK 1 0 LAND MASK (1 FOR LAND, 0 FOR WATER) (-)
>VEGFRA 1 0 VEGETATION FRACTION (-)
>SST 1 0 SEA SURFACE TEMPERATURE (K)
>HGT 1 0 Terrain Height (m)
>TSK 1 0 SURFACE SKIN TEMPERATURE (K)
>RAINC 1 0 ACCUMULATED TOTAL CUMULUS PRECIPITATION (mm)
>RAINNC 1 0 ACCUMULATED TOTAL GRID SCALE PRECIPITATION (mm)
>XLAT 1 0 LATITUDE, SOUTH IS NEGATIVE (degree_north)
>XLONG 1 0 LONGITUDE, WEST IS NEGATIVE (degree_east)
>TMN 1 0 SOIL TEMPERATURE AT LOWER BOUNDARY (K)
>PBLH 1 0 PBL HEIGHT (m)
>ENDVARSAnd the fortran code which I have used to extract the variables is given 
>below.parameter(nx=129,ny=144,nz=171,ns=30,nt=41,nd=5)
>dimension 
>x(nx,ny,nz),dum(nx,ny)open(12,file='./dom3.dat',form='unformatted',access='direct',recl=nx*ny*4)do
> 5 it = 1,nt
>do 10 iz = 1,nz
>read (12) dum
>do 15 ix = 1,nx
>do 16 iy = 1,ny
>x(ix,iy,iz) = dum(ix,iy)
>16 continue
>15 continue
>10 continuedo 99 k=1,10
>p1(k,it) = x(ix,iy,166)+x(ix,iy,167)write (*,*) p1(k,it)
>99 continue
>stop
>end
>
>I am using a 64 bit Linux Machine with pgf77 compiler.
>Can any one clarify what options byteswapped means?
>Kindly Help Me
>
>Regards
>Sravan
> ############################
>M.SRAVAN KUMAR
>SCIENTIST
>RDC SASE
>CHANDIGARH - 160036
>INDIA
>############################ 
>
>_______________________________________________
>gradsusr mailing list
>gradsusr at gradsusr.org
>http://gradsusr.org/mailman/listinfo/gradsusr
>

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20101028/9a782da1/attachment-0003.html 