[gradsusr] Problem to plot a figure in GrAds from gaussian latitude grid

Ricardo Hallak hallak at model.iag.usp.br
Fri Dec 23 10:53:05 EST 2011


José,

every time I face "strange" problems like that, I find the source of error is some kind of typo in the ctl file. Maybe an "invisible" Crtl^M in some place of the .ctl file. It happens when a .ctl file is transfered from a Windows system to a Linux system with the ftp or similar service.

So, I suggest you to rewrite the string 73.9475 since the very first blank space until the last one blank space (that is, the entire field of this number). Remove the field and rewrite it (using vi would be better).

Actually, the number 73.9475 in your .ctl seems to be the only one with 2 blank spaces after the number, when compared to the others.

If it doesn't work, the hardiest other suggestion is to rewrite the .ctl again in a new and separated file. Don't use copy/paste.

Regards,
Ricardo

On Fri, 23 Dec 2011 08:18:02 -0400, Jose Augusto Paixão Veiga wrote
> Ok, thank you Jeff. I'll wait for Arlindo or Jennifer possible contact. I believe they know how to fix it.
> 
>  
> 
> 
> 
> 
> 
> 
> José Augusto P. Veiga, 
> 
> ======================================
> Universidade do Estado do Amazonas 
> Departamento de Meteorologia
> Escola Superior de Tecnologia (EST)
> -----------------------------------------------------------------------------
> Av. Darcy Vargas, 1200, Manaus-AM Brasil
> Work phone:    (92) 3878 4333
> Skype: veiga_j.a.p.
> 
> -----------------------------------------------------------------------------
> CV: http://lattes.cnpq.br/4027612512091565
> URL:http://scientificmet.wordpress.com/
> http://nmet.wordpress.com/
> -----------------------------------------------------------------------------
> 
> On Dec 22, 2011, at 11:41 PM, Jeff Duda wrote:Jose,
> No, I don't know any other tricks. Like I said, this sounds like a parsing error when Grads reads the control file.  Arlindo or Jennifer will probably have the answer.
> 
> Jeff
> 
> On Thu, Dec 22, 2011 at 5:20 PM, Jose Augusto Paixão Veiga <veiga.uea at gmail.com> wrote:
> 
> Yes Jeff.
> I tried use a different number like 73.9480 and GrAds split it into 73 and .9480. So, I used 73.0 instead of the original value of 73.9475 (after your hint) and GrAds produced the values 73 and .0
> 
> I really don't know what could be happening.
> 
> Any other hint Jeff? 
> 
> Thanks,
> 
> 
> 
> 
> 
> 
> José Augusto P. Veiga, 
> 
> ======================================
> Universidade do Estado do Amazonas 
> Departamento de Meteorologia
> Escola Superior de Tecnologia (EST)
> -----------------------------------------------------------------------------
> Av. Darcy Vargas, 1200, Manaus-AM Brasil
> Work phone:    (92) 3878 4333
> Skype: veiga_j.a.p.
> 
> -----------------------------------------------------------------------------
> CV: http://lattes.cnpq.br/4027612512091565
> URL:http://scientificmet.wordpress.com/
> http://nmet.wordpress.com/
> -----------------------------------------------------------------------------
> 
> On Dec 22, 2011, at 6:59 PM, Jeff Duda wrote:
> Jose,
> Have you tried using other numbers that are similar to 73.9475 but have fewer digits?  It sounds like some odd parsing error.
> 
> Jeff Duda
> 
> On Thu, Dec 22, 2011 at 4:11 PM, Jose Augusto Paixão Veiga <veiga.uea at gmail.com> wrote:
> 
> 
> Dear all,
> 
> 
> I am facing the following problem in GrAds:
> 
> 
> After producing my CTL file and plotted the results I get surprised with plot result (Figure in attach).   Note that the fifth point from the last latitude is shifted for the equatorial latitude. So, I decided to use the command q ctlinfo to get some information about the latitudes ingested by GrAds after I've been opened the CTL file. Them, I realize that GrAds "make" some confusing with one specific latitude point.   Pay attention,
> 
> 
> 
> This is a part of the result from the command q ctlinfo :
> 
> 
> 
> xdef 1 linear 1 1
> ydef 64 levels -87.8638 -85.0965 -82.3129 -79.5256 -76.7369 -73.9475 -71.1577 -68.3678
>  -65.5776 -62.7873 -59.997 -57.2066 -54.4162 -51.6257 -48.8352 -46.0447 -43.2542 -40.4636
>  -37.6731 -34.8825 -32.0919 -29.3014 -26.5108 -23.7202 -20.9296 -18.139 -15.3484 -12.5578
>  -9.76715 -6.97653 -4.18592 -1.39531 1.39531 4.18592 6.97653 9.76715 12.5578 15.3484
>  18.139 20.9296 23.7202 26.5108 29.3014 32.0919 34.8825 37.6731 40.4636 43.2542
>  46.0447 48.8352 51.6257 54.4162 57.2066 59.997 62.7873 65.5776 68.3678 71.1577
>  73 0.9475 76.7369 79.5256 82.3129 85.0965
> zdef 1 linear 1 1
> 
> 
> Thus, GrAds separates one latitude value into two others. I told this because the latitude point 73.9475 is duplicated into the points 73 and 0.9475. When I change the value for the integer 73 everything goes well and the value is not shifted for the "center of figure". However I wouldn't like to to this all the time I face this.
> 
> 
> I would appreciate any hint and help.  I don't know what is producing this split in the latitude value.
> 
> 
> 
> 
> This is part of my ctl file 
> 
> 
> 
> title CNRM  model output prepared for IPCC Fourth Assessment climate of the 20th Century experiment (20C3M)
> xdef 1 linear 1 1
> ydef  64 levels -87.8638 -85.0965 -82.3129 -79.5256 -76.7369 -73.9475 -71.1577 -68.3678 -65.5776 -62.7873 -59.997 -57.2066 -54.4162 -51.6257 -48.8352 -46.0447 -43.2542 -40.4636 -37.6731 -34.8825 -32.0919 -29.3014 -26.5108 -23.7202 -20.9296 -18.139 -15.3484 -12.5578 -9.76715 -6.97653 -4.18592 -1.39531 1.39531 4.18592 6.97653 9.76715 12.5578 15.3484 18.139 20.9296 23.7202 26.5108 29.3014 32.0919 34.8825 37.6731 40.4636 43.2542 46.0447 48.8352 51.6257 54.4162 57.2066 59.997 62.7873 65.5776 68.3678 71.1577 73.9475  76.7369 79.5256 82.3129 85.0965 87.8638 
> zdef 1 linear 1 1
> tdef 252 linear JAN1979 1MO
> vars  17
> gz 0 33,100 generation of available potential energy
> 
> 
> 
> 
> Cheers,
> 
> 
> 
> 
> 
> José Augusto P. Veiga, 
> 
> ======================================
> Universidade do Estado do Amazonas 
> Departamento de Meteorologia
> Escola Superior de Tecnologia (EST)
> -----------------------------------------------------------------------------
> Av. Darcy Vargas, 1200, Manaus-AM Brasil
> Work phone:    (92) 3878 4333
> Skype: veiga_j.a.p.
> 
> -----------------------------------------------------------------------------
> CV: http://lattes.cnpq.br/4027612512091565
> URL:http://scientificmet.wordpress.com/
> http://nmet.wordpress.com/
> -----------------------------------------------------------------------------
> 
> 
> _______________________________________________
> 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
> 
> _______________________________________________
> 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
> 
> 
> 
> -- 
> Jeff Duda
> Graduate research assistant
> University of Oklahoma School of Meteorology
> Center for Analysis and Prediction of Storms
> 
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>

Ricardo Hallak

Department of Atmospheric Sciences (DCA) 
Institute of Astronomy, Geophysics and Atmospheric Sciences (IAG) 
University of São Paulo (USP) 
Rua do Matão, 1226 - Cidade Universitária - São Paulo SP 05508-900 
Phone.: +55 (11) 3091-2850; 3091-4703 
Fax: +55 (11) 3091-4714

 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20111223/3e293e1d/attachment-0003.html 


More information about the gradsusr mailing list