[RegCNET] Grid problem!

[Ho Thi Minh Ha]

Dear Monica,

Maybe there's resolution difference between two runs. So you can try

set lat "lat1" lat2"
set lon "lon1' "lon2"
define a1=ave(ave(aext8,z=1,z=18)*18,t=1,t=4)
define a2=ave(ave(aext8,z=1,z=18)*18,t=1,t=4)
a1new=regrid2(a1,0.5)      ! 0.5 is resolution of new grid
a2new=regrid2(a2,0.5)
d a1new-a2new

If it still doesn't work I think you should fwrite the data of the file 1 in binary format and
then interpolate into the grid of the file 2. Because of the slight difference between 2
resolutions you just need do that base on the nearest point interpolation method.

Best wishes,
Ha

> Dear Faisal and Ha,
> thank your for your reply...
>
> I tested setting
>  set lat "lat1" lat2"
> set lon "lon1' "lon2"
>
> before defining a1=ave(ave(aext8,z=1,z=18)*18,t=1,t=4) and
> a2=ave(ave(aext8,z=1,z=18)*18,t=1,t=4) but when I make
>
> d a1-a2 I have:
>
> Operation error:  Incompatable grids
>   Dimension ranges aren't equivalent
>   Dimension = 0
>   1st grid range = 199 570   2nd = 183 554
>   Error ocurred at column 6
> DISPLAY error:  Invalid expression
>   Expression = a2-a1
>
>
> Shoul you suggest me what I have to do?
>
> Thanks again
>
> Monica
>
>
> CMCC - Physics Department - University of Salento
> via Arnesano - 73100 Lecce, Italy
>
> e-mail: monica.santese at le.infn.it
> phone:  +39-0832-297502 (office)
>              +39-0832-297481 (laboratory)
> fax:       +39 -0832-297505
> ******************************
> ----- Original Message -----
> From: "Ho Thi Minh Ha" <hahtm at vnu.vn>
> To: "Monica Santese" <monica.santese at le.infn.it>
> Cc: "regcnet" <regcnet at lists.ictp.it>
> Sent: Monday, July 05, 2010 3:52 PM
> Subject: Re: [RegCNET] Grid problem!
>
>
>> Dear Monica Santese,
>>
>> You should try this:
>>
>> 'open CHE.YYYYMMDD0100'  (file1)
>> 'set lon lonstart lonend'         (e.g., lonstart, lonend are values of
>> start and end longitude,
>> respecitively)
>> 'set lat latstart latend'
>> 'define a1=ave(ave(aext8,z=1,z=18)*18,t=1,t=4)'
>> 'close 1'
>>
>> 'open CHE.YYYYMMDD0100'  (file2)
>> 'set lon lonstart lonend'         (the same region as above)
>> 'set lat latstart latend'
>> 'define a2=ave(ave(aext8,z=1,z=18)*18,t=1,t=4)'
>> 'd a1-a2'
>>
>> 'close 1'
>>
>> Hope it would work
>>
>> Best wishes,
>> Ha
>>
>>
>>> Dear users,
>>>
>>> I have two different files CHE that come from two different runs with
>>> different numbers of pixels:
>>> One for iy=200 and jx = 200    and one for iy=185 and jx = 200.
>>>
>>> Now, if I want to plot the aerosol optical depth map of the differences
>>> between the AOD from the
>>> first run and AOD from the second one using :
>>> a1=ave(ave(aext8,z=1,z=18)*18,t=1,t=4) and
>>> a2=ave(ave(aext8,z=1,z=18)*18,t=1,t=4)
>>> dif=a1-a2
>>>
>>> I have the error message:
>>>
>>>
>>> Operation error:  Incompatable grids
>>>   Dimension ranges aren't equivalent
>>>   Dimension = 0
>>>   1st grid range = -15 756   2nd = 1 771
>>>   Error ocurred at column 6
>>> DEFINE error:  Invalid expression.
>>>
>>>
>>> Is there a way to plot the AOD differences for a lat-lon area that are
>>> common to the two files
>>> CHE?
>>>
>>> Thank you fory your help!!
>>>
>>>
>>> PhD.Monica Santese
>>> CMCC - Physics Department - University of Salento
>>> via Arnesano - 73100 Lecce, Italy
>>>
>>> e-mail: monica.santese at le.infn.it
>>> phone:  +39-0832-297502 (office)
>>>              +39-0832-297481 (laboratory)
>>> fax:       +39 -0832-297505
>>> ******************************_______________________________________________
>>> RegCNET mailing list
>>> RegCNET at lists.ictp.it
>>> https://lists.ictp.it/mailman/listinfo/regcnet
>>
>>
>> Ho Thi Minh Ha, PhD.
>>
>> Department of Meteorology, Hanoi University of Sciences
>> 334 Nguyen Trai street, Thanh Xuan district, Hanoi city, VIETNAM
>> Phone: +84435583811
>> Email: hahtm at vnu.vn
>> or      ha_meteo at yahoo.com
>>
>> _______________________________________________
>> RegCNET mailing list
>> RegCNET at lists.ictp.it
>> https://lists.ictp.it/mailman/listinfo/regcnet
>>
>>
>
> _______________________________________________
> RegCNET mailing list
> RegCNET at lists.ictp.it
> https://lists.ictp.it/mailman/listinfo/regcnet
>


Ho Thi Minh Ha, PhD.

Department of Meteorology, Hanoi University of Sciences
334 Nguyen Trai street, Thanh Xuan district, Hanoi city, VIETNAM
Phone: +84435583811
Email: hahtm at vnu.vn
or      ha_meteo at yahoo.com