[gradsusr] ARW Post for Grads
Sravan
sravan_sase at yahoo.in
Wed Oct 27 01:08:30 EDT 2010
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)
ENDVARSparameter(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)
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
############################ write (*,*) p1(k,it)And the fortran code which I
have used to extract the variables is given below.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20101027/78bb8464/attachment-0003.html
More information about the gradsusr
mailing list