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 9124 for branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC – NEMO

Ignore:
Timestamp:
2017-12-19T09:26:25+01:00 (6 years ago)
Author:
gm
Message:

dev_merge_2017: ln_timing instead of nn_timing + restricted timing to nemo_init and routine called by step in OPA_SRC

Location:
branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC
Files:
60 edited

Legend:

Unmodified
Added
Removed

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/AGE/trcsms_age.F90

    r7068 r9124  
    4545      !! ** Method  : - 
    4646      !!---------------------------------------------------------------------- 
    47       ! 
    4847      INTEGER, INTENT(in) ::   kt   ! ocean time-step index 
    4948      INTEGER ::   jn, jk   ! dummy loop index 
    5049      !!---------------------------------------------------------------------- 
    5150      ! 
    52       IF( nn_timing == 1 )  CALL timing_start('trc_sms_age') 
     51      IF( ln_timing )   CALL timing_start('trc_sms_age') 
    5352      ! 
    5453      IF(lwp) WRITE(numout,*) 
    … …  
    5958      DO jk = 1, nla_age 
    6059         tra(:,:,jk,jp_age) = rn_age_kill_rate * trb(:,:,jk,jp_age) 
    61       ENDDO 
     60      END DO 
    6261      ! 
    6362      tra(:,:,nl_age,jp_age) = frac_kill_age * rn_age_kill_rate * trb(:,:,nl_age,jp_age)  & 
    … …  
    6665      DO jk = nlb_age, jpk 
    6766         tra(:,:,jk,jp_age) = tmask(:,:,jk) * rryear 
    68       ENDDO 
     67      END DO 
    6968      ! 
    7069      IF( l_trdtrc ) CALL trd_trc( tra(:,:,:,jp_age), jn, jptra_sms, kt )   ! save trends 
    7170      ! 
    72       IF( nn_timing == 1 )  CALL timing_stop('trc_sms_age') 
     71      IF( ln_timing )   CALL timing_stop('trc_sms_age') 
    7372      ! 
    7473   END SUBROUTINE trc_sms_age 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/C14/trcatm_c14.F90

    r7192 r9124  
    4949      !!---------------------------------------------------------------------- 
    5050      ! 
    51       IF( nn_timing == 1 )  CALL timing_start('trc_atm_c14_ini') 
     51      IF( ln_timing )   CALL timing_start('trc_atm_c14_ini') 
    5252      ! 
    5353       
    … …  
    192192      ENDIF 
    193193      ! 
    194       IF( nn_timing == 1 )  CALL timing_stop('trc_atm_c14_ini') 
     194      IF( ln_timing )   CALL timing_stop('trc_atm_c14_ini') 
    195195      ! 
    196196   END SUBROUTINE trc_atm_c14_ini 
    … …  
    216216      !!---------------------------------------------------------------------- 
    217217      ! 
    218       IF( nn_timing == 1 )  CALL timing_start('trc_atm_c14') 
     218      IF( ln_timing )   CALL timing_start('trc_atm_c14') 
    219219      ! 
    220220      IF( kc14typ == 0) THEN 
    … …  
    297297      ENDIF 
    298298      ! 
    299       IF( nn_timing == 1 )  CALL timing_stop('trc_atm_c14') 
     299      IF( ln_timing )   CALL timing_stop('trc_atm_c14') 
    300300      ! 
    301301   END SUBROUTINE trc_atm_c14 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/C14/trcsms_c14.F90

    r7068 r9124  
    2828   !!---------------------------------------------------------------------- 
    2929   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    30    !! $                                                                    $  
     30   !! $Id:$  
    3131   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
    3232   !!---------------------------------------------------------------------- 
    … …  
    6060      !!---------------------------------------------------------------------- 
    6161      ! 
    62       IF( nn_timing == 1 )  CALL timing_start('trc_sms_c14') 
     62      IF( ln_timing )   CALL timing_start('trc_sms_c14') 
    6363      ! 
    6464      IF( kt ==  nittrc000 ) THEN 
    … …  
    6868      ENDIF 
    6969      ! 
    70   
    7170      ! Get co2sbc & c14sbc(ji,jj): at 1st iter for all, at each time step for transient 
    7271      IF( kc14typ >= 1 .OR.  kt ==  nittrc000 )   CALL trc_atm_c14( kt, co2sbc, c14sbc )  
    … …  
    8079      DO jj = 1, jpj 
    8180         DO ji = 1, jpi   
    82             IF( tmask(ji,jj,1) >  0.) THEN 
    83  
    84                zt   = MIN( 40., tsn(ji,jj,1,jp_tem)) 
    85  
     81            IF( tmask(ji,jj,1) >  0. ) THEN 
     82               ! 
     83               zt   = MIN( 40. , tsn(ji,jj,1,jp_tem) ) 
     84               ! 
    8685               !  Computation of solubility zsol in [mol/(L * atm)] 
    8786               !   after Wanninkhof (2014) referencing Weiss (1974) 
    … …  
    143142         END DO 
    144143      END DO 
    145  
    146144      ! 
    147145      IF( lrst_trc ) THEN 
    … …  
    149147         IF(lwp) WRITE(numout,*) ' trc_rst_wri_c14 : Write specific variables from c14 model ' 
    150148         IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~~~~~' 
    151       ! 
     149         ! 
    152150         CALL iom_rstput( kt, nitrst, numrtw, 'co2sbc', co2sbc )       ! These five need      & 
    153151         CALL iom_rstput( kt, nitrst, numrtw, 'c14sbc', c14sbc )     ! &    to be written   & 
    … …  
    156154         CALL iom_rstput( kt, nitrst, numrtw, 'qtr_c14', qtr_c14 )   ! &    to be coherent. 
    157155         CALL iom_rstput( kt, nitrst, numrtw, 'qint_c14', qint_c14 ) ! Cumulative 
    158       ! 
     156         ! 
    159157      ENDIF 
    160158 
    161159      IF( l_trdtrc )  CALL trd_trc( tra(:,:,:,jp_c14), 1, jptra_sms, kt )   ! save trends 
    162160      ! 
    163       IF( nn_timing == 1 )  CALL timing_stop('trc_sms_c14') 
     161      IF( ln_timing )   CALL timing_stop('trc_sms_c14') 
    164162      ! 
    165163   END SUBROUTINE trc_sms_c14 
    166164 
    167  
    168165  !!====================================================================== 
    169166END MODULE trcsms_c14 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/CFC/trcsms_cfc.F90

    r8397 r9124  
    8585      !!---------------------------------------------------------------------- 
    8686      ! 
    87       ! 
    88       IF( nn_timing == 1 )  CALL timing_start('trc_sms_cfc') 
     87      IF( ln_timing )   CALL timing_start('trc_sms_cfc') 
    8988      ! 
    9089      ALLOCATE( zpatm(jphem,jp_cfc), STAT=ierr ) 
    … …  
    194193      END IF 
    195194      ! 
    196       IF( nn_timing == 1 )  CALL timing_stop('trc_sms_cfc') 
     195      IF( ln_timing )   CALL timing_stop('trc_sms_cfc') 
    197196      ! 
    198197   END SUBROUTINE trc_sms_cfc 
    … …  
    205204      !! ** Purpose : sets constants for CFC model 
    206205      !!--------------------------------------------------------------------- 
    207       INTEGER :: jn 
     206      INTEGER ::   jn   ! 
    208207      !!---------------------------------------------------------------------- 
    209208      ! 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/MY_TRC/trcsms_my_trc.F90

    r7646 r9124  
    4646      !!---------------------------------------------------------------------- 
    4747      ! 
    48       IF( nn_timing == 1 )  CALL timing_start('trc_sms_my_trc') 
     48      IF( ln_timing )   CALL timing_start('trc_sms_my_trc') 
    4949      ! 
    5050      IF(lwp) WRITE(numout,*) 
    … …  
    6767      END IF 
    6868      ! 
    69       IF( nn_timing == 1 )  CALL timing_stop('trc_sms_my_trc') 
     69      IF( ln_timing )   CALL timing_stop('trc_sms_my_trc') 
    7070      ! 
    7171   END SUBROUTINE trc_sms_my_trc 
     72 
    7273 
    7374   INTEGER FUNCTION trc_sms_my_trc_alloc() 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/P2Z/p2zbio.F90

    r9085 r9124  
    101101      !!--------------------------------------------------------------------- 
    102102      ! 
    103       IF( nn_timing == 1 )  CALL timing_start('p2z_bio') 
     103      IF( ln_timing )   CALL timing_start('p2z_bio') 
    104104      ! 
    105105      IF( lk_iomput ) THEN 
    … …  
    382382      ENDIF 
    383383      ! 
    384       IF( nn_timing == 1 )  CALL timing_stop('p2z_bio') 
     384      IF( ln_timing )  CALL timing_stop('p2z_bio') 
    385385      ! 
    386386   END SUBROUTINE p2z_bio 
     387 
    387388 
    388389   SUBROUTINE p2z_bio_init 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/P2Z/p2zexp.F90

    r7646 r9124  
    6868      !!--------------------------------------------------------------------- 
    6969      ! 
    70       IF( nn_timing == 1 )  CALL timing_start('p2z_exp') 
     70      IF( ln_timing )   CALL timing_start('p2z_exp') 
    7171      ! 
    7272      IF( kt == nittrc000 )   CALL p2z_exp_init 
    … …  
    155155      ENDIF 
    156156      ! 
    157       IF( nn_timing == 1 )  CALL timing_stop('p2z_exp') 
     157      IF( ln_timing )  CALL timing_stop('p2z_exp') 
    158158      ! 
    159159   END SUBROUTINE p2z_exp 
     160 
    160161 
    161162   SUBROUTINE p2z_exp_init 
    … …  
    166167      INTEGER  ::   ji, jj, jk 
    167168      REAL(wp) ::   zmaskt, zfluo, zfluu 
    168       REAL(wp), POINTER, DIMENSION(:,:  ) :: zrro 
    169       REAL(wp), POINTER, DIMENSION(:,:,:) :: zdm0 
    170       !!--------------------------------------------------------------------- 
    171  
     169      REAL(wp), POINTER, DIMENSION(:,:  ) ::   zrro 
     170      REAL(wp), POINTER, DIMENSION(:,:,:) ::   zdm0 
     171      !!--------------------------------------------------------------------- 
     172      ! 
    172173      IF(lwp) THEN 
    173174         WRITE(numout,*) 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/P2Z/p2zopt.F90

    r9019 r9124  
    1111   !!             3.2  !  2009-04  (C. Ethe, G. Madec)  minor optimisation + style 
    1212   !!---------------------------------------------------------------------- 
     13 
     14   !!---------------------------------------------------------------------- 
    1315   !!   p2z_opt        :   Compute the light availability in the water column 
    1416   !!---------------------------------------------------------------------- 
    … …  
    6971      !!--------------------------------------------------------------------- 
    7072      ! 
    71       IF( nn_timing == 1 )  CALL timing_start('p2z_opt') 
     73      IF( ln_timing )   CALL timing_start('p2z_opt') 
    7274      ! 
    7375      ! Allocate temporary workspace 
    … …  
    143145      CALL wrk_dealloc( jpi, jpj, jpk, zparr, zparg    ) 
    144146      ! 
    145       IF( nn_timing == 1 )  CALL timing_stop('p2z_opt') 
     147      IF( ln_timing )   CALL timing_stop('p2z_opt') 
    146148      ! 
    147149   END SUBROUTINE p2z_opt 
    … …  
    157159      !! 
    158160      !!---------------------------------------------------------------------- 
     161      INTEGER ::   ios   ! Local integer 
     162      !! 
    159163      NAMELIST/namlobopt/ xkg0, xkr0, xkgp, xkrp, xlg, xlr, rpig 
    160164      NAMELIST/namlobrat/ rcchl, redf, reddom 
    161       INTEGER :: ios                 ! Local integer output status for namelist read 
    162165      !!---------------------------------------------------------------------- 
    163166 
    … …  
    168171      REWIND( numnatp_cfg )              ! Namelist namlobopt in configuration namelist : Lobster options 
    169172      READ  ( numnatp_cfg, namlobopt, IOSTAT = ios, ERR = 902 ) 
    170 902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namlobopt in configuration namelist', lwp ) 
     173902   IF( ios > 0 ) CALL ctl_nam ( ios , 'namlobopt in configuration namelist', lwp ) 
    171174      IF(lwm) WRITE ( numonp, namlobopt ) 
    172175 
    … …  
    190193      REWIND( numnatp_cfg )              ! Namelist namlobrat in configuration namelist : Lobster ratios 
    191194      READ  ( numnatp_cfg, namlobrat, IOSTAT = ios, ERR = 904 ) 
    192 904   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namlobrat in configuration namelist', lwp ) 
     195904   IF( ios > 0 ) CALL ctl_nam ( ios , 'namlobrat in configuration namelist', lwp ) 
    193196      IF(lwm) WRITE ( numonp, namlobrat ) 
    194197 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/P2Z/p2zsed.F90

    r7646 r9124  
    3535   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    3636   !!---------------------------------------------------------------------- 
    37  
    3837CONTAINS 
    3938 
    … …  
    5655      !!              trend of passive tracers is saved for futher diagnostics. 
    5756      !!--------------------------------------------------------------------- 
    58       !! 
    5957      INTEGER, INTENT( in ) ::   kt      ! ocean time-step index       
    60       !! 
     58      ! 
    6159      INTEGER  ::   ji, jj, jk, jl, ierr 
    6260      CHARACTER (len=25) :: charout 
    … …  
    6563      !!--------------------------------------------------------------------- 
    6664      ! 
    67       IF( nn_timing == 1 )  CALL timing_start('p2z_sed') 
     65      IF( ln_timing )   CALL timing_start('p2z_sed') 
    6866      ! 
    6967      IF( kt == nittrc000 ) THEN 
    … …  
    119117      ENDIF 
    120118      ! 
    121       IF( nn_timing == 1 )  CALL timing_stop('p2z_sed') 
     119      IF( ln_timing )   CALL timing_stop('p2z_sed') 
    122120      ! 
    123121   END SUBROUTINE p2z_sed 
     122 
    124123 
    125124   SUBROUTINE p2z_sed_init 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/P2Z/p2zsms.F90

    r7681 r9124  
    66   !! History :   1.0  !            M. Levy 
    77   !!             2.0  !  2007-12  (C. Ethe, G. Madec)  revised architecture 
     8   !!---------------------------------------------------------------------- 
     9 
    810   !!---------------------------------------------------------------------- 
    911   !!   p2zsms        :  Time loop of passive tracers sms 
    … …  
    4244      !! ** Method  : - ??? 
    4345      !! -------------------------------------------------------------------- 
    44       !! 
    4546      INTEGER, INTENT( in ) ::   kt      ! ocean time-step index       
    4647      ! 
    47       INTEGER :: jn 
     48      INTEGER ::   jn   ! dummy loop index 
    4849      !! -------------------------------------------------------------------- 
    4950      ! 
    50       IF( nn_timing == 1 )  CALL timing_start('p2z_sms') 
     51      IF( ln_timing )   CALL timing_start('p2z_sms') 
    5152      ! 
    5253      CALL p2z_opt( kt )      ! optical model 
    … …  
    5455      CALL p2z_sed( kt )      ! sedimentation model 
    5556      CALL p2z_exp( kt )      ! export 
    56  
     57      ! 
    5758      IF( l_trdtrc ) THEN 
    5859         DO jn = jp_pcs0, jp_pcs1 
    … …  
    6061         END DO 
    6162      END IF 
    62  
    6363      ! 
    6464      IF ( lwm .AND. kt == nittrc000 ) CALL FLUSH    ( numonp )     ! flush output namelist PISCES 
    65       IF( nn_timing == 1 )  CALL timing_stop('p2z_sms') 
     65      IF( ln_timing )   CALL timing_stop('p2z_sms') 
    6666      ! 
    6767   END SUBROUTINE p2z_sms 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zagg.F90

    r7753 r9124  
    1010   !!             3.6  !  2015-05  (O. Aumont) PISCES quota 
    1111   !!---------------------------------------------------------------------- 
     12 
    1213   !!---------------------------------------------------------------------- 
    1314   !!   p4z_agg       :  Compute aggregation of particles 
    … …  
    3839      !! ** Method  : - ??? 
    3940      !!--------------------------------------------------------------------- 
    40       INTEGER, INTENT(in) :: kt, knt 
     41      INTEGER, INTENT(in) ::   kt, knt   ! 
     42      ! 
    4143      INTEGER  ::   ji, jj, jk 
    4244      REAL(wp) ::   zagg, zagg1, zagg2, zagg3, zagg4 
    … …  
    4951      !!--------------------------------------------------------------------- 
    5052      ! 
    51       IF( nn_timing == 1 )  CALL timing_start('p4z_agg') 
    52  
     53      IF( ln_timing )   CALL timing_start('p4z_agg') 
    5354      ! 
    5455      !  Exchange between organic matter compartments due to coagulation/disaggregation 
    … …  
    178179      ENDIF 
    179180      ! 
    180       IF( nn_timing == 1 )  CALL timing_stop('p4z_agg') 
     181      IF( ln_timing )   CALL timing_stop('p4z_agg') 
    181182      ! 
    182183   END SUBROUTINE p4z_agg 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zbio.F90

    r7753 r9124  
    5656      !!--------------------------------------------------------------------- 
    5757      INTEGER, INTENT(in) :: kt, knt 
     58      ! 
    5859      INTEGER             :: ji, jj, jk, jn 
    5960      CHARACTER (len=25) :: charout 
    60  
    6161      !!--------------------------------------------------------------------- 
    6262      ! 
    63       IF( nn_timing == 1 )  CALL timing_start('p4z_bio') 
     63      IF( ln_timing )   CALL timing_start('p4z_bio') 
    6464      ! 
    6565      !     ASSIGN THE SHEAR RATE THAT IS USED FOR AGGREGATION 
    … …  
    112112      ENDIF 
    113113      ! 
    114       IF( nn_timing == 1 )  CALL timing_stop('p4z_bio') 
     114      IF( ln_timing )   CALL timing_stop('p4z_bio') 
    115115      ! 
    116116   END SUBROUTINE p4z_bio 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zche.F90

    r7753 r9124  
    155155      REAL(wp) ::   zst  , zft  , zcks  , zckf  , zaksp1 
    156156      REAL(wp) ::   total2free, free2SWS, total2SWS, SWS2total 
    157  
    158157      !!--------------------------------------------------------------------- 
    159158      ! 
    160       IF( nn_timing == 1 )  CALL timing_start('p4z_che') 
     159      IF( ln_timing )   CALL timing_start('p4z_che') 
    161160      ! 
    162161      ! Computation of chemical constants require practical salinity 
    … …  
    447446      END DO 
    448447      ! 
    449       IF( nn_timing == 1 )  CALL timing_stop('p4z_che') 
     448      IF( ln_timing )  CALL timing_stop('p4z_che') 
    450449      ! 
    451450   END SUBROUTINE p4z_che 
     451 
    452452 
    453453   SUBROUTINE ahini_for_at(p_hini) 
    … …  
    470470      REAL(wp)  ::  za2, za1, za0 
    471471      REAL(wp)  ::  p_dictot, p_bortot, p_alkcb  
    472  
    473       IF( nn_timing == 1 )  CALL timing_start('ahini_for_at') 
     472      !!--------------------------------------------------------------------- 
     473 
     474      IF( ln_timing )  CALL timing_start('ahini_for_at') 
    474475      ! 
    475476      DO jk = 1, jpk 
    … …  
    512513      END DO 
    513514      ! 
    514       IF( nn_timing == 1 )  CALL timing_stop('ahini_for_at') 
     515      IF( ln_timing )  CALL timing_stop('ahini_for_at') 
    515516      ! 
    516517   END SUBROUTINE ahini_for_at 
    517518 
    518519   !=============================================================================== 
     520 
    519521   SUBROUTINE anw_infsup( p_alknw_inf, p_alknw_sup ) 
    520522 
    … …  
    564566   REAL(wp), POINTER, DIMENSION(:,:,:) :: zalknw_inf, zalknw_sup, rmask, zh_min, zh_max, zeqn_absmin 
    565567 
    566    IF( nn_timing == 1 )  CALL timing_start('solve_at_general') 
     568   IF( ln_timing )  CALL timing_start('solve_at_general') 
    567569      !  Allocate temporary workspace 
    568570   CALL wrk_alloc( jpi, jpj, jpk, zalknw_inf, zalknw_sup, rmask ) 
    … …  
    801803 
    802804 
    803    IF( nn_timing == 1 )  CALL timing_stop('solve_at_general') 
    804  
    805  
     805      IF( ln_timing )   CALL timing_stop('solve_at_general') 
     806      ! 
    806807   END SUBROUTINE solve_at_general 
     808 
    807809 
    808810   INTEGER FUNCTION p4z_che_alloc() 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zfechem.F90

    r7753 r9124  
    7979      !!--------------------------------------------------------------------- 
    8080      ! 
    81       IF( nn_timing == 1 )  CALL timing_start('p4z_fechem') 
     81      IF( ln_timing )   CALL timing_start('p4z_fechem') 
    8282      ! 
    8383      ! Allocate temporary workspace 
    … …  
    366366      ENDIF 
    367367      ! 
    368       IF( nn_timing == 1 )  CALL timing_stop('p4z_fechem') 
     368      IF( ln_timing )   CALL timing_stop('p4z_fechem') 
    369369      ! 
    370370   END SUBROUTINE p4z_fechem 
    … …  
    383383      !! 
    384384      !!---------------------------------------------------------------------- 
     385      INTEGER ::   ios   ! Local integer  
     386      !! 
    385387      NAMELIST/nampisfer/ ln_fechem, ln_ligvar, ln_fecolloid, xlam1, xlamdust, ligand, kfep  
    386       INTEGER :: ios                 ! Local integer output status for namelist read 
     388      !!---------------------------------------------------------------------- 
    387389 
    388390      REWIND( numnatp_ref )              ! Namelist nampisfer in reference namelist : Pisces iron chemistry 
    … …  
    392394      REWIND( numnatp_cfg )              ! Namelist nampisfer in configuration namelist : Pisces iron chemistry 
    393395      READ  ( numnatp_cfg, nampisfer, IOSTAT = ios, ERR = 902 ) 
    394 902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampisfer in configuration namelist', lwp ) 
     396902   IF( ios > 0 ) CALL ctl_nam ( ios , 'nampisfer in configuration namelist', lwp ) 
    395397      IF(lwm) WRITE ( numonp, nampisfer ) 
    396398 
    … …  
    434436      ! 
    435437   END SUBROUTINE p4z_fechem_init 
     438    
    436439   !!====================================================================== 
    437440END MODULE p4zfechem 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zflx.F90

    r7753 r9124  
    8383      !!--------------------------------------------------------------------- 
    8484      ! 
    85       IF( nn_timing == 1 )  CALL timing_start('p4z_flx') 
     85      IF( ln_timing )   CALL timing_start('p4z_flx') 
    8686      ! 
    8787      CALL wrk_alloc( jpi, jpj, zkgco2, zkgo2, zh2co3, zoflx, zpco2atm ) 
    … …  
    216216      CALL wrk_dealloc( jpi, jpj, zkgco2, zkgo2, zh2co3, zoflx, zpco2atm ) 
    217217      ! 
    218       IF( nn_timing == 1 )  CALL timing_stop('p4z_flx') 
     218      IF( ln_timing )   CALL timing_stop('p4z_flx') 
    219219      ! 
    220220   END SUBROUTINE p4z_flx 
    … …  
    231231      !! ** input   :   Namelist nampisext 
    232232      !!---------------------------------------------------------------------- 
     233      INTEGER ::   jm 
     234      INTEGER ::   ios   ! Local integer  
     235      ! 
    233236      NAMELIST/nampisext/ln_co2int, atcco2, clname, nn_offset 
    234       INTEGER :: jm 
    235       INTEGER :: ios                 ! Local integer output status for namelist read 
    236237      !!---------------------------------------------------------------------- 
    237238      ! 
    … …  
    243244      REWIND( numnatp_cfg )              ! Namelist nampisext in configuration namelist : Pisces atm. conditions 
    244245      READ  ( numnatp_cfg, nampisext, IOSTAT = ios, ERR = 902 ) 
    245 902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampisext in configuration namelist', lwp ) 
     246902   IF( ios > 0 ) CALL ctl_nam ( ios , 'nampisext in configuration namelist', lwp ) 
    246247      IF(lwm) WRITE ( numonp, nampisext ) 
    247248      ! 
    … …  
    292293         ENDIF 
    293294      ENDIF 
    294  
    295295      ! 
    296296      oce_co2(:,:)  = 0._wp                ! Initialization of Flux of Carbon 
    … …  
    302302   END SUBROUTINE p4z_flx_init 
    303303 
     304 
    304305   SUBROUTINE p4z_patm( kt ) 
    305  
    306306      !!---------------------------------------------------------------------- 
    307307      !!                  ***  ROUTINE p4z_atm  *** 
    … …  
    311311      !! 
    312312      !!---------------------------------------------------------------------- 
    313       !! * arguments 
    314313      INTEGER, INTENT( in  ) ::   kt   ! ocean time step 
    315314      ! 
    … …  
    321320      !! 
    322321      NAMELIST/nampisatm/ ln_presatm, ln_presatmco2, sn_patm, sn_atmco2, cn_dir 
     322      !!---------------------------------------------------------------------- 
    323323 
    324324      !                                         ! ----------------------- ! 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zint.F90

    r7753 r9124  
    3333      !! 
    3434      !!--------------------------------------------------------------------- 
    35       ! 
    3635      INTEGER, INTENT( in ) ::   kt      ! ocean time-step index 
    3736      ! 
    … …  
    4039      !!--------------------------------------------------------------------- 
    4140      ! 
    42       IF( nn_timing == 1 )  CALL timing_start('p4z_int') 
     41      IF( ln_timing )   CALL timing_start('p4z_int') 
    4342      ! 
    4443      ! Computation of phyto and zoo metabolic rate 
    … …  
    6160      ENDIF 
    6261      ! 
    63       IF( nn_timing == 1 )  CALL timing_stop('p4z_int') 
     62      IF( ln_timing )   CALL timing_stop('p4z_int') 
    6463      ! 
    6564   END SUBROUTINE p4z_int 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zligand.F90

    r7646 r9124  
    5050      !!--------------------------------------------------------------------- 
    5151      ! 
    52       IF( nn_timing == 1 )  CALL timing_start('p4z_ligand') 
     52      IF( ln_timing )   CALL timing_start('p4z_ligand') 
    5353      ! 
    5454      ! ------------------------------------------------------------------ 
    … …  
    7575      ! Dissolution of nanoparticle Fe 
    7676      ! ---------------------------------------------------------- 
    77  
    7877      DO jk = 1, jpkm1 
    7978         DO jj = 1, jpj 
    … …  
    9796       ENDIF 
    9897      ! 
    99       IF( nn_timing == 1 )  CALL timing_stop('p4z_ligand') 
     98      IF( ln_timing )   CALL timing_stop('p4z_ligand') 
    10099      ! 
    101100   END SUBROUTINE p4z_ligand 
    … …  
    114113      !! 
    115114      !!---------------------------------------------------------------------- 
     115      INTEGER ::   ios   ! Local integer  
     116      ! 
    116117      NAMELIST/nampislig/ rlgw, prlgw, rlgs, rfep, rlig 
    117       INTEGER :: ios                 ! Local integer output status for namelist read 
     118      !!---------------------------------------------------------------------- 
    118119 
    119120      REWIND( numnatp_ref )              ! Namelist nampislig in reference namelist : Pisces remineralization 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zlim.F90

    r8533 r9124  
    7272   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    7373   !!---------------------------------------------------------------------- 
    74  
    7574CONTAINS 
    7675 
    … …  
    8483      !! ** Method  : - ??? 
    8584      !!--------------------------------------------------------------------- 
    86       ! 
    8785      INTEGER, INTENT(in)  :: kt, knt 
    8886      ! 
    … …  
    9593      !!--------------------------------------------------------------------- 
    9694      ! 
    97       IF( nn_timing == 1 )  CALL timing_start('p4z_lim') 
     95      IF( ln_timing )   CALL timing_start('p4z_lim') 
    9896      ! 
    9997      DO jk = 1, jpkm1 
    … …  
    224222      ENDIF 
    225223      ! 
    226       IF( nn_timing == 1 )  CALL timing_stop('p4z_lim') 
     224      IF( ln_timing )   CALL timing_stop('p4z_lim') 
    227225      ! 
    228226   END SUBROUTINE p4z_lim 
    229227 
     228 
    230229   SUBROUTINE p4z_lim_init 
    231  
    232230      !!---------------------------------------------------------------------- 
    233231      !!                  ***  ROUTINE p4z_lim_init  *** 
    … …  
    241239      !! 
    242240      !!---------------------------------------------------------------------- 
    243  
     241      INTEGER ::   ios   ! Local integer 
     242      ! 
    244243      NAMELIST/namp4zlim/ concnno3, concdno3, concnnh4, concdnh4, concnfer, concdfer, concbfe,   & 
    245244         &                concbno3, concbnh4, xsizedia, xsizephy, xsizern, xsizerd,          &  
    246245         &                xksi1, xksi2, xkdoc, qnfelim, qdfelim, caco3r, oxymin 
    247       INTEGER :: ios                 ! Local integer output status for namelist read 
    248  
     246      !!---------------------------------------------------------------------- 
     247      ! 
    249248      REWIND( numnatp_ref )              ! Namelist nampislim in reference namelist : Pisces nutrient limitation parameters 
    250249      READ  ( numnatp_ref, namp4zlim, IOSTAT = ios, ERR = 901) 
    251250901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namp4zlim in reference namelist', lwp ) 
    252  
     251      ! 
    253252      REWIND( numnatp_cfg )              ! Namelist nampislim in configuration namelist : Pisces nutrient limitation parameters  
    254253      READ  ( numnatp_cfg, namp4zlim, IOSTAT = ios, ERR = 902 ) 
    255 902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namp4zlim in configuration namelist', lwp ) 
     254902   IF( ios > 0 ) CALL ctl_nam ( ios , 'namp4zlim in configuration namelist', lwp ) 
    256255      IF(lwm) WRITE ( numonp, namp4zlim ) 
    257  
     256      ! 
    258257      IF(lwp) THEN                         ! control print 
    259258         WRITE(numout,*) ' ' 
    … …  
    286285   END SUBROUTINE p4z_lim_init 
    287286 
     287 
    288288   INTEGER FUNCTION p4z_lim_alloc() 
    289289      !!---------------------------------------------------------------------- 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zlys.F90

    r7753 r9124  
    6565      !!--------------------------------------------------------------------- 
    6666      ! 
    67       IF( nn_timing == 1 )  CALL timing_start('p4z_lys') 
     67      IF( ln_timing )  CALL timing_start('p4z_lys') 
    6868      ! 
    6969      CALL wrk_alloc( jpi, jpj, jpk, zco3, zcaldiss, zhinit, zhi, zco3sat ) 
    … …  
    141141      CALL wrk_dealloc( jpi, jpj, jpk, zco3, zcaldiss, zhinit, zhi, zco3sat ) 
    142142      ! 
    143       IF( nn_timing == 1 )  CALL timing_stop('p4z_lys') 
     143      IF( ln_timing )   CALL timing_stop('p4z_lys') 
    144144      ! 
    145145   END SUBROUTINE p4z_lys 
    146146 
     147 
    147148   SUBROUTINE p4z_lys_init 
    148  
    149149      !!---------------------------------------------------------------------- 
    150150      !!                  ***  ROUTINE p4z_lys_init  *** 
    … …  
    158158      !! 
    159159      !!---------------------------------------------------------------------- 
    160       INTEGER  ::  ios                 ! Local integer output status for namelist read 
    161  
     160      INTEGER ::   ios   ! Local integer 
     161      ! 
    162162      NAMELIST/nampiscal/ kdca, nca 
    163163      !!---------------------------------------------------------------------- 
    164  
     164      ! 
    165165      REWIND( numnatp_ref )              ! Namelist nampiscal in reference namelist : Pisces CaCO3 dissolution 
    166166      READ  ( numnatp_ref, nampiscal, IOSTAT = ios, ERR = 901) 
    167167901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampiscal in reference namelist', lwp ) 
    168  
     168      ! 
    169169      REWIND( numnatp_cfg )              ! Namelist nampiscal in configuration namelist : Pisces CaCO3 dissolution 
    170170      READ  ( numnatp_cfg, nampiscal, IOSTAT = ios, ERR = 902 ) 
    171 902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampiscal in configuration namelist', lwp ) 
     171902   IF( ios > 0 ) CALL ctl_nam ( ios , 'nampiscal in configuration namelist', lwp ) 
    172172      IF(lwm) WRITE ( numonp, nampiscal ) 
    173  
     173      ! 
    174174      IF(lwp) THEN                         ! control print 
    175175         WRITE(numout,*) ' ' 
    … …  
    179179         WRITE(numout,*) '    order of reaction for calcite dissolution nca       =', nca 
    180180      ENDIF 
    181  
     181      ! 
    182182      ! Number of seconds per month  
    183183      rmtss =  nyear_len(1) * rday / raamo 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zmeso.F90

    r8533 r9124  
    7676      !!--------------------------------------------------------------------- 
    7777      ! 
    78       IF( nn_timing == 1 )  CALL timing_start('p4z_meso') 
     78      IF( ln_timing )   CALL timing_start('p4z_meso') 
    7979      ! 
    8080      CALL wrk_alloc( jpi, jpj, jpk, zgrazing ) 
    … …  
    240240      CALL wrk_dealloc( jpi, jpj, jpk, zgrazing ) 
    241241      ! 
    242       IF( nn_timing == 1 )  CALL timing_stop('p4z_meso') 
     242      IF( ln_timing )   CALL timing_stop('p4z_meso') 
    243243      ! 
    244244   END SUBROUTINE p4z_meso 
    245245 
     246 
    246247   SUBROUTINE p4z_meso_init 
    247  
    248248      !!---------------------------------------------------------------------- 
    249249      !!                  ***  ROUTINE p4z_meso_init  *** 
    … …  
    257257      !! 
    258258      !!---------------------------------------------------------------------- 
    259  
     259      INTEGER ::   ios   ! Local integer 
     260      ! 
    260261      NAMELIST/namp4zmes/ part2, grazrat2, resrat2, mzrat2, xprefc, xprefp, xprefz,   & 
    261262         &                xprefpoc, xthresh2dia, xthresh2phy, xthresh2zoo, xthresh2poc, & 
    262263         &                xthresh2, xkgraz2, epsher2, sigma2, unass2, grazflux 
    263       INTEGER :: ios                 ! Local integer output status for namelist read 
    264  
     264      !!---------------------------------------------------------------------- 
     265      ! 
    265266      REWIND( numnatp_ref )              ! Namelist nampismes in reference namelist : Pisces mesozooplankton 
    266267      READ  ( numnatp_ref, namp4zmes, IOSTAT = ios, ERR = 901) 
    267268901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namp4zmes in reference namelist', lwp ) 
    268  
     269      ! 
    269270      REWIND( numnatp_cfg )              ! Namelist nampismes in configuration namelist : Pisces mesozooplankton 
    270271      READ  ( numnatp_cfg, namp4zmes, IOSTAT = ios, ERR = 902 ) 
    271 902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namp4zmes in configuration namelist', lwp ) 
     272902   IF( ios > 0 ) CALL ctl_nam ( ios , 'namp4zmes in configuration namelist', lwp ) 
    272273      IF(lwm) WRITE ( numonp, namp4zmes ) 
    273  
    274  
     274      ! 
    275275      IF(lwp) THEN                         ! control print 
    276276         WRITE(numout,*) ' '  
    … …  
    296296         WRITE(numout,*) '    half sturation constant for grazing 2          xkgraz2      =', xkgraz2 
    297297      ENDIF 
    298  
    299  
     298      ! 
    300299   END SUBROUTINE p4z_meso_init 
    301300 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zmicro.F90

    r8533 r9124  
    4242   REAL(wp), PUBLIC ::  epsher      !: half sturation constant for grazing 1  
    4343 
    44  
    4544   !!---------------------------------------------------------------------- 
    4645   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    … …  
    4847   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    4948   !!---------------------------------------------------------------------- 
    50  
    5149CONTAINS 
    5250 
    … …  
    7573      !!--------------------------------------------------------------------- 
    7674      ! 
    77       IF( nn_timing == 1 )  CALL timing_start('p4z_micro') 
     75      IF( ln_timing )   CALL timing_start('p4z_micro') 
    7876      ! 
    7977      CALL wrk_alloc( jpi, jpj, jpk, zgrazing ) 
    … …  
    196194      CALL wrk_dealloc( jpi, jpj, jpk, zgrazing ) 
    197195      ! 
    198       IF( nn_timing == 1 )  CALL timing_stop('p4z_micro') 
     196      IF( ln_timing )   CALL timing_stop('p4z_micro') 
    199197      ! 
    200198   END SUBROUTINE p4z_micro 
    … …  
    202200 
    203201   SUBROUTINE p4z_micro_init 
    204  
    205202      !!---------------------------------------------------------------------- 
    206203      !!                  ***  ROUTINE p4z_micro_init  *** 
    … …  
    214211      !! 
    215212      !!---------------------------------------------------------------------- 
    216  
     213      INTEGER ::   ios   ! Local integer 
     214      ! 
    217215      NAMELIST/namp4zzoo/ part, grazrat, resrat, mzrat, xpref2c, xpref2p, & 
    218216         &                xpref2d,  xthreshdia,  xthreshphy,  xthreshpoc, & 
    219217         &                xthresh, xkgraz, epsher, sigma1, unass 
    220       INTEGER :: ios                 ! Local integer output status for namelist read 
    221  
     218      !!---------------------------------------------------------------------- 
     219      ! 
    222220      REWIND( numnatp_ref )              ! Namelist nampiszoo in reference namelist : Pisces microzooplankton 
    223221      READ  ( numnatp_ref, namp4zzoo, IOSTAT = ios, ERR = 901) 
    224222901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namp4zzoo in reference namelist', lwp ) 
    225  
     223      ! 
    226224      REWIND( numnatp_cfg )              ! Namelist nampiszoo in configuration namelist : Pisces microzooplankton 
    227225      READ  ( numnatp_cfg, namp4zzoo, IOSTAT = ios, ERR = 902 ) 
    228 902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namp4zzoo in configuration namelist', lwp ) 
     226902   IF( ios > 0 ) CALL ctl_nam ( ios , 'namp4zzoo in configuration namelist', lwp ) 
    229227      IF(lwm) WRITE ( numonp, namp4zzoo ) 
    230  
     228      ! 
    231229      IF(lwp) THEN                         ! control print 
    232230         WRITE(numout,*) ' ' 
    … …  
    249247         WRITE(numout,*) '    half sturation constant for grazing 1           xkgraz      =', xkgraz 
    250248      ENDIF 
    251  
     249      ! 
    252250   END SUBROUTINE p4z_micro_init 
    253251 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zmort.F90

    r7753 r9124  
    7272      !!--------------------------------------------------------------------- 
    7373      ! 
    74       IF( nn_timing == 1 )  CALL timing_start('p4z_nano') 
     74      IF( ln_timing )   CALL timing_start('p4z_nano') 
    7575      ! 
    7676      prodcal(:,:,:) = 0.  !: calcite production variable set to zero 
    … …  
    126126       ENDIF 
    127127      ! 
    128       IF( nn_timing == 1 )  CALL timing_stop('p4z_nano') 
     128      IF( ln_timing )   CALL timing_stop('p4z_nano') 
    129129      ! 
    130130   END SUBROUTINE p4z_nano 
     131 
    131132 
    132133   SUBROUTINE p4z_diat 
    … …  
    145146      !!--------------------------------------------------------------------- 
    146147      ! 
    147       IF( nn_timing == 1 )  CALL timing_start('p4z_diat') 
     148      IF( ln_timing )   CALL timing_start('p4z_diat') 
    148149      ! 
    149150 
    … …  
    201202      ENDIF 
    202203      ! 
    203       IF( nn_timing == 1 )  CALL timing_stop('p4z_diat') 
     204      IF( ln_timing )   CALL timing_stop('p4z_diat') 
    204205      ! 
    205206   END SUBROUTINE p4z_diat 
    206207 
     208 
    207209   SUBROUTINE p4z_mort_init 
    208  
    209210      !!---------------------------------------------------------------------- 
    210211      !!                  ***  ROUTINE p4z_mort_init  *** 
    … …  
    218219      !! 
    219220      !!---------------------------------------------------------------------- 
    220  
     221      INTEGER ::   ios   ! Local integer 
     222      ! 
    221223      NAMELIST/namp4zmort/ wchl, wchld, wchldm, mprat, mprat2 
    222       INTEGER :: ios                 ! Local integer output status for namelist read 
    223  
     224      !!---------------------------------------------------------------------- 
     225      ! 
    224226      REWIND( numnatp_ref )              ! Namelist nampismort in reference namelist : Pisces phytoplankton 
    225227      READ  ( numnatp_ref, namp4zmort, IOSTAT = ios, ERR = 901) 
    226228901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namp4zmort in reference namelist', lwp ) 
    227  
     229      ! 
    228230      REWIND( numnatp_cfg )              ! Namelist nampismort in configuration namelist : Pisces phytoplankton 
    229231      READ  ( numnatp_cfg, namp4zmort, IOSTAT = ios, ERR = 902 ) 
    230 902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namp4zmort in configuration namelist', lwp ) 
     232902   IF( ios > 0 ) CALL ctl_nam ( ios , 'namp4zmort in configuration namelist', lwp ) 
    231233      IF(lwm) WRITE ( numonp, namp4zmort ) 
    232  
     234      ! 
    233235      IF(lwp) THEN                         ! control print 
    234236         WRITE(numout,*) ' ' 
    … …  
    241243         WRITE(numout,*) '    Diatoms mortality rate                    mprat2    =', mprat2 
    242244      ENDIF 
    243  
     245      ! 
    244246   END SUBROUTINE p4z_mort_init 
    245247 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zopt.F90

    r7970 r9124  
    1818   USE prtctl_trc      !  print control for debugging 
    1919 
    20  
    2120   IMPLICIT NONE 
    2221   PRIVATE 
    … …  
    5958      !! ** Method  : - ??? 
    6059      !!--------------------------------------------------------------------- 
    61       ! 
    6260      INTEGER, INTENT(in) ::   kt, knt   ! ocean time step 
    6361      ! 
    … …  
    7270      !!--------------------------------------------------------------------- 
    7371      ! 
    74       IF( nn_timing == 1 )  CALL timing_start('p4z_opt') 
     72      IF( ln_timing )  CALL timing_start('p4z_opt') 
    7573      ! 
    7674      ! Allocate temporary workspace 
    … …  
    250248                   CALL wrk_dealloc( jpi, jpj, jpk, zpar   ,  ze0, ze1, ze2, ze3, zchl3d ) 
    251249      ! 
    252       IF( nn_timing == 1 )  CALL timing_stop('p4z_opt') 
     250      IF( ln_timing )  CALL timing_stop('p4z_opt') 
    253251      ! 
    254252   END SUBROUTINE p4z_opt 
     253 
    255254 
    256255   SUBROUTINE p4z_opt_par( kt, pqsr, pe1, pe2, pe3, pe0, pqsr100 )  
    … …  
    334333      !! 
    335334      !!---------------------------------------------------------------------- 
    336       !! * arguments 
    337335      INTEGER ,                INTENT(in) ::   kt     ! ocean time step 
    338  
    339       !! * local declarations 
     336      ! 
    340337      INTEGER  :: ji,jj 
    341338      REAL(wp) :: zcoef 
    342339      !!--------------------------------------------------------------------- 
    343340      ! 
    344       IF( nn_timing == 1 )  CALL timing_start('p4z_optsbc') 
     341      IF( ln_timing )  CALL timing_start('p4z_optsbc') 
    345342      ! 
    346343      ! Compute par_varsw at nit000 or only if there is more than 1 time record in par coefficient file 
    … …  
    352349      ENDIF 
    353350      ! 
    354       IF( nn_timing == 1 )  CALL timing_stop('p4z_optsbc') 
     351      IF( ln_timing )  CALL timing_stop('p4z_optsbc') 
    355352      ! 
    356353   END SUBROUTINE p4z_opt_sbc 
     354 
    357355 
    358356   SUBROUTINE p4z_opt_init 
    … …  
    365363      !! ** Input   :   external ascii and netcdf files 
    366364      !!---------------------------------------------------------------------- 
    367       ! 
    368365      INTEGER :: numpar 
    369366      INTEGER :: ierr 
    … …  
    375372      ! 
    376373      NAMELIST/nampisopt/cn_dir, sn_par, ln_varpar, parlux 
    377  
    378       !!---------------------------------------------------------------------- 
    379  
    380       IF( nn_timing == 1 )  CALL timing_start('p4z_opt_init') 
     374      !!---------------------------------------------------------------------- 
    381375 
    382376      REWIND( numnatp_ref )              ! Namelist nampisopt in reference namelist : Pisces attenuation coef. and PAR 
    … …  
    434428      IF( ln_qsr_bio )   etot3    (:,:,:) = 0._wp 
    435429      !  
    436       IF( nn_timing == 1 )  CALL timing_stop('p4z_opt_init') 
    437       ! 
    438430   END SUBROUTINE p4z_opt_init 
    439431 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zpoc.F90

    r7753 r9124  
    1717   USE prtctl_trc      !  print control for debugging 
    1818   USE iom             !  I/O manager 
    19  
    2019 
    2120   IMPLICIT NONE 
    … …  
    6968      !!--------------------------------------------------------------------- 
    7069      ! 
    71       IF( nn_timing == 1 )  CALL timing_start('p4z_poc') 
     70      IF( ln_timing )  CALL timing_start('p4z_poc') 
    7271      ! 
    7372      ! Allocate temporary workspace 
    … …  
    470469      DEALLOCATE( alphag ) 
    471470      ! 
    472       IF( nn_timing == 1 )  CALL timing_stop('p4z_poc') 
     471      IF( ln_timing )   CALL timing_stop('p4z_poc') 
    473472      ! 
    474473   END SUBROUTINE p4z_poc 
    … …  
    487486      !! 
    488487      !!---------------------------------------------------------------------- 
    489       INTEGER :: jn 
     488      INTEGER ::   ios, ifault   ! Local integer 
     489      INTEGER ::   jn 
    490490      REAL(wp) :: remindelta, reminup, remindown 
    491       INTEGER  :: ifault 
    492  
    493       NAMELIST/nampispoc/ xremip, jcpoc, rshape,  & 
    494          &     xremipc, xremipn, xremipp 
    495  
    496  
    497       INTEGER :: ios                 ! Local integer output status for namelist read 
     491      !! 
     492      NAMELIST/nampispoc/ xremip , jcpoc  , rshape,  & 
     493         &                xremipc, xremipn, xremipp 
     494      !!---------------------------------------------------------------------- 
    498495 
    499496      REWIND( numnatp_ref )              ! Namelist nampisrem in reference namelist : Pisces remineralization 
    … …  
    503500      REWIND( numnatp_cfg )              ! Namelist nampisrem in configuration namelist : Pisces remineralization 
    504501      READ  ( numnatp_cfg, nampispoc, IOSTAT = ios, ERR = 902 ) 
    505 902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampispoc in configuration namelist', lwp ) 
     502902   IF( ios > 0 ) CALL ctl_nam ( ios , 'nampispoc in configuration namelist', lwp ) 
    506503      IF(lwm) WRITE ( numonp, nampispoc ) 
    507504 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zprod.F90

    r7753 r9124  
    6565      !! ** Method  : - ??? 
    6666      !!--------------------------------------------------------------------- 
    67       ! 
    6867      INTEGER, INTENT(in) :: kt, knt 
    6968      ! 
    … …  
    8584      !!--------------------------------------------------------------------- 
    8685      ! 
    87       IF( nn_timing == 1 )  CALL timing_start('p4z_prod') 
     86      IF( ln_timing )   CALL timing_start('p4z_prod') 
    8887      ! 
    8988      !  Allocate temporary workspace 
    … …  
    486485     CALL wrk_dealloc( jpi, jpj, jpk, zprorcan, zprorcad, zprofed, zprofen, zpronewn, zpronewd ) 
    487486     ! 
    488      IF( nn_timing == 1 )  CALL timing_stop('p4z_prod') 
     487     IF( ln_timing )  CALL timing_stop('p4z_prod') 
    489488     ! 
    490489   END SUBROUTINE p4z_prod 
    … …  
    502501      !! ** input   :   Namelist nampisprod 
    503502      !!---------------------------------------------------------------------- 
     503      INTEGER :: ios                 ! Local integer output status for namelist read 
    504504      ! 
    505505      NAMELIST/namp4zprod/ pislopen, pisloped, xadap, ln_newprod, bresp, excretn, excretd,  & 
    506506         &                 chlcnm, chlcdm, chlcmin, fecnm, fecdm, grosip 
    507       INTEGER :: ios                 ! Local integer output status for namelist read 
    508507      !!---------------------------------------------------------------------- 
    509508 
    … …  
    514513      REWIND( numnatp_cfg )              ! Namelist nampisprod in configuration namelist : Pisces phytoplankton production 
    515514      READ  ( numnatp_cfg, namp4zprod, IOSTAT = ios, ERR = 902 ) 
    516 902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namp4zprod in configuration namelist', lwp ) 
     515902   IF( ios > 0 ) CALL ctl_nam ( ios , 'namp4zprod in configuration namelist', lwp ) 
    517516      IF(lwm) WRITE ( numonp, namp4zprod ) 
    518517 
    … …  
    555554      ! 
    556555   END FUNCTION p4z_prod_alloc 
     556 
    557557   !!====================================================================== 
    558558END MODULE p4zprod 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zrem.F90

    r8533 r9124  
    5858      !! ** Method  : - ??? 
    5959      !!--------------------------------------------------------------------- 
    60       ! 
    6160      INTEGER, INTENT(in) ::   kt, knt ! ocean time step 
    6261      ! 
    … …  
    7271      !!--------------------------------------------------------------------- 
    7372      ! 
    74       IF( nn_timing == 1 )  CALL timing_start('p4z_rem') 
     73      IF( ln_timing )   CALL timing_start('p4z_rem') 
    7574      ! 
    7675      ! Allocate temporary workspace 
    … …  
    289288      CALL wrk_dealloc( jpi, jpj, jpk, zdepbac, zdepprod, zolimi, zfacsi, zfacsib ) 
    290289      ! 
    291       IF( nn_timing == 1 )  CALL timing_stop('p4z_rem') 
     290      IF( ln_timing )   CALL timing_stop('p4z_rem') 
    292291      ! 
    293292   END SUBROUTINE p4z_rem 
    … …  
    309308         &                xremikc, xremikn, xremikp 
    310309      INTEGER :: ios                 ! Local integer output status for namelist read 
     310      !!---------------------------------------------------------------------- 
    311311 
    312312      REWIND( numnatp_ref )              ! Namelist nampisrem in reference namelist : Pisces remineralization 
    … …  
    316316      REWIND( numnatp_cfg )              ! Namelist nampisrem in configuration namelist : Pisces remineralization 
    317317      READ  ( numnatp_cfg, nampisrem, IOSTAT = ios, ERR = 902 ) 
    318 902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampisrem in configuration namelist', lwp ) 
     318902   IF( ios > 0 ) CALL ctl_nam ( ios , 'nampisrem in configuration namelist', lwp ) 
    319319      IF(lwm) WRITE ( numonp, nampisrem ) 
    320320 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsbc.F90

    r7753 r9124  
    108108      !!--------------------------------------------------------------------- 
    109109      ! 
    110       IF( nn_timing == 1 )  CALL timing_start('p4z_sbc') 
    111  
     110      IF( ln_timing )  CALL timing_start('p4z_sbc') 
    112111      ! 
    113112      ! Compute dust at nit000 or only if there is more than 1 time record in dust file 
    … …  
    189188      ENDIF 
    190189      ! 
    191       IF( nn_timing == 1 )  CALL timing_stop('p4z_sbc') 
     190      IF( ln_timing )  CALL timing_stop('p4z_sbc') 
    192191      ! 
    193192   END SUBROUTINE p4z_sbc 
    194193 
     194 
    195195   SUBROUTINE p4z_sbc_init 
    196  
    197196      !!---------------------------------------------------------------------- 
    198197      !!                  ***  routine p4z_sbc_init  *** 
    … …  
    232231        &                hratio, fep_rats, fep_rath, lgw_rath 
    233232      !!---------------------------------------------------------------------- 
    234       ! 
    235       IF( nn_timing == 1 )  CALL timing_start('p4z_sbc_init') 
    236233      ! 
    237234      !                            !* set file information 
    … …  
    535532      ENDIF 
    536533      ! 
    537       IF( nn_timing == 1 )  CALL timing_stop('p4z_sbc_init') 
    538       ! 
    539534   END SUBROUTINE p4z_sbc_init 
    540535 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsed.F90

    r8541 r9124  
    6868      REAL(wp), POINTER, DIMENSION(:,:,:) :: ztrpo4, ztrdop, zirondep, zsoufer, zpdep, zlight 
    6969      REAL(wp), POINTER, DIMENSION(:,:  ) :: zwsfep 
    70  
    7170      !!--------------------------------------------------------------------- 
    7271      ! 
    73       IF( nn_timing == 1 )  CALL timing_start('p4z_sed') 
     72      IF( ln_timing )  CALL timing_start('p4z_sed') 
    7473      ! 
    7574      IF( kt == nittrc000 .AND. knt == 1 )   r1_rday  = 1. / rday 
    … …  
    522521      IF( ln_ligand ) CALL wrk_dealloc( jpi, jpj, zwsfep ) 
    523522      ! 
    524       IF( nn_timing == 1 )  CALL timing_stop('p4z_sed') 
     523      IF( ln_timing )  CALL timing_stop('p4z_sed') 
    525524      ! 
    526525   END SUBROUTINE p4z_sed 
    … …  
    537536   END FUNCTION p4z_sed_alloc 
    538537 
    539  
    540538   !!====================================================================== 
    541539END MODULE p4zsed 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsink.F90

    r7753 r9124  
    6969      !!--------------------------------------------------------------------- 
    7070      ! 
    71       IF( nn_timing == 1 )  CALL timing_start('p4z_sink') 
     71      IF( ln_timing )   CALL timing_start('p4z_sink') 
    7272 
    7373 
    … …  
    261261      ENDIF 
    262262      ! 
    263       IF( nn_timing == 1 )  CALL timing_stop('p4z_sink') 
     263      IF( ln_timing )   CALL timing_stop('p4z_sink') 
    264264      ! 
    265265   END SUBROUTINE p4z_sink 
    266266 
     267 
    267268   SUBROUTINE p4z_sink_init 
    268269      !!---------------------------------------------------------------------- 
    … …  
    270271      !!---------------------------------------------------------------------- 
    271272      INTEGER :: jk 
    272  
     273      !!---------------------------------------------------------------------- 
     274      ! 
    273275      ik100 = 10        !  last level where depth less than 100 m 
    274276      DO jk = jpkm1, 1, -1 
    … …  
    282284      ! 
    283285   END SUBROUTINE p4z_sink_init 
     286 
    284287 
    285288   SUBROUTINE p4z_sink2( pwsink, psinkflx, jp_tra, kiter ) 
    … …  
    294297      !!      transport term, i.e.  div(u*tra). 
    295298      !!--------------------------------------------------------------------- 
    296       ! 
    297299      INTEGER , INTENT(in   )                         ::   jp_tra    ! tracer index index       
    298300      INTEGER , INTENT(in   )                         ::   kiter     ! number of iterations for time-splitting  
    299301      REAL(wp), INTENT(in   ), DIMENSION(jpi,jpj,jpk) ::   pwsink    ! sinking speed 
    300302      REAL(wp), INTENT(inout), DIMENSION(jpi,jpj,jpk) ::   psinkflx  ! sinking fluxe 
    301       !! 
     303      ! 
    302304      INTEGER  ::   ji, jj, jk, jn 
    303305      REAL(wp) ::   zigma,zew,zign, zflx, zstep 
    … …  
    305307      !!--------------------------------------------------------------------- 
    306308      ! 
    307       IF( nn_timing == 1 )  CALL timing_start('p4z_sink2') 
     309      IF( ln_timing )   CALL timing_start('p4z_sink2') 
    308310      ! 
    309311      ! Allocate temporary workspace 
    … …  
    391393      CALL wrk_dealloc( jpi, jpj, jpk, ztraz, zakz, zwsink2, ztrb ) 
    392394      ! 
    393       IF( nn_timing == 1 )  CALL timing_stop('p4z_sink2') 
     395      IF( ln_timing )  CALL timing_stop('p4z_sink2') 
    394396      ! 
    395397   END SUBROUTINE p4z_sink2 
    … …  
    401403      !!---------------------------------------------------------------------- 
    402404      INTEGER :: ierr(3) 
    403  
     405      !!---------------------------------------------------------------------- 
     406      ! 
    404407      ierr(:) = 0 
    405408      ! 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsms.F90

    r9081 r9124  
    4040   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   xnegtr     !: ??? 
    4141 
    42  
    4342   !!---------------------------------------------------------------------- 
    4443   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    … …  
    6867      !!--------------------------------------------------------------------- 
    6968      ! 
    70       IF( nn_timing == 1 )  CALL timing_start('p4z_sms') 
     69      IF( ln_timing )   CALL timing_start('p4z_sms') 
    7170      ! 
    7271      IF( kt == nittrc000 ) THEN 
    … …  
    171170      ! 
    172171 
    173       IF( lk_iomput .OR. ln_check_mass )  CALL p4z_chk_mass( kt ) ! Mass conservation checking 
    174  
    175       IF ( lwm .AND. kt == nittrc000 ) CALL FLUSH    ( numonp )     ! flush output namelist PISCES 
    176       IF( nn_timing == 1 )  CALL timing_stop('p4z_sms') 
    177       ! 
     172      IF( lk_iomput .OR. ln_check_mass )  CALL p4z_chk_mass( kt )    ! Mass conservation checking 
     173 
     174      IF( lwm .AND. kt == nittrc000    )  CALL FLUSH( numonp )       ! flush output namelist PISCES 
     175      ! 
     176      IF( ln_timing )  CALL timing_stop('p4z_sms') 
    178177      ! 
    179178   END SUBROUTINE p4z_sms 
     179 
    180180 
    181181   SUBROUTINE p4z_sms_init 
    … …  
    188188      !!             namelist: natext, natbio, natsms 
    189189      !!---------------------------------------------------------------------- 
     190      INTEGER :: ios                 ! Local integer output status for namelist read 
     191      !! 
    190192      NAMELIST/nampisbio/ nrdttrc, wsbio, xkmort, ferat3, wsbio2, wsbio2max, wsbio2scale,    & 
    191193         &                   niter1max, niter2max, wfep, ldocp, ldocz, lthet,  & 
    192194         &                   no3rat3, po4rat3 
    193  
     195         ! 
    194196      NAMELIST/nampisdmp/ ln_pisdmp, nn_pisdmp 
    195197      NAMELIST/nampismass/ ln_check_mass 
    196       INTEGER :: ios                 ! Local integer output status for namelist read 
    197198      !!---------------------------------------------------------------------- 
    198199 
    … …  
    203204      REWIND( numnatp_cfg )              ! Namelist nampisbio in configuration namelist : Pisces variables 
    204205      READ  ( numnatp_cfg, nampisbio, IOSTAT = ios, ERR = 902 ) 
    205 902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampisbio in configuration namelist', lwp ) 
     206902   IF( ios > 0 ) CALL ctl_nam ( ios , 'nampisbio in configuration namelist', lwp ) 
    206207      IF(lwm) WRITE ( numonp, nampisbio ) 
    207208 
    … …  
    238239      REWIND( numnatp_cfg )              ! Namelist nampisdmp in configuration namelist : Pisces damping 
    239240      READ  ( numnatp_cfg, nampisdmp, IOSTAT = ios, ERR = 906 ) 
    240 906   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampisdmp in configuration namelist', lwp ) 
     241906   IF( ios > 0 ) CALL ctl_nam ( ios , 'nampisdmp in configuration namelist', lwp ) 
    241242      IF(lwm) WRITE ( numonp, nampisdmp ) 
    242243 
    … …  
    255256      REWIND( numnatp_cfg )              ! Namelist nampismass in configuration namelist : Pisces mass conservation check  
    256257      READ  ( numnatp_cfg, nampismass, IOSTAT = ios, ERR = 908 ) 
    257 908   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampismass in configuration namelist', lwp ) 
     258908   IF( ios > 0 ) CALL ctl_nam ( ios , 'nampismass in configuration namelist', lwp ) 
    258259      IF(lwm) WRITE ( numonp, nampismass ) 
    259260 
    … …  
    264265         WRITE(numout,*) '    Flag to check mass conservation of NO3/Si/TALK ln_check_mass = ', ln_check_mass 
    265266      ENDIF 
    266  
     267      ! 
    267268   END SUBROUTINE p4z_sms_init 
     269 
    268270 
    269271   SUBROUTINE p4z_rst( kt, cdrw ) 
    … …  
    279281      INTEGER         , INTENT(in) ::   kt         ! ocean time-step 
    280282      CHARACTER(len=*), INTENT(in) ::   cdrw       ! "READ"/"WRITE" flag 
    281       ! 
    282       !!--------------------------------------------------------------------- 
    283  
     283      !!--------------------------------------------------------------------- 
     284      ! 
    284285      IF( TRIM(cdrw) == 'READ' ) THEN 
    285286         ! 
    … …  
    337338   END SUBROUTINE p4z_rst 
    338339 
     340 
    339341   SUBROUTINE p4z_dmp( kt ) 
    340342      !!---------------------------------------------------------------------- 
    … …  
    354356      REAL(wp) :: zalksumb, zpo4sumb, zno3sumb, zsilsumb 
    355357      !!--------------------------------------------------------------------- 
    356  
    357358 
    358359      IF(lwp)  WRITE(numout,*) 
    … …  
    415416      !! 
    416417      !!--------------------------------------------------------------------- 
    417       ! 
    418418      INTEGER, INTENT( in ) ::   kt      ! ocean time-step index       
    419419      REAL(wp)             ::  zrdenittot, zsdenittot, znitrpottot 
    … …  
    422422      REAL(wp), POINTER, DIMENSION(:,:,:) :: zwork 
    423423      !!---------------------------------------------------------------------- 
    424  
    425       ! 
    426       !!--------------------------------------------------------------------- 
    427  
     424      ! 
    428425      IF( kt == nittrc000 ) THEN  
    429426         xfact1 = rfact2r * 12. / 1.e15 * ryyss    ! conversion molC/kt --> PgC/yr 
    … …  
    535532         ! 
    536533         IF( lwp ) WRITE(numnit,9200) ndastp, znitrpottot * xfact2  , & 
    537          &                             zrdenittot  * xfact2  , & 
    538          &                             zsdenittot  * xfact2 
    539  
     534            &                             zrdenittot  * xfact2  , & 
     535            &                             zsdenittot  * xfact2 
    540536      ENDIF 
    541537      ! 
    … …  
    543539 9100  FORMAT(i8,5e18.10) 
    544540 9200  FORMAT(i8,3f10.5) 
    545  
    546541       ! 
    547542   END SUBROUTINE p4z_chk_mass 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p5zlim.F90

    r7646 r9124  
    154154      !!--------------------------------------------------------------------- 
    155155      ! 
    156       IF( nn_timing == 1 )  CALL timing_start('p5z_lim') 
     156      IF( ln_timing )   CALL timing_start('p5z_lim') 
    157157      ! 
    158158      zratchl = 6.0 
    … …  
    448448      ENDIF 
    449449      ! 
    450       IF( nn_timing == 1 )  CALL timing_stop('p5z_lim') 
     450      IF( ln_timing )  CALL timing_stop('p5z_lim') 
    451451      ! 
    452452   END SUBROUTINE p5z_lim 
    453453 
     454 
    454455   SUBROUTINE p5z_lim_init 
    455  
    456456      !!---------------------------------------------------------------------- 
    457457      !!                  ***  ROUTINE p5z_lim_init  *** 
    … …  
    465465      !! 
    466466      !!---------------------------------------------------------------------- 
    467  
     467      INTEGER :: ios                 ! Local integer output status for namelist read 
     468      !! 
    468469      NAMELIST/namp5zlim/ concnno3, concpno3, concdno3, concnnh4, concpnh4, concdnh4,  & 
    469470         &                concnfer, concpfer, concdfer, concbfe, concnpo4, concppo4,   & 
    … …  
    471472         &                xsizephy, xsizern, xsizerp, xsizerd, xksi1, xksi2, xkdoc,    & 
    472473         &                caco3r, oxymin 
    473  
     474         ! 
    474475      NAMELIST/namp5zquota/ qnnmin, qnnmax, qpnmin, qpnmax, qnpmin, qnpmax, qppmin,      & 
    475476         &                  qppmax, qndmin, qndmax, qpdmin, qpdmax, qfnmax, qfpmax, qfdmax,  & 
    476477         &                  qfnopt, qfpopt, qfdopt 
    477   
    478  
    479       INTEGER :: ios                 ! Local integer output status for namelist read 
    480  
     478      !!---------------------------------------------------------------------- 
     479      ! 
    481480      REWIND( numnatp_ref )              ! Namelist nampislim in reference namelist : Pisces nutrient limitation parameters 
    482481      READ  ( numnatp_ref, namp5zlim, IOSTAT = ios, ERR = 901) 
    483482901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampislim in reference namelist', lwp ) 
    484  
     483      ! 
    485484      REWIND( numnatp_cfg )              ! Namelist nampislim in configuration namelist : Pisces nutrient limitation parameters  
    486485      READ  ( numnatp_cfg, namp5zlim, IOSTAT = ios, ERR = 902 ) 
    487 902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampislim in configuration namelist', lwp ) 
     486902   IF( ios > 0 ) CALL ctl_nam ( ios , 'nampislim in configuration namelist', lwp ) 
    488487      IF(lwm) WRITE ( numonp, namp5zlim ) 
    489  
     488      ! 
    490489      IF(lwp) THEN                         ! control print 
    491490         WRITE(numout,*) ' ' 
    … …  
    523522      READ  ( numnatp_ref, namp5zquota, IOSTAT = ios, ERR = 903) 
    524523903   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampisquota in reference namelist', lwp ) 
    525  
     524      ! 
    526525      REWIND( numnatp_cfg )              ! Namelist nampislim in configuration namelist : Pisces nutrient limitation parameters  
    527526      READ  ( numnatp_cfg, namp5zquota, IOSTAT = ios, ERR = 904 ) 
    528 904   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampisquota in configuration namelist', lwp ) 
     527904   IF( ios > 0 ) CALL ctl_nam ( ios , 'nampisquota in configuration namelist', lwp ) 
    529528      IF(lwm) WRITE ( numonp, namp5zquota ) 
    530  
     529      ! 
    531530      IF(lwp) THEN                         ! control print 
    532531         WRITE(numout,*) ' ' 
    … …  
    552551         WRITE(numout,*) '    Maximal Fe quota for diatoms             qfdmax    = ', qfdmax 
    553552      ENDIF 
    554  
     553      ! 
    555554      zpsino3 = 2.3 * rno3 
    556555      zpsinh4 = 1.8 * rno3 
    … …  
    560559      ! 
    561560   END SUBROUTINE p5z_lim_init 
     561 
    562562 
    563563   INTEGER FUNCTION p5z_lim_alloc() 
    … …  
    569569      !!---------------------------------------------------------------------- 
    570570      ierr(:) = 0 
    571  
     571      ! 
    572572      !*  Biological arrays for phytoplankton growth 
    573573      ALLOCATE( xpicono3(jpi,jpj,jpk), xpiconh4(jpi,jpj,jpk),       & 
    … …  
    585585         &      fvpuptk (jpi,jpj,jpk), xlimpic (jpi,jpj,jpk),       & 
    586586         &      xlimsi  (jpi,jpj,jpk), STAT=ierr(1) ) 
    587  
     587         ! 
    588588      !*  Minimum/maximum quotas of phytoplankton 
    589589      ALLOCATE( xqnnmin (jpi,jpj,jpk), xqnnmax(jpi,jpj,jpk),       & 
    … …  
    593593         &      xqndmin (jpi,jpj,jpk), xqndmax(jpi,jpj,jpk),       & 
    594594         &      xqpdmin (jpi,jpj,jpk), xqpdmax(jpi,jpj,jpk),     STAT=ierr(2) ) 
    595       ! 
     595         ! 
    596596      p5z_lim_alloc = MAXVAL( ierr ) 
    597597      ! 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p5zmeso.F90

    r7646 r9124  
    8989      !!--------------------------------------------------------------------- 
    9090      ! 
    91       IF( nn_timing == 1 )  CALL timing_start('p5z_meso') 
     91      IF( ln_timing )   CALL timing_start('p5z_meso') 
    9292      ! 
    9393      CALL wrk_alloc( jpi, jpj, jpk, zgrazing ) 
    … …  
    362362      CALL wrk_dealloc( jpi, jpj, jpk, zgrazing ) 
    363363      ! 
    364       IF( nn_timing == 1 )  CALL timing_stop('p5z_meso') 
     364      IF( ln_timing )   CALL timing_stop('p5z_meso') 
    365365      ! 
    366366   END SUBROUTINE p5z_meso 
    367367 
     368 
    368369   SUBROUTINE p5z_meso_init 
    369  
    370370      !!---------------------------------------------------------------------- 
    371371      !!                  ***  ROUTINE p5z_meso_init  *** 
    … …  
    379379      !! 
    380380      !!---------------------------------------------------------------------- 
    381  
     381      INTEGER :: ios                 ! Local integer output status for namelist read 
     382      !! 
    382383      NAMELIST/namp5zmes/part2, bmetexc2, grazrat2, resrat2, mzrat2, xpref2c, xpref2p, xpref2z, & 
    383384         &                xpref2m, xpref2d, xthresh2dia, xthresh2phy, xthresh2zoo, xthresh2poc, & 
    384385         &                xthresh2mes, xthresh2, xkgraz2, epsher2, ssigma2, unass2c, & 
    385       &                   unass2n, unass2p, srespir2, grazflux 
    386  
    387       INTEGER :: ios                 ! Local integer output status for namelist read 
    388  
     386         &                unass2n, unass2p, srespir2, grazflux 
     387      !!---------------------------------------------------------------------- 
     388      ! 
    389389      REWIND( numnatp_ref )              ! Namelist nampismes in reference namelist : Pisces mesozooplankton 
    390390      READ  ( numnatp_ref, namp5zmes, IOSTAT = ios, ERR = 901) 
    391391901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampismes in reference namelist', lwp ) 
    392  
     392      ! 
    393393      REWIND( numnatp_cfg )              ! Namelist nampismes in configuration namelist : Pisces mesozooplankton 
    394394      READ  ( numnatp_cfg, namp5zmes, IOSTAT = ios, ERR = 902 ) 
    395 902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampismes in configuration namelist', lwp ) 
     395902   IF( ios > 0 ) CALL ctl_nam ( ios , 'nampismes in configuration namelist', lwp ) 
    396396      IF(lwm) WRITE ( numonp, namp5zmes ) 
    397  
    398  
     397      ! 
    399398      IF(lwp) THEN                         ! control print 
    400399         WRITE(numout,*) ' '  
    … …  
    426425         WRITE(numout,*) '    Use excess carbon for respiration              bmetexc2    = ', bmetexc2 
    427426      ENDIF 
    428  
    429  
     427      ! 
    430428   END SUBROUTINE p5z_meso_init 
    431429 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p5zmicro.F90

    r7646 r9124  
    8686      !!--------------------------------------------------------------------- 
    8787      ! 
    88       IF( nn_timing == 1 )  CALL timing_start('p5z_micro') 
     88      IF( ln_timing )   CALL timing_start('p5z_micro') 
    8989      ! 
    9090      CALL wrk_alloc( jpi, jpj, jpk, zgrazing ) 
    … …  
    306306      CALL wrk_dealloc( jpi, jpj, jpk, zgrazing ) 
    307307      ! 
    308       IF( nn_timing == 1 )  CALL timing_stop('p5z_micro') 
     308      IF( ln_timing )   CALL timing_stop('p5z_micro') 
    309309      ! 
    310310   END SUBROUTINE p5z_micro 
    … …  
    312312 
    313313   SUBROUTINE p5z_micro_init 
    314  
    315314      !!---------------------------------------------------------------------- 
    316315      !!                  ***  ROUTINE p5z_micro_init  *** 
    … …  
    324323      !! 
    325324      !!---------------------------------------------------------------------- 
    326  
     325      INTEGER ::   ios   ! Local integer 
     326      !! 
    327327      NAMELIST/namp5zzoo/ part, grazrat, bmetexc, resrat, mzrat, xprefc, xprefn, & 
    328328         &                xprefp, xprefd, xprefz, xthreshdia, xthreshphy, & 
    329329         &                xthreshpic, xthreshpoc, xthreshzoo, xthresh, xkgraz, & 
    330330         &                epsher, ssigma, srespir, unassc, unassn, unassp 
    331  
    332       INTEGER :: ios                 ! Local integer output status for namelist read 
    333  
     331      !!---------------------------------------------------------------------- 
     332      ! 
    334333      REWIND( numnatp_ref )              ! Namelist nampiszoo in reference namelist : Pisces microzooplankton 
    335334      READ  ( numnatp_ref, namp5zzoo, IOSTAT = ios, ERR = 901) 
    336335901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namp5zzoo in reference namelist', lwp ) 
    337  
     336      ! 
    338337      REWIND( numnatp_cfg )              ! Namelist nampiszoo in configuration namelist : Pisces microzooplankton 
    339338      READ  ( numnatp_cfg, namp5zzoo, IOSTAT = ios, ERR = 902 ) 
    340 902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namp5zzoo in configuration namelist', lwp ) 
     339902   IF( ios > 0 ) CALL ctl_nam ( ios , 'namp5zzoo in configuration namelist', lwp ) 
    341340      IF(lwm) WRITE ( numonp, namp5zzoo ) 
    342  
     341      ! 
    343342      IF(lwp) THEN                         ! control print 
    344343         WRITE(numout,*) ' ' 
    … …  
    369368         WRITE(numout,*) '    Use of excess carbon for respiration            bmetexc     =', bmetexc 
    370369      ENDIF 
    371  
     370      ! 
    372371   END SUBROUTINE p5z_micro_init 
    373372 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p5zmort.F90

    r7646 r9124  
    7474      !!--------------------------------------------------------------------- 
    7575      ! 
    76       IF( nn_timing == 1 )  CALL timing_start('p5z_nano') 
     76      IF( ln_timing )   CALL timing_start('p5z_nano') 
    7777      ! 
    7878      prodcal(:,:,:) = 0.  !: calcite production variable set to zero 
    … …  
    124124       ENDIF 
    125125      ! 
    126       IF( nn_timing == 1 )  CALL timing_stop('p5z_nano') 
     126      IF( ln_timing )   CALL timing_stop('p5z_nano') 
    127127      ! 
    128128   END SUBROUTINE p5z_nano 
     129 
    129130 
    130131   SUBROUTINE p5z_pico 
    … …  
    143144      !!--------------------------------------------------------------------- 
    144145      ! 
    145       IF( nn_timing == 1 )  CALL timing_start('p5z_pico') 
     146      IF( ln_timing )   CALL timing_start('p5z_pico') 
    146147      ! 
    147148      DO jk = 1, jpkm1 
    … …  
    184185       ENDIF 
    185186      ! 
    186       IF( nn_timing == 1 )  CALL timing_stop('p5z_pico') 
     187      IF( ln_timing )   CALL timing_stop('p5z_pico') 
    187188      ! 
    188189   END SUBROUTINE p5z_pico 
    … …  
    204205      !!--------------------------------------------------------------------- 
    205206      ! 
    206       IF( nn_timing == 1 )  CALL timing_start('p5z_diat') 
     207      IF( ln_timing )   CALL timing_start('p5z_diat') 
    207208      ! 
    208209 
    … …  
    262263      ENDIF 
    263264      ! 
    264       IF( nn_timing == 1 )  CALL timing_stop('p5z_diat') 
     265      IF( ln_timing )   CALL timing_stop('p5z_diat') 
    265266      ! 
    266267   END SUBROUTINE p5z_diat 
    267268 
     269 
    268270   SUBROUTINE p5z_mort_init 
    269  
    270271      !!---------------------------------------------------------------------- 
    271272      !!                  ***  ROUTINE p5z_mort_init  *** 
    … …  
    279280      !! 
    280281      !!---------------------------------------------------------------------- 
    281  
     282      INTEGER :: ios                 ! Local integer output status for namelist read 
     283      !! 
    282284      NAMELIST/namp5zmort/ wchln, wchlp, wchld, wchldm, mpratn, mpratp, mpratd 
    283       INTEGER :: ios                 ! Local integer output status for namelist read 
     285      !!---------------------------------------------------------------------- 
    284286 
    285287      REWIND( numnatp_ref )              ! Namelist nampismort in reference namelist : Pisces phytoplankton 
    … …  
    289291      REWIND( numnatp_cfg )              ! Namelist nampismort in configuration namelist : Pisces phytoplankton 
    290292      READ  ( numnatp_cfg, namp5zmort, IOSTAT = ios, ERR = 902 ) 
    291 902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namp5zmort in configuration namelist', lwp ) 
     293902   IF( ios > 0 ) CALL ctl_nam ( ios , 'namp5zmort in configuration namelist', lwp ) 
    292294      IF(lwm) WRITE ( numonp, namp5zmort ) 
    293295 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p5zprod.F90

    r7646 r9124  
    9898      !!--------------------------------------------------------------------- 
    9999      ! 
    100       IF( nn_timing == 1 )  CALL timing_start('p5z_prod') 
     100      IF( ln_timing )   CALL timing_start('p5z_prod') 
    101101      ! 
    102102      !  Allocate temporary workspace 
    … …  
    571571      CALL wrk_dealloc( jpi, jpj, jpk, zcroissp, zcroissn, zcroissd ) 
    572572      ! 
    573       IF( nn_timing == 1 )  CALL timing_stop('p5z_prod') 
     573      IF( ln_timing )   CALL timing_stop('p5z_prod') 
    574574      ! 
    575575   END SUBROUTINE p5z_prod 
    … …  
    587587      !! ** input   :   Namelist nampisprod 
    588588      !!---------------------------------------------------------------------- 
    589       ! 
     589      INTEGER :: ios                 ! Local integer output status for namelist read 
     590      !! 
    590591      NAMELIST/namp5zprod/ pislopen, pislopep, pisloped, excretn, excretp, excretd,     & 
    591592         &                 thetannm, thetanpm, thetandm, chlcmin, grosip, bresp, xadap 
    592  
    593       INTEGER :: ios                 ! Local integer output status for namelist read 
    594593      !!---------------------------------------------------------------------- 
    595594 
    … …  
    600599      REWIND( numnatp_cfg )              ! Namelist nampisprod in configuration namelist : Pisces phytoplankton production 
    601600      READ  ( numnatp_cfg, namp5zprod, IOSTAT = ios, ERR = 902 ) 
    602 902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namp5zprod in configuration namelist', lwp ) 
     601902   IF( ios > 0 ) CALL ctl_nam ( ios , 'namp5zprod in configuration namelist', lwp ) 
    603602      IF(lwm) WRITE ( numonp, namp5zprod ) 
    604603 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/SED/par_sed.F90

    r7646 r9124  
    77   !!        !  06-12  (C. Ethe)  Orignal 
    88   !!---------------------------------------------------------------------- 
    9    !! $Id$ 
    109#if defined key_sed 
     10   !!---------------------------------------------------------------------- 
     11   !!   'key_sed'                                                  Sediment  
     12   !!---------------------------------------------------------------------- 
     13 
    1114   !! Domain characteristics 
    1215   USE par_kind 
    … …  
    3740      jpwat = 7,         &   !: number of pore water component 
    3841      jpwatp1 = jpwat + 1 
    39  
    4042    
    4143   ! pore water components        
    … …  
    6163      jpdia2dsed = 7 
    6264#endif 
     65 
     66   !!---------------------------------------------------------------------- 
     67   !! NEMO/TOP 4.0 , NEMO Consortium (2017) 
     68   !! $Id$  
     69   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
     70   !!====================================================================== 
    6371END MODULE par_sed 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/SED/sed.F90

    r7646 r9124  
    44   !! Sediment :   set sediment global variables 
    55   !!====================================================================== 
    6 #if defined key_sed 
    76   !! History : 
    87   !!        !  06-12  (C. Ethe)  Orignal 
     8   !!---------------------------------------------------------------------- 
     9#if defined key_sed 
     10   !!---------------------------------------------------------------------- 
     11   !!   'key_sed'                                                  Sediment  
    912   !!---------------------------------------------------------------------- 
    1013   USE par_sed 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/SED/sedadv.F90

    r5215 r9124  
    11MODULE sedadv 
    2 #if defined key_sed 
    32   !!====================================================================== 
    43   !!              ***  MODULE  sedadv  *** 
    54   !!    Sediment : vertical advection and burial 
    65   !!===================================================================== 
    7    !! * Modules used 
     6#if defined key_sed 
     7   !!---------------------------------------------------------------------- 
     8   !!   'key_sed'                                                  Sediment  
    89   !!---------------------------------------------------------------------- 
    910   !!   sed_adv : 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/par_pisces.F90

    r7646 r9124  
    55   !!====================================================================== 
    66   !! History :   2.0  !  2007-12  (C. Ethe, G. Madec)  revised architecture 
    7    !!---------------------------------------------------------------------- 
    8    !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
    9    !! $Id$  
    10    !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    117   !!---------------------------------------------------------------------- 
    128 
    … …  
    6864   INTEGER, PUBLIC  ::   jp_pcs1  !: Last  index of PISCES tracers 
    6965 
     66   !!---------------------------------------------------------------------- 
     67   !! NEMO/TOP 3.3 , NEMO Consortium (2010) 
     68   !! $Id$  
     69   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    7070   !!====================================================================== 
    7171END MODULE par_pisces 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/sms_pisces.F90

    r8533 r9124  
    137137CONTAINS 
    138138 
    139  
    140139   INTEGER FUNCTION sms_pisces_alloc() 
    141140      !!---------------------------------------------------------------------- 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/trcice_pisces.F90

    r7646 r9124  
    2020   PUBLIC   trc_ice_ini_pisces ! called by trcini.F90 module 
    2121 
     22   !!---------------------------------------------------------------------- 
     23   !! NEMO/TOP 4.0 , NEMO Consortium (2017) 
     24   !! $Id: trcsub.F90 9019 2017-12-13 14:58:53Z timgraham $  
     25   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
     26   !!---------------------------------------------------------------------- 
    2227CONTAINS 
    2328 
    24  
    2529   SUBROUTINE trc_ice_ini_pisces 
    2630      !!---------------------------------------------------------------------- 
    … …  
    2933      !! ** Purpose :   Initialisation of the PISCES biochemical model 
    3034      !!---------------------------------------------------------------------- 
    31  
     35      ! 
    3236      IF( ln_p4z .OR. ln_p5z ) THEN  ;   CALL p4z_ice_ini   !  PISCES 
    3337      ELSE                           ;   CALL p2z_ice_ini   !  LOBSTER 
    3438      ENDIF 
    35  
     39      ! 
    3640   END SUBROUTINE trc_ice_ini_pisces 
    3741 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/trcnam_pisces.F90

    r7646 r9124  
    1818   USE iom             ! I/O manager 
    1919 
    20  
    2120   IMPLICIT NONE 
    2221   PRIVATE 
    2322 
    2423   PUBLIC   trc_nam_pisces   ! called by trcnam.F90 module 
    25  
    2624 
    2725   !!---------------------------------------------------------------------- 
    … …  
    3028   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    3129   !!---------------------------------------------------------------------- 
    32  
    3330CONTAINS 
    3431 
    … …  
    4239      !!             namelist: natext, natbio, natsms 
    4340      !!---------------------------------------------------------------------- 
    44       !! 
    4541      INTEGER :: jl, jn 
    4642      INTEGER :: ios, ioptio                 ! Local integer output status for namelist read 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/trcsms_pisces.F90

    r7646 r9124  
    2323   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 
    2424   !!---------------------------------------------------------------------- 
    25  
    2625CONTAINS 
    27  
    28       !!---------------------------------------------------------------------- 
    29       !!                   ***  ROUTINE trc_ini_pisces *** 
    30       !! 
    31       !! ** Purpose :   Initialisation of the PISCES biochemical model 
    32       !!---------------------------------------------------------------------- 
    33  
    3426 
    3527   SUBROUTINE trc_sms_pisces( kt ) 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/TRP/trcbbl.F90

    r9019 r9124  
    5050      !!---------------------------------------------------------------------- 
    5151      ! 
    52       IF( nn_timing == 1 )  CALL timing_start('trc_bbl') 
     52      IF( ln_timing )   CALL timing_start('trc_bbl') 
    5353      ! 
    5454      IF( .NOT. l_offline .AND. nn_dttrc == 1 ) THEN 
    … …  
    9292      ENDIF 
    9393      ! 
    94       IF( nn_timing == 1 ) CALL timing_stop('trc_bbl') 
     94      IF( ln_timing )  CALL timing_stop('trc_bbl') 
    9595      ! 
    9696   END SUBROUTINE trc_bbl 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/TRP/trcdmp.F90

    r9019 r9124  
    8989      !!---------------------------------------------------------------------- 
    9090      ! 
    91       IF( nn_timing == 1 )  CALL timing_start('trc_dmp') 
     91      IF( ln_timing )   CALL timing_start('trc_dmp') 
    9292      ! 
    9393      IF( l_trdtrc )   CALL wrk_alloc( jpi, jpj, jpk, ztrtrd )   ! temporary save of trends 
    … …  
    161161      ENDIF 
    162162      ! 
    163       IF( nn_timing == 1 )  CALL timing_stop('trc_dmp') 
     163      IF( ln_timing )   CALL timing_stop('trc_dmp') 
    164164      ! 
    165165   END SUBROUTINE trc_dmp 
    … …  
    179179      NAMELIST/namtrc_dmp/ nn_zdmp_tr , cn_resto_tr 
    180180      !!---------------------------------------------------------------------- 
    181       ! 
    182       IF( nn_timing == 1 )  CALL timing_start('trc_dmp_ini') 
    183181      ! 
    184182      REWIND( numnat_ref )              ! Namelist namtrc_dmp in reference namelist : Passive tracers newtonian damping 
    … …  
    222220         ! 
    223221      ENDIF 
    224       IF( nn_timing == 1 )  CALL timing_stop('trc_dmp_ini') 
    225222      ! 
    226223   END SUBROUTINE trc_dmp_ini 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/TRP/trcldf.F90

    r9019 r9124  
    7474      !!---------------------------------------------------------------------- 
    7575      ! 
    76       IF( nn_timing == 1 )   CALL timing_start('trc_ldf') 
     76      IF( ln_timing )   CALL timing_start('trc_ldf') 
    7777      ! 
    7878      IF( l_trdtrc )  THEN 
    … …  
    124124      CALL wrk_dealloc( jpi,jpj,jpk,   zahu, zahv ) 
    125125      ! 
    126       IF( nn_timing == 1 )   CALL timing_stop('trc_ldf') 
     126      IF( ln_timing )   CALL timing_stop('trc_ldf') 
    127127      ! 
    128128   END SUBROUTINE trc_ldf 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/TRP/trcnxt.F90

    r9081 r9124  
    8686      !!---------------------------------------------------------------------- 
    8787      ! 
    88       IF( nn_timing == 1 )  CALL timing_start('trc_nxt') 
     88      IF( ln_timing )   CALL timing_start('trc_nxt') 
    8989      ! 
    9090      IF( kt == nittrc000 .AND. lwp ) THEN 
    … …  
    143143      ENDIF 
    144144      ! 
    145       IF( nn_timing == 1 )  CALL timing_stop('trc_nxt') 
     145      IF( ln_timing )   CALL timing_stop('trc_nxt') 
    146146      ! 
    147147   END SUBROUTINE trc_nxt 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/TRP/trcrad.F90

    r7753 r9124  
    5454      !!---------------------------------------------------------------------- 
    5555      ! 
    56       IF( nn_timing == 1 )  CALL timing_start('trc_rad') 
     56      IF( ln_timing )   CALL timing_start('trc_rad') 
    5757      ! 
    5858      IF( kt == nittrc000 ) THEN 
    … …  
    7575      ENDIF 
    7676      ! 
    77       IF( nn_timing == 1 )  CALL timing_stop('trc_rad') 
     77      IF( ln_timing )   CALL timing_stop('trc_rad') 
    7878      ! 
    7979   END SUBROUTINE trc_rad 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/TRP/trcsbc.F90

    r9019 r9124  
    6868      !!--------------------------------------------------------------------- 
    6969      ! 
    70       IF( nn_timing == 1 )  CALL timing_start('trc_sbc') 
     70      IF( ln_timing )   CALL timing_start('trc_sbc') 
    7171      ! 
    7272      ! Allocate temporary workspace 
    … …  
    182182      IF( l_trdtrc )  CALL wrk_dealloc( jpi,jpj,jpk,   ztrtrd ) 
    183183      ! 
    184       IF( nn_timing == 1 )  CALL timing_stop('trc_sbc') 
     184      IF( ln_timing )   CALL timing_stop('trc_sbc') 
    185185      ! 
    186186   END SUBROUTINE trc_sbc 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/TRP/trctrp.F90

    r9033 r9124  
    5656      !! --------------------------------------------------------------------- 
    5757      ! 
    58       IF( nn_timing == 1 )   CALL timing_start('trc_trp') 
     58      IF( ln_timing )   CALL timing_start('trc_trp') 
    5959      ! 
    6060      IF( .NOT. lk_c1d ) THEN 
    … …  
    9191      END IF 
    9292      ! 
    93       IF( nn_timing == 1 )   CALL timing_stop('trc_trp') 
     93      IF( ln_timing )   CALL timing_stop('trc_trp') 
    9494      ! 
    9595   END SUBROUTINE trc_trp 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/TRP/trczdf.F90

    r9019 r9124  
    5050      !!--------------------------------------------------------------------- 
    5151      ! 
    52       IF( nn_timing == 1 )  CALL timing_start('trc_zdf') 
     52      IF( ln_timing )   CALL timing_start('trc_zdf') 
    5353      ! 
    5454      IF( l_trdtrc )   ztrtrd(:,:,:,:)  = tra(:,:,:,:) 
    … …  
    7171      END IF 
    7272      ! 
    73       IF( nn_timing == 1 )  CALL timing_stop('trc_zdf') 
     73      IF( ln_timing )  CALL timing_stop('trc_zdf') 
    7474      ! 
    7575   END SUBROUTINE trc_zdf 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/trcbc.F90

    r9119 r9124  
    6060      !!              - allocates passive tracer BC data structure  
    6161      !!---------------------------------------------------------------------- 
    62       ! 
    6362      INTEGER,INTENT(IN) :: ntrc                           ! number of tracers 
    6463      INTEGER            :: jl, jn , ib, ibd, ii, ij, ik   ! dummy loop indices 
    … …  
    8180                        & cn_dir_cbc, sn_trccbc, rn_trcfac, ln_rnf_ctl, rn_bc_time 
    8281      NAMELIST/namtrc_bdy/ cn_trc_dflt, cn_trc, nn_trcdmp_bdy 
    83  
    8482      !!---------------------------------------------------------------------- 
    85       IF( nn_timing == 1 )  CALL timing_start('trc_bc_ini') 
    8683      ! 
    8784      IF( lwp ) THEN 
    … …  
    332329      ! 
    333330      DEALLOCATE( slf_i )          ! deallocate local field structure 
    334       IF( nn_timing == 1 )  CALL timing_stop('trc_bc_ini') 
    335331      ! 
    336332   END SUBROUTINE trc_bc_ini 
    … …  
    348344      !!---------------------------------------------------------------------- 
    349345      USE fldread 
    350        
    351       !! * Arguments 
     346      !!       
    352347      INTEGER, INTENT( in ) ::   kt              ! ocean time-step index 
    353348      INTEGER, INTENT( in ), OPTIONAL ::   jit   ! subcycle time-step index (for timesplitting option) 
    … …  
    357352      !!--------------------------------------------------------------------- 
    358353      ! 
    359       IF( nn_timing == 1 )  CALL timing_start('trc_bc') 
     354      IF( ln_timing )   CALL timing_start('trc_bc') 
    360355 
    361356      IF( kt == nit000 .AND. lwp) THEN 
    … …  
    452447      END DO                                                     ! tracer loop 
    453448      !                                                          ! =========== 
    454       ! 
    455       IF( nn_timing == 1 )  CALL timing_stop('trc_bc') 
     449      IF( ln_timing )   CALL timing_stop('trc_bc') 
    456450      ! 
    457451   END SUBROUTINE trc_bc 
    … …  
    462456   !!---------------------------------------------------------------------- 
    463457CONTAINS 
    464  
    465458   SUBROUTINE trc_bc_ini( ntrc )        ! Empty routine 
    466459      INTEGER,INTENT(IN) :: ntrc                           ! number of tracers 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/trcbdy.F90

    r7646 r9124  
    5151      !!---------------------------------------------------------------------- 
    5252      ! 
    53       IF( nn_timing == 1 ) CALL timing_start('trc_bdy') 
     53      IF( ln_timing )  CALL timing_start('trc_bdy') 
    5454      ! 
    5555      igrd = 1  
    … …  
    7676      END DO 
    7777      ! 
    78       IF( nn_timing == 1 ) CALL timing_stop('trc_bdy') 
     78      IF( ln_timing )  CALL timing_stop('trc_bdy') 
    7979 
    8080   END SUBROUTINE trc_bdy 
     81 
    8182 
    8283   SUBROUTINE trc_bdy_dmp( kt ) 
    … …  
    9899      !!---------------------------------------------------------------------- 
    99100      ! 
    100       IF( nn_timing == 1 ) CALL timing_start('trc_bdy_dmp') 
     101      IF( ln_timing )  CALL timing_start('trc_bdy_dmp') 
    101102      ! 
    102103      DO jn = 1, jptra 
    103104         DO ib_bdy=1, nb_bdy 
    104             IF ( trcdta_bdy(jn, ib_bdy)%dmp ) THEN 
     105            IF( trcdta_bdy(jn, ib_bdy)%dmp ) THEN 
    105106               igrd = 1                       ! Everything is at T-points here 
    106107               DO ib = 1, idx_bdy(ib_bdy)%nblen(igrd) 
    … …  
    114115               END DO 
    115116            ENDIF 
    116          ENDDO 
    117       ENDDO 
     117         END DO 
     118      END DO 
    118119      ! 
    119       IF( nn_timing == 1 ) CALL timing_stop('trc_bdy_dmp') 
     120      IF( ln_timing )  CALL timing_stop('trc_bdy_dmp') 
    120121      ! 
    121122   END SUBROUTINE trc_bdy_dmp 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/trcdta.F90

    r8542 r9124  
    2020   USE oce_trc       !  shared variables between ocean and passive tracers 
    2121   USE trc           !  passive tracers common variables 
     22   ! 
    2223   USE iom           !  I/O manager 
    2324   USE lib_mpp       !  MPP library 
    … …  
    6970      NAMELIST/namtrc_dta/ sn_trcdta, cn_dir, rn_trfac  
    7071      !!---------------------------------------------------------------------- 
    71       ! 
    72       IF( nn_timing == 1 )  CALL timing_start('trc_dta_ini') 
    7372      ! 
    7473      IF( lwp ) THEN 
    … …  
    154153      ! 
    155154      DEALLOCATE( slf_i )          ! deallocate local field structure 
    156       IF( nn_timing == 1 )  CALL timing_stop('trc_dta_ini') 
    157155      ! 
    158156   END SUBROUTINE trc_dta_ini 
    … …  
    182180      !!---------------------------------------------------------------------- 
    183181      ! 
    184       IF( nn_timing == 1 )  CALL timing_start('trc_dta') 
     182      IF( ln_timing )   CALL timing_start('trc_dta') 
    185183      ! 
    186184      IF( nb_trcdta > 0 ) THEN 
    … …  
    246244      ENDIF 
    247245      ! 
    248       IF( nn_timing == 1 )  CALL timing_stop('trc_dta') 
     246      IF( ln_timing )  CALL timing_stop('trc_dta') 
    249247      ! 
    250248   END SUBROUTINE trc_dta 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/trcice.F90

    r7753 r9124  
    1212   !!   trc_ice   :  Call the appropriate sea ice tracer subroutine 
    1313   !!---------------------------------------------------------------------- 
    14  
    1514   USE oce_trc         ! shared variables between ocean and passive tracers 
    1615   USE trc             ! passive tracers common variables 
    … …  
    2423   PRIVATE 
    2524    
    26    PUBLIC   trc_ice_ini ! called by trc_nam 
     25   PUBLIC   trc_ice_ini   ! called by trc_nam 
    2726 
     27   !!---------------------------------------------------------------------- 
     28   !! NEMO/TOP 4.0 , NEMO Consortium (2011) 
     29   !! $Id: trcini.F90 9019 2017-12-13 14:58:53Z timgraham $  
     30   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt) 
     31   !!---------------------------------------------------------------------- 
    2832CONTAINS 
    2933    
    … …  
    3741      !!             
    3842      !!--------------------------------------------------------------------- 
    39       ! --- Variable declarations --- ! 
    40  
     43      ! 
    4144      IF(lwp) THEN 
    4245         WRITE(numout,*) 
    … …  
    4447         WRITE(numout,*) '~~~~~~~~~~~~~' 
    4548      ENDIF 
    46  
    47       IF( nn_timing == 1 )  CALL timing_start('trc_ice_ini') 
    4849      ! 
    4950      CALL trc_nam_ice 
    5051      ! 
    51       trc_i(:,:,:) = 0.0d0 ! by default 
    52       trc_o(:,:,:) = 0.0d0 ! by default 
    53  
     52      trc_i(:,:,:) = 0._wp  ! by default 
     53      trc_o(:,:,:) = 0._wp  ! by default 
     54      ! 
    5455      IF ( nn_ice_tr == 1 ) THEN 
    5556         IF( ln_pisces  )    CALL trc_ice_ini_pisces       ! PISCES  bio-model 
    … …  
    5960         IF( ln_my_trc  )    CALL trc_ice_ini_my_trc       ! MY_TRC  tracers 
    6061      ENDIF 
    61  
    62       IF( nn_timing == 1 )   CALL timing_stop('trc_ice_ini') 
    6362      ! 
    6463   END SUBROUTINE trc_ice_ini 
     64 
    6565 
    6666   SUBROUTINE trc_nam_ice 
    … …  
    8686         WRITE(numout,*) '~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~' 
    8787      ENDIF 
    88  
    89       IF( nn_timing == 1 )  CALL timing_start('trc_nam_ice') 
    90  
    9188      ! 
    9289      REWIND( numnat_ref )              ! Namelist namtrc_ice in reference namelist : Passive tracer input data 
    … …  
    110107         cn_trc_o      (jn) = sn_tri_tracer(jn)%ctrc_o 
    111108      END DO 
    112  
    113       IF( nn_timing == 1 )   CALL timing_stop('trc_nam_ice') 
    114109      ! 
    115110   END SUBROUTINE trc_nam_ice 
    … …  
    122117   SUBROUTINE trc_ice_ini                   ! Dummy routine    
    123118   END SUBROUTINE trc_ice_ini 
    124  
    125119   SUBROUTINE trc_nam_ice 
    126120   END SUBROUTINE trc_nam_ice 
    127  
    128121#endif 
    129122 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/trcini.F90

    r9019 r9124  
    5252      !!--------------------------------------------------------------------- 
    5353      ! 
    54       IF( nn_timing == 1 )   CALL timing_start('trc_init') 
     54      IF( ln_timing )   CALL timing_start('trc_init') 
    5555      ! 
    5656      IF(lwp) WRITE(numout,*) 
    … …  
    8080      CALL trc_ini_inv   ! Inventories 
    8181      ! 
    82       IF( nn_timing == 1 )   CALL timing_stop('trc_init') 
     82      IF( ln_timing )   CALL timing_stop('trc_init') 
    8383      ! 
    8484   END SUBROUTINE trc_init 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/trcsms.F90

    r7646 r9124  
    4747      !!--------------------------------------------------------------------- 
    4848      ! 
    49       IF( nn_timing == 1 )   CALL timing_start('trc_sms') 
     49      IF( ln_timing )   CALL timing_start('trc_sms') 
    5050      ! 
    5151      IF( ln_pisces  )   CALL trc_sms_pisces ( kt )    ! main program of PISCES  
    … …  
    6161      ENDIF 
    6262      ! 
    63       IF( nn_timing == 1 )   CALL timing_stop('trc_sms') 
     63      IF( ln_timing )   CALL timing_stop('trc_sms') 
    6464      ! 
    6565   END SUBROUTINE trc_sms 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/trcstp.F90

    r9019 r9124  
    5959      !!------------------------------------------------------------------- 
    6060      ! 
    61       IF( nn_timing == 1 )   CALL timing_start('trc_stp') 
     61      IF( ln_timing )   CALL timing_start('trc_stp') 
    6262      ! 
    6363      IF( ( neuler == 0 .AND. kt == nittrc000 ) .OR. ln_top_euler ) THEN     ! at nittrc000 
    … …  
    1121129300  FORMAT(i10,D23.16) 
    113113      ! 
    114       IF( nn_timing == 1 )   CALL timing_stop('trc_stp') 
     114      IF( ln_timing )   CALL timing_stop('trc_stp') 
    115115      ! 
    116116   END SUBROUTINE trc_stp 
    … …  
    136136      CHARACTER(len=2) ::   cl2   ! 2 characters 
    137137      !!---------------------------------------------------------------------- 
     138      ! 
     139      IF( ln_timing )   CALL timing_start('trc_mean_qsr') 
    138140      ! 
    139141      IF( kt == nittrc000 ) THEN 
    … …  
    227229      ENDIF 
    228230      ! 
     231      IF( ln_timing )   CALL timing_stop('trc_mean_qsr') 
     232      ! 
    229233   END SUBROUTINE trc_mean_qsr 
    230234 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/trcsub.F90

    r9019 r9124  
    9797      !!------------------------------------------------------------------- 
    9898      ! 
    99       IF( nn_timing == 1 )  CALL timing_start('trc_sub_stp') 
     99      IF( ln_timing )   CALL timing_start('trc_sub_stp') 
    100100      ! 
    101101      IF( kt == nit000 ) THEN 
    … …  
    286286      ENDIF 
    287287      ! 
    288       IF( nn_timing == 1 )  CALL timing_stop('trc_sub_stp') 
     288      IF( ln_timing )   CALL timing_stop('trc_sub_stp') 
    289289      ! 
    290290   END SUBROUTINE trc_sub_stp 
    … …  
    302302      INTEGER ::   ierr 
    303303      !!------------------------------------------------------------------- 
    304       ! 
    305       IF( nn_timing == 1 )  CALL timing_start('trc_sub_ini') 
    306304      ! 
    307305      IF(lwp) WRITE(numout,*) 
    … …  
    348346      ENDIF 
    349347      ! 
    350       IF( nn_timing == 1 )  CALL timing_stop('trc_sub_ini') 
    351       ! 
    352348   END SUBROUTINE trc_sub_ini 
     349 
    353350 
    354351   SUBROUTINE trc_sub_reset( kt ) 
    … …  
    365362      !!------------------------------------------------------------------- 
    366363      ! 
    367       IF( nn_timing == 1 )  CALL timing_start('trc_sub_reset') 
     364      IF( ln_timing )   CALL timing_start('trc_sub_reset') 
    368365      ! 
    369366      !   restore physics variables 
    … …  
    443440      ! 
    444441      ! 
    445       IF( nn_timing == 1 )  CALL timing_stop('trc_sub_reset') 
     442      IF( ln_timing )   CALL timing_stop('trc_sub_reset') 
    446443      ! 
    447444   END SUBROUTINE trc_sub_reset 
    … …  
    473470      !!--------------------------------------------------------------------- 
    474471      ! 
    475       IF( nn_timing == 1 )  CALL timing_start('trc_sub_ssh') 
     472      IF( ln_timing )   CALL timing_start('trc_sub_ssh') 
    476473      ! 
    477474      ! Allocate temporary workspace 
    … …  
    533530      CALL wrk_dealloc( jpi,jpj,   zhdiv ) 
    534531      ! 
    535       IF( nn_timing == 1 )  CALL timing_stop('trc_sub_ssh') 
     532      IF( ln_timing )   CALL timing_stop('trc_sub_ssh') 
    536533      ! 
    537534   END SUBROUTINE trc_sub_ssh 

  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/trcwri.F90

    r7646 r9124  
    4444      !!--------------------------------------------------------------------- 
    4545      ! 
    46       IF( nn_timing == 1 )  CALL timing_start('trc_wri') 
     46      IF( ln_timing )   CALL timing_start('trc_wri') 
    4747      ! 
    4848      IF( l_offline .AND. kt == nittrc000 .AND. lwp ) THEN    ! WRITE root name in date.file for use by postpro 
    … …  
    6060      IF( ln_my_trc  )   CALL trc_wri_my_trc     ! MY_TRC  tracers 
    6161      ! 
    62       IF( nn_timing == 1 )  CALL timing_stop('trc_wri') 
     62      IF( ln_timing )   CALL timing_stop('trc_wri') 
    6363      ! 
    6464   END SUBROUTINE trc_wri 
Note: See TracChangeset for help on using the changeset viewer.