Opened 11 days ago

#2507 new Bug

PAR output averaged on free-of-ice cell area

Reported by: cperruche Owned by: systeam
Priority: low Milestone:
Component: TOP Version:
Severity: minor Keywords:
Cc:

Description

Context

In p4zopt.F90, PAR is computed on the part of the cell that is free of ice. Below sea ice, it is assumed that there is no primary production

Analysis

But when we save it in output "PAR", I think it should be multiplied by the fraction of free-of-ice water to get the PAR mean over the whole cell

Fix

CALL iom_put( "PARDM", zpar(:,:,:) * tmask(:,:,:) * (1-fr_i(:,:)) )  ! Photosynthetically Available Radiation
CALL iom_put( "PAR"  , emoy(:,:,:) * tmask(:,:,:) * (1-fr_i(:,:)) )  ! Photosynthetically Available Radiation

instead of

CALL iom_put( "PARDM", zpar(:,:,:) * tmask(:,:,:) )  ! Photosynthetically Available Radiation
CALL iom_put( "PAR"  , emoy(:,:,:) * tmask(:,:,:) )  ! Photosynthetically Available Radiation

Commit History (0)

(No commits)

Change History (0)

Note: See TracTickets for help on using tickets.