Changeset 2715 for trunk/NEMOGCM/NEMO/OPA_SRC/BDY/bdyini.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/BDY/bdyini.F90 (modified) (4 diffs)

trunk/NEMOGCM/NEMO/OPA_SRC/BDY/bdyini.F90

@@ -21 +21 @@
    USE obc_par         ! ocean open boundary conditions
    USE bdy_oce         ! unstructured open boundary conditions
+   USE bdydta, ONLY: bdy_dta_alloc ! open boundary data
    USE bdytides        ! tides at open boundaries initialization (tide_init routine)
    USE in_out_manager  ! I/O units
@@ -33 +34 @@
 
    !!----------------------------------------------------------------------
-   !! 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)
+   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt)
    !!----------------------------------------------------------------------
 CONTAINS
@@ -44 +45 @@
       !!         
       !! ** Purpose :   Initialization of the dynamics and tracer fields with 
-      !!      unstructured open boundaries.
+      !!              unstructured open boundaries.
       !!
-      !! ** Method  :  Read initialization arrays (mask, indices) to identify 
-      !!               an unstructured open boundary
+      !! ** Method  :   Read initialization arrays (mask, indices) to identify 
+      !!              an unstructured open boundary
       !!
       !! ** Input   :  bdy_init.nc, input file for unstructured open boundaries
       !!----------------------------------------------------------------------      
-      INTEGER ::   ii, ij, ik, igrd, ib, ir   ! dummy loop indices
-      INTEGER ::   icount, icountr
-      INTEGER ::   ib_len, ibr_max
-      INTEGER ::   iw, ie, is, in 
-      INTEGER ::   inum                 ! local logical unit
-      INTEGER ::   id_dummy             ! local integers
-      INTEGER ::   igrd_start, igrd_end ! start and end of loops on igrd
+      INTEGER  ::   ii, ij, ik, igrd, ib, ir   ! dummy loop indices
+      INTEGER  ::   icount, icountr, ib_len, ibr_max   ! local integers
+      INTEGER  ::   iw, ie, is, in, inum, id_dummy     !   -       -
+      INTEGER  ::   igrd_start, igrd_end               !   -       -
+      REAL(wp) ::   zefl, zwfl, znfl, zsfl              ! local scalars
       INTEGER, DIMENSION (2)             ::   kdimsz
       INTEGER, DIMENSION(jpbdta, jpbgrd) ::   nbidta, nbjdta   ! Index arrays: i and j indices of bdy dta
       INTEGER, DIMENSION(jpbdta, jpbgrd) ::   nbrdta           ! Discrete distance from rim points
-      REAL(wp) :: zefl, zwfl, znfl, zsfl                       ! temporary scalars
-      REAL(wp) , DIMENSION(jpidta,jpjdta) ::   zmask           ! global domain mask
-      REAL(wp) , DIMENSION(jpbdta,1)      ::   zdta            ! temporary array 
-      CHARACTER(LEN=80),DIMENSION(6)      ::   clfile
+      REAL(wp), DIMENSION(jpidta,jpjdta) ::   zmask            ! global domain mask
+      REAL(wp), DIMENSION(jpbdta,1)      ::   zdta             ! temporary array 
+      CHARACTER(LEN=80),DIMENSION(6)     ::   clfile
       !!
-      NAMELIST/nambdy/cn_mask, cn_dta_frs_T, cn_dta_frs_U, cn_dta_frs_V,          &
-         &            cn_dta_fla_T, cn_dta_fla_U, cn_dta_fla_V,              &
-         &            ln_tides, ln_clim, ln_vol, ln_mask,                &
-         &            ln_dyn_fla, ln_dyn_frs, ln_tra_frs,ln_ice_frs,     &
+      NAMELIST/nambdy/cn_mask, cn_dta_frs_T, cn_dta_frs_U, cn_dta_frs_V,   &
+         &            cn_dta_fla_T, cn_dta_fla_U, cn_dta_fla_V,            &
+         &            ln_tides, ln_clim, ln_vol, ln_mask,                  &
+         &            ln_dyn_fla, ln_dyn_frs, ln_tra_frs,ln_ice_frs,       &
          &            nn_dtactl, nn_rimwidth, nn_volctl
       !!----------------------------------------------------------------------
@@ -77 +75 @@
       IF(lwp) WRITE(numout,*) '~~~~~~~~'
       !
+      !                                      ! allocate bdy_oce arrays
+      IF( bdy_oce_alloc() /= 0 )   CALL ctl_stop( 'STOP', 'bdy_init : unable to allocate oce arrays' )
+      IF( bdy_dta_alloc() /= 0 )   CALL ctl_stop( 'STOP', 'bdy_init : unable to allocate dta arrays' )
+
       IF( jperio /= 0 )   CALL ctl_stop( 'Cyclic or symmetric,',   &
          &                               ' and unstructured open boundary condition are not compatible' )