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 2910 for branches/2011/dev_r2802_TOP_substepping/NEMOGCM/NEMO/OPA_SRC – NEMO

Ignore:
Timestamp:
2011-10-13T11:41:49+02:00 (13 years ago)
Author:
kpedwards
Message:

Updates from Christian - use kit000 in local TRA code; plus a few style corrections in TRC code.

Location:
branches/2011/dev_r2802_TOP_substepping/NEMOGCM/NEMO/OPA_SRC/TRA
Files:
18 edited

Legend:

Unmodified
Added
Removed

  • branches/2011/dev_r2802_TOP_substepping/NEMOGCM/NEMO/OPA_SRC/TRA/traadv.F90

    r2892 r2910  
    100100 
    101101      SELECT CASE ( nadv )                            !==  compute advection trend and add it to general trend  ==! 
    102       CASE ( 1 )   ;    CALL tra_adv_cen2  ( kt,nit000, 'TRA',         zun, zvn, zwn, tsb, tsn, tsa, jpts )   !  2nd order centered 
    103       CASE ( 2 )   ;    CALL tra_adv_tvd   ( kt,nit000, 'TRA', r2dtra, zun, zvn, zwn, tsb, tsn, tsa, jpts )   !  TVD  
    104       CASE ( 3 )   ;    CALL tra_adv_muscl ( kt,nit000, 'TRA', r2dtra, zun, zvn, zwn, tsb,      tsa, jpts )   !  MUSCL  
    105       CASE ( 4 )   ;    CALL tra_adv_muscl2( kt,nit000, 'TRA', r2dtra, zun, zvn, zwn, tsb, tsn, tsa, jpts )   !  MUSCL2  
    106       CASE ( 5 )   ;    CALL tra_adv_ubs   ( kt,nit000, 'TRA', r2dtra, zun, zvn, zwn, tsb, tsn, tsa, jpts )   !  UBS  
    107       CASE ( 6 )   ;    CALL tra_adv_qck   ( kt,nit000, 'TRA', r2dtra, zun, zvn, zwn, tsb, tsn, tsa, jpts )   !  QUICKEST  
     102      CASE ( 1 )   ;    CALL tra_adv_cen2  ( kt, nit000, 'TRA',         zun, zvn, zwn, tsb, tsn, tsa, jpts )   !  2nd order centered 
     103      CASE ( 2 )   ;    CALL tra_adv_tvd   ( kt, nit000, 'TRA', r2dtra, zun, zvn, zwn, tsb, tsn, tsa, jpts )   !  TVD  
     104      CASE ( 3 )   ;    CALL tra_adv_muscl ( kt, nit000, 'TRA', r2dtra, zun, zvn, zwn, tsb,      tsa, jpts )   !  MUSCL  
     105      CASE ( 4 )   ;    CALL tra_adv_muscl2( kt, nit000, 'TRA', r2dtra, zun, zvn, zwn, tsb, tsn, tsa, jpts )   !  MUSCL2  
     106      CASE ( 5 )   ;    CALL tra_adv_ubs   ( kt, nit000, 'TRA', r2dtra, zun, zvn, zwn, tsb, tsn, tsa, jpts )   !  UBS  
     107      CASE ( 6 )   ;    CALL tra_adv_qck   ( kt, nit000, 'TRA', r2dtra, zun, zvn, zwn, tsb, tsn, tsa, jpts )   !  QUICKEST  
    108108      ! 
    109109      CASE (-1 )                                      !==  esopa: test all possibility with control print  ==! 
    110          CALL tra_adv_cen2  ( kt,nit000, 'TRA',         zun, zvn, zwn, tsb, tsn, tsa, jpts )           
     110         CALL tra_adv_cen2  ( kt, nit000, 'TRA',         zun, zvn, zwn, tsb, tsn, tsa, jpts )           
    111111         CALL prt_ctl( tab3d_1=tsa(:,:,:,jp_tem), clinfo1=' adv0 - Ta: ', mask1=tmask,               & 
    112112            &          tab3d_2=tsa(:,:,:,jp_sal), clinfo2=       ' Sa: ', mask2=tmask, clinfo3='tra' ) 
    113          CALL tra_adv_tvd   ( kt,nit000, 'TRA', r2dtra, zun, zvn, zwn, tsb, tsn, tsa, jpts )           
     113         CALL tra_adv_tvd   ( kt, nit000, 'TRA', r2dtra, zun, zvn, zwn, tsb, tsn, tsa, jpts )           
    114114         CALL prt_ctl( tab3d_1=tsa(:,:,:,jp_tem), clinfo1=' adv1 - Ta: ', mask1=tmask,               & 
    115115            &          tab3d_2=tsa(:,:,:,jp_sal), clinfo2=       ' Sa: ', mask2=tmask, clinfo3='tra' ) 
    116          CALL tra_adv_muscl ( kt,nit000, 'TRA', r2dtra, zun, zvn, zwn, tsb,      tsa, jpts )           
     116         CALL tra_adv_muscl ( kt, nit000, 'TRA', r2dtra, zun, zvn, zwn, tsb,      tsa, jpts )           
    117117         CALL prt_ctl( tab3d_1=tsa(:,:,:,jp_tem), clinfo1=' adv3 - Ta: ', mask1=tmask,               & 
    118118            &          tab3d_2=tsa(:,:,:,jp_sal), clinfo2=       ' Sa: ', mask2=tmask, clinfo3='tra' ) 
    119          CALL tra_adv_muscl2( kt,nit000, 'TRA', r2dtra, zun, zvn, zwn, tsb, tsn, tsa, jpts )           
     119         CALL tra_adv_muscl2( kt, nit000, 'TRA', r2dtra, zun, zvn, zwn, tsb, tsn, tsa, jpts )           
    120120         CALL prt_ctl( tab3d_1=tsa(:,:,:,jp_tem), clinfo1=' adv4 - Ta: ', mask1=tmask,               & 
    121121            &          tab3d_2=tsa(:,:,:,jp_sal), clinfo2=       ' Sa: ', mask2=tmask, clinfo3='tra' ) 
    122          CALL tra_adv_ubs   ( kt,nit000, 'TRA', r2dtra, zun, zvn, zwn, tsb, tsn, tsa, jpts )           
     122         CALL tra_adv_ubs   ( kt, nit000, 'TRA', r2dtra, zun, zvn, zwn, tsb, tsn, tsa, jpts )           
    123123         CALL prt_ctl( tab3d_1=tsa(:,:,:,jp_tem), clinfo1=' adv5 - Ta: ', mask1=tmask,               & 
    124124            &          tab3d_2=tsa(:,:,:,jp_sal), clinfo2=       ' Sa: ', mask2=tmask, clinfo3='tra' ) 
    125          CALL tra_adv_qck   ( kt,nit000, 'TRA', r2dtra, zun, zvn, zwn, tsb, tsn, tsa, jpts )           
     125         CALL tra_adv_qck   ( kt, nit000, 'TRA', r2dtra, zun, zvn, zwn, tsb, tsn, tsa, jpts )           
    126126         CALL prt_ctl( tab3d_1=tsa(:,:,:,jp_tem), clinfo1=' adv6 - Ta: ', mask1=tmask,               & 
    127127            &          tab3d_2=tsa(:,:,:,jp_sal), clinfo2=       ' Sa: ', mask2=tmask, clinfo3='tra' ) 

  • branches/2011/dev_r2802_TOP_substepping/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_cen2.F90

    r2892 r2910  
    5353CONTAINS 
    5454 
    55    SUBROUTINE tra_adv_cen2( kt,nit000, cdtype, pun, pvn, pwn,     & 
     55   SUBROUTINE tra_adv_cen2( kt, kit000, cdtype, pun, pvn, pwn,     & 
    5656      &                                 ptb, ptn, pta, kjpt   )  
    5757      !!---------------------------------------------------------------------- 
    … …  
    116116      ! 
    117117      INTEGER                              , INTENT(in   ) ::   kt              ! ocean time-step index 
    118       INTEGER                              , INTENT(in   ) ::   nit000          ! first time step index 
     118      INTEGER                              , INTENT(in   ) ::   kit000          ! first time step index 
    119119      CHARACTER(len=3)                     , INTENT(in   ) ::   cdtype          ! =TRA or TRC (tracer indicator) 
    120120      INTEGER                              , INTENT(in   ) ::   kjpt            ! number of tracers 
    … …  
    136136      ENDIF 
    137137 
    138       IF( kt == nit000 )  THEN 
     138      IF( kt == kit000 )  THEN 
    139139         IF(lwp) WRITE(numout,*) 
    140140         IF(lwp) WRITE(numout,*) 'tra_adv_cen2 : 2nd order centered advection scheme on ', cdtype 

  • branches/2011/dev_r2802_TOP_substepping/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_eiv.F90

    r2892 r2910  
    4545CONTAINS 
    4646 
    47    SUBROUTINE tra_adv_eiv( kt,nit000, pun, pvn, pwn, cdtype ) 
     47   SUBROUTINE tra_adv_eiv( kt, kit000, pun, pvn, pwn, cdtype ) 
    4848      !!---------------------------------------------------------------------- 
    4949      !!                  ***  ROUTINE tra_adv_eiv  *** 
    … …  
    6969#endif 
    7070      INTEGER                         , INTENT(in   ) ::   kt       ! ocean time-step index 
    71       INTEGER                         , INTENT(in   ) ::   nit000   ! first time step index 
     71      INTEGER                         , INTENT(in   ) ::   kit000   ! first time step index 
    7272      CHARACTER(len=3)                , INTENT(in   ) ::   cdtype   ! =TRA or TRC (tracer indicator) 
    7373      REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(inout) ::   pun      ! in : 3 ocean velocity components  
    … …  
    9191      ENDIF 
    9292 
    93       IF( kt == nit000 )  THEN 
     93      IF( kt == kit000 )  THEN 
    9494         IF(lwp) WRITE(numout,*) 
    9595         IF(lwp) WRITE(numout,*) 'tra_adv_eiv : eddy induced advection on ', cdtype,' :' 
    … …  
    204204   !!---------------------------------------------------------------------- 
    205205CONTAINS 
    206    SUBROUTINE tra_adv_eiv( kt,nit000, pun, pvn, pwn, cdtype )              ! Empty routine 
     206   SUBROUTINE tra_adv_eiv( kt, kit000, pun, pvn, pwn, cdtype )              ! Empty routine 
    207207      INTEGER  ::   kt     
    208       INTEGER  ::   nit000     
     208      INTEGER  ::   kit000     
    209209      CHARACTER(len=3) ::   cdtype 
    210210      REAL, DIMENSION(:,:,:) ::   pun, pvn, pwn 

  • branches/2011/dev_r2802_TOP_substepping/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_muscl.F90

    r2892 r2910  
    4444CONTAINS 
    4545 
    46    SUBROUTINE tra_adv_muscl( kt,nit000,cdtype, p2dt, pun, pvn, pwn, & 
     46   SUBROUTINE tra_adv_muscl( kt, kit000, cdtype, p2dt, pun, pvn, pwn, & 
    4747      &                                        ptb, pta, kjpt ) 
    4848      !!---------------------------------------------------------------------- 
    … …  
    6666      ! 
    6767      INTEGER                              , INTENT(in   ) ::   kt              ! ocean time-step index 
    68       INTEGER                              , INTENT(in   ) ::   nit000          ! first time step index 
     68      INTEGER                              , INTENT(in   ) ::   kit000          ! first time step index 
    6969      CHARACTER(len=3)                     , INTENT(in   ) ::   cdtype          ! =TRA or TRC (tracer indicator) 
    7070      INTEGER                              , INTENT(in   ) ::   kjpt            ! number of tracers 
    … …  
    8484      ENDIF 
    8585 
    86       IF( kt == nit000 )  THEN 
     86      IF( kt == kit000 )  THEN 
    8787         IF(lwp) WRITE(numout,*) 
    8888         IF(lwp) WRITE(numout,*) 'tra_adv : MUSCL advection scheme on ', cdtype 

  • branches/2011/dev_r2802_TOP_substepping/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_muscl2.F90

    r2892 r2910  
    4242CONTAINS 
    4343 
    44    SUBROUTINE tra_adv_muscl2( kt,nit000, cdtype, p2dt, pun, pvn, pwn,      & 
     44   SUBROUTINE tra_adv_muscl2( kt, kit000, cdtype, p2dt, pun, pvn, pwn,      & 
    4545      &                                         ptb, ptn, pta, kjpt ) 
    4646      !!---------------------------------------------------------------------- 
    … …  
    6464      !! 
    6565      INTEGER                              , INTENT(in   ) ::   kt              ! ocean time-step index 
    66       INTEGER                              , INTENT(in   ) ::   nit000          ! first time step index 
     66      INTEGER                              , INTENT(in   ) ::   kit000          ! first time step index 
    6767      CHARACTER(len=3)                     , INTENT(in   ) ::   cdtype          ! =TRA or TRC (tracer indicator) 
    6868      INTEGER                              , INTENT(in   ) ::   kjpt            ! number of tracers 
    … …  
    8282      ENDIF 
    8383 
    84       IF( kt == nit000 )  THEN 
     84      IF( kt == kit000 )  THEN 
    8585         IF(lwp) WRITE(numout,*) 
    8686         IF(lwp) WRITE(numout,*) 'tra_adv_muscl2 : MUSCL2 advection scheme on ', cdtype 

  • branches/2011/dev_r2802_TOP_substepping/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_qck.F90

    r2892 r2910  
    4545CONTAINS 
    4646 
    47    SUBROUTINE tra_adv_qck ( kt, nit000, cdtype, p2dt, pun, pvn, pwn,      & 
     47   SUBROUTINE tra_adv_qck ( kt, kit000, cdtype, p2dt, pun, pvn, pwn,      & 
    4848      &                                       ptb, ptn, pta, kjpt ) 
    4949      !!---------------------------------------------------------------------- 
    … …  
    8282      !!---------------------------------------------------------------------- 
    8383      INTEGER                              , INTENT(in   ) ::   kt              ! ocean time-step index 
    84       INTEGER                              , INTENT(in   ) ::   nit000          ! first time step index 
     84      INTEGER                              , INTENT(in   ) ::   kit000          ! first time step index 
    8585      CHARACTER(len=3)                     , INTENT(in   ) ::   cdtype          ! =TRA or TRC (tracer indicator) 
    8686      INTEGER                              , INTENT(in   ) ::   kjpt            ! number of tracers 
    … …  
    9191      !!---------------------------------------------------------------------- 
    9292 
    93       IF( kt == nit000 )  THEN 
     93      IF( kt == kit000 )  THEN 
    9494         IF(lwp) WRITE(numout,*) 
    9595         IF(lwp) WRITE(numout,*) 'tra_adv_qck : 3rd order quickest advection scheme on ', cdtype 

  • branches/2011/dev_r2802_TOP_substepping/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_tvd.F90

    r2892 r2910  
    5151CONTAINS 
    5252 
    53    SUBROUTINE tra_adv_tvd ( kt,nit000, cdtype, p2dt, pun, pvn, pwn,      & 
     53   SUBROUTINE tra_adv_tvd ( kt, kit000, cdtype, p2dt, pun, pvn, pwn,      & 
    5454      &                                       ptb, ptn, pta, kjpt ) 
    5555      !!---------------------------------------------------------------------- 
    … …  
    7171      ! 
    7272      INTEGER                              , INTENT(in   ) ::   kt              ! ocean time-step index 
    73       INTEGER                              , INTENT(in   ) ::   nit000          ! first time step index 
     73      INTEGER                              , INTENT(in   ) ::   kit000          ! first time step index 
    7474      CHARACTER(len=3)                     , INTENT(in   ) ::   cdtype          ! =TRA or TRC (tracer indicator) 
    7575      INTEGER                              , INTENT(in   ) ::   kjpt            ! number of tracers 
    … …  
    9090      ENDIF 
    9191 
    92       IF( kt == nit000 )  THEN 
     92      IF( kt == kit000 )  THEN 
    9393         IF(lwp) WRITE(numout,*) 
    9494         IF(lwp) WRITE(numout,*) 'tra_adv_tvd : TVD advection scheme on ', cdtype 

  • branches/2011/dev_r2802_TOP_substepping/NEMOGCM/NEMO/OPA_SRC/TRA/traadv_ubs.F90

    r2892 r2910  
    4040CONTAINS 
    4141 
    42    SUBROUTINE tra_adv_ubs ( kt, nit000, cdtype, p2dt, pun, pvn, pwn,      & 
     42   SUBROUTINE tra_adv_ubs ( kt, kit000, cdtype, p2dt, pun, pvn, pwn,      & 
    4343      &                                       ptb, ptn, pta, kjpt ) 
    4444      !!---------------------------------------------------------------------- 
    … …  
    8080      ! 
    8181      INTEGER                              , INTENT(in   ) ::   kt              ! ocean time-step index 
    82       INTEGER                              , INTENT(in   ) ::   nit000          ! first time step index 
     82      INTEGER                              , INTENT(in   ) ::   kit000          ! first time step index 
    8383      CHARACTER(len=3)                     , INTENT(in   ) ::   cdtype          ! =TRA or TRC (tracer indicator) 
    8484      INTEGER                              , INTENT(in   ) ::   kjpt            ! number of tracers 
    … …  
    9898      ENDIF 
    9999 
    100       IF( kt == nit000 )  THEN 
     100      IF( kt == kit000 )  THEN 
    101101         IF(lwp) WRITE(numout,*) 
    102102         IF(lwp) WRITE(numout,*) 'tra_adv_ubs :  horizontal UBS advection scheme on ', cdtype 

  • branches/2011/dev_r2802_TOP_substepping/NEMOGCM/NEMO/OPA_SRC/TRA/trabbl.F90

    r2892 r2910  
    311311 
    312312 
    313    SUBROUTINE bbl( kt, nit000, cdtype ) 
     313   SUBROUTINE bbl( kt, kit000, cdtype ) 
    314314      !!---------------------------------------------------------------------- 
    315315      !!                  ***  ROUTINE bbl  *** 
    … …  
    343343      ! 
    344344      INTEGER         , INTENT(in   ) ::   kt       ! ocean time-step index 
    345       INTEGER         , INTENT(in   ) ::   nit000          ! first time step index 
     345      INTEGER         , INTENT(in   ) ::   kit000          ! first time step index 
    346346      CHARACTER(len=3), INTENT(in   ) ::   cdtype   ! =TRA or TRC (tracer indicator) 
    347347      !! 
    … …  
    390390      ENDIF 
    391391      
    392       IF( kt == nit000 )  THEN 
     392      IF( kt == kit000 )  THEN 
    393393         IF(lwp)  WRITE(numout,*) 
    394394         IF(lwp)  WRITE(numout,*) 'trabbl:bbl : Compute bbl velocities and diffusive coefficients in ', cdtype 
    … …  
    538538      !! 
    539539      !! ** Method  :   Read the nambbl namelist and check the parameters 
    540       !!              called by nemo_init at the first timestep (nit000) 
     540      !!              called by nemo_init at the first timestep (kit000) 
    541541      !!---------------------------------------------------------------------- 
    542542      USE wrk_nemo, ONLY:   wrk_in_use, wrk_not_released 

  • branches/2011/dev_r2802_TOP_substepping/NEMOGCM/NEMO/OPA_SRC/TRA/traldf.F90

    r2892 r2910  
    6868 
    6969      SELECT CASE ( nldf )                       ! compute lateral mixing trend and add it to the general trend 
    70       CASE ( 0 )   ;   CALL tra_ldf_lap     ( kt,nit000, 'TRA', gtsu, gtsv, tsb, tsa, jpts        )  ! iso-level laplacian 
     70      CASE ( 0 )   ;   CALL tra_ldf_lap     ( kt, nit000, 'TRA', gtsu, gtsv, tsb, tsa, jpts        )  ! iso-level laplacian 
    7171      CASE ( 1 )                                                                              ! rotated laplacian 
    7272         IF( ln_traldf_grif ) THEN                                                           
    … …  
    297297      ! Compute the ldf trends 
    298298      ! ---------------------- 
    299       CALL tra_ldf( nit000+1 )      ! horizontal components (+1: no more init) 
    300       CALL tra_zdf( nit000   )      ! vertical component (if necessary nit000 to performed the init) 
     299      CALL tra_ldf( nit000 + 1 )      ! horizontal components (+1: no more init) 
     300      CALL tra_zdf( nit000     )      ! vertical component (if necessary nit000 to performed the init) 
    301301 
    302302      ! finalise the computation and recover all arrays 

  • branches/2011/dev_r2802_TOP_substepping/NEMOGCM/NEMO/OPA_SRC/TRA/traldf_bilap.F90

    r2892 r2910  
    4747CONTAINS 
    4848  
    49    SUBROUTINE tra_ldf_bilap( kt, nit000, cdtype, pgu, pgv,      & 
     49   SUBROUTINE tra_ldf_bilap( kt, kit000, cdtype, pgu, pgv,      & 
    5050      &                                  ptb, pta, kjpt )   
    5151      !!---------------------------------------------------------------------- 
    … …  
    7979      !! 
    8080      INTEGER                              , INTENT(in   ) ::   kt         ! ocean time-step index 
    81       INTEGER                              , INTENT(in   ) ::   nit000          ! first time step index 
     81      INTEGER                              , INTENT(in   ) ::   kit000          ! first time step index 
    8282      CHARACTER(len=3)                     , INTENT(in   ) ::   cdtype     ! =TRA or TRC (tracer indicator) 
    8383      INTEGER                              , INTENT(in   ) ::   kjpt       ! number of tracers 
    … …  
    9494      ENDIF 
    9595 
    96       IF( kt == nit000 )  THEN 
     96      IF( kt == kit000 )  THEN 
    9797         IF(lwp) WRITE(numout,*) 
    9898         IF(lwp) WRITE(numout,*) 'tra_ldf_bilap : iso-level biharmonic operator on ', cdtype 

  • branches/2011/dev_r2802_TOP_substepping/NEMOGCM/NEMO/OPA_SRC/TRA/traldf_bilapg.F90

    r2892 r2910  
    4242CONTAINS 
    4343 
    44    SUBROUTINE tra_ldf_bilapg( kt, nit000, cdtype, ptb, pta, kjpt ) 
     44   SUBROUTINE tra_ldf_bilapg( kt, kit000, cdtype, ptb, pta, kjpt ) 
    4545      !!---------------------------------------------------------------------- 
    4646      !!                 ***  ROUTINE tra_ldf_bilapg  *** 
    … …  
    7070      ! 
    7171      INTEGER         , INTENT(in   )                      ::   kt       ! ocean time-step index 
    72       INTEGER         , INTENT(in   )                      ::   nit000   ! first time step index 
     72      INTEGER         , INTENT(in   )                      ::   kit000   ! first time step index 
    7373      CHARACTER(len=3), INTENT(in   )                      ::   cdtype   ! =TRA or TRC (tracer indicator) 
    7474      INTEGER         , INTENT(in   )                      ::   kjpt     ! number of tracers 
    … …  
    8383      ENDIF 
    8484 
    85       IF( kt == nit000 )  THEN 
     85      IF( kt == kit000 )  THEN 
    8686         IF(lwp) WRITE(numout,*) 
    8787         IF(lwp) WRITE(numout,*) 'tra_ldf_bilapg : horizontal biharmonic operator in s-coordinate on ', cdtype 
    … …  
    346346   !!---------------------------------------------------------------------- 
    347347CONTAINS 
    348    SUBROUTINE tra_ldf_bilapg( kt,nit000, cdtype, ptb, pta, kjpt )      ! Empty routine 
    349       INTEGER :: kt, nit000 
     348   SUBROUTINE tra_ldf_bilapg( kt, kit000, cdtype, ptb, pta, kjpt )      ! Empty routine 
     349      INTEGER :: kt, kit000 
    350350      CHARACTER(len=3) ::   cdtype 
    351351      REAL, DIMENSION(:,:,:,:) ::   ptb, pta 
    352       WRITE(*,*) 'tra_ldf_iso: You should not have seen this print! error?', kt, cdtype, ptb(1,1,1,1), pta(1,1,1,1), kjpt 
     352      WRITE(*,*) 'tra_ldf_iso: You should not have seen this print! error?', & 
     353        &         kt, kit000, cdtype, ptb(1,1,1,1), pta(1,1,1,1), kjpt 
    353354   END SUBROUTINE tra_ldf_bilapg 
    354355#endif 

  • branches/2011/dev_r2802_TOP_substepping/NEMOGCM/NEMO/OPA_SRC/TRA/traldf_iso.F90

    r2892 r2910  
    4949CONTAINS 
    5050 
    51    SUBROUTINE tra_ldf_iso( kt,nit000, cdtype, pgu, pgv,              & 
     51   SUBROUTINE tra_ldf_iso( kt, kit000, cdtype, pgu, pgv,              & 
    5252      &                                ptb, pta, kjpt, pahtb0 ) 
    5353      !!---------------------------------------------------------------------- 
    … …  
    9696      ! 
    9797      INTEGER                              , INTENT(in   ) ::   kt         ! ocean time-step index 
    98       INTEGER                              , INTENT(in   ) ::   nit000          ! first time step index 
     98      INTEGER                              , INTENT(in   ) ::   kit000          ! first time step index 
    9999      CHARACTER(len=3)                     , INTENT(in   ) ::   cdtype     ! =TRA or TRC (tracer indicator) 
    100100      INTEGER                              , INTENT(in   ) ::   kjpt       ! number of tracers 
    … …  
    117117      ENDIF 
    118118 
    119       IF( kt == nit000 )  THEN 
     119      IF( kt == kit000 )  THEN 
    120120         IF(lwp) WRITE(numout,*) 
    121121         IF(lwp) WRITE(numout,*) 'tra_ldf_iso : rotated laplacian diffusion operator on ', cdtype 
    … …  
    302302   !!---------------------------------------------------------------------- 
    303303CONTAINS 
    304    SUBROUTINE tra_ldf_iso( kt, nit000,cdtype, pgu, pgv, ptb, pta, kjpt, pahtb0 )      ! Empty routine 
    305       INTEGER:: kt, nit000 
     304   SUBROUTINE tra_ldf_iso( kt, kit000,cdtype, pgu, pgv, ptb, pta, kjpt, pahtb0 )      ! Empty routine 
     305      INTEGER:: kt, kit000 
    306306      CHARACTER(len=3) ::   cdtype 
    307307      REAL, DIMENSION(:,:,:) ::   pgu, pgv   ! tracer gradient at pstep levels 
    308308      REAL, DIMENSION(:,:,:,:) ::   ptb, pta 
    309       WRITE(*,*) 'tra_ldf_iso: You should not have seen this print! error?', kt, cdtype, pgu(1,1,1), pgv(1,1,1),   & 
    310          &                                                             ptb(1,1,1,1), pta(1,1,1,1), kjpt, pahtb0 
     309      WRITE(*,*) 'tra_ldf_iso: You should not have seen this print! error?', kt, kit000, cdtype,   & 
     310         &                       pgu(1,1,1), pgv(1,1,1), ptb(1,1,1,1), pta(1,1,1,1), kjpt, pahtb0 
    311311   END SUBROUTINE tra_ldf_iso 
    312312#endif 

  • branches/2011/dev_r2802_TOP_substepping/NEMOGCM/NEMO/OPA_SRC/TRA/traldf_iso_grif.F90

    r2892 r2910  
    4848CONTAINS 
    4949 
    50   SUBROUTINE tra_ldf_iso_grif( kt, nit000, cdtype, pgu, pgv,              & 
     50  SUBROUTINE tra_ldf_iso_grif( kt, kit000, cdtype, pgu, pgv,              & 
    5151       &                                   ptb, pta, kjpt, pahtb0 ) 
    5252      !!---------------------------------------------------------------------- 
    … …  
    9595      ! 
    9696      INTEGER                              , INTENT(in   ) ::   kt         ! ocean time-step index 
    97       INTEGER                              , INTENT(in   ) ::   nit000          ! first time step index 
     97      INTEGER                              , INTENT(in   ) ::   kit000          ! first time step index 
    9898      CHARACTER(len=3)                     , INTENT(in   ) ::   cdtype     ! =TRA or TRC (tracer indicator) 
    9999      INTEGER                              , INTENT(in   ) ::   kjpt       ! number of tracers 
    … …  
    128128      !zdkt(1:jpi,1:jpj,0:1) => wrk_3d_9(:,:,1:2) 
    129129 
    130       IF( kt == nit000 )  THEN 
     130      IF( kt == kit000 )  THEN 
    131131         IF(lwp) WRITE(numout,*) 
    132132         IF(lwp) WRITE(numout,*) 'tra_ldf_iso_grif : rotated laplacian diffusion operator on ', cdtype 

  • branches/2011/dev_r2802_TOP_substepping/NEMOGCM/NEMO/OPA_SRC/TRA/traldf_lap.F90

    r2892 r2910  
    4444CONTAINS 
    4545 
    46    SUBROUTINE tra_ldf_lap( kt, nit000, cdtype, pgu, pgv,      & 
     46   SUBROUTINE tra_ldf_lap( kt, kit000, cdtype, pgu, pgv,      & 
    4747      &                                ptb, pta, kjpt )  
    4848      !!---------------------------------------------------------------------- 
    … …  
    6666      ! 
    6767      INTEGER                              , INTENT(in   ) ::   kt         ! ocean time-step index 
    68       INTEGER                              , INTENT(in   ) ::   nit000          ! first time step index 
     68      INTEGER                              , INTENT(in   ) ::   kit000          ! first time step index 
    6969      CHARACTER(len=3)                     , INTENT(in   ) ::   cdtype     ! =TRA or TRC (tracer indicator) 
    7070      INTEGER                              , INTENT(in   ) ::   kjpt       ! number of tracers 
    … …  
    7878      !!---------------------------------------------------------------------- 
    7979       
    80       IF( kt == nit000 )  THEN 
     80      IF( kt == kit000 )  THEN 
    8181         IF(lwp) WRITE(numout,*) 
    8282         IF(lwp) WRITE(numout,*) 'tra_ldf_lap : iso-level laplacian diffusion on ', cdtype 

  • branches/2011/dev_r2802_TOP_substepping/NEMOGCM/NEMO/OPA_SRC/TRA/tranxt.F90

    r2892 r2910  
    148148      ELSE                                            ! Leap-Frog + Asselin filter time stepping 
    149149         ! 
    150          IF( lk_vvl )  THEN   ;   CALL tra_nxt_vvl( kt,nit000, 'TRA', tsb, tsn, tsa, jpts )  ! variable volume level (vvl)      
    151          ELSE                 ;   CALL tra_nxt_fix( kt,nit000, 'TRA', tsb, tsn, tsa, jpts )  ! fixed    volume level  
     150         IF( lk_vvl )  THEN   ;   CALL tra_nxt_vvl( kt, nit000, 'TRA', tsb, tsn, tsa, jpts )  ! variable volume level (vvl)      
     151         ELSE                 ;   CALL tra_nxt_fix( kt, nit000, 'TRA', tsb, tsn, tsa, jpts )  ! fixed    volume level  
    152152         ENDIF 
    153153      ENDIF  
    … …  
    179179 
    180180 
    181    SUBROUTINE tra_nxt_fix( kt, nit000, cdtype, ptb, ptn, pta, kjpt ) 
     181   SUBROUTINE tra_nxt_fix( kt, kit000, cdtype, ptb, ptn, pta, kjpt ) 
    182182      !!---------------------------------------------------------------------- 
    183183      !!                   ***  ROUTINE tra_nxt_fix  *** 
    … …  
    203203      !!---------------------------------------------------------------------- 
    204204      INTEGER         , INTENT(in   )                               ::   kt       ! ocean time-step index 
    205       INTEGER         , INTENT(in   )                               ::   nit000   ! first time step index 
     205      INTEGER         , INTENT(in   )                               ::   kit000   ! first time step index 
    206206      CHARACTER(len=3), INTENT(in   )                               ::   cdtype   ! =TRA or TRC (tracer indicator) 
    207207      INTEGER         , INTENT(in   )                               ::   kjpt     ! number of tracers 
    … …  
    215215      !!---------------------------------------------------------------------- 
    216216 
    217       IF( kt == nit000 )  THEN 
     217      IF( kt == kit000 )  THEN 
    218218         IF(lwp) WRITE(numout,*) 
    219219         IF(lwp) WRITE(numout,*) 'tra_nxt_fix : time stepping', cdtype 
    … …  
    246246 
    247247 
    248    SUBROUTINE tra_nxt_vvl( kt, nit000, cdtype, ptb, ptn, pta, kjpt ) 
     248   SUBROUTINE tra_nxt_vvl( kt, kit000, cdtype, ptb, ptn, pta, kjpt ) 
    249249      !!---------------------------------------------------------------------- 
    250250      !!                   ***  ROUTINE tra_nxt_vvl  *** 
    … …  
    271271      !!---------------------------------------------------------------------- 
    272272      INTEGER         , INTENT(in   )                               ::   kt       ! ocean time-step index 
    273       INTEGER         , INTENT(in   )                               ::   nit000   ! first time step index 
     273      INTEGER         , INTENT(in   )                               ::   kit000   ! first time step index 
    274274      CHARACTER(len=3), INTENT(in   )                               ::   cdtype   ! =TRA or TRC (tracer indicator) 
    275275      INTEGER         , INTENT(in   )                               ::   kjpt     ! number of tracers 
    … …  
    284284      !!---------------------------------------------------------------------- 
    285285 
    286       IF( kt == nit000 )  THEN 
     286      IF( kt == kit000 )  THEN 
    287287         IF(lwp) WRITE(numout,*) 
    288288         IF(lwp) WRITE(numout,*) 'tra_nxt_vvl : time stepping', cdtype 

  • branches/2011/dev_r2802_TOP_substepping/NEMOGCM/NEMO/OPA_SRC/TRA/trazdf_exp.F90

    r2892 r2910  
    4848CONTAINS 
    4949 
    50    SUBROUTINE tra_zdf_exp( kt,nit000, cdtype, p2dt, kn_zdfexp,   & 
     50   SUBROUTINE tra_zdf_exp( kt, kit000, cdtype, p2dt, kn_zdfexp,   & 
    5151      &                                ptb , pta      , kjpt ) 
    5252      !!---------------------------------------------------------------------- 
    … …  
    7777      ! 
    7878      INTEGER                              , INTENT(in   ) ::   kt          ! ocean time-step index 
    79       INTEGER                              , INTENT(in   ) ::   nit000      ! first time step index 
     79      INTEGER                              , INTENT(in   ) ::   kit000      ! first time step index 
    8080      CHARACTER(len=3)                     , INTENT(in   ) ::   cdtype      ! =TRA or TRC (tracer indicator) 
    8181      INTEGER                              , INTENT(in   ) ::   kjpt        ! number of tracers 
    … …  
    9494      ENDIF 
    9595 
    96       IF( kt == nit000 )  THEN 
     96      IF( kt == kit000 )  THEN 
    9797         IF(lwp) WRITE(numout,*) 
    9898         IF(lwp) WRITE(numout,*) 'tra_zdf_exp : explicit vertical mixing on ', cdtype 

  • branches/2011/dev_r2802_TOP_substepping/NEMOGCM/NEMO/OPA_SRC/TRA/trazdf_imp.F90

    r2892 r2910  
    5555CONTAINS 
    5656  
    57    SUBROUTINE tra_zdf_imp( kt,nit000, cdtype, p2dt, ptb, pta, kjpt )  
     57   SUBROUTINE tra_zdf_imp( kt, kit000, cdtype, p2dt, ptb, pta, kjpt )  
    5858      !!---------------------------------------------------------------------- 
    5959      !!                  ***  ROUTINE tra_zdf_imp  *** 
    … …  
    7979      ! 
    8080      INTEGER                              , INTENT(in   ) ::   kt       ! ocean time-step index 
    81       INTEGER                              , INTENT(in   ) ::   nit000          ! first time step index 
     81      INTEGER                              , INTENT(in   ) ::   kit000          ! first time step index 
    8282      CHARACTER(len=3)                     , INTENT(in   ) ::   cdtype   ! =TRA or TRC (tracer indicator) 
    8383      INTEGER                              , INTENT(in   ) ::   kjpt     ! number of tracers 
    … …  
    9494      ENDIF 
    9595 
    96       IF( kt == nit000 )  THEN 
     96      IF( kt == kit000 )  THEN 
    9797         IF(lwp)WRITE(numout,*) 
    9898         IF(lwp)WRITE(numout,*) 'tra_zdf_imp : implicit vertical mixing on ', cdtype 
Note: See TracChangeset for help on using the changeset viewer.