New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 15547 for NEMO/branches/UKMO/NEMO_4.0.4_CO9_shelf_climate/src/OCE/DIA/diaharm_fast.F90 – NEMO

Ignore:
Timestamp:
2021-11-27T11:04:33+01:00 (9 months ago)
Author:
hadjt
Message:

Moving ftide, utide, v0tide output to sbctides, adding tide_t output to diaharm_fast

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/UKMO/NEMO_4.0.4_CO9_shelf_climate/src/OCE/DIA/diaharm_fast.F90

    r15502 r15547  
    77 
    88   !!---------------------------------------------------------------------- 
    9    !!   'key_harm_ana'  :                Calculate harmonic analysis 
     9   !!                   :                Calculate harmonic analysis 
    1010   !!---------------------------------------------------------------------- 
    1111   !!   harm_ana        : 
     
    162162 
    163163         
    164     DO jh = 1, nb_ana        
    165  
    166         tmp_name=TRIM(Wave(ntide_all(jh))%cname_tide)//'_anau_utide' 
    167         IF( iom_use(TRIM(tmp_name)) )  THEN 
    168         !    IF(lwp) WRITE(numout,*) "harm_ana_out: iom_put: ",TRIM(tmp_name),'; shape = ', SHAPE(anau(jh) ) 
    169             CALL iom_put( TRIM(tmp_name), anau(jh) ) 
    170         !ELSE 
    171         !    IF(lwp) WRITE(numout,*) "harm_ana_out: not requested: ",TRIM(tmp_name) 
    172         ENDIF     
    173  
    174         tmp_name=TRIM(Wave(ntide_all(jh))%cname_tide)//'_anav_v0tide' 
    175         IF( iom_use(TRIM(tmp_name)) )  THEN 
    176         !    IF(lwp) WRITE(numout,*) "harm_ana_out: iom_put: ",TRIM(tmp_name),'; shape = ', SHAPE(anav(jh) ) 
    177             CALL iom_put( TRIM(tmp_name), anav(jh) ) 
    178         !ELSE 
    179         !    IF(lwp) WRITE(numout,*) "harm_ana_out: not requested: ",TRIM(tmp_name) 
    180         ENDIF 
    181  
    182         tmp_name=TRIM(Wave(ntide_all(jh))%cname_tide)//'_anaf_ftide' 
    183         IF( iom_use(TRIM(tmp_name)) )  THEN 
    184         !    IF(lwp) WRITE(numout,*) "harm_ana_out: iom_put: ",TRIM(tmp_name),'; shape = ', SHAPE(anaf(jh) ) 
    185             CALL iom_put( TRIM(tmp_name), anaf(jh) ) 
    186         !ELSE 
    187         !    IF(lwp) WRITE(numout,*) "harm_ana_out: not requested: ",TRIM(tmp_name) 
    188         ENDIF 
    189  
    190      END DO 
    191         
     164!    DO jh = 1, nb_ana        
     165 
     166!        tmp_name=TRIM(Wave(ntide_all(jh))%cname_tide)//'_anau_utide' 
     167!        IF( iom_use(TRIM(tmp_name)) )  THEN 
     168!        !    IF(lwp) WRITE(numout,*) "harm_ana_out: iom_put: ",TRIM(tmp_name),'; shape = ', SHAPE(anau(jh) ) 
     169!            CALL iom_put( TRIM(tmp_name), anau(jh) ) 
     170!        !ELSE 
     171!        !    IF(lwp) WRITE(numout,*) "harm_ana_out: not requested: ",TRIM(tmp_name) 
     172!        ENDIF     
     173 
     174!        tmp_name=TRIM(Wave(ntide_all(jh))%cname_tide)//'_anav_v0tide' 
     175!        IF( iom_use(TRIM(tmp_name)) )  THEN 
     176!        !    IF(lwp) WRITE(numout,*) "harm_ana_out: iom_put: ",TRIM(tmp_name),'; shape = ', SHAPE(anav(jh) ) 
     177!            CALL iom_put( TRIM(tmp_name), anav(jh) ) 
     178!        !ELSE 
     179!        !    IF(lwp) WRITE(numout,*) "harm_ana_out: not requested: ",TRIM(tmp_name) 
     180!        ENDIF 
     181 
     182!        tmp_name=TRIM(Wave(ntide_all(jh))%cname_tide)//'_anaf_ftide' 
     183!        IF( iom_use(TRIM(tmp_name)) )  THEN 
     184!        !    IF(lwp) WRITE(numout,*) "harm_ana_out: iom_put: ",TRIM(tmp_name),'; shape = ', SHAPE(anaf(jh) ) 
     185!            CALL iom_put( TRIM(tmp_name), anaf(jh) ) 
     186!        !ELSE 
     187!        !    IF(lwp) WRITE(numout,*) "harm_ana_out: not requested: ",TRIM(tmp_name) 
     188!        ENDIF 
     189 
     190!     END DO 
     191!        
    192192 
    193193     IF ( ln_diaharm_fast .and. ln_diaharm_store .and. ( lk_diaharm_2D .or. lk_diaharm_3D) ) THEN 
     
    199199          sec2start_old = nint( (fjulday-fjulday_startharm)*86400._wp )  
    200200          sec2start = nsec_day - NINT(0.5_wp * rdt) 
     201 
     202 
     203 
     204 
     205          IF( iom_use('tide_t') ) CALL iom_put( 'tide_t', sec2start ) 
     206 
     207 
     208 
     209 
    201210          !IF(lwp) WRITE(numout,*) "ztime NEW", kt, sec2start, fjulday_startharm 
    202211 
Note: See TracChangeset for help on using the changeset viewer.