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 863 – NEMO

Changeset 863


Ignore:
Timestamp:
2008-03-14T14:34:10+01:00 (16 years ago)
Author:
ctlod
Message:

Add control prints

Location:
trunk/NEMO/LIM_SRC_3
Files:
6 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/LIM_SRC_3/limdyn.F90

    r834 r863  
    1818   USE taumod 
    1919   USE ice 
     20   USE par_ice 
    2021   USE ice_oce 
    2122   USE iceini 
     
    6263      !!------------------------------------------------------------------------------------ 
    6364      !! * Local variables 
    64       INTEGER ::   ji, jj            ! dummy loop indices 
     65      INTEGER ::   ji, jj, jl, ja    ! dummy loop indices 
    6566      INTEGER :: i_j1, i_jpj         ! Starting/ending j-indices for rheology 
    6667 
     
    269270      CALL lbc_lnk( tio_v, 'V', -1. )   ! I-point (i.e. ice U-V point) 
    270271 
    271       IF(ln_ctl) THEN  
    272         CALL prt_ctl(tab2d_1=tio_u , clinfo1=' lim_dyn  : tio_u :', tab2d_2=tio_v , clinfo2=' tio_v :') 
    273         CALL prt_ctl(tab2d_1=ust2s , clinfo1=' lim_dyn  : ust2s :') 
    274         CALL prt_ctl(tab2d_1=at_i  , clinfo1=' lim_dyn  : at_i  :') 
     272      IF(ln_ctl) THEN   ! Control print 
     273         CALL prt_ctl(tab2d_1=tio_u     , clinfo1=' lim_dyn  : tio_u     :', tab2d_2=tio_v , clinfo2=' tio_v :') 
     274         CALL prt_ctl(tab2d_1=ust2s     , clinfo1=' lim_dyn  : ust2s     :') 
     275         CALL prt_ctl(tab2d_1=divu_i    , clinfo1=' lim_dyn  : divu_i    :') 
     276         CALL prt_ctl(tab2d_1=delta_i   , clinfo1=' lim_dyn  : delta_i   :') 
     277         CALL prt_ctl(tab2d_1=strength  , clinfo1=' lim_dyn  : strength  :') 
     278         CALL prt_ctl(tab2d_1=area      , clinfo1=' lim_dyn  : cell area :') 
     279         CALL prt_ctl(tab2d_1=at_i      , clinfo1=' lim_dyn  : at_i      :') 
     280         CALL prt_ctl(tab2d_1=vt_i      , clinfo1=' lim_dyn  : vt_i      :') 
     281         CALL prt_ctl(tab2d_1=vt_s      , clinfo1=' lim_dyn  : vt_s      :') 
     282         CALL prt_ctl(tab2d_1=stress1_i , clinfo1=' lim_dyn  : stress1_i :') 
     283         CALL prt_ctl(tab2d_1=stress2_i , clinfo1=' lim_dyn  : stress2_i :') 
     284         CALL prt_ctl(tab2d_1=stress12_i, clinfo1=' lim_dyn  : stress12_i:') 
     285         DO jl = 1, jpl 
     286            CALL prt_ctl_info(' - Category : ', ivar1=jl) 
     287            CALL prt_ctl_info('   ~~~~~~~~~~') 
     288            CALL prt_ctl(tab2d_1=a_i   (:,:,jl)   , clinfo1= ' lim_dyn  : a_i      : ') 
     289            CALL prt_ctl(tab2d_1=ht_i  (:,:,jl)   , clinfo1= ' lim_dyn  : ht_i     : ') 
     290            CALL prt_ctl(tab2d_1=ht_s  (:,:,jl)   , clinfo1= ' lim_dyn  : ht_s     : ') 
     291            CALL prt_ctl(tab2d_1=v_i   (:,:,jl)   , clinfo1= ' lim_dyn  : v_i      : ') 
     292            CALL prt_ctl(tab2d_1=v_s   (:,:,jl)   , clinfo1= ' lim_dyn  : v_s      : ') 
     293            CALL prt_ctl(tab2d_1=e_s   (:,:,1,jl) , clinfo1= ' lim_dyn  : e_s      : ') 
     294            CALL prt_ctl(tab2d_1=t_su  (:,:,jl)   , clinfo1= ' lim_dyn  : t_su     : ') 
     295            CALL prt_ctl(tab2d_1=t_s   (:,:,1,jl) , clinfo1= ' lim_dyn  : t_snow   : ') 
     296            CALL prt_ctl(tab2d_1=sm_i  (:,:,jl)   , clinfo1= ' lim_dyn  : sm_i     : ') 
     297            CALL prt_ctl(tab2d_1=smv_i (:,:,jl)   , clinfo1= ' lim_dyn  : smv_i    : ') 
     298            DO ja = 1, nlay_i 
     299               CALL prt_ctl_info(' - Layer : ', ivar1=ja) 
     300               CALL prt_ctl_info('   ~~~~~~~') 
     301               CALL prt_ctl(tab2d_1=t_i(:,:,ja,jl) , clinfo1= ' lim_dyn  : t_i      : ') 
     302               CALL prt_ctl(tab2d_1=e_i(:,:,ja,jl) , clinfo1= ' lim_dyn  : e_i      : ') 
     303            END DO 
     304         END DO 
    275305      ENDIF 
    276306 
  • trunk/NEMO/LIM_SRC_3/limitd_me.F90

    r862 r863  
    206206         CALL prt_ctl(tab2d_1=ato_i , clinfo1=' lim_itd_me: ato_i  : ', tab2d_2=at_i   , clinfo2=' at_i    : ') 
    207207         CALL prt_ctl(tab2d_1=divu_i, clinfo1=' lim_itd_me: divu_i : ', tab2d_2=delta_i, clinfo2=' delta_i : ') 
    208          CALL prt_ctl(tab3d_1=a_i   , clinfo1=' lim_itd_me: a_i    : ', kdim=jpl) 
    209          CALL prt_ctl(tab3d_1=v_i   , clinfo1=' lim_itd_me: v_i    : ', kdim=jpl) 
    210208      ENDIF 
    211209 
     
    479477      IF ( ( num_sal .EQ. 2 ) .OR. ( num_sal .EQ. 4 ) ) &  
    480478      smv_i(:,:,:)  = old_smv_i(:,:,:) 
     479 
     480      IF(ln_ctl) THEN     ! Control print 
     481         CALL prt_ctl(tab2d_1=area , clinfo1=' lim_itd_me  : cell area :') 
     482         CALL prt_ctl(tab2d_1=at_i , clinfo1=' lim_itd_me  : at_i      :') 
     483         CALL prt_ctl(tab2d_1=vt_i , clinfo1=' lim_itd_me  : vt_i      :') 
     484         CALL prt_ctl(tab2d_1=vt_s , clinfo1=' lim_itd_me  : vt_s      :') 
     485         DO jl = 1, jpl 
     486            CALL prt_ctl_info(' - Category : ', ivar1=jl) 
     487            CALL prt_ctl_info('   ~~~~~~~~~~') 
     488            CALL prt_ctl(tab2d_1=a_i   (:,:,jl)   , clinfo1= ' lim_itd_me  : a_i      : ') 
     489            CALL prt_ctl(tab2d_1=ht_i  (:,:,jl)   , clinfo1= ' lim_itd_me  : ht_i     : ') 
     490            CALL prt_ctl(tab2d_1=ht_s  (:,:,jl)   , clinfo1= ' lim_itd_me  : ht_s     : ') 
     491            CALL prt_ctl(tab2d_1=v_i   (:,:,jl)   , clinfo1= ' lim_itd_me  : v_i      : ') 
     492            CALL prt_ctl(tab2d_1=v_s   (:,:,jl)   , clinfo1= ' lim_itd_me  : v_s      : ') 
     493            CALL prt_ctl(tab2d_1=e_s   (:,:,1,jl) , clinfo1= ' lim_itd_me  : e_s      : ') 
     494            CALL prt_ctl(tab2d_1=t_su  (:,:,jl)   , clinfo1= ' lim_itd_me  : t_su     : ') 
     495            CALL prt_ctl(tab2d_1=t_s   (:,:,1,jl) , clinfo1= ' lim_itd_me  : t_snow   : ') 
     496            CALL prt_ctl(tab2d_1=sm_i  (:,:,jl)   , clinfo1= ' lim_itd_me  : sm_i     : ') 
     497            CALL prt_ctl(tab2d_1=smv_i (:,:,jl)   , clinfo1= ' lim_itd_me  : smv_i    : ') 
     498            DO jk = 1, nlay_i 
     499               CALL prt_ctl_info(' - Layer : ', ivar1=jk) 
     500               CALL prt_ctl_info('   ~~~~~~~') 
     501               CALL prt_ctl(tab2d_1=t_i(:,:,jk,jl) , clinfo1= ' lim_itd_me  : t_i      : ') 
     502               CALL prt_ctl(tab2d_1=e_i(:,:,jk,jl) , clinfo1= ' lim_itd_me  : e_i      : ') 
     503            END DO 
     504         END DO 
     505      ENDIF 
    481506 
    482507 
  • trunk/NEMO/LIM_SRC_3/limitd_th.F90

    r862 r863  
    2626   USE iceini 
    2727   USE limcons 
     28   USE prtctl           ! Print control 
    2829  
    2930   IMPLICIT NONE 
     
    7980 
    8081       !! * Local variables 
    81        INTEGER ::   jm,       &   ! ice types    dummy loop index 
     82       INTEGER ::   jl, ja,   &   ! ice category, layers 
     83                    jm,       &   ! ice types    dummy loop index 
    8284                    jbnd1,    & 
    8385                    jbnd2 
     
    141143      IF ( ( num_sal .EQ. 2 ) .OR. ( num_sal .EQ. 4 ) ) & 
    142144      smv_i(:,:,:)       = old_smv_i (:,:,:) 
     145 
     146      IF(ln_ctl) THEN   ! Control print 
     147         CALL prt_ctl(tab2d_1=area , clinfo1=' lim_itd_th  : cell area :') 
     148         CALL prt_ctl(tab2d_1=at_i , clinfo1=' lim_itd_th  : at_i      :') 
     149         CALL prt_ctl(tab2d_1=vt_i , clinfo1=' lim_itd_th  : vt_i      :') 
     150         CALL prt_ctl(tab2d_1=vt_s , clinfo1=' lim_itd_th  : vt_s      :') 
     151         DO jl = 1, jpl 
     152            CALL prt_ctl_info(' - Category : ', ivar1=jl) 
     153            CALL prt_ctl_info('   ~~~~~~~~~~') 
     154            CALL prt_ctl(tab2d_1=a_i   (:,:,jl)   , clinfo1= ' lim_itd_th  : a_i      : ') 
     155            CALL prt_ctl(tab2d_1=ht_i  (:,:,jl)   , clinfo1= ' lim_itd_th  : ht_i     : ') 
     156            CALL prt_ctl(tab2d_1=ht_s  (:,:,jl)   , clinfo1= ' lim_itd_th  : ht_s     : ') 
     157            CALL prt_ctl(tab2d_1=v_i   (:,:,jl)   , clinfo1= ' lim_itd_th  : v_i      : ') 
     158            CALL prt_ctl(tab2d_1=v_s   (:,:,jl)   , clinfo1= ' lim_itd_th  : v_s      : ') 
     159            CALL prt_ctl(tab2d_1=e_s   (:,:,1,jl) , clinfo1= ' lim_itd_th  : e_s      : ') 
     160            CALL prt_ctl(tab2d_1=t_su  (:,:,jl)   , clinfo1= ' lim_itd_th  : t_su     : ') 
     161            CALL prt_ctl(tab2d_1=t_s   (:,:,1,jl) , clinfo1= ' lim_itd_th  : t_snow   : ') 
     162            CALL prt_ctl(tab2d_1=sm_i  (:,:,jl)   , clinfo1= ' lim_itd_th  : sm_i     : ') 
     163            CALL prt_ctl(tab2d_1=smv_i (:,:,jl)   , clinfo1= ' lim_itd_th  : smv_i    : ') 
     164            DO ja = 1, nlay_i 
     165               CALL prt_ctl_info(' - Layer : ', ivar1=ja) 
     166               CALL prt_ctl_info('   ~~~~~~~') 
     167               CALL prt_ctl(tab2d_1=t_i(:,:,ja,jl) , clinfo1= ' lim_itd_th  : t_i      : ') 
     168               CALL prt_ctl(tab2d_1=e_i(:,:,ja,jl) , clinfo1= ' lim_itd_th  : e_i      : ') 
     169            END DO 
     170         END DO 
     171      ENDIF 
     172       
    143173 
    144174      END SUBROUTINE lim_itd_th 
  • trunk/NEMO/LIM_SRC_3/limthd.F90

    r834 r863  
    2929   USE par_ice 
    3030   USE limvar 
     31   USE prtctl          ! Print control 
    3132 
    3233   IMPLICIT NONE 
     
    4243         epsi06 = 1e-06   ,  & 
    4344         epsi04 = 1e-04   ,  & 
    44          zzero  = 0.e0     ,  & 
     45         zzero  = 0.e0     , & 
    4546         zone   = 1.e0 
    4647 
     
    473474 
    474475      IF ( con_i ) fbif(:,:) = fbif(:,:) + zqlbsbq(:,:) 
     476 
     477      IF(ln_ctl) THEN   ! Control print 
     478         CALL prt_ctl(tab2d_1=area , clinfo1=' lim_thd  : cell area :') 
     479         CALL prt_ctl(tab2d_1=at_i , clinfo1=' lim_thd  : at_i      :') 
     480         CALL prt_ctl(tab2d_1=vt_i , clinfo1=' lim_thd  : vt_i      :') 
     481         CALL prt_ctl(tab2d_1=vt_s , clinfo1=' lim_thd  : vt_s      :') 
     482         DO jl = 1, jpl 
     483            CALL prt_ctl_info(' - Category : ', ivar1=jl) 
     484            CALL prt_ctl_info('   ~~~~~~~~~~') 
     485            CALL prt_ctl(tab2d_1=a_i   (:,:,jl)   , clinfo1= ' lim_thd  : a_i      : ') 
     486            CALL prt_ctl(tab2d_1=ht_i  (:,:,jl)   , clinfo1= ' lim_thd  : ht_i     : ') 
     487            CALL prt_ctl(tab2d_1=ht_s  (:,:,jl)   , clinfo1= ' lim_thd  : ht_s     : ') 
     488            CALL prt_ctl(tab2d_1=v_i   (:,:,jl)   , clinfo1= ' lim_thd  : v_i      : ') 
     489            CALL prt_ctl(tab2d_1=v_s   (:,:,jl)   , clinfo1= ' lim_thd  : v_s      : ') 
     490            CALL prt_ctl(tab2d_1=e_s   (:,:,1,jl) , clinfo1= ' lim_thd  : e_s      : ') 
     491            CALL prt_ctl(tab2d_1=t_su  (:,:,jl)   , clinfo1= ' lim_thd  : t_su     : ') 
     492            CALL prt_ctl(tab2d_1=t_s   (:,:,1,jl) , clinfo1= ' lim_thd  : t_snow   : ') 
     493            CALL prt_ctl(tab2d_1=sm_i  (:,:,jl)   , clinfo1= ' lim_thd  : sm_i     : ') 
     494            CALL prt_ctl(tab2d_1=smv_i (:,:,jl)   , clinfo1= ' lim_thd  : smv_i    : ') 
     495            DO jk = 1, nlay_i 
     496               CALL prt_ctl_info(' - Layer : ', ivar1=jk) 
     497               CALL prt_ctl_info('   ~~~~~~~') 
     498               CALL prt_ctl(tab2d_1=t_i(:,:,jk,jl) , clinfo1= ' lim_thd  : t_i      : ') 
     499               CALL prt_ctl(tab2d_1=e_i(:,:,jk,jl) , clinfo1= ' lim_thd  : e_i      : ') 
     500            END DO 
     501         END DO 
     502 
     503      ENDIF 
    475504 
    476505   END SUBROUTINE lim_thd 
  • trunk/NEMO/LIM_SRC_3/limtrp.F90

    r834 r863  
    2626   USE lib_mpp 
    2727   USE par_ice 
     28   USE prtctl          ! Print control 
    2829 
    2930   IMPLICIT NONE 
     
    552553      ENDIF 
    553554 
     555      IF(ln_ctl) THEN   ! Control print 
     556         CALL prt_ctl(tab2d_1=area , clinfo1=' lim_trp  : cell area :') 
     557         CALL prt_ctl(tab2d_1=at_i , clinfo1=' lim_trp  : at_i      :') 
     558         CALL prt_ctl(tab2d_1=vt_i , clinfo1=' lim_trp  : vt_i      :') 
     559         CALL prt_ctl(tab2d_1=vt_s , clinfo1=' lim_trp  : vt_s      :') 
     560         DO jl = 1, jpl 
     561            CALL prt_ctl_info(' - Category : ', ivar1=jl) 
     562            CALL prt_ctl_info('   ~~~~~~~~~~') 
     563            CALL prt_ctl(tab2d_1=a_i   (:,:,jl)   , clinfo1= ' lim_trp  : a_i      : ') 
     564            CALL prt_ctl(tab2d_1=ht_i  (:,:,jl)   , clinfo1= ' lim_trp  : ht_i     : ') 
     565            CALL prt_ctl(tab2d_1=ht_s  (:,:,jl)   , clinfo1= ' lim_trp  : ht_s     : ') 
     566            CALL prt_ctl(tab2d_1=v_i   (:,:,jl)   , clinfo1= ' lim_trp  : v_i      : ') 
     567            CALL prt_ctl(tab2d_1=v_s   (:,:,jl)   , clinfo1= ' lim_trp  : v_s      : ') 
     568            CALL prt_ctl(tab2d_1=e_s   (:,:,1,jl) , clinfo1= ' lim_trp  : e_s      : ') 
     569            CALL prt_ctl(tab2d_1=t_su  (:,:,jl)   , clinfo1= ' lim_trp  : t_su     : ') 
     570            CALL prt_ctl(tab2d_1=t_s   (:,:,1,jl) , clinfo1= ' lim_trp  : t_snow   : ') 
     571            CALL prt_ctl(tab2d_1=sm_i  (:,:,jl)   , clinfo1= ' lim_trp  : sm_i     : ') 
     572            CALL prt_ctl(tab2d_1=smv_i (:,:,jl)   , clinfo1= ' lim_trp  : smv_i    : ') 
     573            DO jk = 1, nlay_i 
     574               CALL prt_ctl_info(' - Layer : ', ivar1=jk) 
     575               CALL prt_ctl_info('   ~~~~~~~') 
     576               CALL prt_ctl(tab2d_1=t_i(:,:,jk,jl) , clinfo1= ' lim_trp  : t_i      : ') 
     577               CALL prt_ctl(tab2d_1=e_i(:,:,jk,jl) , clinfo1= ' lim_trp  : e_i      : ') 
     578            END DO 
     579         END DO 
     580      ENDIF 
     581 
    554582   END SUBROUTINE lim_trp 
    555583 
  • trunk/NEMO/LIM_SRC_3/limupdate.F90

    r834 r863  
    1818   !! * Modules used 
    1919   USE limistate 
    20    USE limrhg     ! ice rheology 
     20   USE limrhg          ! ice rheology 
    2121   USE lbclnk 
    2222 
    2323   USE dom_oce 
    24    USE oce  ! dynamics and tracers variables 
     24   USE oce             ! dynamics and tracers variables 
    2525   USE in_out_manager 
    2626   USE ice_oce         ! ice variables 
     
    4545   USE limitd_th 
    4646   USE limvar 
     47   USE prtctl          ! Print control 
    4748 
    4849 
     
    996997 
    997998 
     999      IF(ln_ctl) THEN   ! Control print 
     1000         CALL prt_ctl_info(' - Cell values : ') 
     1001         CALL prt_ctl_info('   ~~~~~~~~~~~~~ ') 
     1002         CALL prt_ctl(tab2d_1=area       , clinfo1=' lim_update  : cell area   :') 
     1003         CALL prt_ctl(tab2d_1=at_i       , clinfo1=' lim_update  : at_i        :') 
     1004         CALL prt_ctl(tab2d_1=vt_i       , clinfo1=' lim_update  : vt_i        :') 
     1005         CALL prt_ctl(tab2d_1=vt_s       , clinfo1=' lim_update  : vt_s        :') 
     1006         CALL prt_ctl(tab2d_1=strength   , clinfo1=' lim_update  : strength    :') 
     1007         CALL prt_ctl(tab2d_1=u_ice      , clinfo1=' lim_update  : u_ice       :', tab2d_2=v_ice      , clinfo2=' v_ice       :') 
     1008         CALL prt_ctl(tab2d_1=d_u_ice_dyn, clinfo1=' lim_update  : d_u_ice_dyn :', tab2d_2=d_v_ice_dyn, clinfo2=' d_v_ice_dyn :') 
     1009         CALL prt_ctl(tab2d_1=old_u_ice  , clinfo1=' lim_update  : old_u_ice   :', tab2d_2=old_v_ice  , clinfo2=' old_v_ice   :') 
     1010 
     1011         DO jl = 1, jpl 
     1012            CALL prt_ctl_info(' - Category : ', ivar1=jl) 
     1013            CALL prt_ctl_info('   ~~~~~~~~~~') 
     1014            CALL prt_ctl(tab2d_1=ht_i       (:,:,jl)        , clinfo1= ' lim_update  : ht_i        : ') 
     1015            CALL prt_ctl(tab2d_1=ht_s       (:,:,jl)        , clinfo1= ' lim_update  : ht_s        : ') 
     1016            CALL prt_ctl(tab2d_1=t_su       (:,:,jl)        , clinfo1= ' lim_update  : t_su        : ') 
     1017            CALL prt_ctl(tab2d_1=t_s        (:,:,1,jl)      , clinfo1= ' lim_update  : t_snow      : ') 
     1018            CALL prt_ctl(tab2d_1=sm_i       (:,:,jl)        , clinfo1= ' lim_update  : sm_i        : ') 
     1019            CALL prt_ctl(tab2d_1=o_i        (:,:,jl)        , clinfo1= ' lim_update  : o_i         : ') 
     1020            CALL prt_ctl(tab2d_1=a_i        (:,:,jl)        , clinfo1= ' lim_update  : a_i         : ') 
     1021            CALL prt_ctl(tab2d_1=old_a_i    (:,:,jl)        , clinfo1= ' lim_update  : old_a_i     : ') 
     1022            CALL prt_ctl(tab2d_1=d_a_i_trp  (:,:,jl)        , clinfo1= ' lim_update  : d_a_i_trp   : ') 
     1023            CALL prt_ctl(tab2d_1=d_a_i_thd  (:,:,jl)        , clinfo1= ' lim_update  : d_a_i_thd   : ') 
     1024            CALL prt_ctl(tab2d_1=v_i        (:,:,jl)        , clinfo1= ' lim_update  : v_i         : ') 
     1025            CALL prt_ctl(tab2d_1=old_v_i    (:,:,jl)        , clinfo1= ' lim_update  : old_v_i     : ') 
     1026            CALL prt_ctl(tab2d_1=d_v_i_trp  (:,:,jl)        , clinfo1= ' lim_update  : d_v_i_trp   : ') 
     1027            CALL prt_ctl(tab2d_1=d_v_i_thd  (:,:,jl)        , clinfo1= ' lim_update  : d_v_i_thd   : ') 
     1028            CALL prt_ctl(tab2d_1=v_s        (:,:,jl)        , clinfo1= ' lim_update  : v_s         : ') 
     1029            CALL prt_ctl(tab2d_1=old_v_s    (:,:,jl)        , clinfo1= ' lim_update  : old_v_s     : ') 
     1030            CALL prt_ctl(tab2d_1=d_v_s_trp  (:,:,jl)        , clinfo1= ' lim_update  : d_v_s_trp   : ') 
     1031            CALL prt_ctl(tab2d_1=d_v_s_thd  (:,:,jl)        , clinfo1= ' lim_update  : d_v_s_thd   : ') 
     1032            CALL prt_ctl(tab2d_1=e_i        (:,:,1,jl)/1.0e9, clinfo1= ' lim_update  : e_i1        : ') 
     1033            CALL prt_ctl(tab2d_1=old_e_i    (:,:,1,jl)/1.0e9, clinfo1= ' lim_update  : old_e_i1    : ') 
     1034            CALL prt_ctl(tab2d_1=d_e_i_trp  (:,:,1,jl)/1.0e9, clinfo1= ' lim_update  : de_i1_trp   : ') 
     1035            CALL prt_ctl(tab2d_1=d_e_i_thd  (:,:,1,jl)/1.0e9, clinfo1= ' lim_update  : de_i1_thd   : ') 
     1036            CALL prt_ctl(tab2d_1=e_i        (:,:,2,jl)/1.0e9, clinfo1= ' lim_update  : e_i2        : ') 
     1037            CALL prt_ctl(tab2d_1=old_e_i    (:,:,2,jl)/1.0e9, clinfo1= ' lim_update  : old_e_i2    : ') 
     1038            CALL prt_ctl(tab2d_1=d_e_i_trp  (:,:,2,jl)/1.0e9, clinfo1= ' lim_update  : de_i2_trp   : ') 
     1039            CALL prt_ctl(tab2d_1=d_e_i_thd  (:,:,2,jl)/1.0e9, clinfo1= ' lim_update  : de_i2_thd   : ') 
     1040            CALL prt_ctl(tab2d_1=e_s        (:,:,1,jl)      , clinfo1= ' lim_update  : e_snow      : ') 
     1041            CALL prt_ctl(tab2d_1=old_e_s    (:,:,1,jl)      , clinfo1= ' lim_update  : old_e_snow  : ') 
     1042            CALL prt_ctl(tab2d_1=d_e_s_trp  (:,:,1,jl)/1.0e9, clinfo1= ' lim_update  : d_e_s_trp   : ') 
     1043            CALL prt_ctl(tab2d_1=d_e_s_thd  (:,:,1,jl)/1.0e9, clinfo1= ' lim_update  : d_e_s_thd   : ') 
     1044            CALL prt_ctl(tab2d_1=smv_i      (:,:,jl)        , clinfo1= ' lim_update  : smv_i       : ') 
     1045            CALL prt_ctl(tab2d_1=old_smv_i  (:,:,jl)        , clinfo1= ' lim_update  : old_smv_i   : ') 
     1046            CALL prt_ctl(tab2d_1=d_smv_i_trp(:,:,jl)        , clinfo1= ' lim_update  : d_smv_i_trp : ') 
     1047            CALL prt_ctl(tab2d_1=d_smv_i_thd(:,:,jl)        , clinfo1= ' lim_update  : d_smv_i_thd : ') 
     1048            CALL prt_ctl(tab2d_1=oa_i       (:,:,jl)        , clinfo1= ' lim_update  : oa_i        : ') 
     1049            CALL prt_ctl(tab2d_1=old_oa_i   (:,:,jl)        , clinfo1= ' lim_update  : old_oa_i    : ') 
     1050            CALL prt_ctl(tab2d_1=d_oa_i_trp (:,:,jl)        , clinfo1= ' lim_update  : d_oa_i_trp  : ') 
     1051            CALL prt_ctl(tab2d_1=d_oa_i_thd (:,:,jl)        , clinfo1= ' lim_update  : d_oa_i_thd  : ') 
     1052            CALL prt_ctl(tab2d_1=REAL(patho_case(:,:,jl))   , clinfo1= ' lim_update  : Path. case  : ') 
     1053 
     1054            DO jk = 1, nlay_i 
     1055               CALL prt_ctl_info(' - Layer : ', ivar1=jk) 
     1056               CALL prt_ctl(tab2d_1=t_i(:,:,jk,jl) , clinfo1= ' lim_dyn  : t_i       : ') 
     1057            END DO 
     1058         END DO 
     1059 
     1060         CALL prt_ctl_info(' - Heat / FW fluxes : ') 
     1061         CALL prt_ctl_info('   ~~~~~~~~~~~~~~~~~~ ') 
     1062         CALL prt_ctl(tab2d_1=fmmec  , clinfo1= ' lim_update : fmmec : ', tab2d_2=fhmec     , clinfo2= ' fhmec     : ') 
     1063         CALL prt_ctl(tab2d_1=sst_io , clinfo1= ' lim_update : sst   : ', tab2d_2=sss_io    , clinfo2= ' sss       : ') 
     1064         CALL prt_ctl(tab2d_1=fhbri  , clinfo1= ' lim_update : fhbri : ', tab2d_2=fheat_rpo , clinfo2= ' fheat_rpo : ') 
     1065 
     1066         CALL prt_ctl_info(' - Stresses : ') 
     1067         CALL prt_ctl_info('   ~~~~~~~~~~ ') 
     1068         CALL prt_ctl(tab2d_1=tauxw , clinfo1= ' lim_update : tauxw : ', tab2d_2=tauyw , clinfo2= ' tauyw : ') 
     1069         CALL prt_ctl(tab2d_1=taux  , clinfo1= ' lim_update : taux  : ', tab2d_2=tauy  , clinfo2= ' tauy  : ') 
     1070         CALL prt_ctl(tab2d_1=ftaux , clinfo1= ' lim_update : ftaux : ', tab2d_2=ftauy , clinfo2= ' ftauy : ') 
     1071         CALL prt_ctl(tab2d_1=gtaux , clinfo1= ' lim_update : gtaux : ', tab2d_2=gtauy , clinfo2= ' gtauy : ') 
     1072         CALL prt_ctl(tab2d_1=u_io  , clinfo1= ' lim_update : u_io  : ', tab2d_2=v_io  , clinfo2= ' v_io  : ') 
     1073      ENDIF 
     1074 
     1075     !--------------------- 
     1076 
    9981077   END SUBROUTINE lim_update 
    9991078#else 
Note: See TracChangeset for help on using the changeset viewer.