 Timestamp:
 20191023T14:02:22+02:00 (5 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

NEMO/branches/2019/dev_r10742_ENHANCE12_SimonMTides/src/OCE/TDE/tide_mod.F90
r11768 r11770 570 570 INTEGER, INTENT(in) :: kformula 571 571 ! 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 574 575 !! 575 576 ! … … 583 584 zf = zf * zf 584 585 ! 585 CASE ( 2 ) !== formule 2, compound waves (78 x 0) === (78)586 zf1= nodal_factort(78)587 zf = nodal_factort( 0)588 zf = zf1 * zf589 !590 586 CASE ( 4 ) !== formule 4, compound waves (78 x 235) 591 587 zf1 = nodal_factort( 78) … … 593 589 zf = zf1 * zf 594 590 ! 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 * zf1599 !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 * zf1604 !605 CASE( 7 ) !== formule 7, compound waves (75 x 75)606 zf = nodal_factort(75)607 zf = zf * zf608 !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 * zf2614 !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 * zf2620 !621 CASE( 10 ) !== formule 10, compound waves (78 x 227)622 zf = nodal_factort( 78)623 zf1 = nodal_factort(227)624 zf = zf * zf1625 !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 * zf1631 !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 * zf1636 !637 CASE( 13 ) !== formule 13, compound waves (78 x 75)638 zf1 = nodal_factort(78)639 zf = nodal_factort(75)640 zf = zf * zf1641 !642 CASE( 14 ) !== formule 14, compound waves (235 x 0) === (235)643 zf = nodal_factort(235)644 zf1 = nodal_factort( 0)645 zf = zf * zf1646 !647 CASE( 15 ) !== formule 15, compound waves (235 x 75)648 zf = nodal_factort(235)649 zf1 = nodal_factort( 75)650 zf = zf * zf1651 !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 * zf1656 !657 CASE( 17 ) !== formule 17, compound waves (227 x 0)658 zf1 = nodal_factort(227)659 zf = nodal_factort( 0)660 zf = zf * zf1661 !662 591 CASE( 18 ) !== formule 18, compound waves (78 x 78 x 78 ) 663 592 zf1 = nodal_factort(78) 664 593 zf = zf1 * zf1 * zf1 665 594 ! 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 missing668 zf = nodal_factort(78)669 zf1 = nodal_factort( 0)670 zf = zf * zf1 * zf1671 !672 595 CASE( 20 ) !== formula 20, compound waves ( 78 x 78 x 78 x 78 ) 673 596 zf1 = nodal_factort(78) … … 689 612 zf = sin(2*sh_I) / 0.7214 690 613 ! 691 CASE( 77 ) !== formule 77692 zs = sin(sh_I/2)693 zf = sin(sh_I) * zs * zs / 0.0164694 !695 614 CASE( 78 ) !== formule 78 696 615 zs = cos(sh_I/2) 697 616 zf = zs * zs * zs * zs / 0.9154 698 617 ! 699 CASE( 79 ) !== formule 79700 zs = sin(sh_I)701 zf = zs * zs / 0.1565702 !703 CASE( 144 ) !== formule 144704 zs = sin(sh_I/2)705 zf = ( 110*zs*zs+15*zs*zs*zs*zs ) * cos(sh_I/2) / 0.5873706 !707 618 CASE( 149 ) !== formule 149 708 619 zs = cos(sh_I/2) … … 721 632 zf = sqrt( 19.0444*zs*zs*zs*zs + 2.7702*zs*zs*cos(2*sh_nu) + .0981 ) 722 633 ! 634 CASE DEFAULT 635 WRITE( clformula, '(I3)' ) kformula 636 CALL ctl_stop('nodal_factort: formula ' // clformula // ' is not available') 723 637 END SELECT 724 638 !
Note: See TracChangeset
for help on using the changeset viewer.