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 1317 for trunk/NEMO/OPA_SRC/TRD/trdmld.F90 – NEMO

Ignore:
Timestamp:
2009-02-17T17:31:42+01:00 (15 years ago)
Author:
smasson
Message:

nwrite = modulo referenced to nit000 in all ouputs, see ticket:339

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/OPA_SRC/TRD/trdmld.F90

    r1312 r1317  
    389389      smltrd(:,:,:) = smltrd(:,:,:) * ucf   !  is no longer used, and is reset to 0. at next time step) 
    390390       
    391       MODULO_NTRD : IF( MOD( kt, ntrd ) == 0 ) THEN        ! nitend MUST be multiple of ntrd 
     391      it = kt - nit000 + 1 
     392 
     393      MODULO_NTRD : IF( MOD( it, ntrd ) == 0 ) THEN        ! nitend MUST be multiple of ntrd 
    392394         ! 
    393395         ztmltot (:,:) = 0.e0   ;   zsmltot (:,:) = 0.e0   ! reset arrays to zero 
     
    574576#if defined key_dimgout 
    575577 
    576       IF( MOD( kt, ntrd ) == 0 ) THEN 
     578      IF( MOD( it, ntrd ) == 0 ) THEN 
    577579         iyear =  ndastp/10000 
    578580         imon  = (ndastp-iyear*10000)/100 
     
    591593      ! ---------------------------------- 
    592594 
    593       IF( lwp .AND. MOD( kt , ntrd ) == 0 ) THEN 
     595      IF( lwp .AND. MOD( it , ntrd ) == 0 ) THEN 
    594596         WRITE(numout,*) ' ' 
    595597         WRITE(numout,*) 'trd_mld : write trends in the NetCDF file :' 
     
    600602      END IF 
    601603          
    602       it = kt - nit000 + 1 
    603  
    604604      !-- Write the trends for T/S instantaneous diagnostics  
    605605      IF( ln_trdmld_instant ) THEN            
     
    683683#endif 
    684684 
    685       IF( MOD( kt , ntrd ) == 0 ) THEN 
     685      IF( MOD( it, ntrd ) == 0 ) THEN 
    686686         ! 
    687687         ! III.5 Reset cumulative arrays to zero 
Note: See TracChangeset for help on using the changeset viewer.