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 8542 – NEMO

Changeset 8542


Ignore:
Timestamp:
2017-09-19T12:19:04+02:00 (7 years ago)
Author:
cetlod
Message:

trunk : minor bugfixes on initialisation sequence of passive tracers

Location:
trunk/NEMOGCM/NEMO/TOP_SRC
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/TOP_SRC/trcdta.F90

    r7646 r8542  
    227227                     ik = mbkt(ji,jj)  
    228228                     IF( ik > 1 ) THEN 
    229                         zl = ( gdept_1d(ik) - gdept_n(ji,jj,ik) ) / ( gdept_1d(ik) - gdept_1d(ik-1) ) 
     229                        zl = ( gdept_1d(ik) - gdept_0(ji,jj,ik) ) / ( gdept_1d(ik) - gdept_1d(ik-1) ) 
    230230                        ptrcdta(ji,jj,ik) = (1.-zl) * ptrcdta(ji,jj,ik) + zl * ptrcdta(ji,jj,ik-1) 
    231231                     ENDIF 
    232232                     ik = mikt(ji,jj) 
    233233                     IF( ik > 1 ) THEN 
    234                         zl = ( gdept_n(ji,jj,ik) - gdept_1d(ik) ) / ( gdept_1d(ik+1) - gdept_1d(ik) ) 
     234                        zl = ( gdept_0(ji,jj,ik) - gdept_1d(ik) ) / ( gdept_1d(ik+1) - gdept_1d(ik) ) 
    235235                        ptrcdta(ji,jj,ik) = (1.-zl) * ptrcdta(ji,jj,ik) + zl * ptrcdta(ji,jj,ik+1) 
    236236                     ENDIF 
  • trunk/NEMOGCM/NEMO/TOP_SRC/trcini.F90

    r7753 r8542  
    6969      ! 
    7070      CALL trc_ini_sms   ! SMS 
    71       CALL trc_ini_inv   ! Inventories 
    7271      CALL trc_ini_trp   ! passive tracers transport 
    7372      CALL trc_ice_ini   ! Tracers in sea ice 
     
    7877      IF( nn_dttrc /= 1 ) & 
    7978      CALL trc_sub_ini    ! Initialize variables for substepping passive tracers 
     79      ! 
     80      CALL trc_ini_inv   ! Inventories 
    8081      ! 
    8182      IF( nn_timing == 1 )   CALL timing_stop('trc_init') 
Note: See TracChangeset for help on using the changeset viewer.