Changeset 14276 for NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/P4Z/p4zint.F90

Timestamp:

2021-01-07T23:09:56+01:00 (3 years ago)

Author:

aumont

Message:

numerous updates to PISCES, PISCES-QUOTA and the sediment module

File:

• NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/P4Z/p4zint.F90 (modified) (3 diffs)

NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/P4Z/p4zint.F90

@@ -36 +36 @@
       !
       INTEGER  :: ji, jj                 ! dummy loop indices
-      REAL(wp) :: zvar                   ! local variable
+      REAL(wp) :: zrum, zcodel, zargu, zvar
       !!---------------------------------------------------------------------
       !
@@ -44 +44 @@
       ! -------------------------------------------
       ! Generic temperature dependence (Eppley, 1972)
-      tgfunc (:,:,:) = EXP( 0.063913 * tsn(:,:,:,jp_tem) )
+      tgfunc (:,:,:) = EXP( 0.0631 * tsn(:,:,:,jp_tem) )
       ! Temperature dependence of mesozooplankton (Buitenhuis et al. (2005))
-      tgfunc2(:,:,:) = EXP( 0.07608  * tsn(:,:,:,jp_tem) )
-      ! Temperature dependence of picophytoplankton (Stawiarsky et al., 2016)
-      tgfunc3(:,:,:) = EXP( 0.0825   * tsn(:,:,:,jp_tem) )
+      tgfunc2(:,:,:) = EXP( 0.0761  * tsn(:,:,:,jp_tem) )
 
       ! Computation of the silicon dependant half saturation  constant for silica uptake
@@ -69 +67 @@
       ENDIF
       !
+      ! compute the day length depending on latitude and the day
+      ! Astronomical parameterization taken from HAMOCC3
+      zrum = REAL( nday_year - 80, wp ) / REAL( nyear_len(1), wp )
+      zcodel = ASIN(  SIN( zrum * rpi * 2._wp ) * SIN( rad * 23.5_wp )  )
+
+      ! 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
+      !
       IF( ln_timing )   CALL timing_stop('p4z_int')
       !