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 7158 for branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/OPA_SRC/SBC/fldread.F90 – NEMO

Ignore:
Timestamp:
2016-10-29T01:21:05+02:00 (7 years ago)
Author:
clem
Message:

debug branch

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/OPA_SRC/SBC/fldread.F90

    r6204 r7158  
    166166         isecsbc = nsec_year + nsec1jan000 + (kit+it_offset)*NINT( rdt/REAL(nn_baro,wp) ) 
    167167      ELSE                      ! middle of sbc time step 
    168          isecsbc = nsec_year + nsec1jan000 + NINT(0.5 * REAL(kn_fsbc - 1,wp) * rdttra(1)) + it_offset * NINT(rdttra(1)) 
     168         isecsbc = nsec_year + nsec1jan000 + NINT(0.5 * REAL(kn_fsbc - 1,wp) * rdt) + it_offset * NINT(rdt) 
    169169      ENDIF 
    170170      imf = SIZE( sd ) 
     
    193193               CALL fld_rec( kn_fsbc, sd(jf), kt_offset = it_offset, kit = kit )    ! update after record informations 
    194194 
    195                ! if kn_fsbc*rdttra is larger than nfreqh (which is kind of odd), 
     195               ! if kn_fsbc*rdt is larger than nfreqh (which is kind of odd), 
    196196               ! it is possible that the before value is no more the good one... we have to re-read it 
    197197               ! if before is not the last record of the file currently opened and after is the first record to be read 
     
    214214               IF( sd(jf)%ln_tint ) THEN 
    215215                   
    216                   ! if kn_fsbc*rdttra is larger than nfreqh (which is kind of odd), 
     216                  ! if kn_fsbc*rdt is larger than nfreqh (which is kind of odd), 
    217217                  ! it is possible that the before value is no more the good one... we have to re-read it 
    218218                  ! if before record is not just just before the after record... 
     
    245245                        ! year/month/week/day file to be not present. If the run continue further than the current 
    246246                        ! year/month/week/day, next year/month/week/day file must exist 
    247                         isecend = nsec_year + nsec1jan000 + (nitend - kt) * NINT(rdttra(1))   ! second at the end of the run  
     247                        isecend = nsec_year + nsec1jan000 + (nitend - kt) * NINT(rdt)   ! second at the end of the run  
    248248                        llstop = isecend > sd(jf)%nrec_a(2)                                   ! read more than 1 record of next year 
    249249                        ! we suppose that the date of next file is next day (should be ok even for weekly files...) 
     
    460460      IF( PRESENT(kt_offset) )   it_offset = kt_offset 
    461461      IF( PRESENT(kit) ) THEN   ;   it_offset = ( kit + it_offset ) * NINT( rdt/REAL(nn_baro,wp) ) 
    462       ELSE                      ;   it_offset =         it_offset   * NINT(       rdttra(1)      ) 
     462      ELSE                      ;   it_offset =         it_offset   * NINT(       rdt            ) 
    463463      ENDIF 
    464464      ! 
     
    537537         ELSE                                           ;   ztmp = REAL(nsec_year ,wp)  ! since 00h on Jan 1 of the current year 
    538538         ENDIF 
    539          ztmp = ztmp + 0.5 * REAL(kn_fsbc - 1, wp) * rdttra(1) + REAL( it_offset, wp )  ! centrered in the middle of sbc time step 
    540          ztmp = ztmp + 0.01 * rdttra(1)                                                 ! avoid truncation error  
     539         ztmp = ztmp + 0.5 * REAL(kn_fsbc - 1, wp) * rdt + REAL( it_offset, wp )  ! centrered in the middle of sbc time step 
     540         ztmp = ztmp + 0.01 * rdt                                                 ! avoid truncation error  
    541541         IF( sdjf%ln_tint ) THEN                ! time interpolation, shift by 1/2 record 
    542542            ! 
Note: See TracChangeset for help on using the changeset viewer.