Changeset 6686


Ignore:
Timestamp:
2016-06-13T09:38:23+02:00 (4 years ago)
Author:
deazer
Message:

Added extra option to compute final harmonic analysis or not

Location:
branches/UKMO/dev_5518_tide_analysis_restart/NEMOGCM/NEMO/OPA_SRC
Files:
3 edited

Legend:

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

    r6370 r6686  
    5151   END TYPE TIDES_DATA 
    5252   INTEGER, PUBLIC, PARAMETER                  ::   jptides_max = 15      !: Max number of tidal contituents 
    53       LOGICAL, PUBLIC                           ::   ln_harm_ana             !: =T  Compute harmonic Analysis 
     53      LOGICAL, PUBLIC                           ::   ln_harm_ana_store    !: =T Stores data for  harmonic Analysis 
     54      LOGICAL, PUBLIC                           ::   ln_harm_ana_compute     !: =T  Compute harmonic Analysis 
    5455      LOGICAL, PUBLIC                           ::   ln_harmana_read         !: =T  Decide to do the analysis  
    5556                                                                             !from scratch or continue previous run 
     
    9495      TYPE(MAP_POINTER), DIMENSION(jpbgrd)      ::   ibmap_ptr           !: array of pointers to nbmap 
    9596      !! 
    96       NAMELIST/nambdy_tide/filtide, ln_bdytide_2ddta, ln_bdytide_conj, ln_harm_ana, ln_harmana_read 
     97      NAMELIST/nambdy_tide/filtide, ln_bdytide_2ddta, ln_bdytide_conj, ln_harm_ana_store, ln_harm_ana_compute, ln_harmana_read 
    9798      !!---------------------------------------------------------------------- 
    9899 
     
    129130            IF(lwp) WRITE(numout,*) '             assume complex conjugate   : ', ln_bdytide_conj 
    130131            IF(lwp) WRITE(numout,*) '             Number of tidal components to read: ', nb_harmo 
    131             IF(lwp) WRITE(numout,*) '             Use PCOMS harmonic ananalysis or not: ', ln_harm_ana 
     132            IF(lwp) WRITE(numout,*) '             Use PCOMS harmonic ananalysis or not: ', ln_harm_ana_store 
     133            IF(lwp) WRITE(numout,*) '             Compute Final  harmonic ananalysis or not: ', ln_harm_ana_compute 
    132134            IF(lwp) WRITE(numout,*) '             Read in previous days harmonic data or start afresh: ', ln_harmana_read 
    133135            IF(lwp) THEN  
  • branches/UKMO/dev_5518_tide_analysis_restart/NEMOGCM/NEMO/OPA_SRC/DIA/diaharmana.F90

    r6370 r6686  
    153153 
    154154 
    155  
     155      if(ln_harm_ana_compute) then 
     156               WRITE(numout,*) "Computing harmonics at last step" 
    156157! 
    157158        cosampu=0.0_wp 
     
    227228! 
    228229        CALL harm_ana_out     ! output analysis (last time step) 
     230     ELSE ! ln_harmana_compute  
     231         WRITE(numout,*) "Skipping Computing harmonics at last step" 
     232     ENDIF 
    229233        
    230234    ENDIF 
  • branches/UKMO/dev_5518_tide_analysis_restart/NEMOGCM/NEMO/OPA_SRC/step.F90

    r6370 r6686  
    347347                               CALL dia_wri_state( 'output.abort', kstp ) 
    348348      ENDIF 
    349       IF( ln_harm_ana   )   CALL harm_ana( kstp )        ! Harmonic analysis of tides  
     349      IF( ln_harm_ana_store   )   CALL harm_ana( kstp )        ! Harmonic analysis of tides  
    350350      IF( kstp == nit000   )   THEN 
    351351                 CALL iom_close( numror )     ! close input  ocean restart file 
Note: See TracChangeset for help on using the changeset viewer.