Re: ECMWF operational analysis

Hi,
Of course there's a problem of definition of lat and long coordinates,
**********************
Variable: T
Type: float
Total Size: 2560000 bytes
           640000 values
Number of Dimensions: 2
Dimensions and sizes: [g50_lat_2 | 800] x [g50_lon_3 | 800]
Coordinates:
Number Of Attributes: 10
 lv_ISBL0 : 500
 center : European Center for Medium-Range Weather Forecasts - Reading
 long_name : Temperature
 units : K
 _FillValue : -999
 level_indicator : 100
 grid_number : 255
 parameter_number : 130
 forecast_time : 0
 initial_time : 03/19/2006 (12:00)
**********************

anyway, further to this, I used a new wxMaximumSize , and I get the same
message:

************
fatal:ContourPlotDraw: Workspace reallocation would exceed maximum size
33554432
fatal:ContourPlotDraw: draw error
************
or
************
fatal:ContourPlotDraw: ARSCAM - MCS TOO SMALL
fatal:ContourPlotDraw: draw error
************
if I increase again the wsMaximumSize.

Thanks,
Andrea

Saji Njarackalazhikam Hameed mi ha detto il 03/23/2006 11:47 PM:

> Andrea,
>
> 1. It will be helpful to print out a summary of the variable you
> are plotting, PrintVarSummary(T), so that one knows how
> (well) the variable T is defined.
>
> 2. Try plotting with
> plot = gsn_contour(wks,T(:,:),res)
>
> If this works, the co-ordinate axes of T are not properly
> described. For instance, longitude should have units of
> degrees_east
> longitude_at_units="degrees_east"
>
> gsn_csm_contour_map?? expects the variable to be properly described
> see:
>
>
> http://www.ncl.ucar.edu/Document/Graphics/Interfaces/gsn_csm_contour_map.shtml
>
>
> I guess the following code snippet may make the plotting with
> gsn_csm_contour_map work
> T!0="latitude"
> T!1="longitude"
> T&longitude=g4_lon_5
> T&longitude_at_units="degrees_east"
> I don't know the name of the latitude co-ordinate
> variable, but follow the above example
>
> 3. To reallocate workspace, get a sample .hluresfile from the
> following link and modify the last entry appropriately:
> http://www.ncl.ucar.edu/Document/Graphics/hlures.shtml
>
> For instance:
> *wsMaximumSize: 32556688
>
> Hope these hints help,
>
> saji
> ....
>
> * Andrea Storto <andrea.storto_at_uniroma1.it> [2006-03-23 13:15:26 +0100]:
>
>
>
>> Hi all,
>>
>
>
>
>
>> maybe it's an old topic, but I couldn't find anything to solve my
>> problems, anyway,
>> I got an ECMWF operational analysis via the MARS archive,
>> (only Z,U,V,Q-rel.hum.-,T for 4 press levels),
>> In my ncl script I print out the variables and the isobaric levels:
>>
>
>
>
>
>> *********************
>> Variable: names
>> Type: string
>> Total Size: 32 bytes
>> 8 values
>> Number of Dimensions: 1
>> Dimensions and sizes: [8]
>> Coordinates:
>> (0) lv_ISBL0
>> (1) g4_lon_5
>> (2) g4_lat_4
>> (3) Q_GDS4_ISBL
>> (4) V_GDS50_ISBL
>> (5) U_GDS50_ISBL
>> (6) T_GDS50_ISBL
>> (7) Z_GDS50_ISBL
>> Variable: lv_ISBL0 (file variable)
>> Type: integer
>> Total Size: 16 bytes
>> 4 values
>> Number of Dimensions: 1
>> Dimensions and sizes: [lv_ISBL0 | 4]
>> Coordinates:
>> lv_ISBL0: [500..1000]
>> Number Of Attributes: 2
>> long_name : isobaric level
>> units : hPa
>> (0) 500
>> (1) 700
>> (2) 850
>> (3) 1000
>> ***********************
>>
>
>
>
>
>> When I try to store a variable and plot it through the functions:
>>
>
>
>
>
>> T = in->T_GDS50_ISBL(0,0,:,:)
>> or
>> T = in->T_GDS50_ISBL
>>
>
>
>
>
>> and (after very simple definition of "res" attributes)
>>
>
>
>
>
>> plot = gsn_csm_contour_map_ce(wks,T(:,:),res)
>> or
>> plot = gsn_csm_contour_map_ce(wks,T(:,:,:,:),res)
>>
>
>
>
>
>> I found two problems:
>>
>
>
>
>
>> 1 - Definition of long and lat:
>>
>
>
>
>
>> *************
>> (0) check_for_y_lat_coord: Warning: Data either does not contain
>> a valid latitude coordinate array or doesn't contain one at all.
>> (0) A valid latitude coordinate array should have a 'units'
>> attribute equal to one of the following values:
>> (0) 'degrees_north' 'degrees-north' 'degree_north' 'degrees
>> north' 'degrees_N' 'Degrees_north' 'degree_N' 'degreeN' 'degreesN'
>> (0) check_for_lon_coord: Warning: Data either does not contain a
>> valid longitude coordinate array or doesn't contain one at all.
>> (0) A valid longitude coordinate array should have a 'units'
>> attribute equal to one of the following values:
>> (0) 'degrees_east' 'degrees-east' 'degree_east' 'degrees
>> east' 'degrees_E' 'Degrees_east' 'degree_E' 'degreeE' 'degreesE'
>> *************
>>
>
>
>
>
>> ...Perhaps... if I print the longitude variable (g4_lon_5) it's
>> ranging from 0 to 360 (all positives)!
>> Something strange or not?
>>
>
>
>
>
>> 2 - Reallocation:
>>
>
>
>
>
>> ****************
>> fatal:ContourPlotDraw: Workspace reallocation would exceed maximum
>> size 16777216
>> fatal:ContourPlotDraw: draw error
>> warning:WorkstationDeactivate: workstation not active or not opened
>> ****************
>>
>
>
>
>
>> I don't know how to solve it!
>>
>
>
>
>
>> Thanks in advance, if the problem was already dealt with, please link
>> the topic!
>>
>
>
>
>
>
>
>
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk_at_ucar.edu
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>
>
>
>

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri Mar 24 2006 - 02:43:33 MST