Changeset 2715 for trunk/NEMOGCM/NEMO/OPA_SRC/DOM/domngb.F90

Timestamp:

2011-03-30T17:58:35+02:00 (10 years ago)

Author:

rblod

Message:

First attempt to put dynamic allocation on the trunk

File:

• trunk/NEMOGCM/NEMO/OPA_SRC/DOM/domngb.F90 (modified) (3 diffs)

trunk/NEMOGCM/NEMO/OPA_SRC/DOM/domngb.F90

@@ -8 +8 @@
 
    !!----------------------------------------------------------------------
-   !!   dom_ngb       : find the closest grid point from a given on/lat position
+   !!   dom_ngb       : find the closest grid point from a given lon/lat position
    !!----------------------------------------------------------------------
-   USE dom_oce         ! ocean space and time domain
-   USE lib_mpp         ! for mppsum
+   USE dom_oce        ! ocean space and time domain
+   USE lib_mpp        ! for mppsum
 
    IMPLICIT NONE
    PRIVATE
 
-   PUBLIC   dom_ngb    ! routine called in iom.F90 module
+   PUBLIC   dom_ngb   ! routine called in iom.F90 module
 
    !!----------------------------------------------------------------------
    !! NEMO/OPA 3.3 , NEMO Consortium (2010)
    !! $Id$ 
-   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
+   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt)
    !!----------------------------------------------------------------------
-
 CONTAINS
 
@@ -30 +29 @@
       !!                    ***  ROUTINE dom_ngb  ***
       !!
-      !! ** Purpose :   find the closest grid point from a given on/lat position
+      !! ** Purpose :   find the closest grid point from a given lon/lat position
       !!
       !! ** Method  :   look for minimum distance in cylindrical projection 
       !!                -> not good if located at too high latitude...
-      !!
       !!----------------------------------------------------------------------
+      USE wrk_nemo, ONLY:   wrk_in_use, wrk_not_released
+      USE wrk_nemo, ONLY:   zglam => wrk_2d_2 , zgphi => wrk_2d_3 , zmask => wrk_2d_4 , zdist => wrk_2d_5
+      !
       REAL(wp)        , INTENT(in   ) ::   plon, plat   ! longitude,latitude of the point
       INTEGER         , INTENT(  out) ::   kii, kjj     ! i-,j-index of the closes grid point
       CHARACTER(len=1), INTENT(in   ) ::   cdgrid       ! grid name 'T', 'U', 'V', 'W'
-      !!
-      INTEGER , DIMENSION(2)        ::   iloc
-      REAL(wp), DIMENSION(jpi,jpj)  ::   zglam, zgphi, zmask, zdist
-      REAL(wp)                      ::   zlon
-      REAL(wp)                      ::   zmini
+      !
+      INTEGER , DIMENSION(2) ::   iloc
+      REAL(wp)               ::   zlon, zmini
       !!--------------------------------------------------------------------
-      
-      zmask(:,:) = 0.
+      !
+      IF( wrk_in_use(2, 2,3,4,5) )   CALL ctl_stop('dom_ngb: Requested workspaces already in use')
+      !
+      zmask(:,:) = 0._wp
       SELECT CASE( cdgrid )
       CASE( 'U' )  ; zglam(:,:) = glamu(:,:) ; zgphi(:,:) = gphiu(:,:) ; zmask(nldi:nlei,nldj:nlej) = umask(nldi:nlei,nldj:nlej,1)
@@ -70 +71 @@
          kjj = iloc(2) + njmpp - 1
       ENDIF
-
+      !
+      IF( wrk_not_released(2, 2,3,4,5) )   CALL ctl_stop('dom_ngb: error releasing workspaces')
+      !
    END SUBROUTINE dom_ngb