Changeset 2715 for trunk/NEMOGCM/NEMO/LIM_SRC_2/limwri_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/limwri_2.F90 (modified) (8 diffs)

trunk/NEMOGCM/NEMO/LIM_SRC_2/limwri_2.F90

@@ -26 +26 @@
    USE ice_2
 
+   USE dianam          ! build name of file (routine)
    USE lbclnk
-   USE dianam          ! build name of file (routine)
    USE in_out_manager
+   USE lib_mpp         ! MPP library
    USE iom
    USE ioipsl
@@ -39 +40 @@
 #endif
    PUBLIC   lim_wri_state_2   ! called by dia_wri_state 
+   PUBLIC   lim_wri_alloc_2   ! called by nemogcm.F90
 
    INTEGER, PARAMETER                       ::   jpnoumax = 40   ! maximum number of variable for ice output
@@ -50 +52 @@
 
    INTEGER ::   nice, nhorid, ndim, niter, ndepid       ! ????
-   INTEGER , DIMENSION( jpij ) ::   ndex51              ! ????
-
-   REAL(wp)  ::            &  ! constant values
-      epsi16 = 1.e-16   ,  &
-      zzero  = 0.e0     ,  &
-      zone   = 1.e0
+   INTEGER, ALLOCATABLE, SAVE, DIMENSION(:) :: ndex51   ! ????
+
+   REAL(wp) ::   epsi16 = 1.e-16_wp   ! constant values
+   REAL(wp) ::   zzero  = 0._wp       !     -      -
+   REAL(wp) ::   zone   = 1._wp       !     -      -
+
+   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   zcmo      ! Workspace array for netcdf writer. 
+
 
    !! * Substitutions
@@ -64 +68 @@
    !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
    !!----------------------------------------------------------------------
-
 CONTAINS
+
+   INTEGER FUNCTION lim_wri_alloc_2()
+      !!-------------------------------------------------------------------
+      !!                  ***   ROUTINE lim_wri_alloc_2  ***
+      !!-------------------------------------------------------------------
+      ALLOCATE( ndex51(jpij), zcmo(jpi,jpj,jpnoumax), STAT=lim_wri_alloc_2)
+      !
+      IF( lk_mpp               )   CALL mpp_sum ( lim_wri_alloc_2 )
+      IF( lim_wri_alloc_2 /= 0 )   CALL ctl_warn('lim_wri_alloc_2: failed to allocate array ndex51')
+      !
+   END FUNCTION lim_wri_alloc_2
+
 
 #if ! defined key_iomput
@@ -85 +100 @@
       !!      of a day
       !!-------------------------------------------------------------------
+      USE wrk_nemo, ONLY: wrk_in_use, wrk_not_released
+      USE wrk_nemo, ONLY: zfield => wrk_2d_1
+      !!
       INTEGER, INTENT(in) ::   kt     ! number of iteration
       !!
@@ -92 +110 @@
          &          zindh, zinda, zindb, ztmu
       REAL(wp), DIMENSION(1)                ::   zdept
-      REAL(wp), DIMENSION(jpi,jpj)          ::   zfield
-      REAL(wp), DIMENSION(jpi,jpj,jpnoumax) ::   zcmo
-      !!-------------------------------------------------------------------
+      !!-------------------------------------------------------------------
+
+      IF( wrk_in_use(2, 1) ) THEN
+         CALL ctl_stop('lim_wri_2 : requested workspace array unavailable')   ;   RETURN
+      ENDIF
                                                  !--------------------!
       IF( kt == nit000 ) THEN                    !   Initialisation   !
          !                                       !--------------------!
+
          CALL lim_wri_init_2 
                            
@@ -186 +207 @@
       IF( ( nn_fsbc * niter ) >= nitend )   CALL histclo( nice ) 
 
+      IF( wrk_not_released(2, 1) )   CALL ctl_stop('lim_wri_2 : failed to release workspace array.')
+      !
    END SUBROUTINE lim_wri_2
      
@@ -222 +245 @@
          field_19
       !!-------------------------------------------------------------------
+      !
+      IF( lim_wri_alloc_2() /= 0 ) THEN      ! allocate lim_wri arrrays
+         CALL ctl_stop( 'STOP', 'lim_wri_init_2 : unable to allocate standard arrays' )   ;   RETURN
+      ENDIF
 
       REWIND ( numnam_ice )                ! Read Namelist namicewri