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 258 for trunk/NEMO/LIM_SRC – NEMO

Changeset 258 for trunk/NEMO/LIM_SRC


Ignore:
Timestamp:
2005-09-02T17:48:36+02:00 (19 years ago)
Author:
opalod
Message:

nemo_v1_update_004 : CT : Integration of the control print option for debugging work

Location:
trunk/NEMO/LIM_SRC
Files:
8 edited

Legend:

Unmodified
Added
Removed

  • trunk/NEMO/LIM_SRC/icestp.F90

    r247 r258  
    3131   USE limwri 
    3232   USE limrst 
     33   USE prtctl          ! Print control 
    3334 
    3435   IMPLICIT NONE 
    … …  
    186187          
    187188 
    188          IF(l_ctl) THEN         ! print mean trends (used for debugging) 
    189             WRITE(numout,*) 'Ice Forcings ' 
    190             WRITE(numout,*) ' qsr_oce  : ', SUM( qsr_oce (2:nictl,2:njctl) ), ' qsr_ice  : ', SUM( qsr_ice (2:nictl,2:njctl) ) 
    191             WRITE(numout,*) ' qnsr_oce : ', SUM( qnsr_oce(2:nictl,2:njctl) ), ' qnsr_ice : ', SUM( qnsr_ice(2:nictl,2:njctl) ) 
    192             WRITE(numout,*) ' evap     : ', SUM( evap    (2:nictl,2:njctl) ) 
    193             WRITE(numout,*) ' precip   : ', SUM( tprecip(2:nictl,2:njctl)  ), ' Snow     : ', SUM( sprecip (2:nictl,2:njctl) ) 
    194             WRITE(numout,*) ' u-stress : ', SUM( gtaux  (2:nictl,2:njctl)  ), ' v-stress : ', SUM( gtauy  (2:nictl,2:njctl) ) 
    195             WRITE(numout,*) ' sst      : ', SUM( sst_io (2:nictl,2:njctl)  ), ' sss      : ', SUM( sss_io (2:nictl,2:njctl) ) 
    196             WRITE(numout,*) ' u_io     : ', SUM( u_io   (2:nictl,2:njctl)  ), ' v_io     : ', SUM( v_io   (2:nictl,2:njctl) ) 
    197             WRITE(numout,*) ' hsnif  1 : ', SUM( hsnif  (2:nictl,2:njctl)  ), ' hicnif   : ', SUM( hicif  (2:nictl,2:njctl) ) 
    198             WRITE(numout,*) ' frld   1 : ', SUM( frld   (2:nictl,2:njctl)  ), ' sist     : ', SUM( sist   (2:nictl,2:njctl) ) 
     189         IF(ln_ctl) THEN         ! print mean trends (used for debugging) 
     190            CALL prt_ctl_info('Ice Forcings ') 
     191            CALL prt_ctl(tab2d_1=qsr_oce ,clinfo1=' qsr_oce  : ', tab2d_2=qsr_ice , clinfo2=' qsr_ice   : ') 
     192            CALL prt_ctl(tab2d_1=qnsr_oce,clinfo1=' qnsr_oce : ', tab2d_2=qnsr_ice, clinfo2=' qnsr_ice  : ') 
     193            CALL prt_ctl(tab2d_1=evap    ,clinfo1=' evap     : ') 
     194            CALL prt_ctl(tab2d_1=tprecip ,clinfo1=' precip   : ', tab2d_2=sprecip , clinfo2=' Snow      : ') 
     195            CALL prt_ctl(tab2d_1=gtaux   ,clinfo1=' u-stress : ', tab2d_2=gtauy   , clinfo2=' v-stress  : ') 
     196            CALL prt_ctl(tab2d_1=sst_io  ,clinfo1=' sst      : ', tab2d_2=sss_io  , clinfo2=' sss       : ') 
     197            CALL prt_ctl(tab2d_1=u_io    ,clinfo1=' u_io     : ', tab2d_2=v_io    , clinfo2=' v_io      : ') 
     198            CALL prt_ctl(tab2d_1=hsnif   ,clinfo1=' hsnif  1 : ', tab2d_2=hicif   , clinfo2=' hicif     : ') 
     199            CALL prt_ctl(tab2d_1=frld    ,clinfo1=' frld   1 : ', tab2d_2=sist    , clinfo2=' sist      : ') 
    199200         ENDIF 
    200201 
    … …  
    205206         CALL lim_dyn                                                ! Ice dynamics !   ( rheology/dynamics ) 
    206207         !                                                           !--------------! 
    207          IF(l_ctl) THEN 
    208             WRITE(numout,*) ' hsnif  2 : ', SUM( hsnif  (2:nictl,2:njctl) ), ' hicnif   : ', SUM( hicif  (2:nictl,2:njctl) ) 
    209             WRITE(numout,*) ' frld   2 : ', SUM( frld   (2:nictl,2:njctl) ), ' sist     : ', SUM( sist   (2:nictl,2:njctl) ) 
     208         IF(ln_ctl) THEN 
     209            CALL prt_ctl(tab2d_1=hsnif ,clinfo1=' hsnif  2 : ', tab2d_2=hicif , clinfo2=' hicif     : ') 
     210            CALL prt_ctl(tab2d_1=frld  ,clinfo1=' frld   2 : ', tab2d_2=sist  , clinfo2=' sist      : ') 
    210211         ENDIF 
    211212 
    … …  
    214215         CALL lim_trp                                                ! Ice transport !  ( Advection/diffusion ) 
    215216         !                                                           !---------------! 
    216          IF(l_ctl) THEN 
    217             WRITE(numout,*) ' hsnif  3 : ', SUM( hsnif  (2:nictl,2:njctl) ), ' hicnif   : ', SUM( hicif  (2:nictl,2:njctl) ) 
    218             WRITE(numout,*) ' frld   3 : ', SUM( frld   (2:nictl,2:njctl) ), ' sist     : ', SUM( sist   (2:nictl,2:njctl) ) 
     217         IF(ln_ctl) THEN 
     218            CALL prt_ctl(tab2d_1=hsnif ,clinfo1=' hsnif  3 : ', tab2d_2=hicif , clinfo2=' hicif     : ') 
     219            CALL prt_ctl(tab2d_1=frld  ,clinfo1=' frld   3 : ', tab2d_2=sist  , clinfo2=' sist      : ') 
    219220         ENDIF 
    220221 
    … …  
    223224         CALL lim_thd                                                ! Ice thermodynamics ! 
    224225         !                                                           !--------------------! 
    225          IF(l_ctl) THEN 
    226             WRITE(numout,*) ' hsnif  4 : ', SUM( hsnif  (2:nictl,2:njctl) ), ' hicnif   : ', SUM( hicif  (2:nictl,2:njctl) ) 
    227             WRITE(numout,*) ' frld   4 : ', SUM( frld   (2:nictl,2:njctl) ), ' sist     : ', SUM( sist   (2:nictl,2:njctl) ) 
     226         IF(ln_ctl) THEN 
     227            CALL prt_ctl(tab2d_1=hsnif ,clinfo1=' hsnif  4 : ', tab2d_2=hicif , clinfo2=' hicif     : ') 
     228            CALL prt_ctl(tab2d_1=frld  ,clinfo1=' frld   4 : ', tab2d_2=sist  , clinfo2=' sist      : ') 
    228229         ENDIF 
    229230 

  • trunk/NEMO/LIM_SRC/limadv.F90

    r247 r258  
    1717   USE in_out_manager  ! I/O manager 
    1818   USE lbclnk 
     19   USE prtctl          ! Print control 
    1920 
    2021   IMPLICIT NONE 
    … …  
    224225      CALL lbc_lnk( psxy, 'T', 1. ) 
    225226 
    226       IF(l_ctl) THEN 
    227          WRITE(numout,*) ' lim_adv_x: psm  ', SUM( psm (2:nictl,2:njctl) ), ' ps0  ', SUM( ps0 (2:nictl,2:njctl) ) 
    228          WRITE(numout,*) ' lim_adv_x: psx  ', SUM( psx (2:nictl,2:njctl) ), ' psxx ', SUM( psxx(2:nictl,2:njctl) ) 
    229          WRITE(numout,*) ' lim_adv_x: psy  ', SUM( psy (2:nictl,2:njctl) ), ' psyy ', SUM( psyy(2:nictl,2:njctl) ) 
    230          WRITE(numout,*) ' lim_adv_x: psxy ', SUM( psxy(2:nictl,2:njctl) ) 
     227      IF(ln_ctl)  THEN 
     228         CALL prt_ctl(tab2d_1=psm  , clinfo1=' lim_adv_x: psm  :', tab2d_2=ps0 , clinfo2=' ps0  : ') 
     229         CALL prt_ctl(tab2d_1=psx  , clinfo1=' lim_adv_x: psx  :', tab2d_2=psxx, clinfo2=' psxx : ') 
     230         CALL prt_ctl(tab2d_1=psy  , clinfo1=' lim_adv_x: psy  :', tab2d_2=psyy, clinfo2=' psyy : ') 
     231         CALL prt_ctl(tab2d_1=psxy , clinfo1=' lim_adv_x: psxy :') 
    231232      ENDIF 
    232233 
    … …  
    424425      CALL lbc_lnk( psxy, 'T', 1. ) 
    425426 
    426       IF(l_ctl) THEN 
    427          WRITE(numout,*) ' lim_adv_y: psm  ', SUM( psm (2:nictl,2:njctl) ), ' ps0  ', SUM( ps0 (2:nictl,2:njctl) ) 
    428          WRITE(numout,*) ' lim_adv_y: psx  ', SUM( psx (2:nictl,2:njctl) ), ' psxx ', SUM( psxx(2:nictl,2:njctl) ) 
    429          WRITE(numout,*) ' lim_adv_y: psy  ', SUM( psy (2:nictl,2:njctl) ), ' psyy ', SUM( psyy(2:nictl,2:njctl) ) 
    430          WRITE(numout,*) ' lim_adv_y: psxy ', SUM( psxy(2:nictl,2:njctl) ) 
     427      IF(ln_ctl) THEN 
     428         CALL prt_ctl(tab2d_1=psm  , clinfo1=' lim_adv_y: psm  :', tab2d_2=ps0 , clinfo2=' ps0  : ') 
     429         CALL prt_ctl(tab2d_1=psx  , clinfo1=' lim_adv_y: psx  :', tab2d_2=psxx, clinfo2=' psxx : ') 
     430         CALL prt_ctl(tab2d_1=psy  , clinfo1=' lim_adv_y: psy  :', tab2d_2=psyy, clinfo2=' psyy : ') 
     431         CALL prt_ctl(tab2d_1=psxy , clinfo1=' lim_adv_y: psxy :') 
    431432      ENDIF 
    432433 

  • trunk/NEMO/LIM_SRC/limdyn.F90

    r247 r258  
    2323   USE lbclnk 
    2424   USE lib_mpp 
     25   USE prtctl          ! Print control 
    2526 
    2627   IMPLICIT NONE 
    … …  
    9293            i_j1 = 1    
    9394            i_jpj = jpj 
    94             IF(l_ctl)   WRITE(numout,*) 'lim_dyn :    i_j1 = ', i_j1, ' ij_jpj = ', i_jpj 
     95            IF(ln_ctl)    THEN 
     96               CALL prt_ctl_info('lim_dyn  :    i_j1 = ', ivar1=i_j1, clinfo2=' ij_jpj = ', ivar2=i_jpj) 
     97            ENDIF 
    9598            CALL lim_rhg( i_j1, i_jpj ) 
    9699 
    … …  
    113116               END DO 
    114117               i_j1 = MAX( 1, i_j1-1 ) 
    115                IF(l_ctl)   WRITE(numout,*) 'lim_dyn : NH i_j1 = ', i_j1, ' ij_jpj = ', i_jpj 
     118               IF(ln_ctl)   WRITE(numout,*) 'lim_dyn : NH i_j1 = ', i_j1, ' ij_jpj = ', i_jpj 
    116119     
    117120               CALL lim_rhg( i_j1, i_jpj ) 
    … …  
    124127               END DO 
    125128               i_jpj = MIN( jpj, i_jpj+2 ) 
    126                IF(l_ctl)   WRITE(numout,*) 'lim_dyn : SH i_j1 = ', i_j1, ' ij_jpj = ', i_jpj 
     129               IF(ln_ctl)   WRITE(numout,*) 'lim_dyn : SH i_j1 = ', i_j1, ' ij_jpj = ', i_jpj 
    127130     
    128131               CALL lim_rhg( i_j1, i_jpj ) 
    … …  
    143146               i_jpj = MIN( jpj, i_jpj+2) 
    144147     
    145                IF(l_ctl)   WRITE(numout,*) 'lim_dyn : one hemisphere: i_j1 = ', i_j1, ' ij_jpj = ', i_jpj 
     148               IF(ln_ctl)   WRITE(numout,*) 'lim_dyn : one hemisphere: i_j1 = ', i_j1, ' ij_jpj = ', i_jpj 
    146149     
    147150               CALL lim_rhg( i_j1, i_jpj ) 
    … …  
    154157         v_ice(:,1) = 0.e0 
    155158 
    156          IF(l_ctl) THEN  
    157             WRITE(numout,*) ' lim_dyn  : u_oce ', SUM( u_oce(2:nictl,2:njctl) ), ' v_oce ', SUM( v_oce(2:nictl,2:njctl) ) 
    158             WRITE(numout,*) ' lim_dyn  : u_ice ', SUM( u_ice(2:nictl,2:njctl) ), ' v_ice ', SUM( v_ice(2:nictl,2:njctl) ) 
     159         IF(ln_ctl)  THEN  
     160            CALL prt_ctl(tab2d_1=u_oce , clinfo1=' lim_dyn  : u_oce :', tab2d_2=v_oce , clinfo2=' v_oce :') 
     161            CALL prt_ctl(tab2d_1=u_ice , clinfo1=' lim_dyn  : u_ice :', tab2d_2=v_ice , clinfo2=' v_ice :') 
    159162         ENDIF 
    160163          
    … …  
    244247      CALL lbc_lnk( tio_v, 'I', -1. )   ! I-point (i.e. ice U-V point) 
    245248 
    246       IF(l_ctl) THEN  
    247          WRITE(numout,*) ' lim_dyn  : tio_u ', SUM( tio_u(2:nictl,2:njctl) ), ' tio_v ', SUM( tio_v(2:nictl,2:njctl) ) 
    248          WRITE(numout,*) ' lim_dyn  : ust2s ', SUM( ust2s(2:nictl,2:njctl) ) 
     249      IF(ln_ctl) THEN  
     250            CALL prt_ctl(tab2d_1=tio_u , clinfo1=' lim_dyn  : tio_u :', tab2d_2=tio_v , clinfo2=' tio_v :') 
     251            CALL prt_ctl(tab2d_1=ust2s , clinfo1=' lim_dyn  : ust2s :') 
    249252      ENDIF 
    250253 

  • trunk/NEMO/LIM_SRC/limflx.F90

    r247 r258  
    1111   !!---------------------------------------------------------------------- 
    1212   !! * Modules used 
    13    USE par_oce 
    14    USE phycst 
    15    USE ocfzpt 
    16    USE ice_oce 
    17    USE flx_oce 
    18    USE ice 
    19    USE flxblk 
    20    USE lbclnk 
    21    USE in_out_manager 
    22    USE albedo 
     13   USE par_oce          ! ocean parameters 
     14   USE phycst           ! physical constants 
     15   USE ocfzpt           ! surface ocean freezing point 
     16   USE ice_oce          ! sea-ice variable 
     17   USE flx_oce          ! sea-ice/ocean forcings variables 
     18   USE ice              ! LIM sea-ice variables 
     19   USE flxblk           ! bulk formulea 
     20   USE lbclnk           ! ocean lateral boundary condition 
     21   USE in_out_manager   ! I/O manager 
     22   USE albedo           ! albedo parameters 
     23   USE prtctl           ! Print control 
    2324 
    2425   IMPLICIT NONE 
    … …  
    195196#endif 
    196197 
    197       IF(l_ctl) THEN 
    198          WRITE(numout,*) ' lim_flx  ' 
    199          WRITE(numout,*) ' fsolar ', SUM( fsolar(2:nictl,2:njctl) ), ' fnsolar', SUM( fnsolar(2:nictl,2:njctl) ) 
    200          WRITE(numout,*) ' fmass  ', SUM( fmass (2:nictl,2:njctl) ), ' fsalt  ', SUM( fsalt  (2:nictl,2:njctl) ) 
    201          WRITE(numout,*) ' ftaux  ', SUM( ftaux (2:nictl,2:njctl) ), ' ftauy  ', SUM( ftauy  (2:nictl,2:njctl) ) 
    202          WRITE(numout,*) ' freeze ', SUM( freeze(2:nictl,2:njctl) ), ' tn_ice ', SUM( tn_ice (2:nictl,2:njctl) ) 
     198      IF(ln_ctl) THEN 
     199         CALL prt_ctl(tab2d_1=fsolar, clinfo1=' lim_flx: fsolar : ', tab2d_2=fnsolar, clinfo2=' fnsolar : ') 
     200         CALL prt_ctl(tab2d_1=fmass , clinfo1=' lim_flx: fmass  : ', tab2d_2=fsalt  , clinfo2=' fsalt   : ') 
     201         CALL prt_ctl(tab2d_1=ftaux , clinfo1=' lim_flx: ftaux  : ', tab2d_2=ftauy  , clinfo2=' ftauy   : ') 
     202         CALL prt_ctl(tab2d_1=freeze, clinfo1=' lim_flx: freeze : ', tab2d_2=tn_ice , clinfo2=' tn_ice  : ') 
    203203      ENDIF  
    204204    

  • trunk/NEMO/LIM_SRC/limhdf.F90

    r247 r258  
    1818   USE lbclnk 
    1919   USE lib_mpp 
     20   USE prtctl          ! Print control 
    2021 
    2122   IMPLICIT NONE 
    … …  
    6768      INTEGER ::  & 
    6869         its, iter            ! temporary integers 
     70      CHARACTER (len=55) :: charout 
    6971      REAL(wp) ::  & 
    7072         zalfa, zrlxint, zconv, zeps   ! temporary scalars 
    … …  
    168170 
    169171      ptab(:,:) = ptab(:,:) 
    170  
    171       IF(l_ctl)   WRITE(numout,*) ' lim_hdf  : ', SUM( ptab(2:nictl,2:njctl)-ptab0(2:nictl,2:njctl) ), & 
    172       &                                        ' zconv= ', zconv, ' iter= ', iter 
     172      IF(ln_ctl)   THEN 
     173         zrlx(:,:) = ptab(:,:) - ptab0(:,:) 
     174         WRITE(charout,FMT="(' lim_hdf  : zconv =',D23.16, ' iter =',I4,2X)") zconv, iter 
     175         CALL prt_ctl(tab2d_1=zrlx, clinfo1=charout) 
     176      ENDIF 
    173177 
    174178   END SUBROUTINE lim_hdf 

  • trunk/NEMO/LIM_SRC/limrhg.F90

    r247 r258  
    1919   USE lib_mpp 
    2020   USE in_out_manager  ! I/O manager 
     21   USE prtctl          ! Print control 
    2122 
    2223   IMPLICIT NONE 
    … …  
    6465         iim1, ijm1, iip1 , ijp1   , & ! temporary integers 
    6566         iter, jter                    !    "          " 
     67 
     68      CHARACTER (len=50) :: charout 
    6669 
    6770      REAL(wp) :: & 
    … …  
    632635      !                                                   ! ==================== ! 
    633636 
    634       IF(l_ctl) THEN 
    635          WRITE(numout,*) ' lim_rhg  : res= ', zresm, 'iter= ', jter,' u_ice ', SUM( u_ice(2:nictl,2:njctl) ) , & 
    636          &                                                          ' v_ice ', SUM( v_ice(2:nictl,2:njctl) ) 
     637      IF(ln_ctl) THEN 
     638         WRITE(charout,FMT="('lim_rhg  : res =',D23.16, ' iter =',I4)") zresm, jter 
     639         CALL prt_ctl_info(charout) 
     640         CALL prt_ctl(tab2d_1=u_ice, clinfo1=' lim_rhg  : u_ice :', tab2d_2=v_ice, clinfo2=' v_ice :') 
    637641      ENDIF 
    638642 

  • trunk/NEMO/LIM_SRC/limthd.F90

    r247 r258  
    2525   USE limthd_lac 
    2626   USE limtab 
     27   USE prtctl          ! Print control 
    2728       
    2829   IMPLICIT NONE 
    … …  
    7879         nbpb  ,               &   ! nb of icy pts for thermo. cal. 
    7980         nbpac                     ! nb of pts for lateral accretion  
     81      CHARACTER (len=22) :: charout 
    8082      REAL(wp) ::  & 
    8183         zfric_umin = 5e-03 ,  &   ! lower bound for the friction velocity 
    … …  
    9294         zhicifp            ,  &   ! ice thickness for outputs 
    9395         zqlbsbq                   ! link with lead energy budget qldif 
     96      REAL(wp), DIMENSION(jpi,jpj,2) :: & 
     97         zmsk                      ! working array 
    9498      !!------------------------------------------------------------------- 
    9599 
    … …  
    112116      rdmicif(:,:) = 0.e0   ! variation of ice mass per unit area 
    113117      hicifp (:,:) = 0.e0   ! daily thermodynamic ice production.  
     118      zmsk (:,:,:) = 0.e0 
    114119 
    115120      DO jj = 1, jpj 
    … …  
    118123         END DO 
    119124      END DO 
    120       IF(l_ctl)   WRITE(numout,*) 'lim_thd  : ', SUM( hsnif(2:nictl,2:njctl) ) , ' hsnif' 
    121        
     125 
     126      IF(ln_ctl)   CALL prt_ctl(tab2d_1=hsnif     , clinfo1=' lim_thd: hsnif   : ') 
    122127       
    123128      !-----------------------------------! 
    … …  
    150155         END DO 
    151156      END DO 
    152       IF(l_ctl) THEN 
    153          WRITE(numout,*) 'lim_thd: hicif : ', SUM( hicif (2:nictl,2:njctl) ) 
    154          WRITE(numout,*) 'lim_thd: hsnif : ', SUM( hsnif (2:nictl,2:njctl) ) 
    155          WRITE(numout,*) 'lim_thd: dmgwi : ', SUM( dmgwi (2:nictl,2:njctl) ) 
    156          WRITE(numout,*) 'lim_thd: qstoif: ', SUM( qstoif(2:nictl,2:njctl) ) 
    157          WRITE(numout,*) 'lim_thd: frld  : ', SUM( frld  (2:nictl,2:njctl) ) 
     157 
     158      IF(ln_ctl) THEN 
     159         CALL prt_ctl(tab2d_1=hicif  , clinfo1=' lim_thd: hicif   : ') 
     160         CALL prt_ctl(tab2d_1=hsnif  , clinfo1=' lim_thd: hsnif   : ') 
     161         CALL prt_ctl(tab2d_1=dmgwi  , clinfo1=' lim_thd: dmgwi   : ') 
     162         CALL prt_ctl(tab2d_1=qstoif , clinfo1=' lim_thd: qstoif  : ') 
     163         CALL prt_ctl(tab2d_1=frld   , clinfo1=' lim_thd: frld    : ') 
    158164      ENDIF 
    159165 
    … …  
    221227         END DO 
    222228      END DO 
    223       IF(l_ctl) THEN 
    224          WRITE(numout,*) 'lim_thd: pfrld ' , SUM( pfrld  (2:nictl,2:njctl) ), ' thcm    ', SUM( thcm   (2:nictl,2:njctl) ) 
    225          WRITE(numout,*) 'lim_thd: fdtcn ' , SUM( fdtcn  (2:nictl,2:njctl) ), ' qdtcn   ', SUM( qdtcn  (2:nictl,2:njctl) ) 
    226          WRITE(numout,*) 'lim_thd: qldif ' , SUM( qldif  (2:nictl,2:njctl) ), ' zqlbsbq ', SUM( zqlbsbq(2:nictl,2:njctl) ) 
    227          WRITE(numout,*) 'lim_thd: qcmif ' , SUM( qcmif  (2:nictl,2:njctl) ), ' fbif    ', SUM( fbif   (2:nictl,2:njctl) ) 
    228          WRITE(numout,*) 'lim_thd: qcmif  ', SUM( qcmif  (2:nictl,2:njctl)*tms(2:nictl,2:njctl) ) 
    229          WRITE(numout,*) 'lim_thd: zhicifp', SUM( zhicifp(2:nictl,2:njctl) ) 
    230          WRITE(numout,*) 'limthd : nbpb = ', nbpb 
     229 
     230      IF(ln_ctl) THEN 
     231         CALL prt_ctl(tab2d_1=pfrld, clinfo1=' lim_thd: pfrld   : ', tab2d_2=thcm   , clinfo2='  thcm    : ') 
     232         CALL prt_ctl(tab2d_1=fdtcn, clinfo1=' lim_thd: fdtcn   : ', tab2d_2=qdtcn  , clinfo2='  qdtcn   : ') 
     233         CALL prt_ctl(tab2d_1=qldif, clinfo1=' lim_thd: qldif   : ', tab2d_2=zqlbsbq, clinfo2='  zqlbsbq : ') 
     234         CALL prt_ctl(tab2d_1=qcmif, clinfo1=' lim_thd: qcmif   : ', tab2d_2=fbif   , clinfo2='  fbif    : ') 
     235         zmsk(:,:,1) = tms(:,:) 
     236         CALL prt_ctl(tab2d_1=qcmif  , clinfo1=' lim_thd: qcmif   : ', mask1=zmsk) 
     237         CALL prt_ctl(tab2d_1=zhicifp, clinfo1=' lim_thd: zhicifp : ') 
     238         WRITE(charout, FMT="('lim_thd: nbpb = ',I4)") nbpb 
     239         CALL prt_ctl_info(charout) 
    231240      ENDIF 
    232241       
    … …  
    327336      END DO 
    328337       
    329       IF(l_ctl) THEN 
    330          WRITE(numout,*) 'lim_thd : phicif ', SUM( phicif(2:nictl,2:njctl) ), ' hicif ', SUM( hicif(2:nictl,2:njctl) ) 
    331          WRITE(numout,*) 'lim_thd : nbpac = ', nbpac 
     338      IF(ln_ctl) THEN 
     339         CALL prt_ctl(tab2d_1=phicif, clinfo1=' lim_thd: phicif  : ', tab2d_2=hicif, clinfo2=' hicif : ') 
     340         WRITE(charout, FMT="('lim_thd: nbpac = ',I4)") nbpac 
     341         CALL prt_ctl_info(charout) 
    332342      ENDIF 
    333343 
    … …  
    381391      END DO 
    382392 
    383       IF(l_ctl) THEN 
    384          WRITE(numout,*) ' lim_thd  end  ' 
    385          WRITE(numout,*) '  hicif ', SUM( hicif (2:nictl,2:njctl)   ), '  hsnif ', SUM( hsnif (2:nictl,2:njctl)   ) 
    386          WRITE(numout,*) '  frld  ', SUM( frld  (2:nictl,2:njctl)   ), '  hicifp', SUM( hicifp(2:nictl,2:njctl)   ) 
    387          WRITE(numout,*) '  phicif', SUM( phicif(2:nictl,2:njctl)   ), '  pfrld ', SUM( pfrld (2:nictl,2:njctl)   ) 
    388          WRITE(numout,*) '  sist  ', SUM( sist  (2:nictl,2:njctl)   ), '  tbif 1', SUM( tbif  (2:nictl,2:njctl,1) ) 
    389          WRITE(numout,*) '  tbif 2', SUM( tbif  (2:nictl,2:njctl,2) ), '  tbif 3', SUM( tbif  (2:nictl,2:njctl,3) ) 
    390          WRITE(numout,*) '  fdtcn ', SUM( fdtcn (2:nictl,2:njctl)   ), '  qdtcn ', SUM( qdtcn (2:nictl,2:njctl)   ) 
    391          WRITE(numout,*) '  qstoif', SUM( qstoif(2:nictl,2:njctl)   ), '  fsbbq ', SUM( fsbbq (2:nictl,2:njctl)   ) 
     393      IF(ln_ctl) THEN 
     394         CALL prt_ctl_info(' lim_thd  end  ') 
     395         CALL prt_ctl(tab2d_1=hicif , clinfo1=' lim_thd: hicif   : ', tab2d_2=hsnif , clinfo2=' hsnif  : ') 
     396         CALL prt_ctl(tab2d_1=frld  , clinfo1=' lim_thd: frld    : ', tab2d_2=hicifp, clinfo2=' hicifp : ') 
     397         CALL prt_ctl(tab2d_1=phicif, clinfo1=' lim_thd: phicif  : ', tab2d_2=pfrld , clinfo2=' pfrld  : ') 
     398         CALL prt_ctl(tab2d_1=sist  , clinfo1=' lim_thd: sist    : ') 
     399         CALL prt_ctl(tab2d_1=tbif(:,:,1), clinfo1=' lim_thd: tbif 1  : ') 
     400         CALL prt_ctl(tab2d_1=tbif(:,:,2), clinfo1=' lim_thd: tbif 2  : ') 
     401         CALL prt_ctl(tab2d_1=tbif(:,:,3), clinfo1=' lim_thd: tbif 3  : ') 
     402         CALL prt_ctl(tab2d_1=fdtcn , clinfo1=' lim_thd: fdtcn   : ', tab2d_2=qdtcn , clinfo2=' qdtcn  : ') 
     403         CALL prt_ctl(tab2d_1=qstoif, clinfo1=' lim_thd: qstoif  : ', tab2d_2=fsbbq , clinfo2=' fsbbq  : ') 
    392404      ENDIF 
    393405 

  • trunk/NEMO/LIM_SRC/limwri.F90

    r247 r258  
    136136      !----------------------------------------------------------------------- 
    137137 
    138       IF(lwp .AND. L_ctl) THEN 
     138      IF(lwp) THEN 
    139139         WRITE(numout,*) 
    140140         WRITE(numout,*) 'lim_wri : write ice outputs in NetCDF files at time : ', nyear, nmonth, nday, numit 
Note: See TracChangeset for help on using the changeset viewer.