Changeset 1180


Ignore:
Timestamp:
2008-09-11T18:41:47+02:00 (12 years ago)
Author:
cetlod
Message:

update PISCES modules to couple with the sediment model, see ticket:249

Location:
trunk/NEMO/TOP_SRC/PISCES
Files:
6 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/TOP_SRC/PISCES/p4zche.F90

    r1152 r1180  
    105105      kf4  =-0.001005 
    106106 
    107  
    108107   REAL(wp) :: &              ! coeff. for 1. dissoc. of boric acid (Dickson and Goyet, 1994) 
    109108      cb0  = -8966.90, & 
     
    129128      cw6 = -0.01615 
    130129  
    131    REAL(wp) :: &             ! coeff. for dissoc. of phosphate (Millero (1974) 
    132       cp10 = 115.54    , & 
    133       cp11 = -4576.752 , & 
    134       cp12 = -18.453   , & 
    135       cp13 = -106.736  , & 
    136       cp14 = 0.69171   , & 
    137       cp15 = -0.65643  , & 
    138       cp16 = -0.01844  , & 
    139       cp20 = 172.1033  , & 
    140       cp21 = -8814.715 , & 
    141       cp22 = -27.927   , & 
    142       cp23 = -160.340  , & 
    143       cp24 = 1.3566    , & 
    144       cp25 = 0.37335   , & 
    145       cp26 = -0.05778  , & 
    146       cp30 = -18.126   , & 
    147       cp31 = -3070.75  , & 
    148       cp32 = 17.27039  , & 
    149       cp33 = 2.81197   , & 
    150       cp34 = -44.99486 , & 
    151       cp35 = -0.09984 
    152  
    153    REAL(wp) :: &             ! coeff. for dissoc. of silicates (Millero (1974)   
    154       cs10 = 117.385   , &   
    155       cs11 = -8904.2   , &  
    156       cs12 = -19.334   , &  
    157       cs13 = -458.79   , &  
    158       cs14 = 3.5913    , &  
    159       cs15 = 188.74    , &  
    160       cs16 = -1.5998   , &  
    161       cs17 = -12.1652  , &  
    162       cs18 = 0.07871   , &  
    163       cs19 = 0.        , &  
    164       cs20 = 1.        , &  
    165       cs21 = -0.001005 
    166  
    167130   REAL(wp) :: &              ! volumetric solubility constants for o2 in ml/l (Weiss, 1974) 
    168131      ox0 = -58.3877   , & 
  • trunk/NEMO/TOP_SRC/PISCES/p4zopt.F90

    r1152 r1180  
    2727      etot, enano, ediat,       &  !: PAR for phyto, nano and diat  
    2828      emoy                         !: averaged PAR in the mixed layer 
    29    REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   &   !: 
    30       heup                         !: Depth of the euphotic zone 
    3129 
    3230   !! * Module variables 
     
    210208!     Computation of the euphotic depth 
    211209!     --------------------------------- 
    212      
    213       heup(:,:) = 300.e0 
     210      ! Euphotic layer bottom level 
     211      neln(:,:) = 1                                           ! initialisation of EL level 
     212      heup(:,:) = 300. 
    214213 
    215214      DO jk = 2, jpkm1 
    216215         DO jj = 1, jpj 
    217             DO ji = 1, jpi 
    218                IF( etot(ji,jj,jk) >= 0.0043 * qsr(ji,jj) )   heup(ji,jj) = fsdepw(ji,jj,jk+1) 
    219             END DO 
    220          END DO 
    221       END DO 
    222  
    223       heup(:,:) = MIN( 300., heup(:,:) ) 
     216           DO ji = 1, jpi 
     217              IF( etot(ji,jj,jk) >= 0.0043 * qsr(ji,jj) )  THEN 
     218                 neln(ji,jj) = jk+1 ! 1rst T-level strictly below EL bottom 
     219              !                                                  ! nb. this is to ensure compatibility with 
     220              !                                                  ! nmld_trc definition in trd_mld_trc_zint 
     221                heup(ji,jj) = fsdepw(ji,jj,jk+1)                 ! Euphotic layer depth 
     222             ENDIF 
     223          END DO 
     224        END DO 
     225     ENDDO 
     226  
     227     heup(:,:) = MIN( 300., heup(:,:) ) 
    224228 
    225229!    Computation of the mean light over the mixed layer depth 
  • trunk/NEMO/TOP_SRC/PISCES/p4zprod.F90

    r1152 r1180  
    1717   USE sms_pisces      !  
    1818   USE prtctl_trc 
     19   USE p4zopt 
    1920   USE p4zint 
    2021   USE p4zlim 
    21    USE p4zopt 
    2222 
    2323   IMPLICIT NONE 
  • trunk/NEMO/TOP_SRC/PISCES/p4zsed.F90

    r1152 r1180  
    6363   !!---------------------------------------------------------------------- 
    6464   !! NEMO/TOP 2.0 , LOCEAN-IPSL (2007)  
    65    !! $Id$  
     65   !! $Header:$  
    6666   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) 
    6767   !!---------------------------------------------------------------------- 
     
    8282      INTEGER  ::   ji, jj, jk 
    8383      INTEGER  ::   ikt 
     84#if ! defined key_sed 
    8485      REAL(wp) ::   zsumsedsi, zsumsedpo4, zsumsedcal 
     86#endif 
    8587      REAL(wp) ::   zconctmp , zdenitot  , znitrpottot 
    8688      REAL(wp) ::   zlim, zconctmp2, zstep, zfact 
     
    137139         &       + zirondep(:,:,jk) + ironsed(:,:,jk) * rfact2 
    138140      END DO 
     141 
     142 
     143#if ! defined key_sed 
    139144 
    140145      ! Initialisation of variables used to compute Sinking Speed 
     
    174179         CALL mpp_sum( zsumsedpo4 )   ! sums over the global domain 
    175180      ENDIF 
     181 
     182#endif 
    176183 
    177184      ! Then this loss is scaled at each bottom grid cell for 
     
    192199 
    193200            trn(ji,jj,ikt,jpdsi) = trn(ji,jj,ikt,jpdsi) - zconctmp 
     201 
     202#if ! defined key_sed 
    194203            trn(ji,jj,ikt,jpsil) = trn(ji,jj,ikt,jpsil) + zconctmp   & 
    195204            &      * ( 1.- ( sumdepsi + rivalkinput / raass / 6. ) / zsumsedsi ) 
     205#endif 
    196206         END DO 
    197207      END DO 
     
    202212            zconctmp = trn(ji,jj,ikt,jpcal) * wscal(ji,jj,ikt) * zstep / fse3t(ji,jj,ikt) 
    203213            trn(ji,jj,ikt,jpcal) = trn(ji,jj,ikt,jpcal) - zconctmp 
     214 
     215#if ! defined key_sed 
    204216            trn(ji,jj,ikt,jptal) = trn(ji,jj,ikt,jptal) + zconctmp   & 
    205217               &   * ( 1.- ( rivalkinput / raass ) / zsumsedcal ) * 2.e0 
    206218            trn(ji,jj,ikt,jpdic) = trn(ji,jj,ikt,jpdic) + zconctmp   & 
    207219               &   * ( 1.- ( rivalkinput / raass ) / zsumsedcal ) 
     220#endif 
    208221         END DO 
    209222      END DO 
     
    219232            trn(ji,jj,ikt,jpgoc) = trn(ji,jj,ikt,jpgoc) - zconctmp  * wsbio4(ji,jj,ikt) * zfact 
    220233            trn(ji,jj,ikt,jppoc) = trn(ji,jj,ikt,jppoc) - zconctmp2 * wsbio3(ji,jj,ikt) * zfact 
     234#if ! defined key_sed 
    221235            trn(ji,jj,ikt,jpdoc) = trn(ji,jj,ikt,jpdoc)    & 
    222236            &      + ( zconctmp  * wsbio4(ji,jj,ikt) + zconctmp2 * wsbio3(ji,jj,ikt) ) * zfact   & 
    223237            &      * ( 1.- rivpo4input / (raass * zsumsedpo4 ) ) 
     238#endif 
    224239            trn(ji,jj,ikt,jpbfe) = trn(ji,jj,ikt,jpbfe) - trn(ji,jj,ikt,jpbfe) * wsbio4(ji,jj,ikt) * zfact 
    225240            trn(ji,jj,ikt,jpsfe) = trn(ji,jj,ikt,jpsfe) - trn(ji,jj,ikt,jpsfe) * wsbio3(ji,jj,ikt) * zfact 
     
    233248            trn(ji,jj,ikt,jppoc) = trn(ji,jj,ikt,jppoc)   & 
    234249            &      - zconctmp2 * wsbio3(ji,jj,ikt) * zfact 
     250#if ! defined key_sed 
    235251            trn(ji,jj,ikt,jpdoc) = trn(ji,jj,ikt,jpdoc)    & 
    236252            &      + ( zconctmp2 * wsbio3(ji,jj,ikt) )   & 
    237253            &      * zfact * ( 1.- rivpo4input / ( raass * zsumsedpo4 ) ) 
     254#endif 
    238255            trn(ji,jj,ikt,jpsfe) = trn(ji,jj,ikt,jpsfe)   & 
    239256            &      - trn(ji,jj,ikt,jpsfe) * wsbio3(ji,jj,ikt) * zfact 
  • trunk/NEMO/TOP_SRC/PISCES/p4zsink.F90

    r1152 r1180  
    2727 
    2828   !! * Module variables 
    29    REAL(wp), DIMENSION(jpi,jpj,jpk) ::   &   !: 
     29   REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpk) ::   &   !: 
    3030     sinking, sinking2,   &    !: POC sinking fluxes (different meanings depending on the parameterization 
    3131     sinkcal, sinksil,    &    !: CaCO3 and BSi sinking fluxes 
  • trunk/NEMO/TOP_SRC/PISCES/sms_pisces.F90

    r1152 r1180  
    4747   !!--------------------------------------------- 
    4848   REAL(wp), DIMENSION(jpi,jpj,jpk) ::   etot3                 !: ??? 
     49   INTEGER , DIMENSION(jpi,jpj)     ::   neln    !: number of levels in the euphotic layer 
     50   REAL(wp), DIMENSION(jpi,jpj)     ::   heup    !: euphotic layer depth 
    4951 
    5052!!---------------------------------------------------------- 
Note: See TracChangeset for help on using the changeset viewer.