Changeset 2715 for trunk/NEMOGCM/NEMO/OPA_SRC/DTA/dtatem.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/OPA_SRC/DTA/dtatem.F90 (modified) (5 diffs)

trunk/NEMOGCM/NEMO/OPA_SRC/DTA/dtatem.F90

@@ -18 +18 @@
    USE oce             ! ocean dynamics and tracers
    USE dom_oce         ! ocean space and time domain
+   USE phycst          ! physical constants
    USE fldread         ! read input fields
    USE in_out_manager  ! I/O manager
-   USE phycst          ! physical constants
+   USE lib_mpp         ! MPP library
 
    IMPLICIT NONE
@@ -27 +28 @@
    PUBLIC   dta_tem    ! called by step.F90 and inidta.F90
 
-   LOGICAL , PUBLIC, PARAMETER ::   lk_dtatem = .TRUE.   !: temperature data flag
-   REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpk) ::  t_dta    !: temperature data at given time-step
+   LOGICAL , PUBLIC, PARAMETER                     ::   lk_dtatem = .TRUE. !: temperature data flag
+   REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:,:) ::   t_dta              !: temperature data at given time-step
 
    TYPE(FLD), ALLOCATABLE, DIMENSION(:) ::   sf_tem      ! structure of input SST (file informations, fields read)
@@ -35 +36 @@
 #  include "domzgr_substitute.h90"
    !!----------------------------------------------------------------------
-   !! NEMO/OPA 3.3 , NEMO Consortium (2010)
+   !! NEMO/OPA 4.0 , NEMO Consortium (2011)
    !! $Id$ 
    !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt)
@@ -58 +59 @@
       !! ** Action  :   define t_dta array at time-step kt
       !!----------------------------------------------------------------------
-      INTEGER, INTENT( in ) ::   kt     ! ocean time-step
+      INTEGER, INTENT( in ) ::   kt   ! ocean time-step
       !
-      INTEGER ::   ji, jj, jk, jl, jkk            ! dummy loop indicies
-      INTEGER ::   ik, ierror                     ! temporary integers
+      INTEGER ::   ji, jj, jk, jl, jkk       ! dummy loop indicies
+      INTEGER ::   ik, ierr, ierr0, ierr1, ierr2   ! local integers
 #if defined key_tradmp
-      INTEGER ::   il0, il1, ii0, ii1, ij0, ij1   ! temporary integers
+      INTEGER ::   il0, il1, ii0, ii1, ij0, ij1   ! local integers
 #endif
       REAL(wp)::   zl
@@ -95 +96 @@
             WRITE(numout,*) '~~~~~~~ '
          ENDIF
-         ALLOCATE( sf_tem(1), STAT=ierror )
-         IF( ierror > 0 ) THEN
-             CALL ctl_stop( 'dta_tem: unable to allocate sf_tem structure' )   ;   RETURN
-         ENDIF
-                                ALLOCATE( sf_tem(1)%fnow(jpi,jpj,jpk)   )
-         IF( sn_tem%ln_tint )   ALLOCATE( sf_tem(1)%fdta(jpi,jpj,jpk,2) )
+
+                                   ! Allocate temperature data array 
+                                ALLOCATE( t_dta(jpi,jpj,jpk)           , STAT=ierr  )
+         IF( ierr > 0              )   CALL ctl_stop( 'STOP', 'dta_tem: unable to allocate t_dta array' )
+                                   ! Allocate sf_tem structure
+                                ierr2 = 0
+                                ALLOCATE( sf_tem(1)                    , STAT=ierr0 )
+                                ALLOCATE( sf_tem(1)%fnow(jpi,jpj,jpk)  , STAT=ierr1 )
+         IF( sn_tem%ln_tint )   ALLOCATE( sf_tem(1)%fdta(jpi,jpj,jpk,2), STAT=ierr2 )
+         IF( ierr0+ierr1+ierr2 > 0 )   CALL ctl_stop( 'STOP', 'dta_tem: unable to allocate sf_tem structure' )
          !                         ! fill sf_tem with sn_tem and control print
          CALL fld_fill( sf_tem, (/ sn_tem /), cn_dir, 'dta_tem', 'Temperature data', 'namdta_tem' )