Re: Coordinate variables must be the same dimension as their dimension

Be careful. SWCF and v_swf_cloud are correct the way they are. The
problem is the way you are subscripting into them. You have the lat
and lon subscripts reversed in many places. nx (144) is associated
with longitudes, and ny (96) is associated with latitudes.

For example, on line 17 you have this:

    v_cloud_yr(ilat,jlon) = avg(v_swf_cloud(ilat,jlon,:))

But I think you mean this:

    v_cloud_yr(ilat,jlon) = avg(v_swf_cloud(jlon,ilat,:))

Now back on line 14, you might want new((/96,144/),double) instead of
new((/144,96/),double). The dim sizes in lines 15 and 16 should be
switched. And so on. Please start at line 14 and make sure every
following line is using subscripts consistently.

--Dave

On Thu, Dec 20, 2012 at 6:19 PM, H.Dang <danghy@gmail.com> wrote:
> Hi Dave,
>
> Thanks for your message! The head of the nc file is short:
>
> dimensions:
> nx = 144 ;
> ny = 96 ;
> time = 12 ;
> variables:
> double lon(nx) ;
> lon:long_name = "longitude" ;
> lon:units = "degrees_east" ;
> double lat(ny) ;
> lat:long_name = "latitude" ;
> lat:units = "degrees_north" ;
> double gw(ny) ;
> gw:long_name = "gauss weights" ;
> double time(time) ;
> time:long_name = "month" ;
> double SWCF(nx, ny, time) ;
> SWCF:long_name = "Short wave cloud forcing" ;
> SWCF:units = "W/m2" ;
> double LWCF(nx, ny, time) ;
> LWCF:long_name = "Long wave cloud forcing" ;
> LWCF:units = "W/m2" ;
>
> So I want to change the "nx" and "ny" in v_swf_cloud (SWCF) to latitude and
> longitude. I don't know if nx and ny influences. Thanks again.
>
> Hongyan
>
> On Thu, Dec 20, 2012 at 7:20 PM, Dave Allured - NOAA Affiliate
> <dave.allured@noaa.gov> wrote:
>>
>> I think you might have reversed the correct meaning of lat and lon
>> subscripts in line 17. What is the result of printVarSummary
>> (v_swf_cloud)?
>>
>> --Dave
>>
>> On Thu, Dec 20, 2012 at 4:35 PM, H.Dang <danghy@gmail.com> wrote:
>> > Hello,
>> >
>> > I met some problem when defining coordinate variables, the error message
>> > is:
>> >
>> > fatal:Coordinate variables must be the same dimension as their dimension
>> > fatal:No coordinate variable exists for dimension (V_lat) in variable
>> > (v_cloud_yr)
>> > fatal:["Execute.c":7556]:Execute: Error occurred at or near line 33 in
>> > file
>> > cloud_cor.ncl
>> >
>> > My code is:
>> >
>> > 1 ;*******************************************
>> > 2 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
>> > 3 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
>> > 4 begin
>> > 5
>> > 6 f1 = addfile("cld_int.nc","r")
>> > 7 V_swf_cloud = f1->SWCF
>> > 8 V_lat = f1->lat
>> > 9 V_lon = f1->lon
>> > 10 v_swf_cloud = V_swf_cloud(:,:,:)
>> > 11 lat_cloud = V_lat(:)
>> > 12 lon_cloud = V_lon(:)
>> > 13
>> > 14 v_cloud_yr = new((/144,96/),double)
>> > 15 do ilat = 0,144-1
>> > 16 do jlon = 0,96-1
>> > 17 v_cloud_yr(ilat,jlon) = avg(v_swf_cloud(ilat,jlon,:))
>> > 18 end do ; jlon
>> > 19 end do ; ilat
>> > 20
>> > 21 lat_cloud!0 = "V_lat"
>> > 22 lat_cloud&V_lat = lat_cloud
>> > 23 lon_cloud!0 = "V_lon"
>> > 24 lon_cloud&V_lon = lon_cloud
>> > 25 printVarSummary(lat_cloud)
>> > 26 printVarSummary(lon_cloud)
>> > 27
>> > 28 printVarSummary(v_cloud_yr)
>> > 29 v_cloud_yr!0 = "V_lat" ;"lat"
>> > 30 v_cloud_yr!1 = "V_lon" ;"lon"
>> > 31 printVarSummary(v_cloud_yr)
>> > 32
>> > 33 v_cloud_yr&V_lat= lat_cloud
>> > 34 v_cloud_yr&V_lon= lon_cloud
>> > 35
>> > 36 end
>> >
>> > Could you please give some suggestion about the line 33? Thanks.
>> >
>> > Hongyan
>> >
>> > P.S. Following is all the information on screen:
>> >
>> > Variable: lat_cloud
>> > Type: double
>> > Total Size: 768 bytes
>> > 96 values
>> > Number of Dimensions: 1
>> > Dimensions and sizes: [V_lat | 96]
>> > Coordinates:
>> > V_lat: [ -90..89.99999999999999]
>> > Number Of Attributes: 2
>> > units : degrees_north
>> > long_name : latitude
>> >
>> > Variable: lon_cloud
>> > Type: double
>> > Total Size: 1152 bytes
>> > 144 values
>> > Number of Dimensions: 1
>> > Dimensions and sizes: [V_lon | 144]
>> > Coordinates:
>> > V_lon: [ 0..357.5]
>> > Number Of Attributes: 2
>> > units : degrees_east
>> > long_name : longitude
>> >
>> > Variable: v_cloud_yr
>> > Type: double
>> > Total Size: 110592 bytes
>> > 13824 values
>> > Number of Dimensions: 2
>> > Dimensions and sizes: [144] x [96]
>> > Coordinates:
>> > Number Of Attributes: 1
>> > _FillValue : 9.969209968386869e+36
>> >
>> > Variable: v_cloud_yr
>> > Type: double
>> > Total Size: 110592 bytes
>> > 13824 values
>> > Number of Dimensions: 2
>> > Dimensions and sizes: [V_lat | 144] x [V_lon | 96]
>> > Coordinates:
>> > Number Of Attributes: 1
>> > _FillValue : 9.969209968386869e+36
>> > fatal:Coordinate variables must be the same dimension as their dimension
>> > fatal:No coordinate variable exists for dimension (V_lat) in variable
>> > (v_cloud_yr)
>> > fatal:["Execute.c":7556]:Execute: Error occurred at or near line 33 in
>> > file
>> > cloud_cor.ncl
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Dec 20 18:51:53 2012