[ncl-talk] Fw: read and plot multiple files with different sizes

Dennis Shea shea at ucar.edu
Tue Dec 12 17:21:45 MST 2017


NCL is a strongly typed language. There are certain operations that are not
allowed.
---------------------
By default, NCL is a strongly typed language. It will not allow certain
operations.

 You can use the NCL *reassignment operator  *  *:=*

Consider:

float  a(10)
double b(10)
float c(30)

By default:
    a  = b     ; NCL will issue a FATAL error; it will not allow a double=>
float because information may be lost
    a := b     ; NCL will truncate

    a  = c     ; FATAL ... it will not allow an array of size 30 to be
reassigned to an array that is size 10
    a := c     ; NCL will overwrite

So I suggest going over your script and anyplace you this a size and or
type might change ... use the := syntax.I

On Tue, Dec 12, 2017 at 3:41 PM, Ehsan Taghizadeh <ehsantaghizadeh at yahoo.com
> wrote:

>
> Hi,
> First I want to say thank you for your help about my previous problem. I
> think this committee is really fast in helping and answering.
> However I have another problem and I'll be thankful if I could have your
> help again.
> I want to draw SMAP (Soil Moisture Active Passive) data and I use a
> provided script at following link:
> https://www.ncl.ucar.edu/Applications/Scripts/smap_l3_2.ncl
> which the parent directory is:
> NCL Applications: HDF, HDF-EOS
> <https://www.ncl.ucar.edu/Applications/HDF.shtml>
>
> NCL Applications: HDF, HDF-EOS
>
> Use of NCL to visualize datasets in HDF format.
> <https://www.ncl.ucar.edu/Applications/HDF.shtml>
>
>
> This script helps to draw multiple files, but after one file I got
> following error:
>       fatal:Dimension sizes on right hand side of assignment do not match
> dimension sizes of left hand side
>       fatal:["Execute.c":8640]:Execute: Error occurred at or near line 80
> in file smap_l3_2.ncl
>
> which mentioned different sizes of my files. For example for one file I
> have:
>       float latitude(281024);
> and for another file I have:
>       float latitude(281711);
> and the error arises when the script wants to do following command:
>      lat2d     = f[nf]->$lat_path$
> also for
>      lon2d     = f[nf]->$lat_path$
>
> Attachment is the script and following lines are output of
> printVarSummary(var):
>
>
> Variable: var
> Type: float
> Total Size: 1124096 bytes
>             281024 values
> Number of Dimensions: 1
> Dimensions and sizes:   [DIM_000 | 281024]
> Coordinates:
> Number Of Attributes: 6
>   _FillValue :  -9999
>   coordinates : /Soil_Moisture_Retrieval_Data/latitude
> /Soil_Moisture_Retrieval_Data/longitude
>   long_name :   Representative soil moisture measurement for the Earth
> based grid cell.
>   units :       cm**3/cm**3
>   valid_max :   0.5
>   valid_min :   0.02
> (0)     Representative soil moisture measurement for the Earth based grid
> cell. (cm**3/cm**3) : min=0.02   max=0.85494
> (0)     nf=1; SMAP_L2_SM_P_E_11167_D_20170305T030914_R14010_001.h5
> fatal:Dimension sizes of left hand side and right hand side of assignment
> do not match
> fatal:["Execute.c":8640]:Execute: Error occurred at or near line 80 in
> file smap_l3_2.ncl
>
>
> So how could I read and plot multiple files with different sizes?
>
> *Sincerely*
> *Ehsan*
>
>
>
>
>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20171212/14ff907b/attachment.html>


More information about the ncl-talk mailing list