Changeset 14416 for NEMO/branches/2021/dev_r14383_PISCES_NEWDEV_PISCO/src/TOP/PISCES/P4Z/p4zint.F90

Timestamp:

2021-02-07T16:08:48+01:00 (3 years ago)

Author:

cetlod

Message:

dev_r14383_PISCES_NEWDEV_PISCO : minor improvments

File:

• NEMO/branches/2021/dev_r14383_PISCES_NEWDEV_PISCO/src/TOP/PISCES/P4Z/p4zint.F90 (modified) (2 diffs)

NEMO/branches/2021/dev_r14383_PISCES_NEWDEV_PISCO/src/TOP/PISCES/P4Z/p4zint.F90

@@ -37 +37 @@
       INTEGER, INTENT( in ) ::   Kbb, Kmm ! time level indices
       !
-      INTEGER  :: ji, jj                 ! dummy loop indices
+      INTEGER  :: ji, jj, jk              ! dummy loop indices
       REAL(wp) :: zrum, zcodel, zargu, zvar
       !!---------------------------------------------------------------------
@@ -75 +75 @@
       ! day length in hours
       strn(:,:) = 0.
-      DO jj = 1, jpj
-         DO ji = 1, jpi
-            zargu = TAN( zcodel ) * TAN( gphit(ji,jj) * rad )
-            zargu = MAX( -1., MIN(  1., zargu ) )
-            strn(ji,jj) = MAX( 0.0, 24. - 2. * ACOS( zargu ) / rad / 15. )
-         END DO
-      END DO
+      DO_2D( 1, 1, 1, 1 )
+         zargu = TAN( zcodel ) * TAN( gphit(ji,jj) * rad )
+         zargu = MAX( -1., MIN(  1., zargu ) )
+         strn(ji,jj) = MAX( 0.0, 24. - 2. * ACOS( zargu ) / rad / 15. )
+      END_2D
+      !
+      DO_3D( 1, 1, 1, 1, 1, jpkm1 )
+        ! denitrification factor computed from O2 levels
+         ! This factor diagnoses below which level of O2 denitrification
+         ! is active
+         nitrfac(ji,jj,jk) = MAX(  0.e0, 0.4 * ( 6.e-6  - tr(ji,jj,jk,jpoxy,Kbb) )    &
+            &                                / ( oxymin + tr(ji,jj,jk,jpoxy,Kbb) )  )
+         nitrfac(ji,jj,jk) = MIN( 1., nitrfac(ji,jj,jk) )
+         !
+         ! redox factor computed from NO3 levels
+         ! This factor diagnoses below which level of NO3 additional redox
+         ! reactions are taking place.
+         nitrfac2(ji,jj,jk) = MAX( 0.e0,       ( 1.E-6 - tr(ji,jj,jk,jpno3,Kbb) )  &
+            &                                / ( 1.E-6 + tr(ji,jj,jk,jpno3,Kbb) ) )
+         nitrfac2(ji,jj,jk) = MIN( 1., nitrfac2(ji,jj,jk) )
+      END_3D
       !
       IF( ln_timing )   CALL timing_stop('p4z_int')