- Timestamp:
- 2016-08-25T04:16:18+02:00 (8 years ago)
- Location:
- branches/UKMO/CO6_Restartable_Tidal_Analysis/NEMOGCM/NEMO/OPA_SRC
- Files:
-
- 1 added
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/CO6_Restartable_Tidal_Analysis/NEMOGCM/NEMO/OPA_SRC/BDY/bdytides.F90
r6822 r6885 50 50 REAL(wp), POINTER, DIMENSION(:,:,:) :: v !: Tidal constituents : V (after nodal cor.) 51 51 END TYPE TIDES_DATA 52 INTEGER, PUBLIC, PARAMETER :: jptides_max = 15 !: Max number of tidal contituents 53 LOGICAL, PUBLIC :: ln_harm_ana_store !: =T Stores data for harmonic Analysis 54 LOGICAL, PUBLIC :: ln_harm_ana_compute !: =T Compute harmonic Analysis 55 LOGICAL, PUBLIC :: ln_harmana_read !: =T Decide to do the analysis 56 !from scratch or continue previous run 52 57 53 58 !$AGRIF_DO_NOT_TREAT … … 90 95 TYPE(MAP_POINTER), DIMENSION(jpbgrd) :: ibmap_ptr !: array of pointers to nbmap 91 96 !! 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 93 98 !!---------------------------------------------------------------------- 94 99 … … 125 130 IF(lwp) WRITE(numout,*) ' assume complex conjugate : ', ln_bdytide_conj 126 131 IF(lwp) WRITE(numout,*) ' Number of tidal components to read: ', nb_harmo 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 134 IF(lwp) WRITE(numout,*) ' Read in previous days harmonic data or start afresh: ', ln_harmana_read 127 135 IF(lwp) THEN 128 136 WRITE(numout,*) ' Tidal components: ' -
branches/UKMO/CO6_Restartable_Tidal_Analysis/NEMOGCM/NEMO/OPA_SRC/step.F90
r6822 r6885 337 337 IF( lk_vvl ) CALL dom_vvl_sf_swp( kstp ) ! swap of vertical scale factors 338 338 ! 339 IF( lrst_oce ) CALL rst_write( kstp ) ! write output ocean restart file340 339 341 340 #if defined key_agrif … … 361 360 CALL dia_wri_state( 'output.abort', kstp ) 362 361 ENDIF 362 IF( ln_harm_ana_store ) CALL harm_ana( kstp ) ! Harmonic analysis of tides 363 363 IF( kstp == nit000 ) THEN 364 364 CALL iom_close( numror ) ! close input ocean restart file … … 366 366 IF( lwm.AND.numoni /= -1 ) CALL FLUSH ( numoni ) ! flush output namelist ice 367 367 ENDIF 368 IF( lrst_oce ) CALL rst_write( kstp ) ! write output ocean restart file 368 369 369 370 !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> -
branches/UKMO/CO6_Restartable_Tidal_Analysis/NEMOGCM/NEMO/OPA_SRC/step_oce.F90
r6822 r6885 106 106 USE prtctl ! Print control (prt_ctl routine) 107 107 108 USE harmonic_analysis ! harmonic analysis of tides (harm_ana routine) 109 USE bdytides ! harmonic analysis of tides (harm_ana routine) 108 110 USE diaobs ! Observation operator 109 111
Note: See TracChangeset
for help on using the changeset viewer.