Changeset 3294 for trunk/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_mfs.F90

Timestamp:

2012-01-28T17:44:18+01:00 (12 years ago)

Author:

rblod

Message:

Merge of 3.4beta into the trunk

File:

• trunk/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_mfs.F90 (copied) (copied from branches/2011/dev_LOCEAN_CMCC_INGV_MERCATOR_2011/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_mfs.F90) (7 diffs)

trunk/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_mfs.F90

@@ -20 +20 @@
    USE in_out_manager  ! I/O manager
    USE lib_mpp         ! distribued memory computing library
+   USE wrk_nemo        ! work arrays
+   USE timing          ! Timing
    USE lbclnk          ! ocean lateral boundary conditions (or mpp link)
    USE prtctl          ! Print control
@@ -119 +121 @@
          &                  sn_tair , sn_rhm , sn_prec 
       !!---------------------------------------------------------------------
-
+      !
+      IF( nn_timing == 1 )  CALL timing_start('sbc_blk_mfs')
+      !
       !                                         ! ====================== !
       IF( kt == nit000 ) THEN                   !  First call kt=nit000  !
@@ -262 +266 @@
 
       ENDIF
-
+      !
+      IF( nn_timing == 1 )  CALL timing_stop('sbc_blk_mfs')
+      !
    END SUBROUTINE sbc_blk_mfs
   
@@ -283 +289 @@
       !!
       !!----------------------------------------------------------------------
-      !!
-
       USE sbcblk_core, ONLY: turb_core_2z ! For wave coupling and Tair/rh from 2 to 10m
-      USE wrk_nemo, ONLY:  wrk_in_use, wrk_not_released
-      USE wrk_nemo, ONLY:  rspeed     => wrk_2d_1
-      USE wrk_nemo, ONLY:  sh10now    => wrk_2d_2
-      USE wrk_nemo, ONLY:  t10now     => wrk_2d_3
-      USE wrk_nemo, ONLY:  cdx        => wrk_2d_4   ! --- drag coeff.
-      USE wrk_nemo, ONLY:  ce         => wrk_2d_5   ! --- turbulent exchange coefficients
-      USE wrk_nemo, ONLY:  shms       => wrk_2d_6
-      USE wrk_nemo, ONLY:  rhom       => wrk_2d_7
-      USE wrk_nemo, ONLY:  sstk       => wrk_2d_8
-      USE wrk_nemo, ONLY:  ch         => wrk_2d_10
-      USE wrk_nemo, ONLY:  rel_windu  => wrk_2d_11
-      USE wrk_nemo, ONLY:  rel_windv  => wrk_2d_12
 
       REAL(wp), INTENT(in   ) :: hour
@@ -310 +302 @@
       REAL(wp)  :: esre, cseep
 
+      REAL(wp), DIMENSION (:,:), POINTER ::   rspeed, sh10now, t10now, cdx, ce, shms
+      REAL(wp), DIMENSION (:,:), POINTER ::   rhom, sstk, ch, rel_windu, rel_windv
       !!----------------------------------------------------------------------
       !!     coefficients ( in MKS )  :
@@ -335 +329 @@
       REAL(wp), DIMENSION(5) :: p_e = (/-0.16,1.0,1.0,1.0,1.0/)
       INTEGER :: kku                        !index varing with wind speed
-
-      ! Set-up access to workspace arrays
-      IF( wrk_in_use(2, 1,2,3,4,5,6,7,8,10,11,12) ) THEN
-         CALL ctl_stop('blk_mfs: requested workspace arrays unavailable')   ;   RETURN
-      END IF
+      !
+      IF( nn_timing == 1 )  CALL timing_start('fluxes_mfs')
+      !
+      CALL wrk_alloc( jpi,jpj, rspeed, sh10now, t10now, cdx, ce, shms )
+      CALL wrk_alloc( jpi,jpj, rhom, sstk, ch, rel_windu, rel_windv )
 
       !!----------------------------------------------------------------------
@@ -501 +495 @@
        tauy(:,:)= rhom(:,:) * cdx(:,:) * rspeed(:,:) * rel_windv(:,:)
 
-
-      IF( wrk_not_released(2, 1,2,3,4,5,6,7,8,10,11,12) )    &
-         CALL ctl_stop('fluxes_mfs: failed to release workspace arrays')
-
-
+      CALL wrk_dealloc( jpi,jpj, rspeed, sh10now, t10now, cdx, ce, shms )
+      CALL wrk_dealloc( jpi,jpj, rhom, sstk, ch, rel_windu, rel_windv )
+      !
+      IF( nn_timing == 1 )  CALL timing_stop('fluxes_mfs')
+      !
    END SUBROUTINE fluxes_mfs