Changeset 2715 for trunk/NEMOGCM/NEMO/LIM_SRC_2/limrhg_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/limrhg_2.F90 (modified) (4 diffs)

trunk/NEMOGCM/NEMO/LIM_SRC_2/limrhg_2.F90

@@ -33 +33 @@
    PRIVATE
 
-   PUBLIC   lim_rhg_2 ! routine called by lim_dyn
+   PUBLIC   lim_rhg_2         ! routine called by lim_dyn
+   PUBLIC   lim_rhg_alloc_2   ! routine called by lim_dyn_alloc_2
 
    REAL(wp) ::   rzero   = 0._wp   ! constant value: zero
    REAL(wp) ::   rone    = 1._wp   !            and  one
+
+   ! 2D workspaces for lim_rhg_2. Can't use wrk_nemo module for them because
+   ! extent in 2nd dimension is > jpj.
+   REAL(wp), SAVE, ALLOCATABLE, DIMENSION(:,:) ::   zu0, zv0
+   REAL(wp), SAVE, ALLOCATABLE, DIMENSION(:,:) ::   zu_n, zv_n
+   REAL(wp), SAVE, ALLOCATABLE, DIMENSION(:,:) ::   zu_a, zv_a
+   REAL(wp), SAVE, ALLOCATABLE, DIMENSION(:,:) ::   zviszeta, zviseta
+   REAL(wp), SAVE, ALLOCATABLE, DIMENSION(:,:) ::   zzfrld, zztms
+   REAL(wp), SAVE, ALLOCATABLE, DIMENSION(:,:) ::   zi1, zi2, zmasst, zpresh
 
    !! * Substitutions
@@ -46 +56 @@
    !!----------------------------------------------------------------------
 CONTAINS
+
+   INTEGER FUNCTION lim_rhg_alloc_2()
+      !!-------------------------------------------------------------------
+      !!               ***  FUNCTION lim_rhg_alloc_2  ***
+      !!-------------------------------------------------------------------
+      ALLOCATE( zu0(jpi,0:jpj+1),      zv0(jpi,0:jpj+1),     &
+         &      zu_n(jpi,0:jpj+1),     zv_n(jpi,0:jpj+1),    &
+         &      zu_a(jpi,0:jpj+1),     zv_a(jpi,0:jpj+1),    &
+         &      zviszeta(jpi,0:jpj+1), zviseta(jpi,0:jpj+1), &
+         &      zzfrld(jpi,0:jpj+1),   zztms(jpi,0:jpj+1),   &
+         &      zi1(jpi,0:jpj+1),      zi2(jpi,0:jpj+1),     &
+         &      zmasst(jpi,0:jpj+1),   zpresh(jpi,0:jpj+1),  &
+         &      Stat=lim_rhg_alloc_2)
+         !
+      IF( lim_rhg_alloc_2 /= 0 )   CALL ctl_warn('lim_rhg_alloc_2 : failed to allocate arrays')
+      !
+   END FUNCTION lim_rhg_alloc_2
+
 
    SUBROUTINE lim_rhg_2( k_j1, k_jpj )
@@ -59 +87 @@
       !!              at I-point
       !!-------------------------------------------------------------------
+      USE wrk_nemo, ONLY: wrk_in_use, wrk_not_released
+      USE wrk_nemo, ONLY: zfrld => wrk_2d_1, zmass => wrk_2d_2, zcorl => wrk_2d_3
+      USE wrk_nemo, ONLY: za1ct => wrk_2d_4, za2ct => wrk_2d_5, zresr => wrk_2d_6
+      USE wrk_nemo, ONLY: zc1u  => wrk_2d_7, zc1v  => wrk_2d_8, zc2u => wrk_2d_9
+      USE wrk_nemo, ONLY: zc2v  => wrk_2d_10, zsang => wrk_2d_11
+      !!
       INTEGER, INTENT(in) ::   k_j1    ! southern j-index for ice computation
       INTEGER, INTENT(in) ::   k_jpj   ! northern j-index for ice computation
@@ -79 +113 @@
       REAL(wp) ::   zs21_11, zs21_12, zs21_21, zs21_22
       REAL(wp) ::   zs22_11, zs22_12, zs22_21, zs22_22
-      REAL(wp), DIMENSION(jpi,  jpj  ) ::   zfrld, zmass, zcorl
-      REAL(wp), DIMENSION(jpi,  jpj  ) ::   za1ct, za2ct, zresr
-      REAL(wp), DIMENSION(jpi,  jpj  ) ::   zc1u, zc1v, zc2u, zc2v
-      REAL(wp), DIMENSION(jpi,  jpj  ) ::   zsang
-      REAL(wp), DIMENSION(jpi,0:jpj+1) ::   zu0, zv0
-      REAL(wp), DIMENSION(jpi,0:jpj+1) ::   zu_n, zv_n
-      REAL(wp), DIMENSION(jpi,0:jpj+1) ::   zu_a, zv_a
-      REAL(wp), DIMENSION(jpi,0:jpj+1) ::   zviszeta, zviseta
-      REAL(wp), DIMENSION(jpi,0:jpj+1) ::   zzfrld, zztms
-      REAL(wp), DIMENSION(jpi,0:jpj+1) ::   zi1, zi2, zmasst, zpresh
       !!-------------------------------------------------------------------