[gradsusr] finding differences between two runs with ctl files having different number of variables

[Jennifer Adams]

Do you get the same message if you do not use define, but just display mcape.2-mcape.1? 

On Jan 23, 2013, at 10:46 PM, Jagabandhu Panda wrote:

> Hello Jennifer
>  
> Version of GrADS: v2.0.a9 (little-endian readline printim grib2 netcdf hdf4-sds hdf5 opendap-grids,stn geotiff shapefile)
>  
> YDEF for both files are same
>  
> I am using 'open'
>  
> The dimension variables and other environments set before executing:
> --------------------------------------------------
> ga-> open sing08040112d05nurbhgtah97p5ic.ctl 
> Scanning description file:  sing08040112d05nurbhgtah97p5ic.ctl
> Data file sing08040112d05nurbhgtah97p5ic.dat is open as file 1
> LON set to 103.552 104.116 
> LAT set to 1.1876 1.53311 
> LEV set to 0.01 0.01 
> Time values set: 2008:4:1:12 2008:4:1:12 
> E set to 1 1 
> ga-> open sing08040112d05nurbhgtxxl.ctl 
> Scanning description file:  sing08040112d05nurbhgtxxl.ctl
> Data file sing08040112d05nurbhgtxxl.dat is open as file 2
> ga-> set dfile 1
> Default file set to: sing08040112d05nurbhgtah97p5ic.dat 
> ga-> define mcape1=mcape
> Define memory allocation size = 216720 bytes
> ga-> set dfile 2
> Default file set to: sing08040112d05nurbhgtxxl.dat 
> ga-> define mcape2=mcape
> Define memory allocation size = 216720 bytes
> ga-> d mcape2-mcape1
> Operation error:  Incompatable grids 
>   Dimension ranges aren't equivalent
>   Dimension = 1
>   1st grid range = 1 129   2nd = 1 129 
>   Error ocurred at column 14
> DISPLAY error:  Invalid expression 
>   Expression = mcape2-mcape1
> ----------------------------------------------------------
> As I mentioned earlier, the number of variables are different in both ctl files because they are output from two different model runs and use different model environments.........
>  
> Best regards
> ~Jagabandhu
>  
> On Wed, Jan 23, 2013 at 10:27 PM, Jennifer Adams <jma at cola.iges.org> wrote:
> What version of GrADS are you using?
> What does the YDEF statement look like in both descriptors? 
> Are you using open or xdfopen?
> What is the dimension environment set to before you execute your display of mcape.2-mcape.1?
> --Jennifer
> 
> 
> On Jan 23, 2013, at 5:30 AM, Jagabandhu Panda wrote:
> 
>> Dear users:
>>  
>>      Hope you guys are doing fine. I need your help or suggestions to overcome a problem. I have WRF output from two model runs with same horizontal and vertical resolutions. I have two ctl files for this purpose with corresponding dat files for them. For both of these outputs, I have xdef, ydef, zdef and tdef same for both the files. However, the number of variables are different. One contains 166 variables and other 123. Since the model was tuned to produce more variables in the output file in the earlier case, the number of variables are more. However, in the later case it was not tuned to do so. On the other hand, they do produce same diagnostics after post processing.
>>  
>> When I want to find the difference between say 2dcape field from these two runs, I am unable to do so as it says:
>>  
>> Operation error:  Incompatable grids 
>>   Dimension ranges aren't equivalent
>>   Dimension = 1
>>   1st grid range = 5 117   2nd = 5 117 
>>   Error ocurred at column 16
>> DISPLAY error:  Invalid expression 
>>   Expression = mcape.2-mcape.1
>> EXEC error:  error in metvar.sh.  EXEC stopped.
>> I am unable to point out how to proceed to find the differences each hour for a meteorological variable from the two runs. Could you please suggest some way to overcome it if possible. I will really be thankful to you for this purpose.
>>  
>> Thanks and regards
>> -- Jagabandhu
>>  
>> _______________________________________________
>> 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
> 
> 
> 
> 
> _______________________________________________
> 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

--
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/20130124/6c7d3e7d/attachment-0003.html