Changeset 14209

Timestamp:

2020-12-18T10:53:36+01:00 (3 years ago)

Author:

cetlod

Message:

trunk : bugfix to diagnose PAR in PISCES as a mean over the whole cell, see ticket #2507

Location:

NEMO/trunk

Files:

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

NEMO/trunk/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/trunk/src/TOP/PISCES/P4Z/p4zopt.F90

@@ -236 +236 @@
       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
       !