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 728 for trunk/NEMO/TOP_SRC/SMS/trcini.pisces.h90 – NEMO

Ignore:
Timestamp:
2007-10-22T12:01:39+02:00 (17 years ago)
Author:
cetlod
Message:

add tmask_i in global sum, see ticket:17

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/TOP_SRC/SMS/trcini.pisces.h90

    r719 r728  
    3737      !!---------------------------------------------------------------------- 
    3838      !!  TOP 1.0 , LOCEAN-IPSL (2005)  
    39       !! $Header$  
     39      !! $Header: /home/opalod/NEMOCVSROOT/NEMO/TOP_SRC/SMS/trcini.pisces.h90,v 1.9 2007/10/12 09:35:04 opalod Exp $  
    4040      !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt  
    4141      !!---------------------------------------------------------------------- 
     
    169169            DO ji = 2, jpim1 
    170170               sumdepsi = sumdepsi + dustmo(ji,jj,jm)/(12.*rmoss)*8.8        & 
    171                   *0.075/28.1*e1t(ji,jj)*e2t(ji,jj)*tmask(ji,jj,1) 
     171                  *0.075/28.1*e1t(ji,jj)*e2t(ji,jj)*tmask(ji,jj,1)*tmask_i(ji,jj) 
    172172            END DO 
    173173         END DO 
     
    200200         DO ji=2,jpim1 
    201201            rivpo4input=rivpo4input+rivinp(ji,jj)*(e1t(ji,jj)*e2t(ji,jj)    & 
    202                *fse3t(ji,jj,1))*tmask(ji,jj,1)*raass 
     202               *fse3t(ji,jj,1))*tmask(ji,jj,1)*tmask_i(ji,jj)*raass 
    203203            rivalkinput=rivalkinput+cotdep(ji,jj)*(e1t(ji,jj)*e2t(ji,jj)    & 
    204                *fse3t(ji,jj,1))*tmask(ji,jj,1)*raass 
     204               *fse3t(ji,jj,1))*tmask(ji,jj,1)*tmask_i(ji,jj)*raass 
    205205            nitdepinput=nitdepinput+nitdep(ji,jj)*(e1t(ji,jj)*e2t(ji,jj)    & 
    206                *fse3t(ji,jj,1))*tmask(ji,jj,1)*raass 
     206               *fse3t(ji,jj,1))*tmask(ji,jj,1)*tmask_i(ji,jj)*raass 
    207207         END DO 
    208208      END DO 
Note: See TracChangeset for help on using the changeset viewer.