Changeset 2715 for trunk/NEMOGCM/NEMO/LIM_SRC_2/limdyn_2.F90

Timestamp:

2011-03-30T17:58:35+02:00 (13 years ago)

Author:

rblod

Message:

First attempt to put dynamic allocation on the trunk

File:

• trunk/NEMOGCM/NEMO/LIM_SRC_2/limdyn_2.F90 (modified) (4 diffs)

trunk/NEMOGCM/NEMO/LIM_SRC_2/limdyn_2.F90

@@ -58 +58 @@
       !!              - treatment of the case if no ice dynamic
       !!---------------------------------------------------------------------
+      USE wrk_nemo, ONLY:   wrk_in_use, wrk_not_released
+      USE wrk_nemo, ONLY:   wrk_1d_1, wrk_1d_2
+      USE wrk_nemo, ONLY:   zu_io => wrk_2d_1, zv_io => wrk_2d_2  ! ice-ocean velocity
+      !
       INTEGER, INTENT(in) ::   kt     ! number of iteration
       !!
@@ -63 +67 @@
       INTEGER  ::   i_j1, i_jpj        ! Starting/ending j-indices for rheology
       REAL(wp) ::   zcoef              ! temporary scalar
-      REAL(wp), DIMENSION(jpj)     ::   zind           ! i-averaged indicator of sea-ice
-      REAL(wp), DIMENSION(jpj)     ::   zmsk           ! i-averaged of tmask
-      REAL(wp), DIMENSION(jpi,jpj) ::   zu_io, zv_io   ! ice-ocean velocity
+      REAL(wp), POINTER, DIMENSION(:) ::   zind     ! i-averaged indicator of sea-ice
+      REAL(wp), POINTER, DIMENSION(:) ::   zmsk     ! i-averaged of tmask
       !!---------------------------------------------------------------------
+
+      IF(  wrk_in_use(1, 1,2)  .OR.  wrk_in_use(2, 1,2)  ) THEN
+         CALL ctl_stop( 'lim_dyn_2 : requested workspace arrays unavailable' )   ;   RETURN
+      ENDIF
+      zind => wrk_1d_1(1:jpj)      ! Set-up pointers to sub-arrays of workspaces
+      zmsk => wrk_1d_2(1:jpj)
 
       IF( kt == nit000 )   CALL lim_dyn_init_2   ! Initialization (first time-step only)
@@ -93 +102 @@
             !
             DO jj = 1, jpj
-               zind(jj) = SUM( frld (:,jj  ) )   ! = FLOAT(jpj) if ocean everywhere on a j-line
-               zmsk(jj) = SUM( tmask(:,jj,1) )   ! = 0          if land  everywhere on a j-line
+               zind(jj) = SUM( frld (:,jj  ) )   ! = REAL(jpj) if ocean everywhere on a j-line
+               zmsk(jj) = SUM( tmask(:,jj,1) )   ! = 0         if land  everywhere on a j-line
             END DO
             !
@@ -200 +209 @@
       !
       IF(ln_ctl)   CALL prt_ctl(tab2d_1=ust2s , clinfo1=' lim_dyn  : ust2s :')
+      !
+      IF( wrk_not_released(1, 1,2)     .OR.   &
+          wrk_not_released(2, 1,2) )   CALL ctl_stop('lim_dyn_2 : failed to release workspace arrays')
       !
    END SUBROUTINE lim_dyn_2