[gradsusr] in netcdf file unable to display 4-D variables

Wed Jul 27 08:22:05 EDT 2011

You did not include the ncdump output that shows the metadata for the  
coordinate axis "lev" -- it may have been insufficient. If that is the  
case, then GrADS won't recognize variables that use that coordinate  
axis. Furthermore, you need to include the version of GrADS you are  
using, since the sdfopen interface has changed over time and this  
information may be quite relevant. After you have sdfopen'ed your  
file, try a 'q file' and see if it has found a z-axis with levels that  
match the ncdump output. If not, then you will need to write a  
descriptor file, either a short version for the xdfopen command, or a  
complete one for the open command. The documentation is at http://iges.org/grads/gadoc/SDFdescriptorfile.html
--Jennifer

On Jul 21, 2011, at 7:37 PM, Charles Seman wrote:

> Jonathan,
>
> This sounds familiar... GrADS may not be able to display both 3D and  
> 4D
> fields if they are in the same netCDF file (don't know why)... One  
> thing
> to try is to extract the 4D fields into a separate netCDF file (e.g.
> using "ncks" if you have it).
>
> Hope this helps,
> Chuck
>
> SMITH, JONATHAN WYNN. (GSFC-6133) wrote:
>> Hello:
>>
>> I have a netcdf file with many variables
>> /    float Q(time, lev, lat, lon) ;
>>        Q:units = "KG/KG" ;
>>    float SA1_VMR_inst(time, lev, lat, lon) ;
>>        SA1_VMR_inst:units = "VMR" ;
>>    float SOA_VMR_inst(time, lev, lat, lon) ;
>>        SOA_VMR_inst:units = "VMR" ;
>>    float T(time, lev, lat, lon) ;
>>        T:units = "K" ;
>>    float TOLUENE_VMR_inst(time, lev, lat, lon) ;
>>        TOLUENE_VMR_inst:units = "VMR" ;
>>    float TROPLEV(time, lat, lon) ;
>>        TROPLEV:units = "KM" ;
>> float U(time, lev, lat, lon) ;
>>        U:units = "M/S" ;
>>    float V(time, lev, lat, lon) ;
>>        V:units = "M/S" ;
>>    char current_mss(nchar) ;
>>        current_mss:long_name = "MSS pathname of this file" ;
>>    int date(time) ;
>>        date:long_name = "current date as 6 digit integer (YYMMDD)" ;
>>    int datesec(time) ;
>>        datesec:long_name = "seconds to complete current date" ;
>>        datesec:units = "s" ;
>>    int days(time) ;
>>        days:long_name = "elapsed simulation days for this case" ;
>>    char first_mss(nchar) ;
>>        first_mss:long_name = "MSS pathname of first file for this  
>> case"/
>>
>> I use the following script:
>>
>> /* B Vs -VTL
>> 'reinit'
>> 'sdfopen
>> /Volumes/JWS/MOZART_OUTPUT/LARGE_DOMAIN/ 
>> mozart4-20110721100736853957.nc
>> <http://mozart4-20110721100736853957.nc>'
>> 'set mpdset hires'
>> 'set display color white'
>> 'run rgbset.gs <http://rgbset.gs>'
>>
>> 'c'
>> "set parea 1 10 1 7.5"
>> "set LON -30 59"
>> "set LAT -12 30"
>> "set LEV 1"
>> 'set grads off'
>> 'set gxout shaded'
>> "set rgb 24 255 120 120"
>> "set clevs 20 40 60 80 100 120 140 160 180 200 220 240"
>> "set ccols 9 14 4 11 5 13 3 10 7 12 8 2 24"
>> 'd days'
>> 'set string 1 l 3'
>> 'draw string 7.5 8.0' _analysis
>> "run cbar.gs <http://cbar.gs>"
>> 'draw title 22 May 2006 MOZART L3\OzoneMixing Ratio (ppbv)'
>> 'printim /Users/jonathanwsmith/MOZART/mozart_22may.gif'/
>>
>> For what ever reason I am only able to plot the one and three
>> dimensional variable like days (1-D) or PS (3-D.  I am am unable to  
>> plot
>> Q or SOA_VMR_inst which are 4-D?
>>
>> What could be going wrong?
>>
>>
>>
>> *Jonathan Wynn Smith*
>>
>> /Student Trainee/
>>
>> _jonathan.smith at nasa.gov_
>>
>> NASA-Goddard Space Flight Center
>>
>> Code 613.3-Atmospheric Chemistry and Dynamics Branch
>>
>> Greenbelt, MD 20771
>>
>> 301-614-5956
>>
>>
>>
>> /Doctoral Candidate/
>>
>> jonathan.smith at bison.howard.edu
>>
>> Howard University Program in Atmospheric Sciences
>>
>> 2355 6th St NW/408 Thirkield Bldg.
>>
>> Washington, DC 20059
>>
>>
>>
>>
>> ------------------------------------------------------------------------
>>
>> _______________________________________________
>> gradsusr mailing list
>> gradsusr at gradsusr.org
>> http://gradsusr.org/mailman/listinfo/gradsusr
>
> -- 
>
> Please note that Charles.Seman at noaa.gov should be considered my NOAA
> email address, not cjs at gfdl.noaa.gov.
>
> ********************************************************************
>  Charles Seman                                Charles.Seman at noaa.gov
>  U.S. Department of Commerce / NOAA / OAR
>  Geophysical Fluid Dynamics Laboratory         voice: (609) 452-6547
>  201 Forrestal Road                              fax: (609) 987-5063
>  Princeton, NJ  08540-6649            http://www.gfdl.noaa.gov/~cjs/
> ********************************************************************
>
> "The contents of this message are mine personally and do not reflect  
> any
> official or unofficial position of the United States Federal  
> Government,
> the United States Department of Commerce, or NOAA."
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr

--
Jennifer M. Adams
IGES/COLA
4041 Powder Mill Road, Suite 302
Calverton, MD 20705
jma at cola.iges.org

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20110727/16dd01a2/attachment-0003.html 