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 2528 for trunk/NEMOGCM/NEMO/TOP_SRC/LOBSTER/trcbio.F90 – NEMO

Ignore:
Timestamp:
2010-12-27T18:33:53+01:00 (13 years ago)
Author:
rblod
Message:

Update NEMOGCM from branch nemo_v3_3_beta

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/NEMOGCM/NEMO/TOP_SRC/LOBSTER/trcbio.F90

    r1795 r2528  
    2020   USE lbclnk          !  
    2121   USE prtctl_trc      ! Print control for debbuging 
    22    USE trdmld_trc 
    23    USE trdmld_trc_oce 
     22   USE trdmod_oce 
     23   USE trdmod_trc 
    2424   USE iom 
    2525    
    … …  
    3232#  include "top_substitute.h90" 
    3333   !!---------------------------------------------------------------------- 
    34    !! NEMO/TOP 2.0 , LOCEAN-IPSL (2007)  
     34   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    3535   !! $Id$  
    36    !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) 
     36   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    3737   !!---------------------------------------------------------------------- 
    3838 
    … …  
    5757      !!                                  source      sink 
    5858      !!         
    59       !!              IF 'key_trc_diabio' defined , the biogeochemical trends 
     59      !!              IF 'key_diabio' defined , the biogeochemical trends 
    6060      !!              for passive tracers are saved for futher diagnostics. 
    6161      !!--------------------------------------------------------------------- 
    … …  
    7171      REAL(wp) ::   zfilpz, zfildz, zphya, zzooa, zno3a 
    7272      REAL(wp) ::   znh4a, zdeta, zdoma, zzoobod, zboddet, zdomaju 
    73 #if defined key_trc_diaadd 
     73#if defined key_diatrc 
    7474      REAL(wp) ::   ze3t 
    7575#endif 
    76 #if defined key_trc_diaadd && defined key_iomput 
     76#if defined key_diatrc && defined key_iomput 
    7777      REAL(wp), DIMENSION(jpi,jpj,17)     :: zw2d 
    78 # if defined key_trc_dia3d 
    7978      REAL(wp), DIMENSION(jpi,jpj,jpk,3) :: zw3d 
    80 # endif 
    8179#endif 
    8280      REAL(wp), DIMENSION(:,:,:,:), ALLOCATABLE ::   ztrbio 
    … …  
    9189 
    9290      fbod(:,:) = 0.e0 
    93 #if defined key_trc_diaadd && ! defined key_iomput 
     91#if defined key_diatrc && ! defined key_iomput 
    9492      DO jl = jp_lob0_2d, jp_lob1_2d 
    9593         trc2d(:,:,jl) = 0.e0 
    9694      END DO  
    9795#endif 
    98 #if defined key_trc_diaadd && defined key_iomput 
     96#if defined key_diatrc && defined key_iomput 
    9997      zw2d(:,:,:) = 0.e0 
    100 # if defined key_trc_dia3d 
    10198      zw3d(:,:,:,:) = 0.e0 
    102 # endif 
    10399#endif 
    104100 
    … …  
    117113 
    118114               ! negative trophic variables DO not contribute to the fluxes 
    119                zdet = MAX( 0.e0, trn(ji,jj,jk,jpdet) ) 
    120                zzoo = MAX( 0.e0, trn(ji,jj,jk,jpzoo) ) 
    121                zphy = MAX( 0.e0, trn(ji,jj,jk,jpphy) ) 
    122                zno3 = MAX( 0.e0, trn(ji,jj,jk,jpno3) ) 
    123                znh4 = MAX( 0.e0, trn(ji,jj,jk,jpnh4) ) 
    124                zdom = MAX( 0.e0, trn(ji,jj,jk,jpdom) ) 
     115               zdet = MAX( 0.e0, trn(ji,jj,jk,jp_lob_det) ) 
     116               zzoo = MAX( 0.e0, trn(ji,jj,jk,jp_lob_zoo) ) 
     117               zphy = MAX( 0.e0, trn(ji,jj,jk,jp_lob_phy) ) 
     118               zno3 = MAX( 0.e0, trn(ji,jj,jk,jp_lob_no3) ) 
     119               znh4 = MAX( 0.e0, trn(ji,jj,jk,jp_lob_nh4) ) 
     120               zdom = MAX( 0.e0, trn(ji,jj,jk,jp_lob_dom) ) 
    125121 
    126122               ! Limitations 
    … …  
    194190 
    195191               ! tracer flux at totox-point added to the general trend 
    196                tra(ji,jj,jk,jpdet) = tra(ji,jj,jk,jpdet) + zdeta 
    197                tra(ji,jj,jk,jpzoo) = tra(ji,jj,jk,jpzoo) + zzooa 
    198                tra(ji,jj,jk,jpphy) = tra(ji,jj,jk,jpphy) + zphya 
    199                tra(ji,jj,jk,jpno3) = tra(ji,jj,jk,jpno3) + zno3a 
    200                tra(ji,jj,jk,jpnh4) = tra(ji,jj,jk,jpnh4) + znh4a 
    201                tra(ji,jj,jk,jpdom) = tra(ji,jj,jk,jpdom) + zdoma 
    202  
    203 #if defined key_trc_diabio 
     192               tra(ji,jj,jk,jp_lob_det) = tra(ji,jj,jk,jp_lob_det) + zdeta 
     193               tra(ji,jj,jk,jp_lob_zoo) = tra(ji,jj,jk,jp_lob_zoo) + zzooa 
     194               tra(ji,jj,jk,jp_lob_phy) = tra(ji,jj,jk,jp_lob_phy) + zphya 
     195               tra(ji,jj,jk,jp_lob_no3) = tra(ji,jj,jk,jp_lob_no3) + zno3a 
     196               tra(ji,jj,jk,jp_lob_nh4) = tra(ji,jj,jk,jp_lob_nh4) + znh4a 
     197               tra(ji,jj,jk,jp_lob_dom) = tra(ji,jj,jk,jp_lob_dom) + zdoma 
     198 
     199#if defined key_diabio 
    204200               trbio(ji,jj,jk,jp_lob0_trd     ) = zno3phy 
    205201               trbio(ji,jj,jk,jp_lob0_trd +  1) = znh4phy 
    … …  
    238234                ENDIF 
    239235 
    240 #if defined key_trc_diaadd 
     236#if defined key_diatrc 
    241237               ! convert fluxes in per day 
    242238               ze3t = fse3t(ji,jj,jk) * 86400. 
    … …  
    282278               zw2d(ji,jj,17) = zw2d(ji,jj,17) + zdetdom * ze3t 
    283279#endif 
    284 #if defined key_trc_dia3d  
     280#if defined key_diatrc  
    285281# if ! defined key_iomput 
    286282               trc3d(ji,jj,jk,jp_lob0_3d    ) = zno3phy * 86400      
    … …  
    307303               !    trophic variables( det, zoo, phy, no3, nh4, dom) 
    308304               !       negative trophic variables DO not contribute to the fluxes 
    309                zdet = MAX( 0.e0, trn(ji,jj,jk,jpdet) ) 
    310                zzoo = MAX( 0.e0, trn(ji,jj,jk,jpzoo) ) 
    311                zphy = MAX( 0.e0, trn(ji,jj,jk,jpphy) ) 
    312                zno3 = MAX( 0.e0, trn(ji,jj,jk,jpno3) ) 
    313                znh4 = MAX( 0.e0, trn(ji,jj,jk,jpnh4) ) 
    314                zdom = MAX( 0.e0, trn(ji,jj,jk,jpdom) ) 
     305               zdet = MAX( 0.e0, trn(ji,jj,jk,jp_lob_det) ) 
     306               zzoo = MAX( 0.e0, trn(ji,jj,jk,jp_lob_zoo) ) 
     307               zphy = MAX( 0.e0, trn(ji,jj,jk,jp_lob_phy) ) 
     308               zno3 = MAX( 0.e0, trn(ji,jj,jk,jp_lob_no3) ) 
     309               znh4 = MAX( 0.e0, trn(ji,jj,jk,jp_lob_nh4) ) 
     310               zdom = MAX( 0.e0, trn(ji,jj,jk,jp_lob_dom) ) 
    315311 
    316312               !    Limitations 
    … …  
    363359 
    364360               ! tracer flux at totox-point added to the general trend 
    365                tra(ji,jj,jk,jpdet) = tra(ji,jj,jk,jpdet) + zdeta 
    366                tra(ji,jj,jk,jpzoo) = tra(ji,jj,jk,jpzoo) + zzooa 
    367                tra(ji,jj,jk,jpphy) = tra(ji,jj,jk,jpphy) + zphya 
    368                tra(ji,jj,jk,jpno3) = tra(ji,jj,jk,jpno3) + zno3a 
    369                tra(ji,jj,jk,jpnh4) = tra(ji,jj,jk,jpnh4) + znh4a 
    370                tra(ji,jj,jk,jpdom) = tra(ji,jj,jk,jpdom) + zdoma 
     361               tra(ji,jj,jk,jp_lob_det) = tra(ji,jj,jk,jp_lob_det) + zdeta 
     362               tra(ji,jj,jk,jp_lob_zoo) = tra(ji,jj,jk,jp_lob_zoo) + zzooa 
     363               tra(ji,jj,jk,jp_lob_phy) = tra(ji,jj,jk,jp_lob_phy) + zphya 
     364               tra(ji,jj,jk,jp_lob_no3) = tra(ji,jj,jk,jp_lob_no3) + zno3a 
     365               tra(ji,jj,jk,jp_lob_nh4) = tra(ji,jj,jk,jp_lob_nh4) + znh4a 
     366               tra(ji,jj,jk,jp_lob_dom) = tra(ji,jj,jk,jp_lob_dom) + zdoma 
    371367               ! 
    372 #if defined key_trc_diabio 
     368#if defined key_diabio 
    373369               trbio(ji,jj,jk,jp_lob0_trd     ) = zno3phy 
    374370               trbio(ji,jj,jk,jp_lob0_trd +  1) = znh4phy 
    … …  
    406402                  !  trend number 17 in trcexp 
    407403                ENDIF 
    408 #if defined key_trc_diaadd && defined key_trc_dia3d 
     404#if defined key_diatrc 
    409405# if ! defined key_iomput 
    410406               trc3d(ji,jj,jk,jp_lob0_3d    ) =  zno3phy * 86400      
    … …  
    421417      END DO 
    422418 
    423 #if defined key_trc_diaadd 
     419#if defined key_diatrc 
    424420      ! Lateral boundary conditions  
    425421# if ! defined key_iomput 
    … …  
    452448#endif 
    453449 
    454 #if defined key_trc_diaadd && defined key_trc_dia3d 
     450#if defined key_diatrc 
    455451      ! Lateral boundary conditions  
    456452# if ! defined key_iomput 
    … …  
    469465#endif 
    470466 
    471 #if defined key_trc_diabio 
     467#if defined key_diabio 
    472468      ! Lateral boundary conditions on trcbio 
    473469      DO jl = jp_lob0_trd, jp_lob1_trd 
Note: See TracChangeset for help on using the changeset viewer.