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

[Jagabandhu Panda]

Dear Jennifer and Amulya:

        As Amulya suggested to use lterp command, I was just trying to
learn how to use that. And just now with this command which Jennifer
suggested, it looks like working! I think this should serve my purpose
hopefully..........

Best regards
~Jagabandhu
On Thu, Jan 24, 2013 at 8:29 PM, Jennifer Adams <jma at cola.iges.org> wrote:

> Also, you can try 'd mcape.2-lterp(mcape.1,mcape.2)'. --Jennifer
>
>
> On Jan 24, 2013, at 6:39 AM, Jennifer Adams wrote:
>
> 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
>
>
>
> _______________________________________________
> 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
>
>


-- 
###########################################################################
Dr. Jagabandhu Panda, PhD
School of Physical and Mathematical Sciences
Nanyang Technological University
Singapore
Phone: +65 - 65922565
Email: jpanda at ntu.edu.sg, jagabandhu at gmail.com
Personal Web Page: http://jagabandhupanda.wordpress.com/
Group Home Page:
http://www1.spms.ntu.edu.sg/~sunshine/WeatherWeb/html/index.html
###########################################################################
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20130124/9269e28d/attachment-0003.html