[gradsusr] How to read Passive Microwave Water Cycle (PMWC) Dataset

[sushant puranik]

Hello
I am trying to read the Passive Microwave Water Cycle (PMWC) Dataset. Data
file is in .dat format and for that purpose i have made a ctl file. It looks
like

The data files contain six 0.25-degree (1440x720) maps:
1. Water Vapor Transport Speed
2. Water Vapor Transport Direction
3. Water Vapor Transport Divergence
4. Evaporation
5. Precipitation
6. Water Vapor

The maps are one byte with valid geophysical data falling between 0
and 250 with special values for: 252 = Sea Ice, 254 = Insufficient
Data, 255 = Land

The scale factors and offsets are:
1. Speed:         x 2.4,         to range between:  0 to  600 mm m/s
2. Direction:     x 1.5,         to range between:  0 to  360 degrees
3. Divergence:    x 0.024 - 3.0, to range between: -3 to    3 mm/hour
4. Evaporation:   x 0.003,       to range between:  0 to 0.75 mm/hour
5. Precipitation: x 0.012,       to range between:  0 to    3 mm/hour
6. Water Vapor:   x 0.3,         to range between:  0 to   75 mm

The water vapor transport direction is given in the oceanographic convention.

The maps have:
Longitude = 0.25*xdim -  0.125 (degrees east)
Latitdue  = 0.25*ydim - 90.125
where
xdim goes from 1 to 1440 and
ydim goes from 1 to 720.

And my ctl file is as

dset /home/Sushant/Desktop/news_2006_05_pmwc_v01b.dat
title "TESTING"
undef 1e+20
xdef 1440 linear 0.125000 0.25
ydef 720 linear 90.125000 0.25
zdef 1 levels 1
tdef 31 linear 1may2006  24hr
vars 6
WVT-U     WATER VAPOR TRANSPORT ZONAL COMPONENT      (MM M/S)
WVT-V     WATER VAPOR TRANSPORT MERIDIONAL COMPONENT (MM M/S)
WVT-D     WATER VAPOR TRANSPORT DIVERGENCE           (MM/HOUR)
EVAP     EVAPORATION                                (MM/HOUR)
PRCP      PRECIPITATION                              (MM/HOUR)
WVAP      COLUMNAR WATER VAPOR                       (MM)
endvars

but in GrADS it shows following error

ga-> open /home/Sushant/Desktop/wvdata.ctl
Scanning description file:  /home/Sushant/Desktop/wvdata.ctl
Open Error:  Invalid variable record
  --> The invalid description file record is:
  --> wvt-u     water vapor transport zonal component      (mm m/s)

  The data file was not opened.


Can anyone came across such a problem?

suggestions are appreciated



Thank you

Sushant
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100416/c2b24065/attachment-0003.html