Changeset 258 for trunk/NEMO/OPA_SRC/SBC

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/OPA_SRC/SBC

Files:

• bulk.F90 (modified) (2 diffs)
• flxblk.F90 (modified) (2 diffs)
• ocesbc.F90 (modified) (5 diffs)

trunk/NEMO/OPA_SRC/SBC/bulk.F90

@@ -21 +21 @@
    USE phycst          ! physical constants
    USE in_out_manager  ! I/O manager
+   USE prtctl          ! Print control
 
    IMPLICIT NONE
@@ -97 +98 @@
 
 # if ! defined key_ice_lim
-         IF(l_ctl) THEN         ! print mean trends (used for debugging)
-            WRITE(numout,*) ' Forcings '
-            WRITE(numout,*) ' qsr_oce  : ', SUM( qsr_oce (1:nictl+1,1:njctl+1) * tmask(1:nictl+1,1:njctl+1,1) )
-            WRITE(numout,*) ' qsr_ice  : ', SUM( qsr_ice (1:nictl+1,1:njctl+1) * tmask(1:nictl+1,1:njctl+1,1) )
-            WRITE(numout,*) ' qnsr_oce : ', SUM( qnsr_oce(1:nictl+1,1:njctl+1) * tmask(1:nictl+1,1:njctl+1,1) )
-            WRITE(numout,*) ' qnsr_ice : ', SUM( qnsr_ice(1:nictl+1,1:njctl+1) * tmask(1:nictl+1,1:njctl+1,1) )
-            WRITE(numout,*) ' evap     : ', SUM( evap    (1:nictl+1,1:njctl+1) * tmask(1:nictl+1,1:njctl+1,1) )
-            WRITE(numout,*) ' precip   : ', SUM( tprecip (1:nictl+1,1:njctl+1) * tmask(1:nictl+1,1:njctl+1,1) ) / rday
-            WRITE(numout,*) ' Snow     : ', SUM( sprecip (1:nictl+1,1:njctl+1) * tmask(1:nictl+1,1:njctl+1,1) ) / rday
-            WRITE(numout,*) ' u-stress : ', SUM( taux    (1:nictl+1,1:njctl+1) * umask(1:nictl+1,1:njctl+1,1) )
-            WRITE(numout,*) ' v-stress : ', SUM( tauy    (1:nictl+1,1:njctl+1) * vmask(1:nictl+1,1:njctl+1,1) )
-            WRITE(numout,*) ' sst      : ', SUM( zsst    (1:nictl+1,1:njctl+1) * tmask(1:nictl+1,1:njctl+1,1) )
-            WRITE(numout,*) ' sss      : ', SUM( zsss    (1:nictl+1,1:njctl+1) * tmask(1:nictl+1,1:njctl+1,1) )
-            WRITE(numout,*)
+         IF(ln_ctl) THEN         ! print mean trends (used for debugging)
+            CALL prt_ctl_info(' Forcings ')
+            CALL prt_ctl(tab2d_1=qsr_oce , clinfo1=' qsr_oce   : ', mask1=tmask, ovlap=1)
+            CALL prt_ctl(tab2d_1=qsr_ice , clinfo1=' qsr_ice   : ', mask1=tmask, ovlap=1)
+            CALL prt_ctl(tab2d_1=qnsr_oce, clinfo1=' qnsr_oce  : ', mask1=tmask, ovlap=1)
+            CALL prt_ctl(tab2d_1=qnsr_ice, clinfo1=' qnsr_ice  : ', mask1=tmask, ovlap=1)
+            CALL prt_ctl(tab2d_1=evap    , clinfo1=' evap      : ', mask1=tmask, ovlap=1)
+            CALL prt_ctl(tab2d_1=tprecip , clinfo1=' precip    : ', mask1=tmask, ovlap=1)
+            CALL prt_ctl(tab2d_1=sprecip , clinfo1=' Snow      : ', mask1=tmask, ovlap=1)
+            CALL prt_ctl(tab2d_1=taux    , clinfo1=' u-stress  : ', mask1=tmask, ovlap=1)
+            CALL prt_ctl(tab2d_1=tauy    , clinfo1=' v-stress  : ', mask1=tmask, ovlap=1)
+            CALL prt_ctl(tab2d_1=zsst    , clinfo1=' sst       : ', mask1=tmask, ovlap=1)
+            CALL prt_ctl(tab2d_1=zsss    , clinfo1=' sss       : ', mask1=tmask, ovlap=1)
          ENDIF
 # endif   

trunk/NEMO/OPA_SRC/SBC/flxblk.F90

@@ -24 +24 @@
    USE lbclnk
    USE albedo
+   USE prtctl          ! Print control
 
    IMPLICIT NONE
@@ -319 +320 @@
 
       iday   = INT( zxday )
-      IF(l_ctl)   WRITE(numout,*) ' declin : iday ', iday, ' nfbulk= ', nfbulk
+
+      IF(ln_ctl) CALL prt_ctl_info('declin : iday ', ivar1=iday, clinfo2=' nfbulk= ', ivar2=nfbulk)
+
       !   computation of the solar declination, his sine and his cosine
       CALL flx_blk_declin( indaet, iday, zdecl )

trunk/NEMO/OPA_SRC/SBC/ocesbc.F90

@@ -13 +13 @@
    USE dom_oce        ! ocean space domain variables
    USE cpl_oce        ! coupled ocean-atmosphere variables
-   USE ice_oce
-   USE blk_oce
-   USE flx_oce
+   USE ice_oce        ! sea-ice variable
+   USE blk_oce        ! bulk variables
+   USE flx_oce        ! sea-ice/ocean forcings variables
    USE phycst         ! Define parameters for the routines
-   USE taumod
-   USE flxmod
-   USE flxrnf
+   USE taumod         ! surface stress forcing
+   USE flxmod         ! thermohaline fluxes
+   USE flxrnf         ! runoffs forcing
    USE tradmp         ! damping salinity trend
-   USE dtatem
-   USE dtasal
-   USE ocfzpt
-   USE lbclnk
-   USE lib_mpp
+   USE dtatem         ! ocean temperature data
+   USE dtasal         ! ocean salinity data
+   USE ocfzpt         ! surface ocean freezing point
+   USE lbclnk         ! ocean lateral boundary condition
+   USE lib_mpp        ! distribued memory computing library
    USE in_out_manager ! I/O manager
+   USE prtctl         ! Print control
 
    IMPLICIT NONE
@@ -301 +302 @@
       !!----------------------------------------------------------------------
       !! * Modules used
-      USE cpl_oce 
+      USE cpl_oce       ! coupled ocean-atmosphere variables
 
       !! * Arguments
@@ -400 +401 @@
       !!----------------------------------------------------------------------
       !! * Modules used
-      USE daymod
+      USE daymod                 ! calendar
 #if ! defined key_dtasst
-      USE dtasst, ONLY : rclice
+      USE dtasst, ONLY : rclice  ! sea surface temperature data
 #endif
 #if defined key_flx_bulk_monthly || defined key_flx_bulk_daily
-      USE blk_oce
+      USE blk_oce                ! bulk variables
 #endif
 #if defined key_flx_forced_daily
-      USE flx_oce
+      USE flx_oce                ! sea-ice/ocean forcings variables
 #endif
 
@@ -769 +770 @@
       REAL(wp), DIMENSION(jpi,jpj)  :: zsss, zfreeze
       REAL(wp) ::   zerp, ztrp, zsrp
+      CHARACTER (len=71) :: charout
 #if defined key_dynspg_fsc
       REAL(wp) ::   zwei
@@ -852 +854 @@
       IF( lk_mpp )   CALL mpp_sum( aplus  )   ! sums over the global domain
       IF( lk_mpp )   CALL mpp_sum( aminus )
-      IF(l_ctl)   WRITE(numout,*) ' oce_sbc_dmp : a+ = ', aplus, ' a- = ', aminus
+
+      IF(ln_ctl)   THEN
+         WRITE(charout,FMT="('oce_sbc_dmp : a+ = ',D23.16, ' a- = ',D23.16)") aplus, aminus
+         CALL prt_ctl_info(charout)
+      ENDIF
 
       zadefi = MIN( aplus, aminus )