Ignore:
Timestamp:
2021-02-26T11:03:18+01:00 (5 months ago)
Author:
techene
Message:

#2605 small bug correction in ln_linssh=T

Location:
NEMO/branches/2021/dev_r14318_RK3_stage1/src/OCE
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2021/dev_r14318_RK3_stage1/src/OCE/DYN/dynspg_ts.F90

    r14547 r14549  
    247247 
    248248!!gm ==>>> !!ts    ISSUe her on en discute    changer les cas ENS ENE  et triad ? 
    249       IF( kt == nit000 .AND. .NOT. ln_linssh )   CALL dyn_cor_2D_init( Kmm )   ! Set zwz, the barotropic Coriolis force coefficient 
     249      IF( kt == nit000 .OR. .NOT. ln_linssh )   CALL dyn_cor_2D_init( Kmm )   ! Set zwz, the barotropic Coriolis force coefficient 
    250250      !                      ! recompute zwz = f/depth  at every time step for (.NOT.ln_linssh) as the water colomn height changes 
    251251      ! 
     
    287287      zhV(:,:) = pvv_b(:,:,Kmm) * hv(:,:,Kmm) * e1v(:,:)        ! NB: FULL domain : put a value in last row and column 
    288288      ! 
    289       CALL dyn_cor_2d( ht(:,:), hu(:,:,Kmm), hv(:,:,Kmm), puu_b(:,:,Kmm), pvv_b(:,:,Kmm), zhU, zhV,  &   ! <<== in 
     289      CALL dyn_cor_2D( ht(:,:), hu(:,:,Kmm), hv(:,:,Kmm), puu_b(:,:,Kmm), pvv_b(:,:,Kmm), zhU, zhV,  &   ! <<== in 
    290290         &                                                                          zu_trd, zv_trd   )   ! ==>> out 
    291291      ! 
     
    612612         ! at each time step. We however keep them constant here for optimization. 
    613613         ! Recall that zhU and zhV hold fluxes at jn+0.5 (extrapolated not backward interpolated) 
    614          CALL dyn_cor_2d( zhtp2_e, zhup2_e, zhvp2_e, ua_e, va_e, zhU, zhV,    zu_trd, zv_trd   ) 
     614         CALL dyn_cor_2D( zhtp2_e, zhup2_e, zhvp2_e, ua_e, va_e, zhU, zhV,    zu_trd, zv_trd   ) 
    615615         ! 
    616616         ! Add tidal astronomical forcing if defined 
     
    12071207      END SELECT 
    12081208      ! 
    1209    END SUBROUTINE dyn_cor_2d_init 
    1210  
    1211  
    1212    SUBROUTINE dyn_cor_2d( pht, phu, phv, punb, pvnb, zhU, zhV,    zu_trd, zv_trd   ) 
     1209   END SUBROUTINE dyn_cor_2D_init 
     1210 
     1211 
     1212   SUBROUTINE dyn_cor_2D( pht, phu, phv, punb, pvnb, zhU, zhV,    zu_trd, zv_trd   ) 
    12131213      !!--------------------------------------------------------------------- 
    1214       !!                   ***  ROUTINE dyn_cor_2d  *** 
     1214      !!                   ***  ROUTINE dyn_cor_2D  *** 
    12151215      !! 
    12161216      !! ** Purpose : Compute u and v coriolis trends 
  • NEMO/branches/2021/dev_r14318_RK3_stage1/src/OCE/stprk3_stg.F90

    r14547 r14549  
    138138      !                     !==  ssh/h0 ratio at Kaa  ==!  
    139139      ! 
    140       CALL dom_qco_r3c( ssh(:,:,Kaa), r3t(:,:,Kaa), r3u(:,:,Kaa), r3v(:,:,Kaa), r3f(:,:) )   ! "after" ssh/h_0 ratio at t,u,v-column 
     140      IF( .NOT.lk_linssh )   CALL dom_qco_r3c( ssh(:,:,Kaa), r3t(:,:,Kaa), r3u(:,:,Kaa), r3v(:,:,Kaa), r3f(:,:) )   ! "after" ssh/h_0 ratio at t,u,v-column 
    141141      ! 
    142142      ! 
Note: See TracChangeset for help on using the changeset viewer.