[gradsusr] about mask.gs

[luan jia]

Hi Matt,
I just ga->run d:\mask.gs
and then got the erro messages.

Now,I made some small changes to it.

'reinit'
'set rgb 90 100 100 100 '
'set rgb 91  50 50 50 '
'set rgb 92 200 200 200'
_fnum = 0
*  put input file name into the ofile argument
rc = ofile('d:\model.ctl')
*  Select gridpoints for the mask
'set mpdset mres'
rc = newar(_fnum)
'set gxout fwrite'

*  output file name as the fwrite operand
'set fwrite d:\maskarea.dat'
  fn = _fnum
  'set x 0.5 '%(_xnum.fn+0.5)
  'set y 0.5 '%(_ynum.fn+0.5)
  'set z 1'
  'set t 1'
'd maskarea'
say result
'disable fwrite'

It works now.
But I have got some new problems:
1:the size of maskarea.dat is (x+2)*(y+2)=93*48,(input file:model.dat
x*y=91*46),why?
2:The mask data I selected from the interface is not equal to
maskarea.dat,how to write  the .ctl file for maskarea.dat(the XDEF,YDEF
part).

thanks in advance

best regards

luan

On Sun, Aug 1, 2010 at 9:12 AM, <matt.alonso at gmail.com> wrote:

> Hi luan
>
> How are you calling the script?
>
> Matt
> Sent from my Verizon Wireless BlackBerry
>
> -----Original Message-----
> From: luan jia <luanjab at gmail.com>
> Sender: gradsusr-bounces at gradsusr.org
> Date: Sun, 1 Aug 2010 09:08:27
> To: <gradsusr at gradsusr.org>
> Reply-To: GrADS Users Forum <gradsusr at gradsusr.org>
> Subject: [gradsusr]  about mask.gs
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
> _______________________________________________
> 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/20100801/f2440cd5/attachment-0003.html