Changeset 6994


Ignore:
Timestamp:
2016-10-05T15:07:10+02:00 (4 years ago)
Author:
clem
Message:

clean online sea ice diagnostics

Location:
branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM
Files:
19 edited

Legend:

Unmodified
Added
Removed
  • branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/CONFIG/SHARED/namelist_ice_lim3_ref

    r6989 r6994  
    3636&namicediag    !   Diagnostics 
    3737!------------------------------------------------------------------------------ 
    38    ln_limdiahsb   =  .false.        !  check online the heat, mass & salt budgets (T) or not (F) 
    39    ln_limdiaout   =  .false.        !  output the heat, mass & salt budgets (T) or not (F) 
    40    ln_icectl      =  .false.        !  ice points output for debug (T or F) 
     38   ln_limdiachk   =  .false.        !  check online the heat, mass & salt budgets (T) or not (F) 
     39   ln_limdiahsb   =  .false.        !  output the heat, mass & salt budgets (T) or not (F) 
     40   ln_limctl      =  .false.        !  ice points output for debug (T or F) 
    4141   iiceprt        =    10           !  i-index for debug 
    4242   jiceprt        =    10           !  j-index for debug 
  • branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/LIM_SRC_3/ice.F90

    r6989 r6994  
    173173    
    174174   !                                     !!** ice-diagnostics namelist (namicediag) ** 
     175   LOGICAL , PUBLIC ::   ln_limdiachk     !: flag for ice diag (T) or not (F) 
    175176   LOGICAL , PUBLIC ::   ln_limdiahsb     !: flag for ice diag (T) or not (F) 
    176    LOGICAL , PUBLIC ::   ln_limdiaout     !: flag for ice diag (T) or not (F) 
    177    LOGICAL , PUBLIC ::   ln_icectl        !: flag for sea-ice points output (T) or not (F) 
     177   LOGICAL , PUBLIC ::   ln_limctl        !: flag for sea-ice points output (T) or not (F) 
    178178   INTEGER , PUBLIC ::   iiceprt          !: debug i-point 
    179179   INTEGER , PUBLIC ::   jiceprt          !: debug j-point 
  • branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limcons.F90

    r6970 r6994  
    164164      !!                     + test if ice concentration and volume are > 0 
    165165      !! 
    166       !! ** Method  : This is an online diagnostics which can be activated with ln_limdiahsb=true 
     166      !! ** Method  : This is an online diagnostics which can be activated with ln_limdiachk=true 
    167167      !!              It prints in ocean.output if there is a violation of conservation at each time-step 
    168168      !!              The thresholds (zv_sill, zs_sill, zh_sill) which determine violations are set to 
     
    274274      !! ** Purpose : Test the conservation of heat, salt and mass at the end of each ice time-step 
    275275      !! 
    276       !! ** Method  : This is an online diagnostics which can be activated with ln_limdiahsb=true 
     276      !! ** Method  : This is an online diagnostics which can be activated with ln_limdiachk=true 
    277277      !!              It prints in ocean.output if there is a violation of conservation at each time-step 
    278278      !!              The thresholds (zv_sill, zs_sill, zh_sill) which determine the violation are set to 
  • branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limctl.F90

    r6746 r6994  
    55   !!====================================================================== 
    66   !! History :  3.5  !  2015-01  (M. Vancoppenolle) Original code 
     7   !!            3.7  !  2016-10  (C. Rousset)       Add routine lim_prt3D 
    78   !!---------------------------------------------------------------------- 
    89#if defined key_lim3 
     
    1213   !!    lim_ctl   : control prints in case of crash 
    1314   !!    lim_prt   : ice control print at a given grid point 
     15   !!    lim_prt3D : control prints of ice arrays 
    1416   !!---------------------------------------------------------------------- 
    1517   USE oce             ! ocean dynamics and tracers 
     
    3436   PUBLIC   lim_ctl 
    3537   PUBLIC   lim_prt 
     38   PUBLIC   lim_prt3D 
    3639 
    3740   !! * Substitutions 
     
    444447   END SUBROUTINE lim_prt 
    445448 
     449   SUBROUTINE lim_prt3D( cd_routine ) 
     450      !!--------------------------------------------------------------------------------------------------------- 
     451      !!                                   ***  ROUTINE lim_prt3D *** 
     452      !! 
     453      !! ** Purpose : CTL prints of ice arrays in case ln_ctl is activated  
     454      !! 
     455      !!--------------------------------------------------------------------------------------------------------- 
     456      CHARACTER(len=*), INTENT(in)  :: cd_routine    ! name of the routine 
     457      INTEGER                       :: jk, jl        ! dummy loop indices 
     458       
     459      CALL prt_ctl_info(' ========== ') 
     460      CALL prt_ctl_info( cd_routine ) 
     461      CALL prt_ctl_info(' ========== ') 
     462      CALL prt_ctl_info(' - Cell values : ') 
     463      CALL prt_ctl_info('   ~~~~~~~~~~~~~ ') 
     464      CALL prt_ctl(tab2d_1=e12t       , clinfo1=' cell area   :') 
     465      CALL prt_ctl(tab2d_1=at_i       , clinfo1=' at_i        :') 
     466      CALL prt_ctl(tab2d_1=ato_i      , clinfo1=' ato_i       :') 
     467      CALL prt_ctl(tab2d_1=vt_i       , clinfo1=' vt_i        :') 
     468      CALL prt_ctl(tab2d_1=vt_s       , clinfo1=' vt_s        :') 
     469      CALL prt_ctl(tab2d_1=divu_i     , clinfo1=' divu_i      :') 
     470      CALL prt_ctl(tab2d_1=delta_i    , clinfo1=' delta_i     :') 
     471      CALL prt_ctl(tab2d_1=stress1_i  , clinfo1=' stress1_i   :') 
     472      CALL prt_ctl(tab2d_1=stress2_i  , clinfo1=' stress2_i   :') 
     473      CALL prt_ctl(tab2d_1=stress12_i , clinfo1=' stress12_i  :') 
     474      CALL prt_ctl(tab2d_1=strength   , clinfo1=' strength    :') 
     475      CALL prt_ctl(tab2d_1=delta_i    , clinfo1=' delta_i     :') 
     476      CALL prt_ctl(tab2d_1=u_ice      , clinfo1=' u_ice       :', tab2d_2=v_ice      , clinfo2=' v_ice       :') 
     477        
     478      DO jl = 1, jpl 
     479         CALL prt_ctl_info(' ') 
     480         CALL prt_ctl_info(' - Category : ', ivar1=jl) 
     481         CALL prt_ctl_info('   ~~~~~~~~~~') 
     482         CALL prt_ctl(tab2d_1=ht_i       (:,:,jl)        , clinfo1= ' ht_i        : ') 
     483         CALL prt_ctl(tab2d_1=ht_s       (:,:,jl)        , clinfo1= ' ht_s        : ') 
     484         CALL prt_ctl(tab2d_1=t_su       (:,:,jl)        , clinfo1= ' t_su        : ') 
     485         CALL prt_ctl(tab2d_1=t_s        (:,:,1,jl)      , clinfo1= ' t_snow      : ') 
     486         CALL prt_ctl(tab2d_1=sm_i       (:,:,jl)        , clinfo1= ' sm_i        : ') 
     487         CALL prt_ctl(tab2d_1=o_i        (:,:,jl)        , clinfo1= ' o_i         : ') 
     488         CALL prt_ctl(tab2d_1=a_i        (:,:,jl)        , clinfo1= ' a_i         : ') 
     489         CALL prt_ctl(tab2d_1=v_i        (:,:,jl)        , clinfo1= ' v_i         : ') 
     490         CALL prt_ctl(tab2d_1=v_s        (:,:,jl)        , clinfo1= ' v_s         : ') 
     491         CALL prt_ctl(tab2d_1=e_i        (:,:,1,jl)      , clinfo1= ' e_i1        : ') 
     492         CALL prt_ctl(tab2d_1=e_s        (:,:,1,jl)      , clinfo1= ' e_snow      : ') 
     493         CALL prt_ctl(tab2d_1=smv_i      (:,:,jl)        , clinfo1= ' smv_i       : ') 
     494         CALL prt_ctl(tab2d_1=oa_i       (:,:,jl)        , clinfo1= ' oa_i        : ') 
     495          
     496         DO jk = 1, nlay_i 
     497            CALL prt_ctl_info(' - Layer : ', ivar1=jk) 
     498            CALL prt_ctl(tab2d_1=t_i(:,:,jk,jl) , clinfo1= ' t_i       : ') 
     499         END DO 
     500      END DO 
     501       
     502      CALL prt_ctl_info(' ') 
     503      CALL prt_ctl_info(' - Heat / FW fluxes : ') 
     504      CALL prt_ctl_info('   ~~~~~~~~~~~~~~~~~~ ') 
     505      CALL prt_ctl(tab2d_1=sst_m  , clinfo1= ' sst   : ', tab2d_2=sss_m     , clinfo2= ' sss       : ') 
     506      CALL prt_ctl(tab2d_1=qsr    , clinfo1= ' qsr   : ', tab2d_2=qns       , clinfo2= ' qns       : ') 
     507      CALL prt_ctl(tab2d_1=emp    , clinfo1= ' emp   : ', tab2d_2=sfx       , clinfo2= ' sfx       : ') 
     508       
     509      CALL prt_ctl_info(' ') 
     510      CALL prt_ctl_info(' - Stresses : ') 
     511      CALL prt_ctl_info('   ~~~~~~~~~~ ') 
     512      CALL prt_ctl(tab2d_1=utau       , clinfo1= ' utau      : ', tab2d_2=vtau       , clinfo2= ' vtau      : ') 
     513      CALL prt_ctl(tab2d_1=utau_ice   , clinfo1= ' utau_ice  : ', tab2d_2=vtau_ice   , clinfo2= ' vtau_ice  : ') 
     514      CALL prt_ctl(tab2d_1=u_oce      , clinfo1= ' u_oce     : ', tab2d_2=v_oce      , clinfo2= ' v_oce     : ') 
     515       
     516   END SUBROUTINE lim_prt3D 
     517 
    446518#else 
    447519   !!-------------------------------------------------------------------------- 
     
    453525   SUBROUTINE lim_prt     ! Empty routine 
    454526   END SUBROUTINE lim_prt 
     527   SUBROUTINE lim_prt3D   ! Empty routine 
     528   END SUBROUTINE lim_prt3D 
    455529#endif 
    456530   !!====================================================================== 
  • branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limdyn.F90

    r6989 r6994  
    2424   USE wrk_nemo         ! work arrays 
    2525   USE in_out_manager   ! I/O manager 
    26    USE prtctl           ! Print control 
    2726   USE lib_fortran      ! glob_sum 
    2827   USE timing           ! Timing 
    2928   USE limcons          ! conservation tests 
     29   USE limctl           ! control prints 
    3030   USE limvar 
    3131 
     
    6767      ! 
    6868      ! conservation test 
    69       IF( ln_limdiahsb ) CALL lim_cons_hsm(0, 'limdyn', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 
     69      IF( ln_limdiachk ) CALL lim_cons_hsm(0, 'limdyn', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 
    7070       
    7171      ! ice velocities before rheology 
     
    8585      CALL lim_rhg 
    8686      ! 
    87       ! 
     87      ! conservation test 
     88      IF( ln_limdiachk ) CALL lim_cons_hsm(1, 'limdyn', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 
     89 
    8890      ! Control prints 
    89       IF(ln_ctl) THEN 
    90          CALL prt_ctl_info(' ') 
    91          CALL prt_ctl_info(' - Cell values : ') 
    92          CALL prt_ctl_info('   ~~~~~~~~~~~~~ ') 
    93          CALL prt_ctl(tab2d_1=divu_i    , clinfo1=' lim_dyn  : divu_i    :') 
    94          CALL prt_ctl(tab2d_1=delta_i   , clinfo1=' lim_dyn  : delta_i   :') 
    95          CALL prt_ctl(tab2d_1=strength  , clinfo1=' lim_dyn  : strength  :') 
    96          CALL prt_ctl(tab2d_1=e12t      , clinfo1=' lim_dyn  : cell area :') 
    97          CALL prt_ctl(tab2d_1=at_i      , clinfo1=' lim_dyn  : at_i      :') 
    98          CALL prt_ctl(tab2d_1=vt_i      , clinfo1=' lim_dyn  : vt_i      :') 
    99          CALL prt_ctl(tab2d_1=vt_s      , clinfo1=' lim_dyn  : vt_s      :') 
    100          CALL prt_ctl(tab2d_1=stress1_i , clinfo1=' lim_dyn  : stress1_i :') 
    101          CALL prt_ctl(tab2d_1=stress2_i , clinfo1=' lim_dyn  : stress2_i :') 
    102          CALL prt_ctl(tab2d_1=stress12_i, clinfo1=' lim_dyn  : stress12_i:') 
    103          DO jl = 1, jpl 
    104             CALL prt_ctl_info(' ') 
    105             CALL prt_ctl_info(' - Category : ', ivar1=jl) 
    106             CALL prt_ctl_info('   ~~~~~~~~~~') 
    107             CALL prt_ctl(tab2d_1=a_i   (:,:,jl)   , clinfo1= ' lim_dyn  : a_i      : ') 
    108             CALL prt_ctl(tab2d_1=ht_i  (:,:,jl)   , clinfo1= ' lim_dyn  : ht_i     : ') 
    109             CALL prt_ctl(tab2d_1=ht_s  (:,:,jl)   , clinfo1= ' lim_dyn  : ht_s     : ') 
    110             CALL prt_ctl(tab2d_1=v_i   (:,:,jl)   , clinfo1= ' lim_dyn  : v_i      : ') 
    111             CALL prt_ctl(tab2d_1=v_s   (:,:,jl)   , clinfo1= ' lim_dyn  : v_s      : ') 
    112             CALL prt_ctl(tab2d_1=e_s   (:,:,1,jl) , clinfo1= ' lim_dyn  : e_s      : ') 
    113             CALL prt_ctl(tab2d_1=t_su  (:,:,jl)   , clinfo1= ' lim_dyn  : t_su     : ') 
    114             CALL prt_ctl(tab2d_1=t_s   (:,:,1,jl) , clinfo1= ' lim_dyn  : t_snow   : ') 
    115             CALL prt_ctl(tab2d_1=sm_i  (:,:,jl)   , clinfo1= ' lim_dyn  : sm_i     : ') 
    116             CALL prt_ctl(tab2d_1=smv_i (:,:,jl)   , clinfo1= ' lim_dyn  : smv_i    : ') 
    117             DO jk = 1, nlay_i 
    118                CALL prt_ctl_info(' ') 
    119                CALL prt_ctl_info(' - Layer : ', ivar1=jk) 
    120                CALL prt_ctl_info('   ~~~~~~~') 
    121                CALL prt_ctl(tab2d_1=t_i(:,:,jk,jl) , clinfo1= ' lim_dyn  : t_i      : ') 
    122                CALL prt_ctl(tab2d_1=e_i(:,:,jk,jl) , clinfo1= ' lim_dyn  : e_i      : ') 
    123             END DO 
    124          END DO 
    125       ENDIF 
    126       ! 
    127       ! conservation test 
    128       IF( ln_limdiahsb ) CALL lim_cons_hsm(1, 'limdyn', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 
     91      IF( ln_ctl )       CALL lim_prt3D( 'limdyn' ) 
    12992      ! 
    13093      IF( nn_timing == 1 )  CALL timing_stop('limdyn') 
  • branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limitd_me.F90

    r6796 r6994  
    2222   USE lib_mpp          ! MPP library 
    2323   USE wrk_nemo         ! work arrays 
    24    USE prtctl           ! Print control 
    2524 
    2625   USE in_out_manager   ! I/O manager 
    2726   USE iom              ! I/O manager 
    2827   USE lib_fortran      ! glob_sum 
    29    USE limdiahsb 
    3028   USE timing           ! Timing 
    3129   USE limcons          ! conservation tests 
     30   USE limctl           ! control prints 
    3231 
    3332   IMPLICIT NONE 
     
    125124      CALL wrk_alloc( jpi,jpj, closing_net, divu_adv, opning, closing_gross ) 
    126125 
    127       IF(ln_ctl) THEN 
    128          CALL prt_ctl(tab2d_1=ato_i , clinfo1=' lim_itd_me: ato_i  : ', tab2d_2=at_i   , clinfo2=' at_i    : ') 
    129          CALL prt_ctl(tab2d_1=divu_i, clinfo1=' lim_itd_me: divu_i : ', tab2d_2=delta_i, clinfo2=' delta_i : ') 
    130       ENDIF 
    131  
    132126      ! conservation test 
    133       IF( ln_limdiahsb ) CALL lim_cons_hsm(0, 'limitd_me', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 
     127      IF( ln_limdiachk ) CALL lim_cons_hsm(0, 'limitd_me', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 
    134128 
    135129      !-----------------------------------------------------------------------------! 
     
    289283      ! control prints 
    290284      !-----------------------------------------------------------------------------! 
    291       IF(ln_ctl) THEN  
    292          CALL lim_var_glo2eqv 
    293  
    294          CALL prt_ctl_info(' ') 
    295          CALL prt_ctl_info(' - Cell values : ') 
    296          CALL prt_ctl_info('   ~~~~~~~~~~~~~ ') 
    297          CALL prt_ctl(tab2d_1=e12t , clinfo1=' lim_itd_me  : cell area :') 
    298          CALL prt_ctl(tab2d_1=at_i , clinfo1=' lim_itd_me  : at_i      :') 
    299          CALL prt_ctl(tab2d_1=vt_i , clinfo1=' lim_itd_me  : vt_i      :') 
    300          CALL prt_ctl(tab2d_1=vt_s , clinfo1=' lim_itd_me  : vt_s      :') 
    301          DO jl = 1, jpl 
    302             CALL prt_ctl_info(' ') 
    303             CALL prt_ctl_info(' - Category : ', ivar1=jl) 
    304             CALL prt_ctl_info('   ~~~~~~~~~~') 
    305             CALL prt_ctl(tab2d_1=a_i   (:,:,jl)   , clinfo1= ' lim_itd_me  : a_i      : ') 
    306             CALL prt_ctl(tab2d_1=ht_i  (:,:,jl)   , clinfo1= ' lim_itd_me  : ht_i     : ') 
    307             CALL prt_ctl(tab2d_1=ht_s  (:,:,jl)   , clinfo1= ' lim_itd_me  : ht_s     : ') 
    308             CALL prt_ctl(tab2d_1=v_i   (:,:,jl)   , clinfo1= ' lim_itd_me  : v_i      : ') 
    309             CALL prt_ctl(tab2d_1=v_s   (:,:,jl)   , clinfo1= ' lim_itd_me  : v_s      : ') 
    310             CALL prt_ctl(tab2d_1=e_s   (:,:,1,jl) , clinfo1= ' lim_itd_me  : e_s      : ') 
    311             CALL prt_ctl(tab2d_1=t_su  (:,:,jl)   , clinfo1= ' lim_itd_me  : t_su     : ') 
    312             CALL prt_ctl(tab2d_1=t_s   (:,:,1,jl) , clinfo1= ' lim_itd_me  : t_snow   : ') 
    313             CALL prt_ctl(tab2d_1=sm_i  (:,:,jl)   , clinfo1= ' lim_itd_me  : sm_i     : ') 
    314             CALL prt_ctl(tab2d_1=smv_i (:,:,jl)   , clinfo1= ' lim_itd_me  : smv_i    : ') 
    315             DO jk = 1, nlay_i 
    316                CALL prt_ctl_info(' ') 
    317                CALL prt_ctl_info(' - Layer : ', ivar1=jk) 
    318                CALL prt_ctl_info('   ~~~~~~~') 
    319                CALL prt_ctl(tab2d_1=t_i(:,:,jk,jl) , clinfo1= ' lim_itd_me  : t_i      : ') 
    320                CALL prt_ctl(tab2d_1=e_i(:,:,jk,jl) , clinfo1= ' lim_itd_me  : e_i      : ') 
    321             END DO 
    322          END DO 
    323       ENDIF 
    324  
    325285      ! conservation test 
    326       IF( ln_limdiahsb ) CALL lim_cons_hsm(1, 'limitd_me', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 
     286      IF( ln_limdiachk ) CALL lim_cons_hsm(1, 'limitd_me', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 
     287 
     288      ! control prints 
     289      IF( ln_ctl )       CALL lim_prt3D( 'limitd_me' ) 
    327290 
    328291      CALL wrk_dealloc( jpi, jpj, closing_net, divu_adv, opning, closing_gross ) 
  • branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limrhg.F90

    r6989 r6994  
    3838   USE bdyice_lim 
    3939#endif 
    40    !!!clem 
    41 !!   USE iom 
    4240 
    4341   IMPLICIT NONE 
     
    703701         ENDIF 
    704702      ENDIF 
    705       ! 
    706  
    707 !!      ! velocity 
    708 !!      IF ( iom_use( "uice_ipa" ) .OR. iom_use( "vice_ipa" ) .OR. iom_use( "icevel" ) ) THEN  
    709 !!         CALL iom_put( "uice_ipa"     , u_ice      )       ! ice velocity u component 
    710 !!         CALL iom_put( "vice_ipa"     , v_ice      )       ! ice velocity v component 
    711 !!      ENDIF 
    712 !!      IF ( iom_use( "iceconc_cat"  ) )  CALL iom_put( "iceconc_cat"      , a_i    )        ! area for categories 
    713 !!      IF ( iom_use( "icethic_cat"  ) )  CALL iom_put( "icethic_cat"      , ht_i   )        ! thickness for categories 
    714       
     703      !      
    715704       
    716705      CALL wrk_dealloc( jpi,jpj, z1_e1t0, z1_e2t0, zp_delt ) 
  • branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limrst.F90

    r6989 r6994  
    9393      ENDIF 
    9494      ! 
    95       IF( ln_icectl )   CALL lim_prt( kt, iiceprt, jiceprt, 1, ' - Beginning the time step - ' )   ! control print 
     95      IF( ln_limctl )   CALL lim_prt( kt, iiceprt, jiceprt, 1, ' - Beginning the time step - ' )   ! control print 
    9696   END SUBROUTINE lim_rst_opn 
    9797 
  • branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limsbc.F90

    r6989 r6994  
    3636   USE wrk_nemo         ! work arrays 
    3737   USE in_out_manager   ! I/O manager 
    38    USE prtctl           ! Print control 
    3938   USE lib_fortran      ! Fortran utilities (allows no signed zero when 'key_nosignedzero' defined)   
    4039   USE traqsr           ! add penetration of solar flux in the calculation of heat budget 
     
    215214 
    216215      ! conservation test 
    217       IF( ln_limdiahsb ) CALL lim_cons_final( 'limsbc' ) 
     216      IF( ln_limdiachk ) CALL lim_cons_final( 'limsbc' ) 
    218217 
    219218      ! control prints 
    220       IF( ln_icectl )   CALL lim_prt( kt, iiceprt, jiceprt, 3, ' - Final state lim_sbc - ' ) 
    221  
    222       IF(ln_ctl) THEN 
    223          CALL prt_ctl( tab2d_1=qsr   , clinfo1=' lim_sbc: qsr    : ', tab2d_2=qns , clinfo2=' qns     : ' ) 
    224          CALL prt_ctl( tab2d_1=emp   , clinfo1=' lim_sbc: emp    : ', tab2d_2=sfx , clinfo2=' sfx     : ' ) 
    225          CALL prt_ctl( tab2d_1=fr_i  , clinfo1=' lim_sbc: fr_i   : ' ) 
    226          CALL prt_ctl( tab3d_1=tn_ice, clinfo1=' lim_sbc: tn_ice : ', kdim=jpl ) 
    227       ENDIF 
     219      IF( ln_limctl )   CALL lim_prt( kt, iiceprt, jiceprt, 3, ' - Final state lim_sbc - ' ) 
     220      IF( ln_ctl )      CALL lim_prt3D( 'limsbc' ) 
    228221 
    229222   END SUBROUTINE lim_sbc_flx 
     
    300293      CALL lbc_lnk_multi( utau, 'U', -1., vtau, 'V', -1. )   ! lateral boundary condition 
    301294      ! 
    302       IF(ln_ctl)   CALL prt_ctl( tab2d_1=utau, clinfo1=' lim_sbc: utau   : ', mask1=umask,   & 
    303          &                       tab2d_2=vtau, clinfo2=' vtau    : '        , mask2=vmask ) 
    304295      !   
    305296   END SUBROUTINE lim_sbc_tau 
  • branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limthd.F90

    r6746 r6994  
    3939   USE wrk_nemo       ! work arrays 
    4040   USE in_out_manager ! I/O manager 
    41    USE prtctl         ! Print control 
    4241   USE lib_fortran    ! Fortran utilities (allows no signed zero when 'key_nosignedzero' defined)   
    4342   USE timing         ! Timing 
     
    104103       
    105104      ! conservation test 
    106       IF( ln_limdiahsb ) CALL lim_cons_hsm(0, 'limthd', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 
     105      IF( ln_limdiachk ) CALL lim_cons_hsm(0, 'limthd', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 
    107106 
    108107      CALL lim_var_glo2eqv 
     
    244243         ! debug point to follow 
    245244         jiindex_1d = 0 
    246          IF( ln_icectl ) THEN 
     245         IF( ln_limctl ) THEN 
    247246            DO ji = mi0(iiceprt), mi1(iiceprt) 
    248247               DO jj = mj0(jiceprt), mj1(jiceprt) 
     
    316315 
    317316      ! control checks 
    318       IF( ln_icectl )   CALL lim_prt( kt, iiceprt, jiceprt, 1, ' - ice thermodyn. - ' )   ! control print 
    319       ! 
    320       IF( ln_limdiahsb ) CALL lim_cons_hsm(1, 'limthd', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 
     317      IF( ln_limctl )    CALL lim_prt( kt, iiceprt, jiceprt, 1, ' - ice thermodyn. - ' )   ! control print 
     318      ! 
     319      IF( ln_limdiachk ) CALL lim_cons_hsm(1, 'limthd', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 
    321320 
    322321      !------------------------------------------------! 
     
    324323      !------------------------------------------------! 
    325324      ! Given thermodynamic growth rates, transport ice between thickness categories. 
    326       IF( ln_limdiahsb ) CALL lim_cons_hsm(0, 'limitd_th_rem', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 
     325      IF( ln_limdiachk ) CALL lim_cons_hsm(0, 'limitd_th_rem', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 
    327326 
    328327      IF( jpl > 1 )      CALL lim_itd_th_rem( 1, jpl, kt ) 
    329328 
    330       IF( ln_limdiahsb ) CALL lim_cons_hsm(1, 'limitd_th_rem', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 
     329      IF( ln_limdiachk ) CALL lim_cons_hsm(1, 'limitd_th_rem', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 
    331330 
    332331      !------------------------------------------------! 
    333332      !  Add frazil ice growing in leads 
    334333      !------------------------------------------------! 
    335       IF( ln_limdiahsb ) CALL lim_cons_hsm(0, 'limthd_lac', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 
     334      IF( ln_limdiachk ) CALL lim_cons_hsm(0, 'limthd_lac', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 
    336335 
    337336      IF( ln_limdO )     CALL lim_thd_lac 
    338337       
    339       IF( ln_limdiahsb ) CALL lim_cons_hsm(1, 'limthd_lac', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 
     338      ! conservation test 
     339      IF( ln_limdiachk ) CALL lim_cons_hsm(1, 'limthd_lac', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 
    340340 
    341341      ! Control print 
    342       IF(ln_ctl) THEN 
    343          CALL lim_var_glo2eqv 
    344  
    345          CALL prt_ctl_info(' ') 
    346          CALL prt_ctl_info(' - Cell values : ') 
    347          CALL prt_ctl_info('   ~~~~~~~~~~~~~ ') 
    348          CALL prt_ctl(tab2d_1=e12t , clinfo1=' lim_itd_th  : cell area :') 
    349          CALL prt_ctl(tab2d_1=at_i , clinfo1=' lim_itd_th  : at_i      :') 
    350          CALL prt_ctl(tab2d_1=vt_i , clinfo1=' lim_itd_th  : vt_i      :') 
    351          CALL prt_ctl(tab2d_1=vt_s , clinfo1=' lim_itd_th  : vt_s      :') 
    352          DO jl = 1, jpl 
    353             CALL prt_ctl_info(' ') 
    354             CALL prt_ctl_info(' - Category : ', ivar1=jl) 
    355             CALL prt_ctl_info('   ~~~~~~~~~~') 
    356             CALL prt_ctl(tab2d_1=a_i   (:,:,jl)   , clinfo1= ' lim_itd_th  : a_i      : ') 
    357             CALL prt_ctl(tab2d_1=ht_i  (:,:,jl)   , clinfo1= ' lim_itd_th  : ht_i     : ') 
    358             CALL prt_ctl(tab2d_1=ht_s  (:,:,jl)   , clinfo1= ' lim_itd_th  : ht_s     : ') 
    359             CALL prt_ctl(tab2d_1=v_i   (:,:,jl)   , clinfo1= ' lim_itd_th  : v_i      : ') 
    360             CALL prt_ctl(tab2d_1=v_s   (:,:,jl)   , clinfo1= ' lim_itd_th  : v_s      : ') 
    361             CALL prt_ctl(tab2d_1=e_s   (:,:,1,jl) , clinfo1= ' lim_itd_th  : e_s      : ') 
    362             CALL prt_ctl(tab2d_1=t_su  (:,:,jl)   , clinfo1= ' lim_itd_th  : t_su     : ') 
    363             CALL prt_ctl(tab2d_1=t_s   (:,:,1,jl) , clinfo1= ' lim_itd_th  : t_snow   : ') 
    364             CALL prt_ctl(tab2d_1=sm_i  (:,:,jl)   , clinfo1= ' lim_itd_th  : sm_i     : ') 
    365             CALL prt_ctl(tab2d_1=smv_i (:,:,jl)   , clinfo1= ' lim_itd_th  : smv_i    : ') 
    366             DO jk = 1, nlay_i 
    367                CALL prt_ctl_info(' ') 
    368                CALL prt_ctl_info(' - Layer : ', ivar1=jk) 
    369                CALL prt_ctl_info('   ~~~~~~~') 
    370                CALL prt_ctl(tab2d_1=t_i(:,:,jk,jl) , clinfo1= ' lim_itd_th  : t_i      : ') 
    371                CALL prt_ctl(tab2d_1=e_i(:,:,jk,jl) , clinfo1= ' lim_itd_th  : e_i      : ') 
    372             END DO 
    373          END DO 
    374       ENDIF 
     342      IF( ln_ctl )       CALL lim_prt3D( 'limthd' ) 
    375343      ! 
    376344      CALL wrk_dealloc( jpi,jpj, zu_io, zv_io, zfric ) 
  • branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limthd_dh.F90

    r6515 r6994  
    613613         hfx_out(ii,ij)  = hfx_out(ii,ij) + ( zq_rema(ji) * a_i_1d(ji) ) * r1_rdtice 
    614614 
    615          IF( ln_icectl .AND. zq_rema(ji) < 0. .AND. lwp ) WRITE(numout,*) 'ALERTE zq_rema <0 = ', zq_rema(ji) 
     615         IF( ln_limctl .AND. zq_rema(ji) < 0. .AND. lwp ) WRITE(numout,*) 'ALERTE zq_rema <0 = ', zq_rema(ji) 
    616616      END DO 
    617617       
  • branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limthd_dif.F90

    r5512 r6994  
    734734      END DO  ! End of the do while iterative procedure 
    735735 
    736       IF( ln_icectl .AND. lwp ) THEN 
     736      IF( ln_limctl .AND. lwp ) THEN 
    737737         WRITE(numout,*) ' zerritmax : ', zerritmax 
    738738         WRITE(numout,*) ' nconv     : ', nconv 
  • branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limthd_lac.F90

    r6746 r6994  
    254254      ! debug point to follow 
    255255      jiindex_1d = 0 
    256       IF( ln_icectl ) THEN 
     256      IF( ln_limctl ) THEN 
    257257         DO ji = mi0(iiceprt), mi1(iiceprt) 
    258258            DO jj = mj0(jiceprt), mj1(jiceprt) 
     
    264264      ENDIF 
    265265    
    266       IF( ln_icectl ) WRITE(numout,*) 'lim_thd_lac : nbpac = ', nbpac 
     266      IF( ln_limctl ) WRITE(numout,*) 'lim_thd_lac : nbpac = ', nbpac 
    267267 
    268268      !------------------------------ 
  • branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limtrp.F90

    r6989 r6994  
    109109 
    110110      ! conservation test 
    111       IF( ln_limdiahsb )   CALL lim_cons_hsm(0, 'limtrp', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 
     111      IF( ln_limdiachk )   CALL lim_cons_hsm(0, 'limtrp', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 
    112112       
    113113      ! store old values for diag 
     
    517517       
    518518      ! conservation test 
    519       IF( ln_limdiahsb )   CALL lim_cons_hsm(1, 'limtrp', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 
     519      IF( ln_limdiachk )   CALL lim_cons_hsm(1, 'limtrp', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 
    520520         
    521521      ! ------------------------------------------------- 
    522522      ! control prints 
    523523      ! ------------------------------------------------- 
    524       IF( ln_icectl )   CALL lim_prt( kt, iiceprt, jiceprt,-1, ' - ice dyn & trp - ' ) 
     524      IF( ln_limctl )   CALL lim_prt( kt, iiceprt, jiceprt,-1, ' - ice dyn & trp - ' ) 
    525525      ! 
    526526      CALL wrk_dealloc( jpi,jpj,                            zatold, zeiold, zesold, zsmvold ) 
  • branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limupdate1.F90

    r6746 r6994  
    2121   USE limitd_th 
    2222   USE limvar 
    23    USE prtctl          ! Print control 
    2423   USE wrk_nemo        ! work arrays 
    2524   USE timing          ! Timing 
    2625   USE limcons         ! conservation tests 
     26   USE limctl          ! control prints 
    2727   USE lib_mpp         ! MPP library 
    2828   USE lib_fortran     ! Fortran utilities (allows no signed zero when 'key_nosignedzero' defined)   
     
    6565 
    6666      ! conservation test 
    67       IF( ln_limdiahsb ) CALL lim_cons_hsm(0, 'limupdate1', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 
     67      IF( ln_limdiachk ) CALL lim_cons_hsm(0, 'limupdate1', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 
    6868 
    6969      !---------------------------------------------------- 
     
    135135 
    136136      ! conservation test 
    137       IF( ln_limdiahsb ) CALL lim_cons_hsm(1, 'limupdate1', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 
     137      IF( ln_limdiachk ) CALL lim_cons_hsm(1, 'limupdate1', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 
    138138 
    139       ! ------------------------------------------------- 
    140139      ! control prints 
    141       ! ------------------------------------------------- 
    142       IF(ln_ctl) THEN   ! Control print 
    143          CALL prt_ctl_info(' ') 
    144          CALL prt_ctl_info(' - Cell values : ') 
    145          CALL prt_ctl_info('   ~~~~~~~~~~~~~ ') 
    146          CALL prt_ctl(tab2d_1=e12t       , clinfo1=' lim_update1  : cell area   :') 
    147          CALL prt_ctl(tab2d_1=at_i       , clinfo1=' lim_update1  : at_i        :') 
    148          CALL prt_ctl(tab2d_1=vt_i       , clinfo1=' lim_update1  : vt_i        :') 
    149          CALL prt_ctl(tab2d_1=vt_s       , clinfo1=' lim_update1  : vt_s        :') 
    150          CALL prt_ctl(tab2d_1=strength   , clinfo1=' lim_update1  : strength    :') 
    151          CALL prt_ctl(tab2d_1=u_ice      , clinfo1=' lim_update1  : u_ice       :', tab2d_2=v_ice      , clinfo2=' v_ice       :') 
    152          CALL prt_ctl(tab2d_1=u_ice_b    , clinfo1=' lim_update1  : u_ice_b     :', tab2d_2=v_ice_b    , clinfo2=' v_ice_b     :') 
    153  
    154          DO jl = 1, jpl 
    155             CALL prt_ctl_info(' ') 
    156             CALL prt_ctl_info(' - Category : ', ivar1=jl) 
    157             CALL prt_ctl_info('   ~~~~~~~~~~') 
    158             CALL prt_ctl(tab2d_1=ht_i       (:,:,jl)        , clinfo1= ' lim_update1  : ht_i        : ') 
    159             CALL prt_ctl(tab2d_1=ht_s       (:,:,jl)        , clinfo1= ' lim_update1  : ht_s        : ') 
    160             CALL prt_ctl(tab2d_1=t_su       (:,:,jl)        , clinfo1= ' lim_update1  : t_su        : ') 
    161             CALL prt_ctl(tab2d_1=t_s        (:,:,1,jl)      , clinfo1= ' lim_update1  : t_snow      : ') 
    162             CALL prt_ctl(tab2d_1=sm_i       (:,:,jl)        , clinfo1= ' lim_update1  : sm_i        : ') 
    163             CALL prt_ctl(tab2d_1=o_i        (:,:,jl)        , clinfo1= ' lim_update1  : o_i         : ') 
    164             CALL prt_ctl(tab2d_1=a_i        (:,:,jl)        , clinfo1= ' lim_update1  : a_i         : ') 
    165             CALL prt_ctl(tab2d_1=a_i_b      (:,:,jl)        , clinfo1= ' lim_update1  : a_i_b       : ') 
    166             CALL prt_ctl(tab2d_1=v_i        (:,:,jl)        , clinfo1= ' lim_update1  : v_i         : ') 
    167             CALL prt_ctl(tab2d_1=v_i_b      (:,:,jl)        , clinfo1= ' lim_update1  : v_i_b       : ') 
    168             CALL prt_ctl(tab2d_1=v_s        (:,:,jl)        , clinfo1= ' lim_update1  : v_s         : ') 
    169             CALL prt_ctl(tab2d_1=v_s_b      (:,:,jl)        , clinfo1= ' lim_update1  : v_s_b       : ') 
    170             CALL prt_ctl(tab2d_1=e_i        (:,:,1,jl)      , clinfo1= ' lim_update1  : e_i1        : ') 
    171             CALL prt_ctl(tab2d_1=e_i_b      (:,:,1,jl)      , clinfo1= ' lim_update1  : e_i1_b      : ') 
    172             CALL prt_ctl(tab2d_1=e_i        (:,:,2,jl)      , clinfo1= ' lim_update1  : e_i2        : ') 
    173             CALL prt_ctl(tab2d_1=e_i_b      (:,:,2,jl)      , clinfo1= ' lim_update1  : e_i2_b      : ') 
    174             CALL prt_ctl(tab2d_1=e_s        (:,:,1,jl)      , clinfo1= ' lim_update1  : e_snow      : ') 
    175             CALL prt_ctl(tab2d_1=e_s_b      (:,:,1,jl)      , clinfo1= ' lim_update1  : e_snow_b    : ') 
    176             CALL prt_ctl(tab2d_1=smv_i      (:,:,jl)        , clinfo1= ' lim_update1  : smv_i       : ') 
    177             CALL prt_ctl(tab2d_1=smv_i_b    (:,:,jl)        , clinfo1= ' lim_update1  : smv_i_b     : ') 
    178             CALL prt_ctl(tab2d_1=oa_i       (:,:,jl)        , clinfo1= ' lim_update1  : oa_i        : ') 
    179             CALL prt_ctl(tab2d_1=oa_i_b     (:,:,jl)        , clinfo1= ' lim_update1  : oa_i_b      : ') 
    180  
    181             DO jk = 1, nlay_i 
    182                CALL prt_ctl_info(' - Layer : ', ivar1=jk) 
    183                CALL prt_ctl(tab2d_1=t_i(:,:,jk,jl) , clinfo1= ' lim_update1  : t_i       : ') 
    184             END DO 
    185          END DO 
    186  
    187          CALL prt_ctl_info(' ') 
    188          CALL prt_ctl_info(' - Heat / FW fluxes : ') 
    189          CALL prt_ctl_info('   ~~~~~~~~~~~~~~~~~~ ') 
    190          CALL prt_ctl(tab2d_1=sst_m  , clinfo1= ' lim_update1 : sst   : ', tab2d_2=sss_m     , clinfo2= ' sss       : ') 
    191  
    192          CALL prt_ctl_info(' ') 
    193          CALL prt_ctl_info(' - Stresses : ') 
    194          CALL prt_ctl_info('   ~~~~~~~~~~ ') 
    195          CALL prt_ctl(tab2d_1=utau       , clinfo1= ' lim_update1 : utau      : ', tab2d_2=vtau       , clinfo2= ' vtau      : ') 
    196          CALL prt_ctl(tab2d_1=utau_ice   , clinfo1= ' lim_update1 : utau_ice  : ', tab2d_2=vtau_ice   , clinfo2= ' vtau_ice  : ') 
    197          CALL prt_ctl(tab2d_1=u_oce      , clinfo1= ' lim_update1 : u_oce     : ', tab2d_2=v_oce      , clinfo2= ' v_oce     : ') 
    198       ENDIF 
     140      IF( ln_ctl )       CALL lim_prt3D( 'limupdate1' ) 
    199141    
    200142      IF( nn_timing == 1 )  CALL timing_stop('limupdate1') 
     143 
    201144   END SUBROUTINE lim_update1 
     145 
    202146#else 
    203147   !!---------------------------------------------------------------------- 
     
    207151   SUBROUTINE lim_update1     ! Empty routine 
    208152   END SUBROUTINE lim_update1 
    209  
    210153#endif 
    211154 
  • branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/LIM_SRC_3/limupdate2.F90

    r6746 r6994  
    2121   USE limitd_th 
    2222   USE limvar 
    23    USE prtctl          ! Print control 
    2423   USE lbclnk          ! lateral boundary condition - MPP exchanges 
    2524   USE wrk_nemo        ! work arrays 
     
    6766 
    6867      ! conservation test 
    69       IF( ln_limdiahsb ) CALL lim_cons_hsm(0, 'limupdate2', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 
     68      IF( ln_limdiachk ) CALL lim_cons_hsm(0, 'limupdate2', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 
    7069 
    7170      !---------------------------------------------------------------------- 
     
    176175 
    177176      ! conservation test 
    178       IF( ln_limdiahsb ) CALL lim_cons_hsm(1, 'limupdate2', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 
     177      IF( ln_limdiachk ) CALL lim_cons_hsm(1, 'limupdate2', zvi_b, zsmv_b, zei_b, zfw_b, zfs_b, zft_b) 
    179178 
    180       ! ------------------------------------------------- 
    181179      ! control prints 
    182       ! ------------------------------------------------- 
    183       IF( ln_icectl )   CALL lim_prt( kt, iiceprt, jiceprt, 2, ' - Final state - ' )   ! control print 
    184  
    185       IF(ln_ctl) THEN   ! Control print 
    186          CALL prt_ctl_info(' ') 
    187          CALL prt_ctl_info(' - Cell values : ') 
    188          CALL prt_ctl_info('   ~~~~~~~~~~~~~ ') 
    189          CALL prt_ctl(tab2d_1=e12t       , clinfo1=' lim_update2  : cell area   :') 
    190          CALL prt_ctl(tab2d_1=at_i       , clinfo1=' lim_update2  : at_i        :') 
    191          CALL prt_ctl(tab2d_1=vt_i       , clinfo1=' lim_update2  : vt_i        :') 
    192          CALL prt_ctl(tab2d_1=vt_s       , clinfo1=' lim_update2  : vt_s        :') 
    193          CALL prt_ctl(tab2d_1=strength   , clinfo1=' lim_update2  : strength    :') 
    194          CALL prt_ctl(tab2d_1=u_ice      , clinfo1=' lim_update2  : u_ice       :', tab2d_2=v_ice      , clinfo2=' v_ice       :') 
    195          CALL prt_ctl(tab2d_1=u_ice_b    , clinfo1=' lim_update2  : u_ice_b     :', tab2d_2=v_ice_b    , clinfo2=' v_ice_b     :') 
    196  
    197          DO jl = 1, jpl 
    198             CALL prt_ctl_info(' ') 
    199             CALL prt_ctl_info(' - Category : ', ivar1=jl) 
    200             CALL prt_ctl_info('   ~~~~~~~~~~') 
    201             CALL prt_ctl(tab2d_1=ht_i       (:,:,jl)        , clinfo1= ' lim_update2  : ht_i        : ') 
    202             CALL prt_ctl(tab2d_1=ht_s       (:,:,jl)        , clinfo1= ' lim_update2  : ht_s        : ') 
    203             CALL prt_ctl(tab2d_1=t_su       (:,:,jl)        , clinfo1= ' lim_update2  : t_su        : ') 
    204             CALL prt_ctl(tab2d_1=t_s        (:,:,1,jl)      , clinfo1= ' lim_update2  : t_snow      : ') 
    205             CALL prt_ctl(tab2d_1=sm_i       (:,:,jl)        , clinfo1= ' lim_update2  : sm_i        : ') 
    206             CALL prt_ctl(tab2d_1=o_i        (:,:,jl)        , clinfo1= ' lim_update2  : o_i         : ') 
    207             CALL prt_ctl(tab2d_1=a_i        (:,:,jl)        , clinfo1= ' lim_update2  : a_i         : ') 
    208             CALL prt_ctl(tab2d_1=a_i_b      (:,:,jl)        , clinfo1= ' lim_update2  : a_i_b       : ') 
    209             CALL prt_ctl(tab2d_1=v_i        (:,:,jl)        , clinfo1= ' lim_update2  : v_i         : ') 
    210             CALL prt_ctl(tab2d_1=v_i_b      (:,:,jl)        , clinfo1= ' lim_update2  : v_i_b       : ') 
    211             CALL prt_ctl(tab2d_1=v_s        (:,:,jl)        , clinfo1= ' lim_update2  : v_s         : ') 
    212             CALL prt_ctl(tab2d_1=v_s_b      (:,:,jl)        , clinfo1= ' lim_update2  : v_s_b       : ') 
    213             CALL prt_ctl(tab2d_1=e_i        (:,:,1,jl)      , clinfo1= ' lim_update2  : e_i1        : ') 
    214             CALL prt_ctl(tab2d_1=e_i_b      (:,:,1,jl)      , clinfo1= ' lim_update2  : e_i1_b      : ') 
    215             CALL prt_ctl(tab2d_1=e_i        (:,:,2,jl)      , clinfo1= ' lim_update2  : e_i2        : ') 
    216             CALL prt_ctl(tab2d_1=e_i_b      (:,:,2,jl)      , clinfo1= ' lim_update2  : e_i2_b      : ') 
    217             CALL prt_ctl(tab2d_1=e_s        (:,:,1,jl)      , clinfo1= ' lim_update2  : e_snow      : ') 
    218             CALL prt_ctl(tab2d_1=e_s_b      (:,:,1,jl)      , clinfo1= ' lim_update2  : e_snow_b    : ') 
    219             CALL prt_ctl(tab2d_1=smv_i      (:,:,jl)        , clinfo1= ' lim_update2  : smv_i       : ') 
    220             CALL prt_ctl(tab2d_1=smv_i_b    (:,:,jl)        , clinfo1= ' lim_update2  : smv_i_b     : ') 
    221             CALL prt_ctl(tab2d_1=oa_i       (:,:,jl)        , clinfo1= ' lim_update2  : oa_i        : ') 
    222             CALL prt_ctl(tab2d_1=oa_i_b     (:,:,jl)        , clinfo1= ' lim_update2  : oa_i_b      : ') 
    223  
    224             DO jk = 1, nlay_i 
    225                CALL prt_ctl_info(' - Layer : ', ivar1=jk) 
    226                CALL prt_ctl(tab2d_1=t_i(:,:,jk,jl) , clinfo1= ' lim_update2  : t_i       : ') 
    227             END DO 
    228          END DO 
    229  
    230          CALL prt_ctl_info(' ') 
    231          CALL prt_ctl_info(' - Heat / FW fluxes : ') 
    232          CALL prt_ctl_info('   ~~~~~~~~~~~~~~~~~~ ') 
    233          CALL prt_ctl(tab2d_1=sst_m  , clinfo1= ' lim_update2 : sst   : ', tab2d_2=sss_m     , clinfo2= ' sss       : ') 
    234  
    235          CALL prt_ctl_info(' ') 
    236          CALL prt_ctl_info(' - Stresses : ') 
    237          CALL prt_ctl_info('   ~~~~~~~~~~ ') 
    238          CALL prt_ctl(tab2d_1=utau       , clinfo1= ' lim_update2 : utau      : ', tab2d_2=vtau       , clinfo2= ' vtau      : ') 
    239          CALL prt_ctl(tab2d_1=utau_ice   , clinfo1= ' lim_update2 : utau_ice  : ', tab2d_2=vtau_ice   , clinfo2= ' vtau_ice  : ') 
    240          CALL prt_ctl(tab2d_1=u_oce      , clinfo1= ' lim_update2 : u_oce     : ', tab2d_2=v_oce      , clinfo2= ' v_oce     : ') 
    241       ENDIF 
     180      IF( ln_limctl )    CALL lim_prt( kt, iiceprt, jiceprt, 2, ' - Final state - ' ) 
     181      IF( ln_ctl )       CALL lim_prt3D( 'limupdate2' ) 
    242182    
    243183      IF( nn_timing == 1 )  CALL timing_stop('limupdate2') 
    244184 
    245185   END SUBROUTINE lim_update2 
     186 
    246187#else 
    247188   !!---------------------------------------------------------------------- 
     
    251192   SUBROUTINE lim_update2     ! Empty routine 
    252193   END SUBROUTINE lim_update2 
    253  
    254194#endif 
    255195 
  • branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/OPA_SRC/BDY/bdyice_lim.F90

    r6746 r6994  
    2828   USE ice             ! LIM_3 ice variables 
    2929   USE limvar 
     30   USE limctl 
    3031#endif  
    3132   USE par_oce         ! ocean parameters 
     
    7980 
    8081#if defined key_lim3 
    81       CALL lim_var_zapsmall 
    82       CALL lim_var_agg(1) 
     82                        CALL lim_var_zapsmall 
     83                        CALL lim_var_agg(1) 
     84      IF( ln_limctl )   CALL lim_prt( kt, iiceprt, jiceprt, 1, ' - ice thermo bdy - ' ) 
    8385#endif 
    8486 
  • branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/OPA_SRC/DIA/diahsb.F90

    r6970 r6994  
    240240   SUBROUTINE dia_hsb_rst( kt, cdrw ) 
    241241     !!--------------------------------------------------------------------- 
    242      !!                   ***  ROUTINE limdia_rst  *** 
     242     !!                   ***  ROUTINE dia_hsb_rst  *** 
    243243     !!                      
    244244     !! ** Purpose :   Read or write DIA file in restart file 
  • branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_lim.F90

    r6989 r6994  
    195195#if defined key_bdy 
    196196         IF( ln_limthd )              CALL bdy_ice_lim( kt )   ! -- bdy ice thermo  
    197          IF( ln_icectl )              CALL lim_prt( kt, iiceprt, jiceprt, 1, ' - ice thermo bdy - ' ) 
    198197#endif 
    199198 
     
    260259                                      CALL lim_sbc_flx( kt )    ! -- Update surface ocean mass, heat and salt fluxes 
    261260                                      ! 
    262          IF(ln_limdiaout)             CALL lim_diahsb( kt )     ! -- Diagnostics and outputs  
     261         IF( ln_limdiahsb )           CALL lim_diahsb( kt )     ! -- Diagnostics and outputs  
    263262         ! 
    264263                                      CALL lim_wri( 1 )         ! -- Ice outputs  
     
    269268         IF( lrst_ice )               CALL lim_rst_write( kt )  ! -- Ice restart file  
    270269         ! 
    271          IF( ln_icectl )              CALL lim_ctl( kt )        ! alerts in case of model crash 
     270         IF( ln_limctl )              CALL lim_ctl( kt )        ! alerts in case of model crash 
    272271         ! 
    273272      ENDIF   ! End sea-ice time step only 
     
    346345      CALL lim_sbc_init                 ! ice surface boundary condition    
    347346      ! 
    348       IF( ln_limdiaout) CALL lim_diahsb_init  ! initialization for diags 
     347      IF( ln_limdiahsb) CALL lim_diahsb_init  ! initialization for diags 
    349348      ! 
    350349      fr_i(:,:)     = at_i(:,:)         ! initialisation of sea-ice fraction 
     
    386385      NAMELIST/namicerun/ jpl, nlay_i, nlay_s, rn_amax_n, rn_amax_s, cn_icerst_in, cn_icerst_indir,   & 
    387386         &                cn_icerst_out, cn_icerst_outdir, ln_limthd, ln_limdyn, nn_limdyn, rn_uice, rn_vice   
    388       NAMELIST/namicediag/ ln_limdiahsb, ln_limdiaout, ln_icectl, iiceprt, jiceprt   
     387      NAMELIST/namicediag/ ln_limdiachk, ln_limdiahsb, ln_limctl, iiceprt, jiceprt   
    389388      !!------------------------------------------------------------------- 
    390389      !                     
     
    427426         WRITE(numout,*) '...and ice diagnostics' 
    428427         WRITE(numout,*) ' ~~~~~~~~~~~~~~~~~~~~' 
    429          WRITE(numout,*) '   Diagnose heat/mass/salt budget or not     ln_limdiahsb  = ', ln_limdiahsb 
    430          WRITE(numout,*) '   Output   heat/mass/salt budget or not     ln_limdiaout  = ', ln_limdiaout 
    431          WRITE(numout,*) '   control prints in ocean.out for (i,j)=(iiceprt,jiceprt) = ', ln_icectl 
    432          WRITE(numout,*) '   i-index for control prints (ln_icectl=true)             = ', iiceprt 
    433          WRITE(numout,*) '   j-index for control prints (ln_icectl=true)             = ', jiceprt 
     428         WRITE(numout,*) '   Diagnose online heat/mass/salt budget     ln_limdiachk  = ', ln_limdiachk 
     429         WRITE(numout,*) '   Output          heat/mass/salt budget     ln_limdiahsb  = ', ln_limdiahsb 
     430         WRITE(numout,*) '   control prints in ocean.out for (i,j)=(iiceprt,jiceprt) = ', ln_limctl 
     431         WRITE(numout,*) '   i-index for control prints (ln_limctl=true)             = ', iiceprt 
     432         WRITE(numout,*) '   j-index for control prints (ln_limctl=true)             = ', jiceprt 
    434433      ENDIF 
    435434      ! 
     
    443442      ! 
    444443#if defined key_bdy 
    445       IF( lwp .AND. ln_limdiahsb )  CALL ctl_warn('online conservation check activated but it does not work with BDY') 
     444      IF( lwp .AND. ln_limdiachk )  CALL ctl_warn('online conservation check activated but it does not work with BDY') 
    446445#endif 
    447446      ! 
Note: See TracChangeset for help on using the changeset viewer.