Changeset 2715 for trunk/NEMOGCM/NEMO/LIM_SRC_3/limdia.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_3/limdia.F90 (modified) (4 diffs)

trunk/NEMOGCM/NEMO/LIM_SRC_3/limdia.F90

@@ -52 +52 @@
 
    REAL(wp), DIMENSION(jpinfmx) ::   vinfom     ! temporary working space
-   REAL(wp), DIMENSION(jpi,jpj) ::   aire       ! masked grid cell area
+   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) ::   aire       ! masked grid cell area
 
    !! * Substitutions
@@ -67 +67 @@
       !!                  ***  ROUTINE lim_dia  ***
       !!   
-      !! ** Purpose : Computation and outputs on file ice.evolu 
-      !!      the temporal evolution of some key variables
+      !! ** Purpose :   Computation and outputs on file ice.evolu 
+      !!              the temporal evolution of some key variables
       !!-------------------------------------------------------------------
       INTEGER  ::   jv, ji, jj, jl   ! dummy loop indices
@@ -410 +410 @@
       !!-------------------------------------------------------------------
       INTEGER  ::   jv    ! dummy loop indice
-      INTEGER  ::   ntot , ndeb , irecl   ! local integers
+      INTEGER  ::   ierr, ntot , ndeb , irecl   ! local integers
       REAL(wp) ::   zxx0, zxx1    ! local scalars
       CHARACTER(len=jpchinf) ::   titinf
@@ -431 +431 @@
       ENDIF
 
-      aire(:,:) = area(:,:) * tms(:,:) * tmask_i(:,:)      ! masked grid cell area (interior domain only)
+      ALLOCATE( aire(jpi,jpj) , STAT=ierr )      ! masked grid cell area (interior domain only)
+      IF( lk_mpp    )   CALL mpp_sum( ierr )
+      IF( ierr /= 0 )   CALL ctl_stop( 'STOP', 'lim_dia_init_2 : unable to allocate arrays' )
+      aire(:,:) = area(:,:) * tms(:,:) * tmask_i(:,:)
 
       ! Titles of ice key variables :
       titvar(1) = 'NoIt'  ! iteration number
       titvar(2) = 'T yr'  ! time step in years
-      nbvt = 2            ! number of time variables
+      nbvt      = 2       ! number of time variables
 
       titvar(3) = 'AI_N'  ! sea ice area in the northern Hemisp.(10^12 km2)