Changeset 863 for trunk/NEMO/LIM_SRC_3/limdyn.F90

Timestamp:

2008-03-14T14:34:10+01:00 (16 years ago)

Author:

ctlod

Message:

Add control prints

File:

• trunk/NEMO/LIM_SRC_3/limdyn.F90 (modified) (3 diffs)

trunk/NEMO/LIM_SRC_3/limdyn.F90

@@ -18 +18 @@
    USE taumod
    USE ice
+   USE par_ice
    USE ice_oce
    USE iceini
@@ -62 +63 @@
       !!------------------------------------------------------------------------------------
       !! * Local variables
-      INTEGER ::   ji, jj            ! dummy loop indices
+      INTEGER ::   ji, jj, jl, ja    ! dummy loop indices
       INTEGER :: i_j1, i_jpj         ! Starting/ending j-indices for rheology
 
@@ -269 +270 @@
       CALL lbc_lnk( tio_v, 'V', -1. )   ! I-point (i.e. ice U-V point)
 
-      IF(ln_ctl) THEN 
-        CALL prt_ctl(tab2d_1=tio_u , clinfo1=' lim_dyn  : tio_u :', tab2d_2=tio_v , clinfo2=' tio_v :')
-        CALL prt_ctl(tab2d_1=ust2s , clinfo1=' lim_dyn  : ust2s :')
-        CALL prt_ctl(tab2d_1=at_i  , clinfo1=' lim_dyn  : at_i  :')
+      IF(ln_ctl) THEN   ! Control print
+         CALL prt_ctl(tab2d_1=tio_u     , clinfo1=' lim_dyn  : tio_u     :', tab2d_2=tio_v , clinfo2=' tio_v :')
+         CALL prt_ctl(tab2d_1=ust2s     , clinfo1=' lim_dyn  : ust2s     :')
+         CALL prt_ctl(tab2d_1=divu_i    , clinfo1=' lim_dyn  : divu_i    :')
+         CALL prt_ctl(tab2d_1=delta_i   , clinfo1=' lim_dyn  : delta_i   :')
+         CALL prt_ctl(tab2d_1=strength  , clinfo1=' lim_dyn  : strength  :')
+         CALL prt_ctl(tab2d_1=area      , clinfo1=' lim_dyn  : cell area :')
+         CALL prt_ctl(tab2d_1=at_i      , clinfo1=' lim_dyn  : at_i      :')
+         CALL prt_ctl(tab2d_1=vt_i      , clinfo1=' lim_dyn  : vt_i      :')
+         CALL prt_ctl(tab2d_1=vt_s      , clinfo1=' lim_dyn  : vt_s      :')
+         CALL prt_ctl(tab2d_1=stress1_i , clinfo1=' lim_dyn  : stress1_i :')
+         CALL prt_ctl(tab2d_1=stress2_i , clinfo1=' lim_dyn  : stress2_i :')
+         CALL prt_ctl(tab2d_1=stress12_i, clinfo1=' lim_dyn  : stress12_i:')
+         DO jl = 1, jpl
+            CALL prt_ctl_info(' - Category : ', ivar1=jl)
+            CALL prt_ctl_info('   ~~~~~~~~~~')
+            CALL prt_ctl(tab2d_1=a_i   (:,:,jl)   , clinfo1= ' lim_dyn  : a_i      : ')
+            CALL prt_ctl(tab2d_1=ht_i  (:,:,jl)   , clinfo1= ' lim_dyn  : ht_i     : ')
+            CALL prt_ctl(tab2d_1=ht_s  (:,:,jl)   , clinfo1= ' lim_dyn  : ht_s     : ')
+            CALL prt_ctl(tab2d_1=v_i   (:,:,jl)   , clinfo1= ' lim_dyn  : v_i      : ')
+            CALL prt_ctl(tab2d_1=v_s   (:,:,jl)   , clinfo1= ' lim_dyn  : v_s      : ')
+            CALL prt_ctl(tab2d_1=e_s   (:,:,1,jl) , clinfo1= ' lim_dyn  : e_s      : ')
+            CALL prt_ctl(tab2d_1=t_su  (:,:,jl)   , clinfo1= ' lim_dyn  : t_su     : ')
+            CALL prt_ctl(tab2d_1=t_s   (:,:,1,jl) , clinfo1= ' lim_dyn  : t_snow   : ')
+            CALL prt_ctl(tab2d_1=sm_i  (:,:,jl)   , clinfo1= ' lim_dyn  : sm_i     : ')
+            CALL prt_ctl(tab2d_1=smv_i (:,:,jl)   , clinfo1= ' lim_dyn  : smv_i    : ')
+            DO ja = 1, nlay_i
+               CALL prt_ctl_info(' - Layer : ', ivar1=ja)
+               CALL prt_ctl_info('   ~~~~~~~')
+               CALL prt_ctl(tab2d_1=t_i(:,:,ja,jl) , clinfo1= ' lim_dyn  : t_i      : ')
+               CALL prt_ctl(tab2d_1=e_i(:,:,ja,jl) , clinfo1= ' lim_dyn  : e_i      : ')
+            END DO
+         END DO
       ENDIF