Changeset 3152 for branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90

Timestamp:

2011-11-18T10:19:26+01:00 (12 years ago)

Author:

smasson

Message:

dev_NEMO_MERGE_2011: new dynamical allocation in IOM and SBC

File:

• branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90 (modified) (4 diffs)

branches/2011/dev_NEMO_MERGE_2011/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90

@@ -18 +18 @@
    USE lib_mpp         ! distributed memory computing library
    USE lbclnk          ! ocean lateral boundary conditions (or mpp link)
+   USE wrk_nemo_2      ! work arrays
    USE daymod          ! calendar
    USE fldread         ! read input fields
@@ -178 +179 @@
      CALL lbc_lnk ( fr_iv , 'V', 1. )
 
-
    END SUBROUTINE cice_sbc_init
 
+   
    SUBROUTINE cice_sbc_in (kt, nsbc)
       !!---------------------------------------------------------------------
       !!                    ***  ROUTINE cice_sbc_in  ***
       !! ** Purpose: Set coupling fields and pass to CICE
-      !!
-      USE wrk_nemo, ONLY:   wrk_in_use, wrk_not_released
-      USE wrk_nemo, ONLY:   wrk_2d_1, wrk_3d_1                    ! Workspace
-      !!
-      INTEGER, INTENT( in  ) ::   kt   ! ocean time step
-      INTEGER, INTENT( in  ) ::   nsbc ! surface forcing type
-      !!---------------------------------------------------------------------
-     INTEGER  ::   ji, jj, jpl                   ! dummy loop indices
-
-     REAL(wp), DIMENSION(:,:), POINTER :: ztmp
-     REAL(wp), DIMENSION(:,:,:), POINTER :: ztmpn
-
-     IF( wrk_in_use(2, 1) .OR. wrk_in_use(3, 1)  ) THEN
-        CALL ctl_stop('cice_sbc_in: requested workspace arrays are unavailable')   ;   RETURN
-     ENDIF
-     ztmp => wrk_2d_1(:,:)
-     ztmpn => wrk_3d_1(:,:,1:ncat)
+      !!---------------------------------------------------------------------
+      INTEGER, INTENT(in   ) ::   kt   ! ocean time step
+      INTEGER, INTENT(in   ) ::   nsbc ! surface forcing type
+
+      INTEGER  ::   ji, jj, jpl                   ! dummy loop indices      
+      REAL(wp), DIMENSION(:,:), POINTER :: ztmp
+      REAL(wp), DIMENSION(:,:,:), POINTER :: ztmpn
+      !!---------------------------------------------------------------------
+
+      CALL wrk_alloc( jpi,jpj, ztmp )
+      CALL wrk_alloc( jpi,jpj,ncat, ztmpn )
 
      IF( kt == nit000 )  THEN
@@ -368 +363 @@
      CALL nemo2cice(ztmp,vocn,'F', -1. )
 
-     IF( wrk_not_released(2, 1) .OR. wrk_not_released(3, 1)  )   CALL ctl_stop('cice_sbc_in: failed to release workspace arrays')
+      CALL wrk_dealloc( jpi,jpj, ztmp )
+      CALL wrk_dealloc( jpi,jpj,ncat, ztmpn )
      !
    END SUBROUTINE cice_sbc_in
 
+
    SUBROUTINE cice_sbc_out (kt,nsbc)
       !!---------------------------------------------------------------------
       !!                    ***  ROUTINE cice_sbc_out  ***
       !! ** Purpose: Get fields from CICE and set surface fields for NEMO
-      !!
-      USE wrk_nemo, ONLY:   wrk_in_use, wrk_not_released
-      USE wrk_nemo, ONLY:   wrk_2d_1                         ! 2D workspace
-      !!
+      !!---------------------------------------------------------------------
       INTEGER, INTENT( in  ) ::   kt   ! ocean time step
       INTEGER, INTENT( in  ) ::   nsbc ! surface forcing type
-      !!---------------------------------------------------------------------
-
-     INTEGER  ::   ji, jj, jpl                 ! dummy loop indices
-
-     REAL(wp), DIMENSION(:,:), POINTER :: ztmp
-
-     IF( kt == nit000 )  THEN
+      
+      INTEGER  ::   ji, jj, jpl                 ! dummy loop indices
+      REAL(wp), DIMENSION(:,:), POINTER :: ztmp
+      !!---------------------------------------------------------------------
+
+      CALL wrk_alloc( jpi,jpj, ztmp )
+      
+      IF( kt == nit000 )  THEN
          IF(lwp) WRITE(numout,*)'cice_sbc_out'
-     ENDIF
-
-     IF( wrk_in_use(2, 1) ) THEN
-        CALL ctl_stop('cice_sbc_out: requested workspace arrays are unavailable')   ;   RETURN
-     ENDIF
-     ztmp => wrk_2d_1(:,:)
-
+      ENDIF
+      
 ! x comp of ocean-ice stress 
      CALL cice2nemo(strocnx,ztmp,'F', -1. )
@@ -532 +522 @@
 ! Release work space
 
-     IF( wrk_not_released(2, 1) )   CALL ctl_stop('cice_sbc_out: failed to release workspace arrays')
+      CALL wrk_dealloc( jpi,jpj, ztmp )
      !
    END SUBROUTINE cice_sbc_out
+
 
 #if defined key_oasis3 || defined key_oasis4