Changeset 1482 for trunk/NEMO/OPA_SRC/SBC/sbcblk_core.F90

Timestamp:

2009-07-03T17:28:06+02:00 (15 years ago)

Author:

smasson

Message:

distribution of iom_put + cleaning of LIM2 outputs, see ticket:437

File:

• trunk/NEMO/OPA_SRC/SBC/sbcblk_core.F90 (modified) (7 diffs)

trunk/NEMO/OPA_SRC/SBC/sbcblk_core.F90

@@ -4 +4 @@
    !! Ocean forcing:  momentum, heat and freshwater flux formulation
    !!=====================================================================
-   !! History :  1.0   !  04-08  (U. Schweckendiek)  Original code
-   !!            2.0   !  05-04  (L. Brodeau, A.M. Treguier) additions: 
-   !!                            -  new bulk routine for efficiency
-   !!                            -  WINDS ARE NOW ASSUMED TO BE AT T POINTS in input files !!!!
-   !!                            -  file names and file characteristics in namelist 
-   !!                            -  Implement reading of 6-hourly fields   
-   !!            3.0   !  06-06  (G. Madec) sbc rewritting
+   !! History :  1.0  !  2004-08  (U. Schweckendiek)  Original code
+   !!            2.0  !  2005-04  (L. Brodeau, A.M. Treguier) additions: 
+   !!                           -  new bulk routine for efficiency
+   !!                           -  WINDS ARE NOW ASSUMED TO BE AT T POINTS in input files !!!!
+   !!                           -  file names and file characteristics in namelist 
+   !!                           -  Implement reading of 6-hourly fields   
+   !!            3.0  !  2006-06  (G. Madec) sbc rewritting   
+   !!            3.2  !  2009-04  (B. Lemaire)  Introduce iom_put
    !!----------------------------------------------------------------------
 
@@ -68 +69 @@
 #  include "vectopt_loop_substitute.h90"
    !!----------------------------------------------------------------------
-   !!   OPA 9.0 , LOCEAN-IPSL (2006) 
+   !! NEMO/OPA 3.2 , LOCEAN-IPSL (2009) 
    !! $Id$
    !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
@@ -190 +191 @@
       !!              - qns     : Non Solar heat flux over the ocean    (W/m2)
       !!              - evap    : Evaporation over the ocean            (kg/m2/s)
-      !!              - tprecip : Total precipitation                   (Kg/m2/s)
-      !!              - sprecip : Solid precipitation                   (Kg/m2/s)
+      !!              - emp(s)  : evaporation minus precipitation       (kg/m2/s)
       !!
       !!  ** Nota  :   sf has to be a dummy argument for AGRIF on NEC
@@ -329 +329 @@
       emp (:,:) = zevap(:,:) - sf(jp_prec)%fnow(:,:) * alpha_precip * tmask(:,:,1)
 !CDIR COLLAPSE
-      emps(:,:) = zevap(:,:) - sf(jp_prec)%fnow(:,:) * alpha_precip * tmask(:,:,1)
+      emps(:,:) = emp(:,:)
+      !
+      CALL iom_put( "qlw_oce",   zqlw )                 ! output downward longwave heat over the ocean
+      CALL iom_put( "qsb_oce", - zqsb )                 ! output downward sensible heat over the ocean
+      CALL iom_put( "qla_oce", - zqla )                 ! output downward latent   heat over the ocean
+      CALL iom_put( "qns_oce",   qns  )                 ! output downward non solar heat over the ocean
+      !
+      IF(ln_ctl) THEN
+         CALL prt_ctl(tab2d_1=zqsb , clinfo1=' blk_oce_core: zqsb   : ', tab2d_2=zqlw , clinfo2=' zqlw  : ')
+         CALL prt_ctl(tab2d_1=zqla , clinfo1=' blk_oce_core: zqla   : ', tab2d_2=qsr  , clinfo2=' qsr   : ')
+         CALL prt_ctl(tab2d_1=pst  , clinfo1=' blk_oce_core: pst    : ', tab2d_2=emp  , clinfo2=' emp   : ')
+         CALL prt_ctl(tab2d_1=utau , clinfo1=' blk_oce_core: utau   : ', mask1=umask,   &
+            &         tab2d_2=vtau , clinfo2=              ' vtau  : ' , mask2=vmask )
+      ENDIF
       !
    END SUBROUTINE blk_oce_core
@@ -377 +390 @@
       REAL(wp) ::             zwndi_t , zwndj_t                  ! relative wind components at T-point
       REAL(wp), DIMENSION(jpi,jpj) ::   z_wnds_t                 ! wind speed ( = | U10m - U_ice | ) at T-point
-      REAL(wp), DIMENSION(jpi,jpj,pdim) ::   z_qlw                ! long wave heat flux over ice
-      REAL(wp), DIMENSION(jpi,jpj,pdim) ::   z_qsb                ! sensible  heat flux over ice
-      REAL(wp), DIMENSION(jpi,jpj,pdim) ::   z_dqlw               ! sensible  heat flux over ice
-      REAL(wp), DIMENSION(jpi,jpj,pdim) ::   z_dqsb               ! sensible  heat flux over ice
+      REAL(wp), DIMENSION(jpi,jpj,pdim) ::   z_qlw               ! long wave heat flux over ice
+      REAL(wp), DIMENSION(jpi,jpj,pdim) ::   z_qsb               ! sensible  heat flux over ice
+      REAL(wp), DIMENSION(jpi,jpj,pdim) ::   z_dqlw              ! long wave heat sensitivity over ice
+      REAL(wp), DIMENSION(jpi,jpj,pdim) ::   z_dqsb              ! sensible  heat sensitivity over ice
       !!---------------------------------------------------------------------
 
@@ -523 +536 @@
 !CDIR COLLAPSE
       p_spr(:,:) = sf(jp_snow)%fnow(:,:) * alpha_precip      ! solid precipitation [kg/m2/s]
+      CALL iom_put( 'snowpre', p_spr )                       ! Snow precipitation 
       !
       IF(ln_ctl) THEN
@@ -680 +694 @@
       !!   9.0  !  06-12  (L. Brodeau) Original code for 2Z
       !!----------------------------------------------------------------------
-      !! * Arguments
       REAL(wp), INTENT(in)   :: &
          zt,      &     ! height for T_zt and q_zt                   [m]