Skip to content

Spectrum

Sub-commands to support spectrum definition

SPECTRUM

Bases: BaseSubComponent

SWAN spectrum specification.

.. code-block:: text

->CIRCLE|SECTOR ([dir1] [dir2]) [mdc] [flow] [fhigh] [msc]

Notes

Directions in the spectra are defined either as a CIRCLE or as a SECTOR. In the case of a SECTOR, both dir1 and dir2 must be specified. In the case of a CIRCLE, neither dir1 nor dir2 should be specified.

At least two of flow, fhigh and msc must be specified in which case the third parameter will be calculated by SWAN such that the frequency resolution df/f = 0.1 (10% increments).

Examples

.. ipython:: python :okwarning:

from rompy_swan.subcomponents.spectrum import SPECTRUM
spec = SPECTRUM(mdc=36, flow=0.04, fhigh=1.0)
print(spec.render())
spec = SPECTRUM(mdc=36, dir1=0, dir2=180, flow=0.04, msc=31)
print(spec.render())
Source code in src/rompy_swan/subcomponents/spectrum.py 
 13
 14
 15
 16
 17
 18
 19
 20
 21
 22
 23
 24
 25
 26
 27
 28
 29
 30
 31
 32
 33
 34
 35
 36
 37
 38
 39
 40
 41
 42
 43
 44
 45
 46
 47
 48
 49
 50
 51
 52
 53
 54
 55
 56
 57
 58
 59
 60
 61
 62
 63
 64
 65
 66
 67
 68
 69
 70
 71
 72
 73
 74
 75
 76
 77
 78
 79
 80
 81
 82
 83
 84
 85
 86
 87
 88
 89
 90
 91
 92
 93
 94
 95
 96
 97
 98
 99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
class SPECTRUM(BaseSubComponent):
    """SWAN spectrum specification.

    .. code-block:: text

        ->CIRCLE|SECTOR ([dir1] [dir2]) [mdc] [flow] [fhigh] [msc]

    Notes
    -----

    Directions in the spectra are defined either as a CIRCLE or as a SECTOR. In the
    case of a SECTOR, both `dir1` and `dir2` must be specified. In the case of a
    CIRCLE, neither `dir1` nor `dir2` should be specified.

    At least two of `flow`, `fhigh` and `msc` must be specified in which case the
    third parameter will be calculated by SWAN such that the frequency resolution
    `df/f = 0.1` (10% increments).

    Examples
    --------

    .. ipython:: python
        :okwarning:

        from rompy_swan.subcomponents.spectrum import SPECTRUM
        spec = SPECTRUM(mdc=36, flow=0.04, fhigh=1.0)
        print(spec.render())
        spec = SPECTRUM(mdc=36, dir1=0, dir2=180, flow=0.04, msc=31)
        print(spec.render())

    """

    model_type: Literal["spectrum", "SPECTRUM"] = Field(
        default="spectrum", description="Model type discriminator"
    )
    mdc: int = Field(
        description=(
            "Number of meshes in theta-space. In the case of CIRCLE, this is the "
            "number of subdivisions of the 360 degrees of a circle so "
            "`dtheta = [360]/[mdc]` is the spectral directional resolution. In the "
            "case of SECTOR, `dtheta = ([dir2] - [dir1])/[mdc]`. The minimum number "
            "of directional bins is 3 per directional quadrant."
        )
    )
    flow: Optional[float] = Field(
        default=None,
        description=(
            "Lowest discrete frequency that is used in the calculation (in Hz)."
        ),
    )
    fhigh: Optional[float] = Field(
        default=None,
        description=(
            "Highest discrete frequency that is used in the calculation (in Hz)."
        ),
    )
    msc: Optional[int] = Field(
        default=None,
        description=(
            "One less than the number of frequencies. This defines the grid "
            "resolution in frequency-space between the lowest discrete frequency "
            "`flow` and the highest discrete frequency `fhigh`. This resolution is "
            "not constant, since the frequencies are distributed logarithmical: "
            "`fi+1 = yfi` where `y` is a constant. The minimum number of frequencies "
            "is 4"
        ),
        ge=3,
    )
    dir1: Optional[float] = Field(
        default=None,
        description=(
            "The direction of the right-hand boundary of the sector when looking "
            "outward from the sector (required for option SECTOR) in degrees."
        ),
    )
    dir2: Optional[float] = Field(
        default=None,
        description=(
            "The direction of the left-hand boundary of the sector when looking "
            "outward from the sector (required for option SECTOR) in degrees."
        ),
    )

    @model_validator(mode="before")
    @classmethod
    def check_direction_definition(cls, data: Any) -> Any:
        """Check that dir1 and dir2 are specified together."""
        dir1 = data.get("dir1")
        dir2 = data.get("dir2")
        if None in [dir1, dir2] and dir1 != dir2:
            raise ValueError("dir1 and dir2 must be specified together")
        return data

    @model_validator(mode="after")
    def check_frequency_definition(self) -> "SPECTRUM":
        """Check spectral frequencies are prescribed correctly."""
        args = [self.flow, self.fhigh, self.msc]
        if None in args:
            args = [arg for arg in args if arg is not None]
            if len(args) != 2:
                raise ValueError("You must specify at least 2 of [flow, fhigh, msc]")
        if self.flow is not None and self.fhigh is not None and self.flow >= self.fhigh:
            raise ValueError("flow must be less than fhigh")
        return self

    @property
    def dir_sector(self):
        if self.dir1 is None and self.dir2 is None:
            return "CIRCLE"
        else:
            return f"SECTOR {self.dir1} {self.dir2}"

    def cmd(self) -> str:
        repr = f"{self.dir_sector} mdc={self.mdc}"
        if self.flow is not None:
            repr += f" flow={self.flow}"
        if self.fhigh is not None:
            repr += f" fhigh={self.fhigh}"
        if self.msc is not None:
            repr += f" msc={self.msc}"
        return repr

Attributes

model_type class-attribute instance-attribute 

model_type: Literal['spectrum', 'SPECTRUM'] = Field(default='spectrum', description='Model type discriminator')

mdc class-attribute instance-attribute 

mdc: int = Field(description='Number of meshes in theta-space. In the case of CIRCLE, this is the number of subdivisions of the 360 degrees of a circle so `dtheta = [360]/[mdc]` is the spectral directional resolution. In the case of SECTOR, `dtheta = ([dir2] - [dir1])/[mdc]`. The minimum number of directional bins is 3 per directional quadrant.')

flow class-attribute instance-attribute 

flow: Optional[float] = Field(default=None, description='Lowest discrete frequency that is used in the calculation (in Hz).')

fhigh class-attribute instance-attribute 

fhigh: Optional[float] = Field(default=None, description='Highest discrete frequency that is used in the calculation (in Hz).')

msc class-attribute instance-attribute 

msc: Optional[int] = Field(default=None, description='One less than the number of frequencies. This defines the grid resolution in frequency-space between the lowest discrete frequency `flow` and the highest discrete frequency `fhigh`. This resolution is not constant, since the frequencies are distributed logarithmical: `fi+1 = yfi` where `y` is a constant. The minimum number of frequencies is 4', ge=3)

dir1 class-attribute instance-attribute 

dir1: Optional[float] = Field(default=None, description='The direction of the right-hand boundary of the sector when looking outward from the sector (required for option SECTOR) in degrees.')

dir2 class-attribute instance-attribute 

dir2: Optional[float] = Field(default=None, description='The direction of the left-hand boundary of the sector when looking outward from the sector (required for option SECTOR) in degrees.')

dir_sector property 

dir_sector

Functions

check_direction_definition classmethod 

check_direction_definition(data: Any) -> Any

Check that dir1 and dir2 are specified together.

Source code in src/rompy_swan/subcomponents/spectrum.py 
 96
 97
 98
 99
100
101
102
103
104
@model_validator(mode="before")
@classmethod
def check_direction_definition(cls, data: Any) -> Any:
    """Check that dir1 and dir2 are specified together."""
    dir1 = data.get("dir1")
    dir2 = data.get("dir2")
    if None in [dir1, dir2] and dir1 != dir2:
        raise ValueError("dir1 and dir2 must be specified together")
    return data

check_frequency_definition 

check_frequency_definition() -> SPECTRUM

Check spectral frequencies are prescribed correctly.

Source code in src/rompy_swan/subcomponents/spectrum.py 
106
107
108
109
110
111
112
113
114
115
116
@model_validator(mode="after")
def check_frequency_definition(self) -> "SPECTRUM":
    """Check spectral frequencies are prescribed correctly."""
    args = [self.flow, self.fhigh, self.msc]
    if None in args:
        args = [arg for arg in args if arg is not None]
        if len(args) != 2:
            raise ValueError("You must specify at least 2 of [flow, fhigh, msc]")
    if self.flow is not None and self.fhigh is not None and self.flow >= self.fhigh:
        raise ValueError("flow must be less than fhigh")
    return self

cmd 

cmd() -> str
Source code in src/rompy_swan/subcomponents/spectrum.py 
125
126
127
128
129
130
131
132
133
def cmd(self) -> str:
    repr = f"{self.dir_sector} mdc={self.mdc}"
    if self.flow is not None:
        repr += f" flow={self.flow}"
    if self.fhigh is not None:
        repr += f" fhigh={self.fhigh}"
    if self.msc is not None:
        repr += f" msc={self.msc}"
    return repr

SHAPESPEC

Bases: BaseSubComponent

Spectral shape specification.

.. code-block:: text

BOUND SHAPESPEC JONSWAP|PM|GAUSS|BIN|TMA PEAK|MEAN DSPR [POWER|DEGREES]

This command BOUND SHAPESPEC defines the shape of the spectra (both in frequency and direction) at the boundary of the computational grid in case of parametric spectral input.

Notes

While technically a component BOUND SHAPESPEC, this is only intended to be used as a subcomponent of the BOUNDSPEC component.

Examples

.. ipython:: python :okwarning:

from rompy_swan.subcomponents.spectrum import SHAPESPEC
shapespec = SHAPESPEC()
print(shapespec.render())
shapespec = SHAPESPEC(
    shape=dict(model_type="tma", gamma=3.1, d=12),
    per_type="mean",
    dspr_type="degrees",
)
print(shapespec.render())
Source code in src/rompy_swan/subcomponents/spectrum.py 
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
class SHAPESPEC(BaseSubComponent):
    """Spectral shape specification.

    .. code-block:: text

        BOUND SHAPESPEC JONSWAP|PM|GAUSS|BIN|TMA PEAK|MEAN DSPR [POWER|DEGREES]

    This command BOUND SHAPESPEC defines the shape of the spectra (both in frequency
    and direction) at the boundary of the computational grid in case of parametric
    spectral input.

    Notes
    -----
    While technically a component `BOUND SHAPESPEC`, this is only intended to be used
    as a subcomponent of the `BOUNDSPEC` component.

    Examples
    --------

    .. ipython:: python
        :okwarning:

        from rompy_swan.subcomponents.spectrum import SHAPESPEC
        shapespec = SHAPESPEC()
        print(shapespec.render())
        shapespec = SHAPESPEC(
            shape=dict(model_type="tma", gamma=3.1, d=12),
            per_type="mean",
            dspr_type="degrees",
        )
        print(shapespec.render())

    """

    model_type: Literal["shapespec", "SHAPESPEC"] = Field(
        default="shapespec", description="Model type discriminator"
    )
    shape: JONSWAP | PM | GAUSS | BIN | TMA = Field(
        default_factory=JONSWAP,
        description="The spectral shape",
    )
    per_type: Literal["peak", "mean"] = Field(
        default="peak",
        description="The type of characteristic wave period",
    )
    dspr_type: Literal["power", "degrees"] = Field(
        default="power",
        description="The type of directional spreading",
    )

    def cmd(self) -> str:
        repr = (
            f"BOUND SHAPESPEC {self.shape.render()} {self.per_type.upper()} "
            f"DSPR {self.dspr_type.upper()}"
        )
        return repr

Attributes

model_type class-attribute instance-attribute 

model_type: Literal['shapespec', 'SHAPESPEC'] = Field(default='shapespec', description='Model type discriminator')

shape class-attribute instance-attribute 

shape: JONSWAP | PM | GAUSS | BIN | TMA = Field(default_factory=JONSWAP, description='The spectral shape')

per_type class-attribute instance-attribute 

per_type: Literal['peak', 'mean'] = Field(default='peak', description='The type of characteristic wave period')

dspr_type class-attribute instance-attribute 

dspr_type: Literal['power', 'degrees'] = Field(default='power', description='The type of directional spreading')

Functions

cmd 

cmd() -> str
Source code in src/rompy_swan/subcomponents/spectrum.py 
331
332
333
334
335
336
def cmd(self) -> str:
    repr = (
        f"BOUND SHAPESPEC {self.shape.render()} {self.per_type.upper()} "
        f"DSPR {self.dspr_type.upper()}"
    )
    return repr

JONSWAP

Bases: BaseSubComponent

Jonswap spectral shape.

.. code-block:: text

JONSWAP [gamma]

Examples

.. ipython:: python :okwarning:

from rompy_swan.subcomponents.spectrum import JONSWAP
shape = JONSWAP(gamma=3.3)
print(shape.render())
Source code in src/rompy_swan/subcomponents/spectrum.py 
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
class JONSWAP(BaseSubComponent):
    """Jonswap spectral shape.

    .. code-block:: text

        JONSWAP [gamma]

    Examples
    --------

    .. ipython:: python
        :okwarning:

        from rompy_swan.subcomponents.spectrum import JONSWAP
        shape = JONSWAP(gamma=3.3)
        print(shape.render())

    """

    model_type: Literal["jonswap", "JONSWAP"] = Field(
        default="jonswap", description="Model type discriminator"
    )
    gamma: float = Field(
        default=3.3,
        description="Peak enhancement parameter of the JONSWAP spectrum.",
        gt=0.0,
    )

    def cmd(self) -> str:
        return f"{super().cmd()} gamma={self.gamma}"

Attributes

model_type class-attribute instance-attribute 

model_type: Literal['jonswap', 'JONSWAP'] = Field(default='jonswap', description='Model type discriminator')

gamma class-attribute instance-attribute 

gamma: float = Field(default=3.3, description='Peak enhancement parameter of the JONSWAP spectrum.', gt=0.0)

Functions

cmd 

cmd() -> str
Source code in src/rompy_swan/subcomponents/spectrum.py 
164
165
def cmd(self) -> str:
    return f"{super().cmd()} gamma={self.gamma}"

TMA

Bases: JONSWAP

TMA spectral shape.

.. code-block:: text

TMA [gamma] [d]

Examples

.. ipython:: python :okwarning:

from rompy_swan.subcomponents.spectrum import TMA
shape = TMA(gamma=2.0, d=18)
print(shape.render())
Source code in src/rompy_swan/subcomponents/spectrum.py 
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
class TMA(JONSWAP):
    """TMA spectral shape.

    .. code-block:: text

        TMA [gamma] [d]

    Examples
    --------

    .. ipython:: python
        :okwarning:

        from rompy_swan.subcomponents.spectrum import TMA
        shape = TMA(gamma=2.0, d=18)
        print(shape.render())

    """

    model_type: Literal["tma", "TMA"] = Field(
        default="tma", description="Model type discriminator"
    )
    d: float = Field(
        description="The reference depth at the wave maker in meters.",
        gt=0.0,
    )

    def cmd(self) -> str:
        return f"{super().cmd()} d={self.d}"

Attributes

model_type class-attribute instance-attribute 

model_type: Literal['tma', 'TMA'] = Field(default='tma', description='Model type discriminator')

d class-attribute instance-attribute 

d: float = Field(description='The reference depth at the wave maker in meters.', gt=0.0)

Functions

cmd 

cmd() -> str
Source code in src/rompy_swan/subcomponents/spectrum.py 
195
196
def cmd(self) -> str:
    return f"{super().cmd()} d={self.d}"

GAUSS

Bases: BaseSubComponent

Gaussian spectral shape.

.. code-block:: text

GAUSS [sigfr]

Examples

.. ipython:: python :okwarning:

from rompy_swan.subcomponents.spectrum import GAUSS
shape = GAUSS(sigfr=0.02)
print(shape.render())
Source code in src/rompy_swan/subcomponents/spectrum.py 
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
class GAUSS(BaseSubComponent):
    """Gaussian spectral shape.

    .. code-block:: text

        GAUSS [sigfr]

    Examples
    --------

    .. ipython:: python
        :okwarning:

        from rompy_swan.subcomponents.spectrum import GAUSS
        shape = GAUSS(sigfr=0.02)
        print(shape.render())

    """

    model_type: Literal["gauss", "GAUSS"] = Field(
        default="gauss", description="Model type discriminator"
    )
    sigfr: float = Field(
        description=(
            "Width of the Gaussian frequency spectrum expressed "
            "as a standard deviation in Hz."
        ),
        gt=0.0,
    )

    def cmd(self) -> str:
        return f"{super().cmd()} sigfr={self.sigfr}"

Attributes

model_type class-attribute instance-attribute 

model_type: Literal['gauss', 'GAUSS'] = Field(default='gauss', description='Model type discriminator')

sigfr class-attribute instance-attribute 

sigfr: float = Field(description='Width of the Gaussian frequency spectrum expressed as a standard deviation in Hz.', gt=0.0)

Functions

cmd 

cmd() -> str
Source code in src/rompy_swan/subcomponents/spectrum.py 
229
230
def cmd(self) -> str:
    return f"{super().cmd()} sigfr={self.sigfr}"

PM

Bases: BaseSubComponent

Pearson-Moskowitz spectral shape.

.. code-block:: text

PM

Examples

.. ipython:: python :okwarning:

from rompy_swan.subcomponents.spectrum import PM
shape = PM()
print(shape.render())
Source code in src/rompy_swan/subcomponents/spectrum.py 
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
class PM(BaseSubComponent):
    """Pearson-Moskowitz spectral shape.

    .. code-block:: text

        PM

    Examples
    --------

    .. ipython:: python
        :okwarning:

        from rompy_swan.subcomponents.spectrum import PM
        shape = PM()
        print(shape.render())

    """

    model_type: Literal["pm", "PM"] = Field(
        default="pm", description="Model type discriminator"
    )

Attributes

model_type class-attribute instance-attribute 

model_type: Literal['pm', 'PM'] = Field(default='pm', description='Model type discriminator')

BIN

Bases: BaseSubComponent

Single frequency bin spectral shape.

.. code-block:: text

BIN

Examples

.. ipython:: python :okwarning:

from rompy_swan.subcomponents.spectrum import BIN
shape = BIN()
print(shape.render())
Source code in src/rompy_swan/subcomponents/spectrum.py 
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
class BIN(BaseSubComponent):
    """Single frequency bin spectral shape.

    .. code-block:: text

        BIN

    Examples
    --------

    .. ipython:: python
        :okwarning:

        from rompy_swan.subcomponents.spectrum import BIN
        shape = BIN()
        print(shape.render())

    """

    model_type: Literal["bin", "BIN"] = Field(
        default="bin", description="Model type discriminator"
    )

Attributes

model_type class-attribute instance-attribute 

model_type: Literal['bin', 'BIN'] = Field(default='bin', description='Model type discriminator')