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

trunk/NEMOGCM/NEMO/LIM_SRC_2/limdia_2.F90

@@ -12 +12 @@
    !!   'key_lim2' :                                  LIM 2.0 sea-ice model
    !!----------------------------------------------------------------------
-   !!----------------------------------------------------------------------
    !!   lim_dia_2      : computation of the time evolution of keys var.
    !!   lim_dia_init_2 : initialization and namelist read
@@ -24 +23 @@
    USE limistate_2     !
    USE in_out_manager  ! I/O manager
+   USE lib_mpp         ! MPP library
 
    IMPLICIT NONE
@@ -29 +29 @@
 
    PUBLIC               lim_dia_2          ! called by sbc_ice_lim_2
+
    INTEGER, PUBLIC ::   ntmoy   = 1 ,   &  !: instantaneous values of ice evolution or averaging ntmoy
       &                 ninfo   = 1        !: frequency of ouputs on file ice_evolu in case of averaging
@@ -52 +53 @@
    REAL(wp)                     ::   epsi06 = 1.e-06      ! ???
    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
@@ -61 +62 @@
    !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
    !!----------------------------------------------------------------------
-
 CONTAINS
 
@@ -173 +173 @@
        !!-------------------------------------------------------------------
        CHARACTER(len=jpchinf) ::   titinf
-       INTEGER  ::   jv            ! dummy loop indice
-       INTEGER  ::   ntot , ndeb 
-       INTEGER  ::   nv            ! indice of variable 
-       REAL(wp) ::   zxx0, zxx1    ! temporary scalars
+       INTEGER  ::   jv   ! dummy loop indice
+       INTEGER  ::   ntot , ndeb, nv, ierr   ! local integer
+       REAL(wp) ::   zxx0, zxx1              ! local scalars
 
        NAMELIST/namicedia/fmtinf, nfrinf, ninfo, ntmoy
        !!-------------------------------------------------------------------
 
-       ! Read Namelist namicedia
-       REWIND ( numnam_ice )
-       READ   ( numnam_ice  , namicedia )
+       REWIND( numnam_ice )                     ! Read Namelist namicedia
+       READ  ( numnam_ice  , namicedia )
        
-       IF(lwp) THEN
+       IF(lwp) THEN                             ! control print
           WRITE(numout,*)
           WRITE(numout,*) 'lim_dia_init_2 : ice parameters for ice diagnostics '
@@ -195 +193 @@
        ENDIF
 
-       ! masked grid cell area
+       ALLOCATE( aire(jpi,jpj) , STAT=ierr )    ! masked grid cell area
+       IF( lk_mpp    )   CALL mpp_sum( ierr )
+       IF( ierr /= 0 )   CALL ctl_stop( 'STOP', 'lim_dia_init_2 : unable to allocate standard arrays' )
        aire(:,:) = area(:,:) * tms(:,:)
 
-       ! Titles of ice key variables :
-       nv = 1
+       nv = 1                                   ! Titles of ice key variables
        titvar(nv) = 'NoIt'  ! iteration number
        nv = nv + 1
        titvar(nv) = 'T yr'  ! time step in years
-       
        nbvt = nv - 1
-
        nv = nv + 1   ;   titvar(nv) = 'AEFN' ! sea ice area in the northern Hemisp.(10^12 km2)
        nv = nv + 1   ;   titvar(nv) = 'AEFS' ! sea ice area in the southern Hemisp.(10^12 km2)