Skip to content

Age tracers hard coded rate of aging based on a 365 days year #521

@NelsonPoumaere

Description

@NelsonPoumaere

Hi,

The rate of aging of age tracers seems to be based on a 365 days year, independent of the calendar type actually used.

In MOM6-examples/src/MOM6/src/tracer/ideal_age_example.F90, line 371:

Isecs_per_year = 1.0 / (365.0*86400.0*US%s_to_T)

and later in the same file, when the aging is implemented, line 454:

 if (G%mask2dT(i,j) > 0.0) then
               layer_frac = BL_layers(i,j)-nk
               CS%tr(i,j,k,m) = (1.-layer_frac) * (CS%tr(i,j,k,m) + G%mask2dT(i,j)*dt &
                             *Isecs_per_year) + layer_frac * young_val
           else
               CS%tr(i,j,k,m) = CS%land_val(m)
endif

I think the same issue exists in ice_age_tracer.F90 and oil_tracer.F90

I'm using a 360 days calendar and for now I can simply replace 365.0 by 360.0.
I'd be happy to make the necessary changes to make the rate of aging compliant with the calendar, but don't have time to do so at the moment (probably in a month).

All the best,
Nelson

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type
    No fields configured for issues without a type.

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions