[gradsusr] Error setting parea in a loop

[Jeff Duda]

Hey Abayomi,
The problem is with this line:

'define ho=loading.1(z=1,lon=2.375,lon=15.125,lat=2.375,lat=15.125)'

As far as I know, you cannot display a variable using more than one
dimension modification (i.e., you can only alter either z OR lat OR lon,
but not more than one).  Instead, set the dimension environment to fix each
dimension before issuing the define command.

Jeff Duda


On Sat, Apr 13, 2013 at 3:40 PM, Abatan Abayomi <yomiabatan69 at gmail.com>wrote:

> Dear All,
>
> I am trying to make multiple plot (e.g. say 20) of a certain variable
> using parea in a loop. I ran into problem shown below. I need your help.
> =========
> No hardcopy metafile open
> All files closed; all defined objects released;
> All GrADS attributes have been reinitialized
> SET error:  Missing or invalid arguments for PAREA option
> Data Request Error: Invalid dimension expression
>   Attempt to set or modify varying dimension
>   Variable = loading.1, Dimension = 0
>   Error ocurred at column 1
> DEFINE error:  Invalid expression.
> Syntax Error:  Invalid Operand
>   'ho' not a variable or function name
>   Error ocurred at column 1
> DISPLAY error:  Invalid expression
>   Expression = ho
> Hardcopy output file is closed
> ga->
>
> ==========
> My script is shown here:
> 'reinit'
> 'set map 1 1 6'
> 'set font 4'
> 'set mproj scaled'
> *=======================
>
> 'sdfopen trial.nc'
> 'q file'
> 'set display color white'
> 'c'
> 'set grid off'
> 'set strmden 5'
> 'set lat 2.375 15.125'
> 'set lon 2.375 15.125'
> 'set mpdset hires'
> *--------------------
> 'set xlab off'
> *'set xlab on'
> *'set xlevs 4 8 12'
> *--------------------
> 'set ylab on'
> 'set ylevs 4 8 12'
> 'set grads off'
> *-------------------------
> dist_apart = 1.12
> t1 = 0
> row = 1
> ylo = 7.08
> yhi = 8.20
> while(row<=1)
> xlo = 0.42
> xhi = 1.54
> col = 1
> while(col <= 1)
> t1 = t1 + 1
> *'set lat 2.375 15.125'
> *'set lon 2.375 15.125'
> 'set parea 'xlo 'xhi' 'ylo' 'yhi'
> 'set t 't1
> 'define ho=loading.1(z=1,lon=2.375,lon=15.125,lat=2.375,lat=15.125)'
> 'set gxout shaded'
> 'd ho'
> xlo = xhi
> xhi = xhi + dist_apart
> col = col + 1
> endwhile
> yhi = ylo
> ylo = yhi - dist_apart
> row = row + 1
> endwhile
> 'close 1'
> *'draw string 0.25 9.4 m'
> *'run cbarn.gs 0.5 1 0.25 8.7'
> *---------------------------------
>
> Thank you,
> Abayomi
>
> --
>
>
>           Trust in the Lord......
>
> *******************************************************************************************
> *   Abayomi A. Abatan                             *   +1 515 572 4490
> (Home)       *
> *   Geo. and Atmos. Sci.                           *   +1 515 294 1321
> (Office)       *
> *   Iowa State University, Ames                 *   +1 515 203 9020
> (Cell)          *
> *   Iowa. USA.                                          *
>                                *
> *******************************************************************************************
>
>
>
> ****************************************************************************************************
> *   Abayomi A. Abatan                             *   234-802-724-7846
> (out of service)     *
> *   Fed. Univ. of Tech.,                              *   234-803-217-8815
> (out of service)    *
> *   SEMS, PMB 704 Akure                        *
>                                  *
> *   Ondo Nigeria.                                      *
>                                          *
> ****************************************************************************************************
>
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>


-- 
Jeff Duda
Graduate research assistant
University of Oklahoma School of Meteorology
Center for Analysis and Prediction of Storms
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20130413/2d7d28fe/attachment-0003.html