NOC-MSM / SE-NEMO

Shelf Enabled Global NEMO
GNU General Public License v3.0
2 stars 1 forks source link

Vertical interpolation code in 4.2.1 migration fails #138

Closed jdha closed 8 months ago

jdha commented 9 months ago

Problem:

Model run hangs when trying to open the depth/mask file used for interpolating the initial conditions onto the native grid.

Analysis:

dtatsd.F90 was modified in 4.0.4 to read depth/mask information in order to linearly interpolate T/S initial conditions onto the native grid. This involved making fld_clopn a PUBLIC routine in fldread.F90 and calling it within dtatsd.F90 to perform the read. It appears that fld_clopn has been refactored into two subroutines in 4.2.1 fld_def and fld_clopn.

Solution:

Prepend the original call to fld_clopn with a call to fld_def:

            CALL fld_def ( sf_tsd(jp_dep) )
            CALL fld_clopn ( sf_tsd(jp_dep) )