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 6989 for branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limrhg.F90 – NEMO

Ignore:
Timestamp:
2016-10-05T09:43:42+02:00 (7 years ago)
Author:
clem
Message:

use a namelist parameter to choose between the different advection schemes

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limrhg.F90

    r6853 r6989  
    114114      CHARACTER (len=50) ::   charout 
    115115 
     116      REAL(wp) ::   zrhoco                                                   ! rau0 * rn_cio 
    116117      REAL(wp) ::   zdtevp, z1_dtevp                                         ! time step for subcycling 
    117118      REAL(wp) ::   ecc2, z1_ecc2                                            ! square of yield ellipse eccenticity 
     
    203204      ! 1) define some variables and initialize arrays 
    204205      !------------------------------------------------------------------------------! 
     206      zrhoco = rau0 * rn_cio  
     207 
    205208      ! ecc2: square of yield ellipse eccenticrity 
    206209      ecc2    = rn_ecc * rn_ecc 
     
    357360                
    358361               ! delta at T points 
    359                zdelta = SQRT( zdiv2 + ( zdt2 + zds2 ) * usecc2 )   
     362               zdelta = SQRT( zdiv2 + ( zdt2 + zds2 ) * z1_ecc2 )   
    360363 
    361364               ! P/delta at T points 
     
    424427 
    425428                  ! tau_io/(v_oce - v_ice) 
    426                   zTauO = zaV(ji,jj) * rhoco * SQRT( ( v_ice (ji,jj) - v_oce (ji,jj) ) * ( v_ice (ji,jj) - v_oce (ji,jj) )  & 
    427                      &                             + ( u_iceV(ji,jj) - u_oceV(ji,jj) ) * ( u_iceV(ji,jj) - u_oceV(ji,jj) ) ) 
     429                  zTauO = zaV(ji,jj) * zrhoco * SQRT( ( v_ice (ji,jj) - v_oce (ji,jj) ) * ( v_ice (ji,jj) - v_oce (ji,jj) )  & 
     430                     &                              + ( u_iceV(ji,jj) - u_oceV(ji,jj) ) * ( u_iceV(ji,jj) - u_oceV(ji,jj) ) ) 
    428431 
    429432                  ! tau_bottom/v_ice 
     
    469472                                
    470473                  ! tau_io/(u_oce - u_ice) 
    471                   zTauO = zaU(ji,jj) * rhoco * SQRT( ( u_ice (ji,jj) - u_oce (ji,jj) ) * ( u_ice (ji,jj) - u_oce (ji,jj) )  & 
    472                      &                             + ( v_iceU(ji,jj) - v_oceU(ji,jj) ) * ( v_iceU(ji,jj) - v_oceU(ji,jj) ) ) 
     474                  zTauO = zaU(ji,jj) * zrhoco * SQRT( ( u_ice (ji,jj) - u_oce (ji,jj) ) * ( u_ice (ji,jj) - u_oce (ji,jj) )  & 
     475                     &                              + ( v_iceU(ji,jj) - v_oceU(ji,jj) ) * ( v_iceU(ji,jj) - v_oceU(ji,jj) ) ) 
    473476 
    474477                  ! tau_bottom/u_ice 
     
    515518 
    516519                  ! tau_io/(u_oce - u_ice) 
    517                   zTauO = zaU(ji,jj) * rhoco * SQRT( ( u_ice (ji,jj) - u_oce (ji,jj) ) * ( u_ice (ji,jj) - u_oce (ji,jj) )  & 
    518                      &                             + ( v_iceU(ji,jj) - v_oceU(ji,jj) ) * ( v_iceU(ji,jj) - v_oceU(ji,jj) ) ) 
     520                  zTauO = zaU(ji,jj) * zrhoco * SQRT( ( u_ice (ji,jj) - u_oce (ji,jj) ) * ( u_ice (ji,jj) - u_oce (ji,jj) )  & 
     521                     &                              + ( v_iceU(ji,jj) - v_oceU(ji,jj) ) * ( v_iceU(ji,jj) - v_oceU(ji,jj) ) ) 
    519522 
    520523                  ! tau_bottom/u_ice 
     
    559562          
    560563                  ! tau_io/(v_oce - v_ice) 
    561                   zTauO = zaV(ji,jj) * rhoco * SQRT( ( v_ice (ji,jj) - v_oce (ji,jj) ) * ( v_ice (ji,jj) - v_oce (ji,jj) )  & 
    562                      &                             + ( u_iceV(ji,jj) - u_oceV(ji,jj) ) * ( u_iceV(ji,jj) - u_oceV(ji,jj) ) ) 
     564                  zTauO = zaV(ji,jj) * zrhoco * SQRT( ( v_ice (ji,jj) - v_oce (ji,jj) ) * ( v_ice (ji,jj) - v_oce (ji,jj) )  & 
     565                     &                              + ( u_iceV(ji,jj) - u_oceV(ji,jj) ) * ( u_iceV(ji,jj) - u_oceV(ji,jj) ) ) 
    563566 
    564567                  ! tau_bottom/v_ice 
     
    651654             
    652655            ! delta at T points 
    653             zdelta         = SQRT( divu_i(ji,jj) * divu_i(ji,jj) + ( zdt2 + zds2 ) * usecc2 )   
     656            zdelta         = SQRT( divu_i(ji,jj) * divu_i(ji,jj) + ( zdt2 + zds2 ) * z1_ecc2 )   
    654657            rswitch        = 1._wp - MAX( 0._wp, SIGN( 1._wp, -zdelta ) ) ! 0 if delta=0 
    655658            delta_i(ji,jj) = zdelta + rn_creepl * rswitch 
Note: See TracChangeset for help on using the changeset viewer.