[gradsusr] Issue Entire Grid Undefined

[James T. Potemra]

Hi Muhammad,

It looks like your first six files are on the same grid (34 lons, 28 
lats) while the seventh file is on a different grid (326,218). GrADS 
will use the definitions in the first open file by default and apply 
that to the subsequent ones.  I think you need to override this by adding

'set dfile 7'

before defining any dimension ranges in the final plot/variable. Note 
you may have to explicitly give the 'set time' and 'set lev' if file 7 
also had different time and altitude dimensions.

Jim

On 4/24/22 4:27 PM, Muhammad Shoaib wrote:
> Hi
>
> i have 7 NETCDF files and i wanted to plot the variable *eir* on 
> portrait virtual page using the following script.
>
> 'set display color white'
> 'c'
> 'set map 1 1 7'
> 'set vpage 0 8.5 0 11'
> 'set digsize 15'
> 'set mpdset mres'
> 'set grads off'
> 'set grid off'
> 'set lat 23 37'
> 'set lon 60 77'
> 'set gxout shaded'
> 'set rbcols 0 14 4 11 5 13 10 7 12'
> 'set clevs 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9'
> 'set mproj latlon'
> 'draw string 1 10.5 GFDL-ESM4  IPSL-CM6A-LR                    
> MPI-ESM2-HR'
> 'draw string 1 7.2 MRI-ESM2 UKESM1-0-LL                     EWEMBI'
> 'draw string 7 3.9 MMEM                                 '
> 'draw string 1 1 Units: (fraction)'
> 'set parea 0.4 2.8 7.5 10.5'
> 'set xlab on'
> 'set ylab on'
> 'set xlopts 1 5 0.115'
> 'set ylopts 1 5 0.115'
> 'set ylint 5'
> 'set xlint 50'
> 'd eir.1'
> 'set ylint 50'
> 'set xlint 50'
> 'set parea 3.2 5.6 7.5 10.5'
> 'd eir.2'
> 'set ylint 50'
> 'set xlint 50'
> 'set parea 6 8.4 7.5 10.5'
> 'd eir.3'
> 'set parea 0.4 2.8 4.2 7.2'
> 'set xlopts 1 2 0.118'
> 'set ylopts 1 2 0.118'
> 'set xlopts 1 5 0.115'
> 'set ylopts 1 5 0.115'
> 'set ylint 5'
> 'set xlint 5'
> 'd eir.4'
> 'set parea 3.2 5.6 4.2 7.2'
> 'set ylint 50'
> 'set xlint 5'
> 'd eir.5'
> 'set ylint 50'
> 'set xlint 5'
> 'set parea 6 8.4 4.2 7.2'
> 'd eir.6'
> 'define MMEM5=(eir.1+eir.2+eir.3+eir.4+eir.5)'
> 'define MMEMeir=(MMEM5/5)'
> 'set ylint 5'
> 'set xlint 5'
> 'set parea 6 8.4 0.9 3.9'
> 'd MMEMeir'
> 'set vpage off'
> 'set vpage 0 8.5 0 11'
> 'set mpdset mres'
> 'set lat 23 37'
> 'set lon 60 77'
> 'set gxout shaded'
> 'set rbcols 0 14 4 11 5 13 10 7 12'
> 'set clevs 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9'
> 'set parea 3.2 5.6 0.9 3.9'
> 'd fig5_eir_region.7'
> 'cbarn'
> 'printim G:\Historical\eir_all_2001-2014.png'
>
>
> but when i display the the last file (file7) using command 'd 
> fig5_eir_region.7' it gives me the following warning and on display 
> "Entire Grid Undefined shows.
>
> Data Request Warning:  Request is completely outside file limits
>   Entire grid contents are set to missing data
>   Grid limits of file:     X = 1 326  Y = 1 218  Z = 1 1  T = 1 1  E = 1 1
>   Grid limits of request:  X = -4 110  Y = 116 210  Z = 1 1  T = 
> 1055579041 1055579041  E = 1 1
>   Warning issued for variable = fig5_eir_region.7
> Cannot contour grid - all undefined values
>
>
> the files and the script is attached with this email.
> if someone knows about this issue kindly guide me i will be very thankful.
>
> THANK YOU
>
> *Muhammad Shoaib*
> *INTERNEE AT CENTER FOR CLIMATE RESEARCH AND DEVELOPMENT *
> *COMSATS UNIVERSITY ISLAMABAD PAKISTAN*
> *
> *
>
> _______________________________________________
> 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/20220424/8e3b476a/attachment.html>