Changeset 14208

Timestamp:

2020-12-18T10:50:15+01:00 (3 years ago)

Author:

cetlod

Message:

r4.0-HEAD : bugfix to diagnose PAR in PISCES as a mean over the whole cell, see ticket #2507

Location:

NEMO/releases/r4.0/r4.0-HEAD

Files:

• cfgs/SHARED/field_def_nemo-pisces.xml (modified) (1 diff)
• src/TOP/PISCES/P4Z/p4zopt.F90 (modified) (1 diff)

NEMO/releases/r4.0/r4.0-HEAD/cfgs/SHARED/field_def_nemo-pisces.xml

@@ -166 +166 @@
        <field id="CO3sat"      long_name="CO3 saturation"                          unit="mol/m3"     grid_ref="grid_T_3D" />
        <field id="PAR"         long_name="Photosynthetically Available Radiation"  unit="W/m2"       grid_ref="grid_T_3D" />
-       <field id="PARDM"       long_name="Daily mean PAR"                          unit="W/m2"       grid_ref="grid_T_3D" />
        <field id="PPPHYN"      long_name="Primary production of nanophyto"         unit="molC/m3/s"  grid_ref="grid_T_3D" />
        <field id="PPPHYP"      long_name="Primary production of picophyto"         unit="molC/m3/s"  grid_ref="grid_T_3D" />

NEMO/releases/r4.0/r4.0-HEAD/src/TOP/PISCES/P4Z/p4zopt.F90

@@ -264 +264 @@
       IF( lk_iomput .AND.  knt == nrdttrc ) THEN
          CALL iom_put( "Heup" , heup(:,:  ) * tmask(:,:,1) )  ! euphotic layer deptht
-         CALL iom_put( "PARDM", zpar(:,:,:) * tmask(:,:,:) )  ! Photosynthetically Available Radiation
-         CALL iom_put( "PAR"  , emoy(:,:,:) * tmask(:,:,:) )  ! Photosynthetically Available Radiation
+         IF( iom_use( "PAR" ) ) THEN
+            zpar(:.:,1) = zpar(:,:,1) * ( 1._wp - fr_i(:,:) )
+            CALL iom_put( "PAR", zpar(:,:,:) * tmask(:,:,:) )  ! Photosynthetically Available Radiation
+         ENDIF
       ENDIF
       !