Changeset 11770


Ignore:
Timestamp:
2019-10-23T14:02:22+02:00 (10 months ago)
Author:
smueller
Message:

Removal of unused node-factor formulae (ticket #2194)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r10742_ENHANCE-12_SimonM-Tides/src/OCE/TDE/tide_mod.F90

    r11768 r11770  
    570570      INTEGER, INTENT(in) :: kformula 
    571571      ! 
    572       REAL(wp) :: zf 
    573       REAL(wp) :: zs, zf1, zf2 
     572      REAL(wp)            :: zf 
     573      REAL(wp)            :: zs, zf1, zf2 
     574      CHARACTER(LEN=3)    :: clformula 
    574575      !!---------------------------------------------------------------------- 
    575576      ! 
     
    583584         zf = zf * zf 
    584585         ! 
    585       CASE ( 2 )                 !==  formule 2, compound waves (78 x 0)  ===  (78)  
    586        zf1= nodal_factort(78) 
    587        zf = nodal_factort( 0) 
    588        zf = zf1 * zf 
    589        ! 
    590586      CASE ( 4 )                 !==  formule 4,  compound waves (78 x 235)  
    591587         zf1 = nodal_factort( 78) 
     
    593589         zf  = zf1 * zf 
    594590         ! 
    595       CASE ( 5 )                 !==  formule 5,  compound waves (78 *78 x 235) 
    596          zf1 = nodal_factort( 78) 
    597          zf  = nodal_factort(235) 
    598          zf  = zf * zf1 * zf1 
    599          ! 
    600       CASE ( 6 )                 !==  formule 6,  compound waves (78 *78 x 0) 
    601          zf1 = nodal_factort(78) 
    602          zf  = nodal_factort( 0) 
    603          zf  = zf * zf1 * zf1  
    604          ! 
    605       CASE( 7 )                  !==  formule 7, compound waves (75 x 75) 
    606          zf = nodal_factort(75) 
    607          zf = zf * zf 
    608          ! 
    609       CASE( 8 )                  !==  formule 8,  compound waves (78 x 0 x 235) 
    610          zf  = nodal_factort( 78) 
    611          zf1 = nodal_factort(  0) 
    612          zf2 = nodal_factort(235) 
    613          zf  = zf * zf1 * zf2 
    614          ! 
    615       CASE( 9 )                  !==  formule 9,  compound waves (78 x 0 x 227) 
    616          zf  = nodal_factort( 78) 
    617          zf1 = nodal_factort(  0) 
    618          zf2 = nodal_factort(227) 
    619          zf  = zf * zf1 * zf2 
    620          ! 
    621       CASE( 10 )                 !==  formule 10,  compound waves (78 x 227) 
    622          zf  = nodal_factort( 78) 
    623          zf1 = nodal_factort(227) 
    624          zf  = zf * zf1 
    625          ! 
    626       CASE( 11 )                 !==  formule 11,  compound waves (75 x 0) 
    627 !!gm bug???? zf 2 fois ! 
    628          zf = nodal_factort(75) 
    629          zf1 = nodal_factort( 0) 
    630          zf = zf * zf1 
    631          ! 
    632       CASE( 12 )                 !==  formule 12,  compound waves (78 x 78 x 78 x 0)  
    633          zf1 = nodal_factort(78) 
    634          zf  = nodal_factort( 0) 
    635          zf  = zf * zf1 * zf1 * zf1 
    636          ! 
    637       CASE( 13 )                 !==  formule 13, compound waves (78 x 75) 
    638          zf1 = nodal_factort(78) 
    639          zf  = nodal_factort(75) 
    640          zf  = zf * zf1 
    641          ! 
    642       CASE( 14 )                 !==  formule 14, compound waves (235 x 0)  ===  (235) 
    643          zf  = nodal_factort(235) 
    644          zf1 = nodal_factort(  0) 
    645          zf  = zf * zf1 
    646          ! 
    647       CASE( 15 )                 !==  formule 15, compound waves (235 x 75)  
    648          zf  = nodal_factort(235) 
    649          zf1 = nodal_factort( 75) 
    650          zf  = zf * zf1 
    651          ! 
    652       CASE( 16 )                 !==  formule 16, compound waves (78 x 0 x 0)  ===  (78) 
    653          zf  = nodal_factort(78) 
    654          zf1 = nodal_factort( 0) 
    655          zf  = zf * zf1 * zf1 
    656          ! 
    657       CASE( 17 )                 !==  formule 17,  compound waves (227 x 0)  
    658          zf1 = nodal_factort(227) 
    659          zf  = nodal_factort(  0) 
    660          zf  = zf * zf1 
    661          ! 
    662591      CASE( 18 )                 !==  formule 18,  compound waves (78 x 78 x 78 ) 
    663592         zf1 = nodal_factort(78) 
    664593         zf  = zf1 * zf1 * zf1 
    665594         ! 
    666       CASE( 19 )                 !==  formule 19, compound waves (78 x 0 x 0 x 0)  ===  (78) 
    667 !!gm bug2 ==>>>   here identical to formule 16,  a third multiplication by zf1 is missing 
    668          zf  = nodal_factort(78) 
    669          zf1 = nodal_factort( 0) 
    670          zf = zf * zf1 * zf1 
    671          ! 
    672595      CASE( 20 )                 !==  formula 20, compound waves ( 78 x 78 x 78 x 78 ) 
    673596         zf1 = nodal_factort(78) 
     
    689612         zf = sin(2*sh_I) / 0.7214 
    690613         ! 
    691       CASE( 77 )                 !==  formule 77 
    692          zs = sin(sh_I/2) 
    693          zf = sin(sh_I) * zs * zs / 0.0164 
    694          ! 
    695614      CASE( 78 )                 !==  formule 78 
    696615         zs = cos(sh_I/2) 
    697616         zf = zs * zs * zs * zs / 0.9154 
    698617         ! 
    699       CASE( 79 )                 !==  formule 79 
    700          zs = sin(sh_I) 
    701          zf = zs * zs / 0.1565 
    702          ! 
    703       CASE( 144 )                !==  formule 144 
    704          zs = sin(sh_I/2) 
    705          zf = ( 1-10*zs*zs+15*zs*zs*zs*zs ) * cos(sh_I/2) / 0.5873 
    706          ! 
    707618      CASE( 149 )                !==  formule 149 
    708619         zs = cos(sh_I/2) 
     
    721632         zf = sqrt( 19.0444*zs*zs*zs*zs + 2.7702*zs*zs*cos(2*sh_nu) + .0981 ) 
    722633         ! 
     634      CASE DEFAULT 
     635         WRITE( clformula, '(I3)' ) kformula 
     636         CALL ctl_stop('nodal_factort: formula ' // clformula // ' is not available') 
    723637      END SELECT 
    724638      ! 
Note: See TracChangeset for help on using the changeset viewer.