[gradsusr] Problem with maskout

Thu Nov 7 20:39:25 EST 2013

Hello Eric,
Hello Kishore,
Thank you very much for your attention and help.
Now the maskout is working fine.
Best wishes.
Celso
____________________________________________
Prof. Celso Bandeira de Melo Ribeiro
Universidade Federal de Juiz de Fora - UFJF
Faculdade de Engenharia, Plataforma 4
Departamento de Engenharia Ambiental e Sanitaria
Bairro Martelos, Juiz de Fora (MG) CEP. 36036-900
tel. (32) 2102-3400
Visiting Scholar in Texas A&M University

On Thu, 7 Nov 2013 19:10:34 -0500 (EST), Eric Altshuler
<ela at cola.iges.org>
wrote:
> Hello Celso, 
> 
> 
> Does your mask data vary in time? If not, there may be only one time
level
> of data in the mask file, in which case you should have 1 instead of 365
in
> your TDEF line, and always use t=1 for the mask variable: 
> 
> 
> 'set t 1' 
> 'd maskout(prec,mask.2(t=1))' 
> 'set t 2' 
> 'd maskout(prec,mask.2(t=1))' 
> ... 
> 
> 
> Best regards, 
> 
> 
> 
> Eric L. Altshuler 
> Research Scientist 
> Center for Ocean-Land-Atmosphere Studies 
> 253 Research Hall, Mail Stop 6C5 
> George Mason University 
> 4400 University Drive 
> Fairfax, VA 22030 USA 
> 
> E-mail: ela at cola.iges.org 
> Phone: (703) 993-5725 
> Fax: (703) 993-5770 
> 
> ----- Original Message -----
> From: "Celso Bandeira de Melo Ribeiro" <celso.bandeira at ufjf.edu.br> 
> To: "GrADS Users Forum" <gradsusr at gradsusr.org> 
> Sent: Thursday, November 7, 2013 6:12:04 PM 
> Subject: Re: [gradsusr] Problem with maskout 
> 
> Hello Kishore, 
> The command "x=maskout(your_var,mask.2(t=1))" works well for the first 
> time. 
> The file.dat don't have problema as I was thinking. 
> But It doesn’t work for times different of 1. 
> See the message of error (attached file): 
> 
> ga-> d maskout(prec,mask.2(t=2)) 
> Low Level I/O Error: Read error on data file 
> Data file name = am_maskaut.dat 
> Error reading 245 bytes at location 306740 
> Data Request Error: Error for variable 'mask.2' 
> Operation Error: Error from maskout function 
> Error ocurred at column 1 
> DISPLAY error: Invalid expression 
> Expression = maskout(prec,mask.2(t=2)) 
> ga-> 
> 
> Is it possible use run maskout for all time (365 days)? 
> Can you help me again? 
> Thank you. 
> Celso 
> 
> 
> On Thu, 7 Nov 2013 09:59:59 +0530, Kishore Ragi <kishoreragi at gmail.com> 
> wrote: 
>> just use the command like the following according to your need .. 
>> 
>> open your file and then open mask file as second ... then use the 
> following 
>> command ... 
>> 
>> x=maskout(your_var,mask.2(t=1)) 
>> 
>> Cheers, 
>> 
>> Kishore 
>> 
>> 
>> On Thu, Nov 7, 2013 at 9:45 AM, Celso Bandeira de Melo Ribeiro < 
>> celso.bandeira at ufjf.edu.br> wrote: 
>> 
>>> Hello Alejandro Meza and grads users, 
>>> My name is Celso and I am trying extract daily precipitation from 
> Merge 
>>> TRMM by INPE (NATIONAL INSTITUTE OF RESEARCH SPATIAL - Brasil) using 
>>> maskout function, same that you did. 
>>> I did one maskout ctl using MetoInfo software, following the same 
>>> proceedings but I don't know how to use this mask in grads. 
>>> The maskout ctl is: 
>>> 
>>> DSET am_maskaut.dat 
>>> TITLE Mask data 
>>> UNDEF -9999 
>>> XDEF 245 LINEAR -82.8 0.2 
>>> YDEF 313 LINEAR -50.2 0.2 
>>> ZDEF 1 LINEAR 1 1 
>>> TDEF 365 LINEAR 1Jan1998 1dy 
>>> VARS 1 
>>> mask 0 99 background mask data 
>>> ENDVARS 
>>> 
>>> Looking forward for any help. 
>>> Thank you, 
>>> Celso 
>>> 
>>> ____________________________________________ 
>>> Celso Bandeira de Melo Ribeiro 
>>> Universidade Federal de Juiz de Fora - UFJF 
>>> Faculdade de Engenharia, Plataforma 4 
>>> Departamento de Engenharia Ambiental e Sanitaria 
>>> Bairro Martelos, Juiz de Fora (MG) CEP. 36036-900 
>>> tel. (32) 2102-3400 
>>> Visiting Scholar in Texas A&M University 
>>> 
>>> 
>>> On Fri, 25 Oct 2013 12:57:30 -0400, Jennifer Adams <jma at cola.iges.org>

>>> wrote: 
>>> > Try using a local override of the time index in your expression,
like 
>>> this: 
>>> > ga-> sst2 = maskout(sst1,mask.2(t=1)) 
>>> > 
>>> > --Jennifer 
>>> > 
>>> > On Oct 25, 2013, at 12:45 PM, Alejandro Meza Mayorga wrote: 
>>> > 
>>> >> Hello grads users. I have netcdf file containing daily sst data
over 
> a 
>>> >> year, and I need to extract the grid data over an irregular area
for 
>>> each 
>>> >> day, so I used the maskout function. 
>>> >> 
>>> >> First I created a shapefile for the area that I needed using ArcGis

>>> >> software (1.png). 
>>> >> 
>>> >> Then I followed the instructions posted in a previous thread by Dr.

>>> >> Yaqiang Wang 
>>> >> (http://gradsusr.org/pipermail/gradsusr/2013-April/034603.html). So

>>> using 
>>> >> MeteoInfo software I made a grads maskout data with the shapefile 
>>> >> (2.png), and I changed the TDEF in the maskout ctl file just as 
> stated 
>>> >> (3.png). 
>>> >> 
>>> >> The maskout works fine for the first day (t=1) of the data (4.png),

>>> >> but 
>>> >> the problem is that when I try to use another 't' different than 1 
> it 
>>> >> gives an error (5.png). 
>>> >> 
>>> >> I hope you could give some ideas of what the error might be. Is it
a 
>>> >> mistake of how to create the shapefile, or an error in the control 
>>> file? 
>>> >> 
>>> >> Looking forward for your help. Thanks in advance. 
>>> >> 
>>> >> Alejandro. 
>>> >> 
>>> >> 
>>> >> 
>>> 
>>> 
>
<1.png><2.png><3.png><4.png><5.png>_______________________________________________
> 
>>> >> gradsusr mailing list 
>>> >> gradsusr at gradsusr.org 
>>> >> http://gradsusr.org/mailman/listinfo/gradsusr 
>>> > 
>>> > -- 
>>> > Jennifer M. Adams 
>>> > Center for Ocean-Land-Atmosphere Studies (COLA) 
>>> > 111 Research Hall, Mail Stop 2B3 
>>> > George Mason University 
>>> > 4400 University Drive 
>>> > Fairfax, VA 22030 
>>> 
>>> -- 
>>> _______________________________________________ 
>>> gradsusr mailing list 
>>> gradsusr at gradsusr.org 
>>> http://gradsusr.org/mailman/listinfo/gradsusr 
>>>

-- 