Ignore:
Timestamp:
2013-11-15T11:40:36+01:00 (8 years ago)
Author:
clem
Message:

corrected step.F90

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2013/dev_LOCEAN_2013/NEMOGCM/NEMO/OPA_SRC/step.F90

    r4161 r4215  
    115115      ! 
    116116      !  VERTICAL PHYSICS 
    117       ! bg jchanut tschanges 
    118       ! One need bottom friction parameter in ssh_wzv routine with time splitting. 
    119       ! The idea could be to move the call below before ssh_wzv. However, "now" scale factors 
    120       ! at U-V points (which are set thanks to sshu_n, sshv_n) are actually available in sshwzv. 
    121       ! These are needed for log bottom friction... 
    122 #if ! defined key_dynspg_ts 
    123117                         CALL zdf_bfr( kstp )         ! bottom friction 
    124 #endif 
    125       ! end jchanut tschanges 
    126118 
    127119      !                                               ! Vertical eddy viscosity and diffusivity coefficients 
     
    224216            &                                          rhd, gru , grv  )      ! of t, s, rd at the last ocean level 
    225217 
    226       ELSE    
    227                                                ! centered hpg  (eos then time stepping) 
    228       ! bg jchanut tschanges 
    229 #if ! defined key_dynspg_ts 
    230       ! eos already called 
     218      ELSE                                                  ! centered hpg  (eos then time stepping) 
    231219                             CALL eos    ( tsn, rhd, rhop )      ! now in situ density for hpg computation 
    232220         IF( ln_zps      )   CALL zps_hde( kstp, jpts, tsn, gtsu, gtsv,  &    ! zps: now hor. derivative 
    233221            &                                          rhd, gru , grv  )      ! of t, s, rd at the last ocean level 
    234 #endif 
    235       ! end jchanut tschanges 
    236222         IF( ln_zdfnpc   )   CALL tra_npc( kstp )                ! update after fields by non-penetrative convection 
    237223                             CALL tra_nxt( kstp )                ! tracer fields at next time step 
     
    241227      ! Dynamics                                    (tsa used as workspace) 
    242228      !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 
    243       ! bg jchanut tschanges 
    244 #if defined key_dynspg_ts       
    245 ! revert to previously computed tendencies: 
    246 ! (not using ua, va as temporary arrays during tracers' update could avoid that) 
    247                                ua(:,:,:) = ua_bak(:,:,:)             
    248                                va(:,:,:) = va_bak(:,:,:) 
    249                                CALL dyn_bfr( kstp )         ! bottom friction 
    250                                CALL dyn_zdf( kstp )         ! vertical diffusion 
    251 #else 
    252       ! end jchanut tschanges 
    253229                               ua(:,:,:) = 0.e0             ! set dynamics trends to zero 
    254230                               va(:,:,:) = 0.e0 
     
    270246                               CALL dyn_zdf( kstp )         ! vertical diffusion 
    271247                               CALL dyn_spg( kstp, indic )  ! surface pressure gradient 
    272       ! bg jchanut tschanges 
    273 #endif 
    274       ! end jchanut tschanges 
    275248                               CALL dyn_nxt( kstp )         ! lateral velocity at next time step 
    276249 
Note: See TracChangeset for help on using the changeset viewer.