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

Timestamp:

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

Author:

ctlod

Message:

Add control prints

File:

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

trunk/NEMO/LIM_SRC_3/limthd.F90

@@ -29 +29 @@
    USE par_ice
    USE limvar
+   USE prtctl          ! Print control
 
    IMPLICIT NONE
@@ -42 +43 @@
          epsi06 = 1e-06   ,  &
          epsi04 = 1e-04   ,  &
-         zzero  = 0.e0     ,  &
+         zzero  = 0.e0     , &
          zone   = 1.e0
 
@@ -473 +474 @@
 
       IF ( con_i ) fbif(:,:) = fbif(:,:) + zqlbsbq(:,:)
+
+      IF(ln_ctl) THEN   ! Control print
+         CALL prt_ctl(tab2d_1=area , clinfo1=' lim_thd  : cell area :')
+         CALL prt_ctl(tab2d_1=at_i , clinfo1=' lim_thd  : at_i      :')
+         CALL prt_ctl(tab2d_1=vt_i , clinfo1=' lim_thd  : vt_i      :')
+         CALL prt_ctl(tab2d_1=vt_s , clinfo1=' lim_thd  : vt_s      :')
+         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_thd  : a_i      : ')
+            CALL prt_ctl(tab2d_1=ht_i  (:,:,jl)   , clinfo1= ' lim_thd  : ht_i     : ')
+            CALL prt_ctl(tab2d_1=ht_s  (:,:,jl)   , clinfo1= ' lim_thd  : ht_s     : ')
+            CALL prt_ctl(tab2d_1=v_i   (:,:,jl)   , clinfo1= ' lim_thd  : v_i      : ')
+            CALL prt_ctl(tab2d_1=v_s   (:,:,jl)   , clinfo1= ' lim_thd  : v_s      : ')
+            CALL prt_ctl(tab2d_1=e_s   (:,:,1,jl) , clinfo1= ' lim_thd  : e_s      : ')
+            CALL prt_ctl(tab2d_1=t_su  (:,:,jl)   , clinfo1= ' lim_thd  : t_su     : ')
+            CALL prt_ctl(tab2d_1=t_s   (:,:,1,jl) , clinfo1= ' lim_thd  : t_snow   : ')
+            CALL prt_ctl(tab2d_1=sm_i  (:,:,jl)   , clinfo1= ' lim_thd  : sm_i     : ')
+            CALL prt_ctl(tab2d_1=smv_i (:,:,jl)   , clinfo1= ' lim_thd  : smv_i    : ')
+            DO jk = 1, nlay_i
+               CALL prt_ctl_info(' - Layer : ', ivar1=jk)
+               CALL prt_ctl_info('   ~~~~~~~')
+               CALL prt_ctl(tab2d_1=t_i(:,:,jk,jl) , clinfo1= ' lim_thd  : t_i      : ')
+               CALL prt_ctl(tab2d_1=e_i(:,:,jk,jl) , clinfo1= ' lim_thd  : e_i      : ')
+            END DO
+         END DO
+
+      ENDIF
 
    END SUBROUTINE lim_thd