[gradsusr] help with plot

[Hui Wang]

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/
> -----------------------------------------------------------------------------
>
>
>
>
> _______________________________________________
> 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/c78a6c61/attachment-0003.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: image/gif
Size: 27555 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20120113/c78a6c61/attachment-0006.gif 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: image/gif
Size: 40486 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20120113/c78a6c61/attachment-0007.gif