Ignore:
Timestamp:
2009-02-03T14:45:10+01:00 (12 years ago)
Author:
cetlod
Message:

computation of horizontal derivatives at the first time step, see ticket:315

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/TOP_SRC/trcini.F90

    r1254 r1286  
    1818   USE oce_trc 
    1919   USE trc 
     20   USE trp_trc 
    2021   USE trcrst 
    2122   USE trcctl 
     
    2627   USE trcini_c14b     ! C14 bomb initialisation 
    2728   USE trcini_my_trc   ! MY_TRC   initialisation 
    28    USE trcdta          !  
    29    USE trdmld_trc_oce 
    30    USE trdmld_trc 
    31  
     29   USE trcdta     
     30   USE zpshde_trc      ! partial step: hor. derivative  
    3231   USE in_out_manager  ! I/O manager 
    3332   USE prtctl_trc      ! Print control passive tracers (prt_ctl_trc_init routine) 
     
    108107 
    109108      IF( .NOT. lrsttr ) THEN  
    110             
    111109# if defined key_dtatrc 
    112110         ! Initialization of tracer from a file that may also be used for damping 
     
    117115# endif 
    118116         trb(:,:,:,:) = trn(:,:,:,:) 
    119  
    120117      ELSE 
    121  
    122118         CALL trc_rst_read      ! restart from a file 
    123  
    124119      ENDIF 
    125120 
    126       IF( lk_trdmld_trc  )  CALL trd_mld_trc_init     ! trends: Mixed-layer 
     121      tra(:,:,:,:) = 0. 
     122       
     123      IF( ln_zps .AND. .NOT. lk_trc_c1d )   &              ! Partial steps: before horizontal gradient of passive 
     124      &                     CALL zps_hde_trc( nittrc000, trb, gtru, gtrv )       ! tracers at the bottom ocean level 
     125 
    127126 
    128127      !                 ! Computation content of all tracers 
Note: See TracChangeset for help on using the changeset viewer.