[gradsusr] help with plot

Jose Augusto Paixão Veiga veiga.uea at gmail.com
Fri Jan 13 16:01:27 EST 2012


Dear Hui, 

thank you for the hint. Now, it is working.

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/
-----------------------------------------------------------------------------

On Jan 13, 2012, at 6:52 PM, Hui Wang wrote:

> Hi Jose Augusto,
> 
> You need to set a missing value in the CTL file, which is the value over land.
> 
> Best,
> Hui
> 
> On 1/13/2012 2:36 PM, Jose Augusto Paixão Veiga wrote:
>> 
>> Dear all,
>> 
>> I am facing problems when displaying  a variable which was recorded into a binary file from a netCDF one.
>> 
>> I put in attach the two figures I've saved from the display mode in GrAds. One is related to the original netCDF file and the other denotes the field displayed from the CTL constructed (which I would like to use from now on).
>> 
>> Here are a piece of the script I am using to save a binary file from the netCDF.
>> 
>> 'reinit'
>> 
>> 'sdfopen wap_A1_1950_Jan_to_1999_Dec.nc'
>> 
>> 
>> ***********************************************************
>> 
>> 'set fwrite uvw.bin'
>> 
>> 'set gxout fwrite'
>> 'set x 1 96 '
>> 'set y 1 72 '
>> 
>> 
>> i=1
>> 
>> while (i <= 256)
>> 'set t ' i
>> 
>> j=1
>> 
>> 
>> while (j <= 12)
>> 'set z ' j
>> 
>> 'd omega' 
>> 
>> j=j+1
>> endwhile
>> 
>> i=i+1
>> endwhile
>> 
>> 'disable fwrite'
>> 
>> **********************************************************
>> write (uvw.ctl,'dset   ^uvw.bin ')
>> write (uvw.ctl,'undef  1e+20 ')
>> write (uvw.ctl,'title  Met Office  model output prepared for IPCC Fourth Assessment climate of the 20th Century experiment (20C3M) ')
>> write (uvw.ctl,'xdef 96 linear     1.875  3.75 ')
>> write (uvw.ctl,'ydef 72 linear   -88.75   2.50 ')
>> write (uvw.ctl,'zdef 12 levels  100000 95000 85000 70000 60000 50000 40000 30000 25000 20000 15000 10000 ')
>> write (uvw.ctl,'tdef 252 linear  JAN1969 1mo '   )
>> write (uvw.ctl,'vars  1 ')
>> write (uvw.ctl,'w 12 11,100  vertical velocity     (W)  [ Pa/s]  ')
>> write (uvw.ctl,'endvars')
>> **********************************************************   
>> 
>> 
>> What should be wrong ?
>> 
>> Any help is welcome so much.
>> 
>> 
>> 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/
>> -----------------------------------------------------------------------------
>> 
>> 
>> <Mail Attachment.gif>
>> 
>> <Mail Attachment.gif>
>> 
>> _______________________________________________
>> 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/20120113/9b4b4fd1/attachment-0003.html 


More information about the gradsusr mailing list