ncl-metno home

shell scripts:


layersection.sh

Syntax

Example: Meridional cross-section of salinity

The figure below was made by issuing the command:

layersection.sh hycom_031_h.nc hycom_031_TS.nc thknss salin lat 62 145

layersection.sh: syntax

layersection.sh --help

        layersection.sh  /  ncl-metno 1.2

>>>
>>>
>>> NOTE: This script assumes that the variable to depict has been
>>>  stored with dimensions in the order x  -y  -layer(-time)
>>>                                   or lon-lat-layer(-time)
>>>  (w/ Fortran style sequence of dimensions).
>>>
>>> Syntax:
>>> =======
>>>
>>>  layersection.sh (layer) <hfile> <varfile> <hname> <varname>\
>>>                                            <dimname> <node> <time>
>>> where
>>>    (layer)     include an initial argument 'layer' if you want isopleths
>>>                 to be drawn for layer interfaces, ignore otherwise
>>>    <hfile>     name of netCDF file w/thickness results
>>>    <varfile>   name of netCDF file w/requested variable
>>>                 you may use '.' if both variables are on the same file
>>>    <hname>     name of thickness variable
>>>    <varname>   name of variable to depict
>>>             NOTE! This script requires the existence of an
>>>                   attribute 'missing_value' to <varname>
>>>    <dimname>   name of crossection's horizontal dimension
>>>                 e.g.,  lat  for a lat-z (meridional) crossection
>>>    <node>      node of non-depicted dimension
>>>                 if <dimname> is  lat , this is the  lon  grid no.
>>>    <time>      time step no. to depict (use  0  for x-y-z & lon-lat-z fields)
>>>                a negative <time> value is interpreted as a flag that
>>>                 stops the ncl script from being deleted
>>>
>>>  The script will produce an eps-file and a png-file.
>>>
>>>
>>> User specifications:
>>> ====================
>>>
>>>  By copying the default spec.s from