Re: (no subject)

This is offline.

My speculation is that you are encountering:

[1] a memory limit in your computer

and/or

[2] you do not have netCDF's Large File Support [LFS} on your computer.
     This allows the netCDF software to spane files larger than 2GB

====
Also, if you are reading or accessing whole arrays, do not us

    u0=in->U(:,:,:,:)

It is better to use

    u0=in->U

It is faster.

=====
You do not *have to* rename dimensions

    u0!0="time"
    u0!1="lev"
    u0!2="lat"
    u0!3="lon"

the WRF files have dimension names

    U ( Time, bottom_top, south_north, west_east_stag )
    V ( Time, bottom_top, south_north_stag, west_east )

  PH ( Time, bottom_top_stag, south_north, west_east )
    T ( Time, bottom_top, south_north, west_east )

=====
so
    u1=u0(lev|:,lat|:,lon|:,time|:)

becomes
    u1 = u0(bottom_top|:,south_north|:, west_east_stag|:,Time|:)
======
You are reading all of your data into memory .... that is not good
programming style. Also, you read into u0, then you reorder but
you do not delete u0. It is no longer used

  u1 = u0(bottom_top|:,south_north|:, west_east_stag|:,Time|:)
  delete(u0) ; no longer used

Do this after each reorder to reduce memory.

There are many other instances where you use a lot of memory unnecessarily.
Please look at your code.
=====
WRF grids are curvilinear not rectilinear. Averaging over "Yime" is ok but
averaging over latitude and longitude is not.

Gang Liu wrote:
> Hello,
>
>
>
> I have run WRF model and then I use the NCL software for post-processing
> the WRF results output in netCDF format. The following are the command I
> executed and the error message that followed:
>
>
>
> ncl nclpost.ncl
>
> Copyright (C) 1995-2007 - All Rights Reserved
>
> University Corporation for Atmospheric Research
>
> NCAR Command Language Version 5.0.0
>
> The use of this software is governed by a License Agreement.
>
> See http://www.ncl.ucar.edu/ for more details.
>
> Segmentation fault (core dumped)
>
>
>
> The type of the operating system is ‘Linux mnode04 2.6.9-55.ELsmp #1 SMP
> Fri Apr 20 16:36:54 EDT 2007 x86_64 x86_64 x86_64 GNU/Linux’. The
> attached is the NCL script file named ‘nclpost.ncl’.
>
>
>
> The data file named ‘wrfout_d01_0001-01-01_09:00:00’ is too large to be
> uploaded, which has 279421216144 bytes. For the original
> four-dimensional data arrays read in from the WRF result, the subscripts
> and dimensions are ‘time’ (nt = 1800), ‘lev’ (nz = 201), ‘lat’ (ny =
> 121) and ‘lon’ (nx = 121) from left to right.
>
>
>
> Could you please provide me any help in solving the difficulty?
>
>
>
> Thank you very much.
>
>
>
> Gang
>
>
>
>
>
> ************************************
> Gang Liu
> Department of Atmospheric Sciences
> Nanjing University
> Nanjing 210093
> P. R. China
> Phone: 86-25-83595185 (o)
> Email: gangliu_at_nju.edu.cn <mailto:gangliu_at_nju.edu.cn>
> ************************************
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Nov 27 2008 - 10:47:48 MST