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 1291 for trunk/NEMO/OFF_SRC/dtadyn.F90 – NEMO

Ignore:
Timestamp:
2009-02-03T15:22:22+01:00 (15 years ago)
Author:
cetlod
Message:

update modules, see ticket:320

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/OFF_SRC/dtadyn.F90

    r1181 r1291  
    206206         flx  ! auxiliary field for 2-D surface boundary conditions 
    207207 
    208  
    209208      ! 0. Initialization 
    210209      ! ----------------- 
    211210 
    212211      IF (lfirdyn) THEN 
    213       ! 
    214       ! time step MUST BE nint000 
    215       ! 
    216           IF (kt.ne.nit000) THEN 
    217               IF (lwp) THEN 
    218                   WRITE (numout,*) ' kt MUST BE EQUAL to nit000. kt=',kt  & 
    219                      ,' nit000=',nit000 
    220               END IF 
     212         ! first time step MUST BE nit000 
     213         IF( kt /= nit000 ) THEN 
     214            IF (lwp) THEN  
     215               WRITE (numout,*) ' kt MUST BE EQUAL to nit000. kt=',kt ,' nit000=',nit000  
    221216              STOP 'dtadyn' 
    222           END if 
    223       ! Initialize the parameters of the interpolation 
    224       CALL dta_dyn_init 
     217            ENDIF 
     218          ENDIF  
     219          ! Initialize the parameters of the interpolation 
     220          CALL dta_dyn_init 
    225221      ENDIF 
    226222 
     
    234230      zweighm1 = 1. - zweigh 
    235231 
    236       IF (lperdyn) THEN 
     232      IF( lperdyn ) THEN 
    237233         iperm1 = MOD(INT(zt),ndtadyn) 
    238234      ELSE 
     
    240236      ENDIF 
    241237      iper = iperm1 + 1 
    242       IF (iperm1 == 0) THEN 
    243           IF (lperdyn) THEN 
     238      IF( iperm1 == 0 ) THEN 
     239          IF( lperdyn ) THEN 
    244240              iperm1 = ndtadyn 
    245241          ELSE  
    246               IF (lfirdyn) THEN 
    247                   IF (lwp) THEN  
    248                       WRITE (numout,*) ' dynamic file is not periodic ' 
    249                       WRITE (numout,*) ' with or without interpolation, ' 
    250                       WRITE (numout,*) ' we take the first value' 
    251                       WRITE (numout,*) ' for the previous period ' 
    252                       WRITE (numout,*) ' iperm1 = 0  ' 
    253                   END IF  
     242              IF( lfirdyn ) THEN  
     243                  IF (lwp) WRITE (numout,*) &  
     244                      &   ' dynamic file is not periodic with or without interpolation  & 
     245                      &   we take the first value for the previous period iperm1 = 0  ' 
    254246              END IF 
    255247          END IF  
Note: See TracChangeset for help on using the changeset viewer.