Skip to content

HmsUtils

Utility functions for HMS Commander.

hms_commander.HmsUtils

HmsUtils - Utility Functions for HEC-HMS Operations

This module provides utility functions for common HEC-HMS operations including file management, unit conversions, and data validation.

All methods are static and designed to be used without instantiation.

HmsUtils

Utility functions for HEC-HMS operations.

Provides helper methods for file operations, unit conversions, time interval handling, and data validation.

All methods are static - no instantiation required.

Example

from hms_commander import HmsUtils interval_min = HmsUtils.parse_time_interval("15 Minutes") print(interval_min) # 15

Source code in hms_commander/HmsUtils.py
 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
134
135
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
166
167
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
197
198
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
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
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
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
class HmsUtils:
    """
    Utility functions for HEC-HMS operations.

    Provides helper methods for file operations, unit conversions,
    time interval handling, and data validation.

    All methods are static - no instantiation required.

    Example:
        >>> from hms_commander import HmsUtils
        >>> interval_min = HmsUtils.parse_time_interval("15 Minutes")
        >>> print(interval_min)  # 15
    """

    # Unit conversion factors - reference centralized constants
    CONVERSION_FACTORS = {
        # Length
        'in_to_mm': INCHES_TO_MM,
        'mm_to_in': MM_TO_INCHES,
        'ft_to_m': FEET_TO_METERS,
        'm_to_ft': METERS_TO_FEET,

        # Area
        'sqmi_to_sqkm': SQMI_TO_SQKM,
        'sqkm_to_sqmi': SQKM_TO_SQMI,
        'acre_to_sqkm': ACRE_TO_SQKM,
        'sqkm_to_acre': SQKM_TO_ACRE,

        # Flow
        'cfs_to_cms': CFS_TO_CMS,
        'cms_to_cfs': CMS_TO_CFS,

        # Volume
        'acft_to_m3': ACFT_TO_M3,
        'm3_to_acft': M3_TO_ACFT,
    }

    # Time interval mappings - reference centralized constants
    # Note: TIME_INTERVALS imported from _constants

    @staticmethod
    @log_call
    def parse_time_interval(interval_str: str) -> int:
        """
        Parse HMS time interval string to minutes.

        Args:
            interval_str: Time interval string (e.g., "15 Minutes", "1 Hour")

        Returns:
            Interval in minutes

        Example:
            >>> HmsUtils.parse_time_interval("15 Minutes")
            15
            >>> HmsUtils.parse_time_interval("1 Hour")
            60
        """
        if interval_str in TIME_INTERVALS:
            return TIME_INTERVALS[interval_str]

        # Try to parse custom format
        match = re.match(r'(\d+)\s*(\w+)', interval_str)
        if match:
            value = int(match.group(1))
            unit = match.group(2).lower()

            if 'min' in unit:
                return value
            elif 'hour' in unit:
                return value * MINUTES_PER_HOUR
            elif 'day' in unit:
                return value * MINUTES_PER_DAY

        raise ValueError(f"Unknown time interval: {interval_str}")

    @staticmethod
    @log_call
    def minutes_to_interval_string(minutes: int) -> str:
        """
        Convert minutes to HMS interval string.

        Args:
            minutes: Number of minutes

        Returns:
            HMS interval string

        Example:
            >>> HmsUtils.minutes_to_interval_string(15)
            '15 Minutes'
            >>> HmsUtils.minutes_to_interval_string(60)
            '1 Hour'
        """
        # Find exact match
        for interval_str, mins in TIME_INTERVALS.items():
            if mins == minutes:
                return interval_str

        # Generate approximate string
        if minutes < MINUTES_PER_HOUR:
            return f"{minutes} Minutes" if minutes > 1 else "1 Minute"
        elif minutes < MINUTES_PER_DAY:
            hours = minutes // MINUTES_PER_HOUR
            return f"{hours} Hours" if hours > 1 else "1 Hour"
        else:
            days = minutes // MINUTES_PER_DAY
            return f"{days} Days" if days > 1 else "1 Day"

    @staticmethod
    @log_call
    def convert_units(
        value: float,
        from_unit: str,
        to_unit: str
    ) -> float:
        """
        Convert between common hydrologic units.

        Args:
            value: Value to convert
            from_unit: Source unit
            to_unit: Target unit

        Returns:
            Converted value

        Example:
            >>> HmsUtils.convert_units(1.0, "in", "mm")
            25.4
            >>> HmsUtils.convert_units(100, "cfs", "cms")
            2.83...
        """
        # Normalize unit names
        from_unit = from_unit.lower().replace(' ', '')
        to_unit = to_unit.lower().replace(' ', '')

        # Same units
        if from_unit == to_unit:
            return value

        # Build conversion key
        key = f"{from_unit}_to_{to_unit}"

        if key in HmsUtils.CONVERSION_FACTORS:
            return value * HmsUtils.CONVERSION_FACTORS[key]

        # Try reverse
        reverse_key = f"{to_unit}_to_{from_unit}"
        if reverse_key in HmsUtils.CONVERSION_FACTORS:
            return value / HmsUtils.CONVERSION_FACTORS[reverse_key]

        raise ValueError(f"Unknown unit conversion: {from_unit} to {to_unit}")

    @staticmethod
    @log_call
    def parse_hms_date(date_str: str, time_str: str = "00:00") -> datetime:
        """
        Parse HMS date/time strings to datetime.

        Handles multiple HMS date formats:
        - "16 January 1973" (full month, spaces)
        - "16 Jan 1973" (abbreviated month, spaces)
        - "16Jan1973" (compact, no spaces)

        Args:
            date_str: Date string in any HMS format
            time_str: Time string (e.g., "00:00")

        Returns:
            datetime object

        Example:
            >>> dt = HmsUtils.parse_hms_date("15 March 2020", "12:30")
            >>> print(dt)
            2020-03-15 12:30:00
        """
        datetime_str = f"{date_str} {time_str}"
        formats = [
            "%d %B %Y %H:%M",   # "16 January 1973 03:00"
            "%d %b %Y %H:%M",   # "16 Jan 1973 03:00"
            "%d%b%Y %H:%M",     # "16Jan1973 03:00"
        ]
        for fmt in formats:
            try:
                return datetime.strptime(datetime_str, fmt)
            except ValueError:
                continue
        raise ValueError(
            f"Error parsing date/time: time data '{datetime_str}' "
            f"does not match any HMS date format"
        )

    @staticmethod
    @log_call
    def format_hms_date(dt: datetime) -> Tuple[str, str]:
        """
        Format datetime to HMS date/time strings.

        Args:
            dt: datetime object

        Returns:
            Tuple of (date_str, time_str)

        Example:
            >>> dt = datetime(2020, 3, 15, 12, 30)
            >>> date_str, time_str = HmsUtils.format_hms_date(dt)
            >>> print(date_str, time_str)
            15Mar2020 12:30
        """
        date_str = dt.strftime("%d%b%Y")
        time_str = dt.strftime("%H:%M")
        return date_str, time_str

    @staticmethod
    @log_call
    def copy_project(
        source_folder: Union[str, Path],
        dest_folder: Union[str, Path],
        overwrite: bool = False
    ) -> Path:
        """
        Copy an HMS project to a new location.

        Args:
            source_folder: Source project folder
            dest_folder: Destination folder
            overwrite: Whether to overwrite existing destination

        Returns:
            Path to the copied project
        """
        source_folder = Path(source_folder)
        dest_folder = Path(dest_folder)

        if not source_folder.exists():
            raise FileNotFoundError(f"Source folder not found: {source_folder}")

        if dest_folder.exists():
            if overwrite:
                shutil.rmtree(dest_folder)
            else:
                raise FileExistsError(f"Destination exists: {dest_folder}")

        logger.info(f"Copying project from {source_folder} to {dest_folder}")
        shutil.copytree(source_folder, dest_folder)

        return dest_folder

    @staticmethod
    @log_call
    def list_project_files(
        project_folder: Union[str, Path]
    ) -> Dict[str, List[Path]]:
        """
        List all HMS files in a project folder by type.

        Args:
            project_folder: Path to the project folder

        Returns:
            Dictionary mapping file type to list of file paths

        Example:
            >>> files = HmsUtils.list_project_files("MyProject")
            >>> print(files['basin'])  # List of .basin files
        """
        project_folder = Path(project_folder)

        if not project_folder.exists():
            raise FileNotFoundError(f"Project folder not found: {project_folder}")

        file_types = {
            'hms': list(project_folder.glob("*.hms")),
            'basin': list(project_folder.glob("*.basin")),
            'met': list(project_folder.glob("*.met")),
            'control': list(project_folder.glob("*.control")),
            'gage': list(project_folder.glob("*.gage")),
            'run': list(project_folder.glob("*.run")),
            'dss': list(project_folder.glob("*.dss")),
            'log': list(project_folder.glob("*.log")),
            'geo': list(project_folder.glob("*.geo")),
            'map': list(project_folder.glob("*.map")),
            'grid': list(project_folder.glob("*.grid")),
            'sqlite': list(project_folder.glob("*.sqlite")),
        }

        return file_types

    @staticmethod
    @log_call
    def clean_project_outputs(
        project_folder: Union[str, Path],
        delete_dss: bool = True,
        delete_logs: bool = True
    ) -> int:
        """
        Clean output files from a project folder.

        Args:
            project_folder: Path to the project folder
            delete_dss: Whether to delete DSS files
            delete_logs: Whether to delete log files

        Returns:
            Number of files deleted
        """
        project_folder = Path(project_folder)
        deleted_count = 0

        patterns = []
        if delete_dss:
            patterns.extend(['*.dss'])
        if delete_logs:
            patterns.extend(['*.log', '*.out', '*.err'])

        for pattern in patterns:
            for file in project_folder.glob(pattern):
                try:
                    file.unlink()
                    deleted_count += 1
                    logger.debug(f"Deleted: {file}")
                except Exception as e:
                    logger.warning(f"Could not delete {file}: {e}")

        logger.info(f"Deleted {deleted_count} output files")
        return deleted_count

    @staticmethod
    @log_call
    def validate_project(
        project_folder: Union[str, Path]
    ) -> Dict[str, Any]:
        """
        Validate an HMS project folder.

        Checks for required files and common issues.

        Args:
            project_folder: Path to the project folder

        Returns:
            Dictionary with validation results

        Example:
            >>> result = HmsUtils.validate_project("MyProject")
            >>> if result['valid']:
            ...     print("Project is valid")
        """
        project_folder = Path(project_folder)
        result = {
            'valid': True,
            'errors': [],
            'warnings': [],
            'project_file': None,
            'basin_files': [],
            'met_files': [],
            'control_files': [],
            'run_files': []
        }

        if not project_folder.exists():
            result['valid'] = False
            result['errors'].append(f"Project folder not found: {project_folder}")
            return result

        # Check for .hms project file
        hms_files = list(project_folder.glob("*.hms"))
        if not hms_files:
            result['valid'] = False
            result['errors'].append("No .hms project file found")
        else:
            result['project_file'] = str(hms_files[0])
            if len(hms_files) > 1:
                result['warnings'].append(f"Multiple .hms files found: {len(hms_files)}")

        # Check for basin files
        result['basin_files'] = [str(f) for f in project_folder.glob("*.basin")]
        if not result['basin_files']:
            result['warnings'].append("No .basin files found")

        # Check for met files
        result['met_files'] = [str(f) for f in project_folder.glob("*.met")]

        # Check for control files
        result['control_files'] = [str(f) for f in project_folder.glob("*.control")]

        # Check for run files
        result['run_files'] = [str(f) for f in project_folder.glob("*.run")]

        return result

    @staticmethod
    @log_call
    def get_project_summary(
        project_folder: Union[str, Path]
    ) -> str:
        """
        Get a text summary of a project.

        Args:
            project_folder: Path to the project folder

        Returns:
            Formatted summary string
        """
        files = HmsUtils.list_project_files(project_folder)

        summary = []
        summary.append(f"HMS Project: {Path(project_folder).name}")
        summary.append("=" * 50)

        for file_type, file_list in files.items():
            if file_list:
                summary.append(f"\n{file_type.upper()} files ({len(file_list)}):")
                for f in file_list[:5]:  # Limit to first 5
                    summary.append(f"  - {f.name}")
                if len(file_list) > 5:
                    summary.append(f"  ... and {len(file_list) - 5} more")

        return "\n".join(summary)

    @staticmethod
    def calculate_cn_from_ia(
        initial_abstraction: float,
        method: str = "standard"
    ) -> float:
        """
        Calculate SCS Curve Number from Initial Abstraction.

        Args:
            initial_abstraction: Initial abstraction in inches
            method: Calculation method ("standard" uses Ia = 0.2*S)

        Returns:
            Curve Number (0-100)

        Example:
            >>> cn = HmsUtils.calculate_cn_from_ia(1.0)
            >>> print(f"CN = {cn:.1f}")
        """
        # Standard method: Ia = IA_RATIO * S, where S = (1000/CN) - 10
        # Solving for CN: CN = 1000 / (S + 10) where S = Ia / IA_RATIO
        s = initial_abstraction / IA_RATIO
        cn = CN_FORMULA_NUMERATOR / (s + CN_FORMULA_BASE)
        return max(CN_MIN, min(CN_MAX, cn))

    @staticmethod
    def calculate_ia_from_cn(
        curve_number: float,
        method: str = "standard"
    ) -> float:
        """
        Calculate Initial Abstraction from SCS Curve Number.

        Args:
            curve_number: SCS Curve Number (0-100)
            method: Calculation method ("standard" uses Ia = 0.2*S)

        Returns:
            Initial Abstraction in inches

        Example:
            >>> ia = HmsUtils.calculate_ia_from_cn(75)
            >>> print(f"Ia = {ia:.2f} inches")
        """
        if curve_number <= CN_MIN or curve_number > CN_MAX:
            raise ValueError(f"Curve Number must be between {CN_MIN} and {CN_MAX}, got {curve_number}")

        # S = (CN_FORMULA_NUMERATOR/CN) - CN_FORMULA_BASE
        s = (CN_FORMULA_NUMERATOR / curve_number) - CN_FORMULA_BASE
        # Ia = IA_RATIO * S
        ia = IA_RATIO * s
        return ia

    @staticmethod
    @log_call
    def clone_file(
        template_path: Union[str, Path],
        new_path: Union[str, Path],
        modify_func: Optional[Any] = None
    ) -> Path:
        """
        Clone a file with optional modification via callback function.

        This is the core clone utility used by all HMS clone operations
        (clone_run, clone_basin, clone_met). It follows the ras-commander
        pattern of reading, modifying via callback, and writing.

        Args:
            template_path: Path to template file
            new_path: Path for new file
            modify_func: Optional callback function(lines: List[str]) -> List[str]
                        that modifies file content

        Returns:
            Path to the cloned file

        Raises:
            FileNotFoundError: If template doesn't exist
            FileExistsError: If new_path already exists

        Example:
            >>> def update_name(lines):
            ...     return [line.replace("Old", "New") for line in lines]
            >>> HmsUtils.clone_file("old.basin", "new.basin", update_name)

        Note:
            This implements the CLB Engineering LLM Forward Approach:
            - Non-destructive (creates new file, preserves original)
            - Traceable (clear source → destination relationship)
            - Modifiable (callback allows precise updates)
        """
        template_path = Path(template_path)
        new_path = Path(new_path)

        if not template_path.exists():
            raise FileNotFoundError(f"Template file not found: {template_path}")

        if new_path.exists():
            raise FileExistsError(f"Destination file already exists: {new_path}")

        content = HmsFileParser.read_file(template_path)

        # Apply modification function if provided
        if modify_func is not None:
            lines = content.splitlines(keepends=True)
            modified_lines = modify_func(lines)
            content = ''.join(modified_lines)

        # Write new file
        HmsFileParser.write_file(new_path, content)
        logger.info(f"Cloned file: {template_path.name}{new_path.name}")

        return new_path

    @staticmethod
    def _legacy_update_project_file(
        hms_file: Path,
        entry_type: str,
        entry_name: str,
    ) -> bool:
        """Fallback for historical flat project-file entries."""

        content = HmsFileParser.read_file(hms_file)
        file_extensions = {
            'Run': 'run',
            'Gage': 'gage',
        }
        extension = file_extensions.get(entry_type, entry_type.lower())

        entry_pattern = rf'^{entry_type}\s+File:\s*{re.escape(entry_name)}\.{extension}\s*$'
        if re.search(entry_pattern, content, re.MULTILINE | re.IGNORECASE):
            logger.info(f"{entry_type} '{entry_name}' already in project file")
            return True

        new_entry = f"{entry_type} File: {entry_name}.{extension}\n"
        type_pattern = rf'^{entry_type}\s+File:.*$'
        matches = list(re.finditer(type_pattern, content, re.MULTILINE | re.IGNORECASE))

        if matches:
            insert_pos = matches[-1].end()
            content = content[:insert_pos] + '\n' + new_entry + content[insert_pos:]
        else:
            end_match = re.search(r'^End:\s*$', content, re.MULTILINE)
            if end_match:
                content = content[:end_match.start()] + new_entry + content[end_match.start():]
            else:
                content = content.rstrip() + '\n' + new_entry

        HmsFileParser.write_file(hms_file, content)
        logger.info(f"Added legacy {entry_type} file entry '{entry_name}' to project file")
        return True

    @staticmethod
    @log_call
    def update_project_file(
        hms_file: Union[str, Path],
        entry_type: str,
        entry_name: str
    ) -> bool:
        """
        Update HMS project file to register a new component.

        When cloning basins, mets, controls, or runs, the .hms project
        file must be updated to register the new component. This follows
        the ras-commander pattern of updating the project file after
        creating new model components.

        Args:
            hms_file: Path to the .hms project file
            entry_type: Type of entry ('Basin', 'Met', 'Control', 'Run', etc.)
            entry_name: Name of the new component (filename without extension)

        Returns:
            True if successful

        Raises:
            FileNotFoundError: If hms_file doesn't exist

        Example:
            >>> HmsUtils.update_project_file(
            ...     "MyProject.hms",
            ...     "Basin",
            ...     "Updated_Atlas14"
            ... )

        Note:
            Basin, Met/Meteorology, and Control entries are written as real HMS
            registry blocks for new registrations. Existing legacy flat entries
            such as "Basin File:" and "Met File:" are treated as already
            registered to avoid duplicates. Other entry types keep the
            historical flat-line fallback for compatibility.
        """
        hms_file = Path(hms_file)

        if not hms_file.exists():
            raise FileNotFoundError(f"HMS project file not found: {hms_file}")

        try:
            register_project_block(
                hms_file,
                entry_type=entry_type,
                logical_name=entry_name,
                allow_existing=True,
            )
        except ValueError:
            logger.warning(
                f"Unknown entry type '{entry_type}', using legacy flat-line registration"
            )
            return HmsUtils._legacy_update_project_file(hms_file, entry_type, entry_name)

        logger.info(f"Added {entry_type} '{entry_name}' to project file")

        return True

parse_time_interval(interval_str) staticmethod

Parse HMS time interval string to minutes.

Parameters:

Name Type Description Default
interval_str str

Time interval string (e.g., "15 Minutes", "1 Hour")

required

Returns:

Type Description
int

Interval in minutes

Example

HmsUtils.parse_time_interval("15 Minutes") 15 HmsUtils.parse_time_interval("1 Hour") 60

Source code in hms_commander/HmsUtils.py
@staticmethod
@log_call
def parse_time_interval(interval_str: str) -> int:
    """
    Parse HMS time interval string to minutes.

    Args:
        interval_str: Time interval string (e.g., "15 Minutes", "1 Hour")

    Returns:
        Interval in minutes

    Example:
        >>> HmsUtils.parse_time_interval("15 Minutes")
        15
        >>> HmsUtils.parse_time_interval("1 Hour")
        60
    """
    if interval_str in TIME_INTERVALS:
        return TIME_INTERVALS[interval_str]

    # Try to parse custom format
    match = re.match(r'(\d+)\s*(\w+)', interval_str)
    if match:
        value = int(match.group(1))
        unit = match.group(2).lower()

        if 'min' in unit:
            return value
        elif 'hour' in unit:
            return value * MINUTES_PER_HOUR
        elif 'day' in unit:
            return value * MINUTES_PER_DAY

    raise ValueError(f"Unknown time interval: {interval_str}")

minutes_to_interval_string(minutes) staticmethod

Convert minutes to HMS interval string.

Parameters:

Name Type Description Default
minutes int

Number of minutes

required

Returns:

Type Description
str

HMS interval string

Example

HmsUtils.minutes_to_interval_string(15) '15 Minutes' HmsUtils.minutes_to_interval_string(60) '1 Hour'

Source code in hms_commander/HmsUtils.py
@staticmethod
@log_call
def minutes_to_interval_string(minutes: int) -> str:
    """
    Convert minutes to HMS interval string.

    Args:
        minutes: Number of minutes

    Returns:
        HMS interval string

    Example:
        >>> HmsUtils.minutes_to_interval_string(15)
        '15 Minutes'
        >>> HmsUtils.minutes_to_interval_string(60)
        '1 Hour'
    """
    # Find exact match
    for interval_str, mins in TIME_INTERVALS.items():
        if mins == minutes:
            return interval_str

    # Generate approximate string
    if minutes < MINUTES_PER_HOUR:
        return f"{minutes} Minutes" if minutes > 1 else "1 Minute"
    elif minutes < MINUTES_PER_DAY:
        hours = minutes // MINUTES_PER_HOUR
        return f"{hours} Hours" if hours > 1 else "1 Hour"
    else:
        days = minutes // MINUTES_PER_DAY
        return f"{days} Days" if days > 1 else "1 Day"

convert_units(value, from_unit, to_unit) staticmethod

Convert between common hydrologic units.

Parameters:

Name Type Description Default
value float

Value to convert

required
from_unit str

Source unit

required
to_unit str

Target unit

required

Returns:

Type Description
float

Converted value

Example

HmsUtils.convert_units(1.0, "in", "mm") 25.4 HmsUtils.convert_units(100, "cfs", "cms") 2.83...

Source code in hms_commander/HmsUtils.py
@staticmethod
@log_call
def convert_units(
    value: float,
    from_unit: str,
    to_unit: str
) -> float:
    """
    Convert between common hydrologic units.

    Args:
        value: Value to convert
        from_unit: Source unit
        to_unit: Target unit

    Returns:
        Converted value

    Example:
        >>> HmsUtils.convert_units(1.0, "in", "mm")
        25.4
        >>> HmsUtils.convert_units(100, "cfs", "cms")
        2.83...
    """
    # Normalize unit names
    from_unit = from_unit.lower().replace(' ', '')
    to_unit = to_unit.lower().replace(' ', '')

    # Same units
    if from_unit == to_unit:
        return value

    # Build conversion key
    key = f"{from_unit}_to_{to_unit}"

    if key in HmsUtils.CONVERSION_FACTORS:
        return value * HmsUtils.CONVERSION_FACTORS[key]

    # Try reverse
    reverse_key = f"{to_unit}_to_{from_unit}"
    if reverse_key in HmsUtils.CONVERSION_FACTORS:
        return value / HmsUtils.CONVERSION_FACTORS[reverse_key]

    raise ValueError(f"Unknown unit conversion: {from_unit} to {to_unit}")

parse_hms_date(date_str, time_str='00:00') staticmethod

Parse HMS date/time strings to datetime.

Handles multiple HMS date formats: - "16 January 1973" (full month, spaces) - "16 Jan 1973" (abbreviated month, spaces) - "16Jan1973" (compact, no spaces)

Parameters:

Name Type Description Default
date_str str

Date string in any HMS format

required
time_str str

Time string (e.g., "00:00")

'00:00'

Returns:

Type Description
datetime

datetime object

Example

dt = HmsUtils.parse_hms_date("15 March 2020", "12:30") print(dt) 2020-03-15 12:30:00

Source code in hms_commander/HmsUtils.py
@staticmethod
@log_call
def parse_hms_date(date_str: str, time_str: str = "00:00") -> datetime:
    """
    Parse HMS date/time strings to datetime.

    Handles multiple HMS date formats:
    - "16 January 1973" (full month, spaces)
    - "16 Jan 1973" (abbreviated month, spaces)
    - "16Jan1973" (compact, no spaces)

    Args:
        date_str: Date string in any HMS format
        time_str: Time string (e.g., "00:00")

    Returns:
        datetime object

    Example:
        >>> dt = HmsUtils.parse_hms_date("15 March 2020", "12:30")
        >>> print(dt)
        2020-03-15 12:30:00
    """
    datetime_str = f"{date_str} {time_str}"
    formats = [
        "%d %B %Y %H:%M",   # "16 January 1973 03:00"
        "%d %b %Y %H:%M",   # "16 Jan 1973 03:00"
        "%d%b%Y %H:%M",     # "16Jan1973 03:00"
    ]
    for fmt in formats:
        try:
            return datetime.strptime(datetime_str, fmt)
        except ValueError:
            continue
    raise ValueError(
        f"Error parsing date/time: time data '{datetime_str}' "
        f"does not match any HMS date format"
    )

format_hms_date(dt) staticmethod

Format datetime to HMS date/time strings.

Parameters:

Name Type Description Default
dt datetime

datetime object

required

Returns:

Type Description
Tuple[str, str]

Tuple of (date_str, time_str)

Example

dt = datetime(2020, 3, 15, 12, 30) date_str, time_str = HmsUtils.format_hms_date(dt) print(date_str, time_str) 15Mar2020 12:30

Source code in hms_commander/HmsUtils.py
@staticmethod
@log_call
def format_hms_date(dt: datetime) -> Tuple[str, str]:
    """
    Format datetime to HMS date/time strings.

    Args:
        dt: datetime object

    Returns:
        Tuple of (date_str, time_str)

    Example:
        >>> dt = datetime(2020, 3, 15, 12, 30)
        >>> date_str, time_str = HmsUtils.format_hms_date(dt)
        >>> print(date_str, time_str)
        15Mar2020 12:30
    """
    date_str = dt.strftime("%d%b%Y")
    time_str = dt.strftime("%H:%M")
    return date_str, time_str

copy_project(source_folder, dest_folder, overwrite=False) staticmethod

Copy an HMS project to a new location.

Parameters:

Name Type Description Default
source_folder Union[str, Path]

Source project folder

required
dest_folder Union[str, Path]

Destination folder

required
overwrite bool

Whether to overwrite existing destination

False

Returns:

Type Description
Path

Path to the copied project

Source code in hms_commander/HmsUtils.py
@staticmethod
@log_call
def copy_project(
    source_folder: Union[str, Path],
    dest_folder: Union[str, Path],
    overwrite: bool = False
) -> Path:
    """
    Copy an HMS project to a new location.

    Args:
        source_folder: Source project folder
        dest_folder: Destination folder
        overwrite: Whether to overwrite existing destination

    Returns:
        Path to the copied project
    """
    source_folder = Path(source_folder)
    dest_folder = Path(dest_folder)

    if not source_folder.exists():
        raise FileNotFoundError(f"Source folder not found: {source_folder}")

    if dest_folder.exists():
        if overwrite:
            shutil.rmtree(dest_folder)
        else:
            raise FileExistsError(f"Destination exists: {dest_folder}")

    logger.info(f"Copying project from {source_folder} to {dest_folder}")
    shutil.copytree(source_folder, dest_folder)

    return dest_folder

list_project_files(project_folder) staticmethod

List all HMS files in a project folder by type.

Parameters:

Name Type Description Default
project_folder Union[str, Path]

Path to the project folder

required

Returns:

Type Description
Dict[str, List[Path]]

Dictionary mapping file type to list of file paths

Example

files = HmsUtils.list_project_files("MyProject") print(files['basin']) # List of .basin files

Source code in hms_commander/HmsUtils.py
@staticmethod
@log_call
def list_project_files(
    project_folder: Union[str, Path]
) -> Dict[str, List[Path]]:
    """
    List all HMS files in a project folder by type.

    Args:
        project_folder: Path to the project folder

    Returns:
        Dictionary mapping file type to list of file paths

    Example:
        >>> files = HmsUtils.list_project_files("MyProject")
        >>> print(files['basin'])  # List of .basin files
    """
    project_folder = Path(project_folder)

    if not project_folder.exists():
        raise FileNotFoundError(f"Project folder not found: {project_folder}")

    file_types = {
        'hms': list(project_folder.glob("*.hms")),
        'basin': list(project_folder.glob("*.basin")),
        'met': list(project_folder.glob("*.met")),
        'control': list(project_folder.glob("*.control")),
        'gage': list(project_folder.glob("*.gage")),
        'run': list(project_folder.glob("*.run")),
        'dss': list(project_folder.glob("*.dss")),
        'log': list(project_folder.glob("*.log")),
        'geo': list(project_folder.glob("*.geo")),
        'map': list(project_folder.glob("*.map")),
        'grid': list(project_folder.glob("*.grid")),
        'sqlite': list(project_folder.glob("*.sqlite")),
    }

    return file_types

clean_project_outputs(project_folder, delete_dss=True, delete_logs=True) staticmethod

Clean output files from a project folder.

Parameters:

Name Type Description Default
project_folder Union[str, Path]

Path to the project folder

required
delete_dss bool

Whether to delete DSS files

True
delete_logs bool

Whether to delete log files

True

Returns:

Type Description
int

Number of files deleted

Source code in hms_commander/HmsUtils.py
@staticmethod
@log_call
def clean_project_outputs(
    project_folder: Union[str, Path],
    delete_dss: bool = True,
    delete_logs: bool = True
) -> int:
    """
    Clean output files from a project folder.

    Args:
        project_folder: Path to the project folder
        delete_dss: Whether to delete DSS files
        delete_logs: Whether to delete log files

    Returns:
        Number of files deleted
    """
    project_folder = Path(project_folder)
    deleted_count = 0

    patterns = []
    if delete_dss:
        patterns.extend(['*.dss'])
    if delete_logs:
        patterns.extend(['*.log', '*.out', '*.err'])

    for pattern in patterns:
        for file in project_folder.glob(pattern):
            try:
                file.unlink()
                deleted_count += 1
                logger.debug(f"Deleted: {file}")
            except Exception as e:
                logger.warning(f"Could not delete {file}: {e}")

    logger.info(f"Deleted {deleted_count} output files")
    return deleted_count

validate_project(project_folder) staticmethod

Validate an HMS project folder.

Checks for required files and common issues.

Parameters:

Name Type Description Default
project_folder Union[str, Path]

Path to the project folder

required

Returns:

Type Description
Dict[str, Any]

Dictionary with validation results

Example

result = HmsUtils.validate_project("MyProject") if result['valid']: ... print("Project is valid")

Source code in hms_commander/HmsUtils.py
@staticmethod
@log_call
def validate_project(
    project_folder: Union[str, Path]
) -> Dict[str, Any]:
    """
    Validate an HMS project folder.

    Checks for required files and common issues.

    Args:
        project_folder: Path to the project folder

    Returns:
        Dictionary with validation results

    Example:
        >>> result = HmsUtils.validate_project("MyProject")
        >>> if result['valid']:
        ...     print("Project is valid")
    """
    project_folder = Path(project_folder)
    result = {
        'valid': True,
        'errors': [],
        'warnings': [],
        'project_file': None,
        'basin_files': [],
        'met_files': [],
        'control_files': [],
        'run_files': []
    }

    if not project_folder.exists():
        result['valid'] = False
        result['errors'].append(f"Project folder not found: {project_folder}")
        return result

    # Check for .hms project file
    hms_files = list(project_folder.glob("*.hms"))
    if not hms_files:
        result['valid'] = False
        result['errors'].append("No .hms project file found")
    else:
        result['project_file'] = str(hms_files[0])
        if len(hms_files) > 1:
            result['warnings'].append(f"Multiple .hms files found: {len(hms_files)}")

    # Check for basin files
    result['basin_files'] = [str(f) for f in project_folder.glob("*.basin")]
    if not result['basin_files']:
        result['warnings'].append("No .basin files found")

    # Check for met files
    result['met_files'] = [str(f) for f in project_folder.glob("*.met")]

    # Check for control files
    result['control_files'] = [str(f) for f in project_folder.glob("*.control")]

    # Check for run files
    result['run_files'] = [str(f) for f in project_folder.glob("*.run")]

    return result

get_project_summary(project_folder) staticmethod

Get a text summary of a project.

Parameters:

Name Type Description Default
project_folder Union[str, Path]

Path to the project folder

required

Returns:

Type Description
str

Formatted summary string

Source code in hms_commander/HmsUtils.py
@staticmethod
@log_call
def get_project_summary(
    project_folder: Union[str, Path]
) -> str:
    """
    Get a text summary of a project.

    Args:
        project_folder: Path to the project folder

    Returns:
        Formatted summary string
    """
    files = HmsUtils.list_project_files(project_folder)

    summary = []
    summary.append(f"HMS Project: {Path(project_folder).name}")
    summary.append("=" * 50)

    for file_type, file_list in files.items():
        if file_list:
            summary.append(f"\n{file_type.upper()} files ({len(file_list)}):")
            for f in file_list[:5]:  # Limit to first 5
                summary.append(f"  - {f.name}")
            if len(file_list) > 5:
                summary.append(f"  ... and {len(file_list) - 5} more")

    return "\n".join(summary)

calculate_cn_from_ia(initial_abstraction, method='standard') staticmethod

Calculate SCS Curve Number from Initial Abstraction.

Parameters:

Name Type Description Default
initial_abstraction float

Initial abstraction in inches

required
method str

Calculation method ("standard" uses Ia = 0.2*S)

'standard'

Returns:

Type Description
float

Curve Number (0-100)

Example

cn = HmsUtils.calculate_cn_from_ia(1.0) print(f"CN = {cn:.1f}")

Source code in hms_commander/HmsUtils.py
@staticmethod
def calculate_cn_from_ia(
    initial_abstraction: float,
    method: str = "standard"
) -> float:
    """
    Calculate SCS Curve Number from Initial Abstraction.

    Args:
        initial_abstraction: Initial abstraction in inches
        method: Calculation method ("standard" uses Ia = 0.2*S)

    Returns:
        Curve Number (0-100)

    Example:
        >>> cn = HmsUtils.calculate_cn_from_ia(1.0)
        >>> print(f"CN = {cn:.1f}")
    """
    # Standard method: Ia = IA_RATIO * S, where S = (1000/CN) - 10
    # Solving for CN: CN = 1000 / (S + 10) where S = Ia / IA_RATIO
    s = initial_abstraction / IA_RATIO
    cn = CN_FORMULA_NUMERATOR / (s + CN_FORMULA_BASE)
    return max(CN_MIN, min(CN_MAX, cn))

calculate_ia_from_cn(curve_number, method='standard') staticmethod

Calculate Initial Abstraction from SCS Curve Number.

Parameters:

Name Type Description Default
curve_number float

SCS Curve Number (0-100)

required
method str

Calculation method ("standard" uses Ia = 0.2*S)

'standard'

Returns:

Type Description
float

Initial Abstraction in inches

Example

ia = HmsUtils.calculate_ia_from_cn(75) print(f"Ia = {ia:.2f} inches")

Source code in hms_commander/HmsUtils.py
@staticmethod
def calculate_ia_from_cn(
    curve_number: float,
    method: str = "standard"
) -> float:
    """
    Calculate Initial Abstraction from SCS Curve Number.

    Args:
        curve_number: SCS Curve Number (0-100)
        method: Calculation method ("standard" uses Ia = 0.2*S)

    Returns:
        Initial Abstraction in inches

    Example:
        >>> ia = HmsUtils.calculate_ia_from_cn(75)
        >>> print(f"Ia = {ia:.2f} inches")
    """
    if curve_number <= CN_MIN or curve_number > CN_MAX:
        raise ValueError(f"Curve Number must be between {CN_MIN} and {CN_MAX}, got {curve_number}")

    # S = (CN_FORMULA_NUMERATOR/CN) - CN_FORMULA_BASE
    s = (CN_FORMULA_NUMERATOR / curve_number) - CN_FORMULA_BASE
    # Ia = IA_RATIO * S
    ia = IA_RATIO * s
    return ia

clone_file(template_path, new_path, modify_func=None) staticmethod

Clone a file with optional modification via callback function.

This is the core clone utility used by all HMS clone operations (clone_run, clone_basin, clone_met). It follows the ras-commander pattern of reading, modifying via callback, and writing.

Parameters:

Name Type Description Default
template_path Union[str, Path]

Path to template file

required
new_path Union[str, Path]

Path for new file

required
modify_func Optional[Any]

Optional callback function(lines: List[str]) -> List[str] that modifies file content

None

Returns:

Type Description
Path

Path to the cloned file

Raises:

Type Description
FileNotFoundError

If template doesn't exist

FileExistsError

If new_path already exists

Example

def update_name(lines): ... return [line.replace("Old", "New") for line in lines] HmsUtils.clone_file("old.basin", "new.basin", update_name)

Note

This implements the CLB Engineering LLM Forward Approach: - Non-destructive (creates new file, preserves original) - Traceable (clear source → destination relationship) - Modifiable (callback allows precise updates)

Source code in hms_commander/HmsUtils.py
@staticmethod
@log_call
def clone_file(
    template_path: Union[str, Path],
    new_path: Union[str, Path],
    modify_func: Optional[Any] = None
) -> Path:
    """
    Clone a file with optional modification via callback function.

    This is the core clone utility used by all HMS clone operations
    (clone_run, clone_basin, clone_met). It follows the ras-commander
    pattern of reading, modifying via callback, and writing.

    Args:
        template_path: Path to template file
        new_path: Path for new file
        modify_func: Optional callback function(lines: List[str]) -> List[str]
                    that modifies file content

    Returns:
        Path to the cloned file

    Raises:
        FileNotFoundError: If template doesn't exist
        FileExistsError: If new_path already exists

    Example:
        >>> def update_name(lines):
        ...     return [line.replace("Old", "New") for line in lines]
        >>> HmsUtils.clone_file("old.basin", "new.basin", update_name)

    Note:
        This implements the CLB Engineering LLM Forward Approach:
        - Non-destructive (creates new file, preserves original)
        - Traceable (clear source → destination relationship)
        - Modifiable (callback allows precise updates)
    """
    template_path = Path(template_path)
    new_path = Path(new_path)

    if not template_path.exists():
        raise FileNotFoundError(f"Template file not found: {template_path}")

    if new_path.exists():
        raise FileExistsError(f"Destination file already exists: {new_path}")

    content = HmsFileParser.read_file(template_path)

    # Apply modification function if provided
    if modify_func is not None:
        lines = content.splitlines(keepends=True)
        modified_lines = modify_func(lines)
        content = ''.join(modified_lines)

    # Write new file
    HmsFileParser.write_file(new_path, content)
    logger.info(f"Cloned file: {template_path.name}{new_path.name}")

    return new_path

update_project_file(hms_file, entry_type, entry_name) staticmethod

Update HMS project file to register a new component.

When cloning basins, mets, controls, or runs, the .hms project file must be updated to register the new component. This follows the ras-commander pattern of updating the project file after creating new model components.

Parameters:

Name Type Description Default
hms_file Union[str, Path]

Path to the .hms project file

required
entry_type str

Type of entry ('Basin', 'Met', 'Control', 'Run', etc.)

required
entry_name str

Name of the new component (filename without extension)

required

Returns:

Type Description
bool

True if successful

Raises:

Type Description
FileNotFoundError

If hms_file doesn't exist

Example

HmsUtils.update_project_file( ... "MyProject.hms", ... "Basin", ... "Updated_Atlas14" ... )

Note

Basin, Met/Meteorology, and Control entries are written as real HMS registry blocks for new registrations. Existing legacy flat entries such as "Basin File:" and "Met File:" are treated as already registered to avoid duplicates. Other entry types keep the historical flat-line fallback for compatibility.

Source code in hms_commander/HmsUtils.py
@staticmethod
@log_call
def update_project_file(
    hms_file: Union[str, Path],
    entry_type: str,
    entry_name: str
) -> bool:
    """
    Update HMS project file to register a new component.

    When cloning basins, mets, controls, or runs, the .hms project
    file must be updated to register the new component. This follows
    the ras-commander pattern of updating the project file after
    creating new model components.

    Args:
        hms_file: Path to the .hms project file
        entry_type: Type of entry ('Basin', 'Met', 'Control', 'Run', etc.)
        entry_name: Name of the new component (filename without extension)

    Returns:
        True if successful

    Raises:
        FileNotFoundError: If hms_file doesn't exist

    Example:
        >>> HmsUtils.update_project_file(
        ...     "MyProject.hms",
        ...     "Basin",
        ...     "Updated_Atlas14"
        ... )

    Note:
        Basin, Met/Meteorology, and Control entries are written as real HMS
        registry blocks for new registrations. Existing legacy flat entries
        such as "Basin File:" and "Met File:" are treated as already
        registered to avoid duplicates. Other entry types keep the
        historical flat-line fallback for compatibility.
    """
    hms_file = Path(hms_file)

    if not hms_file.exists():
        raise FileNotFoundError(f"HMS project file not found: {hms_file}")

    try:
        register_project_block(
            hms_file,
            entry_type=entry_type,
            logical_name=entry_name,
            allow_existing=True,
        )
    except ValueError:
        logger.warning(
            f"Unknown entry type '{entry_type}', using legacy flat-line registration"
        )
        return HmsUtils._legacy_update_project_file(hms_file, entry_type, entry_name)

    logger.info(f"Added {entry_type} '{entry_name}' to project file")

    return True
CLB Engineering Corporation  ·  LLM Forward Engineering
HMS Commander is a free and open-source project maintained by CLB Engineering Corporation. For agencies and firms seeking to modernize H&H workflows with LLM Forward approaches, contact CLB to partner with the engineers who wrote the automation.