[gradsusr] Query Reg Extracting the data in Vertical Levels

[Geeta Geeta]

Dear All, 
I am trying to extract the SPECIFIC HUMIDITY for all the levels in vertical (19 Levels from 1000 to 100hPa) from the Model. 
So I am Using tvrh2q function of Grads. Following is the ctl file. 

'open test.ctl'
*'set lev 1000 100'
'set t 1'
 * while (lev< 100)
*zlev=1000
'set gxout fwrite'
'set fwrite test.dat'
 zlev =1000
'define tvk=tk(lev='zlev')+0.61*QVAPOR(lev='zlev')*tk(lev='zlev')'
'd tvrh2q(tvk(lev='zlev'),rh(lev='zlev'))'
************
 zlev =950
'define tvk=tk(lev='zlev')+0.61*QVAPOR(lev='zlev')*tk(lev='zlev')'
'd tvrh2q(tvk(lev='zlev'),rh(lev='zlev'))' 
*'disable fwrite'
*****************
 zlev =900
'define tvk=tk(lev='zlev')+0.61*QVAPOR(lev='zlev')*tk(lev='zlev')'
'd tvrh2q(tvk(lev='zlev'),rh(lev='zlev'))' 
'disable fwrite'
***********************

* endwhile

I am Unable to Extract the data Using the While Statement in Grads, so I am defining the Zlev separately . 

Kindly Suggest. 
geeta

 		 	   		  
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20101025/86067578/attachment-0003.html