New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 725 for trunk/NEMO – NEMO

Changeset 725 for trunk/NEMO


Ignore:
Timestamp:
2007-10-22T11:54:12+02:00 (16 years ago)
Author:
cetlod
Message:

update of iron aggregation computation, see ticket:14

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/TOP_SRC/SMS/p4zrem.F

    r719 r725  
    11 
    2 CCC $Header$  
     2CCC $Header: /home/opalod/NEMOCVSROOT/NEMO/TOP_SRC/SMS/p4zrem.F,v 1.8 2007/10/12 09:28:41 opalod Exp $  
    33CCC  TOP 1.0 , LOCEAN-IPSL (2005)  
    44C This software is governed by CeCILL licence see modipsl/doc/NEMO_CeCILL.txt  
     
    5050      REAL zsatur,zsatur2,znusil,zdepbac(jpi,jpj,jpk) 
    5151      REAL zlamfac,zstep,fesatur(jpi,jpj,jpk) 
     52      REAL ztempbac(jpi,jpj) 
    5253C 
    5354C      Time step duration for the biology 
     
    5960C      -------------------------------------------------- 
    6061C 
    61         DO jk=1,12 
    62          zdepbac(:,:,jk)=min(0.7*(trn(:,:,jk,jpzoo)+2*trn(:,:,jk,jpmes)) 
    63      &     ,4E-6) 
     62        DO jk=1,jpkm1 
     63          DO jj = 1, jpj 
     64            DO ji = 1, jpi 
     65         IF (fsdept(ji,jj,jk).lt.120.) THEN 
     66         zdepbac(ji,jj,jk)=min(0.7*(trn(ji,jj,jk,jpzoo) 
     67     &     +2*trn(ji,jj,jk,jpmes)),4E-6) 
     68         ztempbac(ji,jj)=zdepbac(ji,jj,jk) 
     69         ELSE 
     70         zdepbac(ji,jj,jk)=min(1.,120./fsdept(ji,jj,jk)) 
     71     &      *ztempbac(ji,jj) 
     72         ENDIF 
     73            END DO 
     74          END DO 
    6475        END DO 
    65 C 
    66 C      Vertical decay of the bacterial activity 
    67 C      ---------------------------------------- 
    68 C 
    69          do jk=13,jpk 
    70            do jj=1,jpj 
    71              do ji=1,jpi 
    72          zdepbac(ji,jj,jk)=min(1.,fsdept(ji,jj,12)/fsdept(ji,jj,jk)) 
    73      &      *zdepbac(ji,jj,12) 
    74              end do 
    75            end do 
    76          end do 
    7776 
    7877         DO jk = 1,jpkm1 
     
    253252         zlamfac=min(1.,zlamfac) 
    254253#if ! defined key_trc_kriest 
    255          xlam1b=(80.*(trn(ji,jj,jk,jpdoc)+40E-6)+698. 
     254         xlam1b=(80.*(trn(ji,jj,jk,jpdoc)+35E-6)+698. 
    256255     &    *trn(ji,jj,jk,jppoc)+1.05E4*trn(ji,jj,jk,jpgoc)) 
    257      &    *zdiss(ji,jj,jk)+1E-5*(1.-zlamfac)+xlam1*max(0., 
     256     &    *zdiss(ji,jj,jk)+1E-4*(1.-zlamfac)+xlam1*max(0., 
    258257     &    (trn(ji,jj,jk,jpfer)*1E9-1.)) 
    259258#else 
    260          xlam1b=(80.*(trn(ji,jj,jk,jpdoc)+40E-6)+698. 
     259         xlam1b=(80.*(trn(ji,jj,jk,jpdoc)+35E-6)+698. 
    261260     &    *trn(ji,jj,jk,jppoc)) 
    262      &    *zdiss(ji,jj,jk)+1E-5*(1.-zlamfac)+xlam1*max(0., 
     261     &    *zdiss(ji,jj,jk)+1E-4*(1.-zlamfac)+xlam1*max(0., 
    263262     &    (trn(ji,jj,jk,jpfer)*1E9-1.)) 
    264263#endif 
    265264 
    266265 
    267          xaggdfe(ji,jj,jk)=xlam1b*zstep*0.76*(trn(ji,jj,jk,jpfer) 
     266         xaggdfe(ji,jj,jk)=xlam1b*zstep*0.5*(trn(ji,jj,jk,jpfer) 
    268267     &     -xfeequi) 
    269268#    if defined key_off_degrad 
Note: See TracChangeset for help on using the changeset viewer.