Changeset 3294 for trunk/NEMOGCM/NEMO/LIM_SRC_2/limsbc_2.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/LIM_SRC_2/limsbc_2.F90 (modified) (9 diffs)

trunk/NEMOGCM/NEMO/LIM_SRC_2/limsbc_2.F90

@@ -27 +27 @@
    USE sbc_ice          ! surface boundary condition: ice
    USE sbc_oce          ! surface boundary condition: ocean
+   USE sbccpl
 
    USE albedo           ! albedo parameters
    USE lbclnk           ! ocean lateral boundary condition - MPP exchanges
    USE lib_mpp          ! MPP library
+   USE wrk_nemo         ! work arrays
    USE in_out_manager   ! I/O manager
    USE diaar5, ONLY :   lk_diaar5
@@ -50 +52 @@
    !
    REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) ::   soce_0, sice_0   ! constant SSS and ice salinity used in levitating sea-ice case
-
    REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) ::   utau_oce, vtau_oce   ! air-ocean surface i- & j-stress              [N/m2]
    REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) ::   tmod_io              ! modulus of the ice-ocean relative velocity   [m/s]
@@ -100 +101 @@
       !!              Tartinville et al. 2001 Ocean Modelling, 3, 95-108.
       !!---------------------------------------------------------------------
-      USE wrk_nemo, ONLY: wrk_not_released, wrk_in_use
-      USE wrk_nemo, ONLY: zqnsoce => wrk_2d_1 ! 2D workspace
-      USE wrk_nemo, ONLY: wrk_3d_4, wrk_3d_5
       INTEGER, INTENT(in) ::   kt    ! number of iteration
       !!
@@ -111 +109 @@
       REAL(wp) ::   zqsr, zqns, zfm            ! local scalars
       REAL(wp) ::   zinda, zfons, zemp         !   -      -
+      REAL(wp), POINTER, DIMENSION(:,:)   ::   zqnsoce       ! 2D workspace
       REAL(wp), POINTER, DIMENSION(:,:,:) ::   zalb, zalbp   ! 2D/3D workspace
       !!---------------------------------------------------------------------
      
-      IF( wrk_in_use(2, 1) .OR. wrk_in_use(3, 4,5) )THEN
-         CALL ctl_stop('lim_sbc_flx_2 : requested workspace arrays unavailable')   ;   RETURN
-      ENDIF
-      zalb  => wrk_3d_4(:,:,1:1)      ! Set-up pointers to sub-arrays of 3d workspaces
-      zalbp => wrk_3d_5(:,:,1:1)
+      CALL wrk_alloc( jpi, jpj, zqnsoce )
+      CALL wrk_alloc( jpi, jpj, 1, zalb, zalbp )
 
       !------------------------------------------!
@@ -234 +230 @@
       !-----------------------------------------------!
 
-      IF( lk_cpl ) THEN          ! coupled case
-         tn_ice(:,:,1) = sist(:,:)          ! sea-ice surface temperature       
-         !                                  ! Computation of snow/ice and ocean albedo
-         CALL albedo_ice( tn_ice, reshape( hicif, (/jpi,jpj,1/) ), reshape( hsnif, (/jpi,jpj,1/) ), zalbp, zalb )
-         alb_ice(:,:,1) =  0.5 * ( zalbp(:,:,1) + zalb (:,:,1) )   ! Ice albedo (mean clear and overcast skys)
-         CALL iom_put( "icealb_cea", alb_ice(:,:,1) * fr_i(:,:) )  ! ice albedo
-      ENDIF
+#if defined key_coupled
+      tn_ice(:,:,1) = sist(:,:)          ! sea-ice surface temperature       
+      ht_i(:,:,1) = hicif(:,:)
+      ht_s(:,:,1) = hsnif(:,:)
+      a_i(:,:,1) = fr_i(:,:)
+      !                                  ! Computation of snow/ice and ocean albedo
+      CALL albedo_ice( tn_ice, ht_i, ht_s, zalbp, zalb )
+      alb_ice(:,:,1) =  0.5 * ( zalbp(:,:,1) + zalb (:,:,1) )   ! Ice albedo (mean clear and overcast skys)
+      CALL iom_put( "icealb_cea", alb_ice(:,:,1) * fr_i(:,:) )  ! ice albedo
+#endif
 
       IF(ln_ctl) THEN            ! control print
@@ -250 +249 @@
       ENDIF 
       !
-      IF( wrk_not_released(2, 1)     .OR.    &
-          wrk_not_released(3, 4,5) )   CALL ctl_stop('lim_sbc_flx_2 : failed to release workspace arrays')
+      CALL wrk_dealloc( jpi, jpj, zqnsoce )
+      CALL wrk_dealloc( jpi, jpj, 1, zalb, zalbp )
       !
    END SUBROUTINE lim_sbc_flx_2
@@ -282 +281 @@
       !!              - taum       : modulus of the surface ocean stress (T-point) updated with ice-ocean fluxes
       !!---------------------------------------------------------------------
-      USE wrk_nemo, ONLY: wrk_in_use, wrk_not_released
-      USE wrk_nemo, ONLY: ztio_u => wrk_2d_1, ztio_v => wrk_2d_2     ! ocean stress below sea-ice
       INTEGER ,                     INTENT(in) ::   kt               ! ocean time-step index
       REAL(wp), DIMENSION(jpi,jpj), INTENT(in) ::   pu_oce, pv_oce   ! surface ocean currents
@@ -291 +288 @@
       REAL(wp) ::   zfrldv, zat_v, zv_i, zvtau_ice, zv_t, zmodi   !   -      -
       REAL(wp) ::   zsang, zumt                                   !    -         -
+      REAL(wp), POINTER, DIMENSION(:,:) ::   ztio_u, ztio_v   ! ocean stress below sea-ice
       !!---------------------------------------------------------------------
       !
-      IF( wrk_in_use(2, 1,2) ) THEN
-         CALL ctl_stop('lim_sbc_tau_2 : requested workspace arrays unavailable.')   ;   RETURN
-      ENDIF
+      CALL wrk_alloc( jpi, jpj, ztio_u, ztio_v )
       !
       SELECT CASE( cp_ice_msh )     
@@ -409 +405 @@
          &                       tab2d_2=vtau, clinfo2=' vtau    : '        , mask2=vmask )
       !  
-      IF( wrk_not_released(2, 1,2) )   CALL ctl_stop('lim_sbc_tau_2 : failed to release workspace arrays')
+      CALL wrk_dealloc( jpi, jpj, ztio_u, ztio_v )
       !
    END SUBROUTINE lim_sbc_tau_2