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

[Charles Seman]

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."