Changeset 12047


Ignore:
Timestamp:
2019-12-04T12:49:03+01:00 (8 weeks ago)
Author:
smueller
Message:

Replacement of local constants in subroutine tide_pulse of module tide_mod by existing module variables (ticket #2194)
This modification slightly changes many tidal-constituent frequencies computed by subroutine tide_pulse, with the exception of the frequencies of the S1, S2, and S4 tidal constituents.

File:
1 edited

Legend:

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

    r12045 r12047  
    529529      INTEGER  ::   jh 
    530530      REAL(wp) ::   zscale 
    531       REAL(wp) ::   zomega_T =  13149000.0_wp 
    532       REAL(wp) ::   zomega_s =    481267.892_wp 
    533       REAL(wp) ::   zomega_h =     36000.76892_wp 
    534       REAL(wp) ::   zomega_p =      4069.0322056_wp 
    535       REAL(wp) ::   zomega_n =      1934.1423972_wp 
    536       REAL(wp) ::   zomega_p1=         1.719175_wp 
    537       !!---------------------------------------------------------------------- 
    538       ! 
    539       zscale =  rad / ( 36525._wp * 86400._wp )  
     531      !!---------------------------------------------------------------------- 
     532      ! 
     533      zscale =  rad / 3600.0_wp 
    540534      ! 
    541535      DO jh = 1, size(ptide_harmo) 
    542          ptide_harmo(jh)%omega = (  zomega_T * ptide_comp( jh )%nT   & 
    543             &                         + zomega_s * ptide_comp( jh )%ns   & 
    544             &                         + zomega_h * ptide_comp( jh )%nh   & 
    545             &                         + zomega_p * ptide_comp( jh )%np   & 
    546             &                         + zomega_p1* ptide_comp( jh )%np1  ) * zscale 
     536         ptide_harmo(jh)%omega = (  romega_T * ptide_comp( jh )%nT   & 
     537            &                     + romega_s * ptide_comp( jh )%ns   & 
     538            &                     + romega_h * ptide_comp( jh )%nh   & 
     539            &                     + romega_p * ptide_comp( jh )%np   & 
     540            &                     + romega_p1* ptide_comp( jh )%np1  ) * zscale 
    547541      END DO 
    548542      ! 
Note: See TracChangeset for help on using the changeset viewer.