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

trunk/NEMOGCM/NEMO/LIM_SRC_2/limthd_2.F90

@@ -75 +75 @@
       !! References :   Goosse et al. 1996, Bul. Soc. Roy. Sc. Liege, 65, 87-90
       !!---------------------------------------------------------------------
+      USE wrk_nemo, ONLY: wrk_in_use, wrk_not_released
+      USE wrk_nemo, ONLY: ztmp    => wrk_2d_1, & ! 2D workspace
+                          zqlbsbq => wrk_2d_2, & ! link with lead energy budget qldif
+                          zlicegr => wrk_2d_3    ! link with lateral ice growth 
+      USE wrk_nemo, ONLY: zmsk => wrk_3d_4       ! 3D workspace
+      USE wrk_nemo, ONLY: zdvosif => wrk_2d_4, & !: Variation of volume at surface
+                          zdvobif => wrk_2d_5, & !: Variation of ice volume at the bottom ice     (outputs only)
+                          zdvolif => wrk_2d_6, & !: Total variation of ice volume                 (outputs only)
+                          zdvonif => wrk_2d_7, & !: Surface accretion Snow to Ice transformation  (outputs only)
+                          zdvomif => wrk_2d_8, & !: Bottom variation of ice volume due to melting (outputs only)
+                          zu_imasstr =>wrk_2d_9, & !: Sea-ice transport along i-axis at U-point     (outputs only) 
+                          zv_imasstr =>wrk_2d_10   !: Sea-ice transport along j-axis at V-point     (outputs only) 
+      !!
       INTEGER, INTENT(in) ::   kt     ! number of iteration
       !!
@@ -91 +104 @@
       REAL(wp) ::   zfontn               ! heat flux from snow thickness
       REAL(wp) ::   zfntlat, zpareff     ! test. the val. of lead heat budget
-      REAL(wp), DIMENSION(jpi,jpj)     ::   ztmp      ! 2D workspace
-      REAL(wp), DIMENSION(jpi,jpj)     ::   zqlbsbq   ! link with lead energy budget qldif
+
       REAL(wp) ::   zuice_m, zvice_m     ! Sea-ice velocities at U & V-points
       REAL(wp) ::   zhice_u, zhice_v     ! Sea-ice volume at U & V-points
@@ -98 +110 @@
       REAL(wp) ::   zrhoij, zrhoijm1     ! temporary scalars
       REAL(wp) ::   zztmp                ! temporary scalars within a loop
-      REAL(wp), DIMENSION(jpi,jpj)     ::   zlicegr   ! link with lateral ice growth 
-      REAL(wp), DIMENSION(jpi,jpj,jpk) ::   zmsk      ! 3D workspace
 !!$      REAL(wp), DIMENSION(jpi,jpj) ::   firic         !: IR flux over the ice            (outputs only)
 !!$      REAL(wp), DIMENSION(jpi,jpj) ::   fcsic         !: Sensible heat flux over the ice (outputs only)
 !!$      REAL(wp), DIMENSION(jpi,jpj) ::   fleic         !: Latent heat flux over the ice   (outputs only)
 !!$      REAL(wp), DIMENSION(jpi,jpj) ::   qlatic        !: latent flux                     (outputs only)
-      REAL(wp), DIMENSION(jpi,jpj) ::   zdvosif       !: Variation of volume at surface                (outputs only)
-      REAL(wp), DIMENSION(jpi,jpj) ::   zdvobif       !: Variation of ice volume at the bottom ice     (outputs only)
-      REAL(wp), DIMENSION(jpi,jpj) ::   zdvolif       !: Total variation of ice volume                 (outputs only)
-      REAL(wp), DIMENSION(jpi,jpj) ::   zdvonif       !: Surface accretion Snow to Ice transformation  (outputs only)
-      REAL(wp), DIMENSION(jpi,jpj) ::   zdvomif       !: Bottom variation of ice volume due to melting (outputs only)
-      REAL(wp), DIMENSION(jpi,jpj) ::   zu_imasstr    !: Sea-ice transport along i-axis at U-point     (outputs only) 
-      REAL(wp), DIMENSION(jpi,jpj) ::   zv_imasstr    !: Sea-ice transport along j-axis at V-point     (outputs only) 
       !!-------------------------------------------------------------------
+
+      IF( wrk_in_use(2, 1,2,3,4,5,6,7,8,9,10)   .OR.  &
+          wrk_in_use(3, 4)                    ) THEN
+         CALL ctl_stop('lim_thd_2 : requested workspace arrays unavailable')   ;   RETURN
+      ENDIF
 
       IF( kt == nit000 )   CALL lim_thd_init_2  ! Initialization (first time-step only)
@@ -512 +520 @@
       ENDIF
        !
+      IF( wrk_not_released(2, 1,2,3,4,5,6,7,8,9,10)   .OR.  &
+          wrk_not_released(3, 4)                    ) THEN
+         CALL ctl_stop('lim_thd_2 : failed to release workspace arrays')
+      ENDIF
+      !
     END SUBROUTINE lim_thd_2