Changeset 6988


Ignore:
Timestamp:
2016-10-04T18:45:23+02:00 (4 years ago)
Author:
kingr
Message:

Added Pete/Enda?'s tidal diagnostic changes.

Location:
branches/UKMO/CO6_KD490_amm7oper/NEMOGCM/NEMO/OPA_SRC
Files:
1 added
3 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/CO6_KD490_amm7oper/NEMOGCM/NEMO/OPA_SRC/BDY/bdytides.F90

    r6332 r6988  
    5050      REAL(wp), POINTER, DIMENSION(:,:,:)    ::   v          !: Tidal constituents : V    (after nodal cor.) 
    5151   END TYPE TIDES_DATA 
     52 
     53   INTEGER, PUBLIC, PARAMETER                  ::   jptides_max = 15      !: Max number of tidal contituents 
     54      LOGICAL, PUBLIC                           ::   ln_harm_ana_store    !: =T Stores data for  harmonic Analysis 
     55      LOGICAL, PUBLIC                           ::   ln_harm_ana_compute     !: =T  Compute harmonic Analysis 
     56      LOGICAL, PUBLIC                           ::   ln_harmana_read         !: =T  Decide to do the analysis 
    5257 
    5358!$AGRIF_DO_NOT_TREAT 
     
    9095      TYPE(MAP_POINTER), DIMENSION(jpbgrd)      ::   ibmap_ptr           !: array of pointers to nbmap 
    9196      !! 
    92       NAMELIST/nambdy_tide/filtide, ln_bdytide_2ddta, ln_bdytide_conj 
     97      NAMELIST/nambdy_tide/filtide, ln_bdytide_2ddta, ln_bdytide_conj,ln_harm_ana_store, ln_harm_ana_compute, ln_harmana_read 
    9398      !!---------------------------------------------------------------------- 
    9499 
     
    126131            IF(lwp) WRITE(numout,*) '             assume complex conjugate   : ', ln_bdytide_conj 
    127132            IF(lwp) WRITE(numout,*) '             Number of tidal components to read: ', nb_harmo 
     133            IF(lwp) WRITE(numout,*) '             Use PCOMS harmonic ananalysis or not: ', ln_harm_ana_store 
     134            IF(lwp) WRITE(numout,*) '             Compute Final  harmonic ananalysis or not: ', ln_harm_ana_compute 
     135            IF(lwp) WRITE(numout,*) '             Read in previous days harmonic data or start afresh: ', ln_harmana_read 
    128136            IF(lwp) THEN  
    129137                    WRITE(numout,*) '             Tidal components: '  
  • branches/UKMO/CO6_KD490_amm7oper/NEMOGCM/NEMO/OPA_SRC/step.F90

    r6610 r6988  
    338338      IF( lk_vvl           )   CALL dom_vvl_sf_swp( kstp )  ! swap of vertical scale factors 
    339339      ! 
    340       IF( lrst_oce         )   CALL rst_write( kstp )       ! write output ocean restart file 
    341340 
    342341#if defined key_agrif 
     
    362361                               CALL dia_wri_state( 'output.abort', kstp ) 
    363362      ENDIF 
     363      IF( ln_harm_ana_store   )   CALL harm_ana( kstp )        ! Harmonic analysis of tides  
    364364      IF( kstp == nit000   )   THEN 
    365365                 CALL iom_close( numror )     ! close input  ocean restart file 
     
    367367         IF( lwm.AND.numoni /= -1 ) CALL FLUSH    ( numoni )     ! flush output namelist ice 
    368368      ENDIF 
     369      IF( lrst_oce         )   CALL rst_write( kstp )       ! write output ocean restart file 
    369370 
    370371      !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
  • branches/UKMO/CO6_KD490_amm7oper/NEMOGCM/NEMO/OPA_SRC/step_oce.F90

    r6610 r6988  
    108108   USE prtctl           ! Print control                    (prt_ctl routine) 
    109109 
     110   USE harmonic_analysis ! harmonic analysis of tides (harm_ana routine) 
     111   USE bdytides          ! harmonic analysis of tides (harm_ana routine)  
    110112   USE diaobs           ! Observation operator 
    111113 
Note: See TracChangeset for help on using the changeset viewer.