Changeset 2194


Ignore:
Timestamp:
2010-10-09T00:08:28+02:00 (11 years ago)
Author:
cetlod
Message:

ensure restartability of PISCES

Location:
branches/DEV_r2006_merge_TRA_TRC/NEMO/TOP_SRC
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • branches/DEV_r2006_merge_TRA_TRC/NEMO/TOP_SRC/TRP/trcadv.F90

    r2124 r2194  
    6767      IF( kt == nit000 )   CALL trc_adv_ctl          ! initialisation & control of options 
    6868 
     69#if ! defined key_pisces 
    6970      IF( neuler == 0 .AND. kt == nit000 ) THEN     ! at nit000 
    7071         r2dt(:) =  rdttra(:) * FLOAT(nn_dttrc)          ! = rdtra (restarting with Euler time stepping) 
     
    7273         r2dt(:) = 2. * rdttra(:) * FLOAT(nn_dttrc)      ! = 2 rdttra (leapfrog) 
    7374      ENDIF 
     75#else 
     76      r2dt(:) =  rdttra(:) * FLOAT(nn_dttrc)          ! = rdtra (restarting with Euler time stepping) 
     77#endif 
    7478 
    7579      !                                                   ! effective transport 
  • branches/DEV_r2006_merge_TRA_TRC/NEMO/TOP_SRC/TRP/trctrp.F90

    r2082 r2194  
    7474#endif 
    7575                                CALL trc_zdf( kstp )            ! vertical mixing and after tracer fields 
     76                                CALL trc_nxt( kstp )            ! tracer fields at next time step      
     77         IF( ln_trcrad )        CALL trc_rad( kstp )            ! Correct artificial negative concentrations 
    7678         IF( ln_zps    )        CALL zps_hde( kstp, jptra, trn, gtru, gtrv )  ! Partial steps: now horizontal gradient of passive 
    7779                                                                ! tracers at the bottom ocean level 
    78                                 CALL trc_nxt( kstp )            ! tracer fields at next time step      
    79          IF( ln_trcrad )        CALL trc_rad( kstp )            ! Correct artificial negative concentrations 
    8080         ! 
    8181      ELSE                                               ! 1D vertical configuration 
  • branches/DEV_r2006_merge_TRA_TRC/NEMO/TOP_SRC/trcini.F90

    r2104 r2194  
    133133       
    134134      IF( ln_zps .AND. .NOT. lk_trc_c1d )   &              ! Partial steps: before horizontal gradient of passive 
    135       &                     CALL zps_hde( nit000, jptra, trb, gtru, gtrv )       ! tracers at the bottom ocean level 
     135      &                     CALL zps_hde( nit000, jptra, trn, gtru, gtrv )       ! tracers at the bottom ocean level 
    136136 
    137137 
Note: See TracChangeset for help on using the changeset viewer.