MODULE p4zint !!====================================================================== !! *** MODULE p4zint *** !! TOP : PISCES interpolation and computation of various accessory fields !!====================================================================== !! History : 1.0 ! 2004-03 (O. Aumont) Original code !! 2.0 ! 2007-12 (C. Ethe, G. Madec) F90 !!---------------------------------------------------------------------- #if defined key_pisces !!---------------------------------------------------------------------- !! 'key_pisces' PISCES bio-model !!---------------------------------------------------------------------- !! p4z_int : interpolation and computation of various accessory fields !!---------------------------------------------------------------------- USE oce_trc ! USE trp_trc USE sms IMPLICIT NONE PRIVATE PUBLIC p4z_int ! called in p4zprg.F90 REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpk) :: & !: Tgfunc, & !: Temp. dependancy of various biological rates Tgfunc2 !: Temp. dependancy of mesozooplankton rates !! * Module variables REAL(wp) :: & xksilim = 16.5E-6 ! Half-saturation constant for the computation of the Si half-saturation constant !!---------------------------------------------------------------------- !! NEMO/TOP 2.0 , LOCEAN-IPSL (2007) !! $Header:$ !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) !!---------------------------------------------------------------------- CONTAINS SUBROUTINE p4z_int !!--------------------------------------------------------------------- !! *** ROUTINE p4z_int *** !! !! ** Purpose : interpolation and computation of various accessory fields !! !! ** Method : - ??? !!--------------------------------------------------------------------- !! INTEGER :: ji, jj REAL(wp) :: zdum !!--------------------------------------------------------------------- ! Computation of phyto and zoo metabolic rate ! ------------------------------------------- Tgfunc (:,:,:) = EXP( 0.063913 * tn(:,:,:) ) Tgfunc2(:,:,:) = EXP( 0.07608 * tn(:,:,:) ) ! Computation of the silicon dependant half saturation ! constant for silica uptake ! --------------------------------------------------- DO ji = 1, jpi DO jj = 1, jpj zdum = trn(ji,jj,1,jpsil) * trn(ji,jj,1,jpsil) xksimax(ji,jj) = MAX( xksimax(ji,jj), ( 1.+ 7.* zdum / ( xksilim * xksilim + zdum ) ) * 1e-6 ) END DO END DO IF( nday_year == 365 ) THEN xksi = xksimax xksimax = 0.e0 ENDIF ! END SUBROUTINE p4z_int #else !!====================================================================== !! Dummy module : No PISCES bio-model !!====================================================================== CONTAINS SUBROUTINE p4z_int ! Empty routine WRITE(*,*) 'p4z_int: You should not have seen this print! error?' END SUBROUTINE p4z_int #endif !!====================================================================== END MODULE p4zint