Changeset 2444 for branches/nemo_v3_3_beta/NEMOGCM/NEMO/OFF_SRC/domrea.F90

Timestamp:

2010-11-29T15:30:48+01:00 (14 years ago)

Author:

cetlod

Message:

Improvment of OFFLINE in v3.3beta (review done by gm) : clean the style in all routines, suppression of key_zdfddm

File:

• branches/nemo_v3_3_beta/NEMOGCM/NEMO/OFF_SRC/domrea.F90 (modified) (7 diffs)

branches/nemo_v3_3_beta/NEMOGCM/NEMO/OFF_SRC/domrea.F90

@@ -11 +11 @@
    !!                         = 3  :   mesh_hgr, mesh_zgr and mask
    !!----------------------------------------------------------------------
-   !! * Modules used
    USE dom_oce         ! ocean space and time domain
-   USE dommsk  
-   USE in_out_manager
+   USE dommsk          ! domain: masks
+   USE in_out_manager  ! I/O manager
 
    IMPLICIT NONE
    PRIVATE
 
-   !! * Accessibility
-   PUBLIC dom_rea        ! routine called by inidom.F90
+   PUBLIC   dom_rea    ! routine called by inidom.F90
    !!----------------------------------------------------------------------
    !! NEMO/OFF 3.3 , NEMO Consortium (2010)
@@ -29 +27 @@
 CONTAINS
 
-#if ( defined key_mpp_mpi || defined key_mpp_shmem ) && defined key_dimgout
+#if   defined key_mpp_mpi   &&   defined key_dimgout
    !!----------------------------------------------------------------------
    !!   'key_mpp_mpi'     OR
-   !!   'key_mpp_shmem'
    !!   'key_dimgout' :         each processor makes its own direct access file 
    !!                      use build_nc_meshmask off line to retrieve 
@@ -38 +35 @@
    !!----------------------------------------------------------------------
 #  include "domrea_dimg.h90"
-
 
 #else
@@ -67 +63 @@
       !!      meshmask.nc  : domain size, horizontal grid-point position,
       !!                     masks, depth and vertical scale factors
+      !!----------------------------------------------------------------------
+      USE iom
       !!
-      !! History :
-      !!        !  97-02  (G. Madec)  Original code
-      !!        !  99-11  (M. Imbard)  NetCDF FORMAT with IOIPSL
-      !!   9.0  !  02-08  (G. Madec)  F90 and several file
-      !!        !  06-07  (C. Ethe )  Use of iom module
+      INTEGER  ::   ji, jj, jk   ! dummy loop indices
+      INTEGER  ::   ik, inum0 , inum1 , inum2 , inum3 , inum4   ! local integers
+      REAL(wp) ::   zrefdep         ! local real
+      REAL(wp), DIMENSION(jpi,jpj) ::   zprt   ! 2D workspace
       !!----------------------------------------------------------------------
-      !! * Modules used
-      USE iom
-
-      !! * Local declarations
-      INTEGER  ::   ji, jj, jk
-      INTEGER  ::                & !!! * temprary units for :
-         inum0 ,                 &  ! 'mesh_mask.nc' file
-         inum1 ,                 &  ! 'mesh.nc'      file
-         inum2 ,                 &  ! 'mask.nc'      file
-         inum3 ,                 &  ! 'mesh_hgr.nc'  file
-         inum4                      ! 'mesh_zgr.nc'  file
- 
-      REAL(wp), DIMENSION(jpi,jpj) :: zprt
-      REAL(wp) ::   zrefdep         ! depth of the reference level (~10m)
-      INTEGER  :: ik
-      !!----------------------------------------------------------------------
-
-       IF(lwp) WRITE(numout,*)
-       IF(lwp) WRITE(numout,*) 'dom_rea : read NetCDF mesh and mask information file(s)'
-       IF(lwp) WRITE(numout,*) '~~~~~~~'
-
-
-      zprt(:,:) = 0.
+
+      IF(lwp) WRITE(numout,*)
+      IF(lwp) WRITE(numout,*) 'dom_rea : read NetCDF mesh and mask information file(s)'
+      IF(lwp) WRITE(numout,*) '~~~~~~~'
+
+      zprt(:,:) = 0._wp
 
       SELECT CASE (nmsh)
@@ -180 +160 @@
          DO jj = 1, jpj
             DO ji = 1, jpi
-               mbathy(ji,jj) = MAX( zprt(ji,jj) * tmask(ji,jj,1), 1. ) + 1
+               mbathy(ji,jj) = MAX( zprt(ji,jj) * tmask(ji,jj,1), 1._wp ) + 1
             ENDDO
          ENDDO
@@ -262 +242 @@
 !!gm BUG in s-coordinate this does not work!
       ! deepest/shallowest W level Above/Below ~10m
-      zrefdep = 10. - ( 0.1*MINVAL(e3w_0) )                          ! ref. depth with tolerance (10% of minimum layer thickness)
+      zrefdep = 10._wp - ( 0.1_wp * MINVAL(e3w_0) )                  ! ref. depth with tolerance (10% of minimum layer thickness)
       nlb10 = MINLOC( gdepw_0, mask = gdepw_0 > zrefdep, dim = 1 )   ! shallowest W level Below ~10m
       nla10 = nlb10 - 1                                              ! deepest    W level Above ~10m
@@ -308 +288 @@
 
       DO jk = 1, jpk
-         IF( e3w_0(jk) <= 0. .OR. e3t_0(jk) <= 0. ) CALL ctl_stop ( ' e3w_0 or e3t_0 =< 0 ' )
-         IF( gdepw_0(jk) < 0. .OR. gdept_0(jk) < 0.) CALL ctl_stop( ' gdepw_0 or gdept_0 < 0 ' )
+         IF( e3w_0  (jk) <= 0._wp .OR. e3t_0  (jk) <= 0._wp )   CALL ctl_stop( ' e3w_0 or e3t_0 =< 0 ' )
+         IF( gdepw_0(jk) <  0._wp .OR. gdept_0(jk) <  0._wp )   CALL ctl_stop( ' gdepw_0 or gdept_0 < 0 ' )
       END DO
-
-         !                                     ! ============================
-         !                                     !        close the files 
-         !                                     ! ============================
-         SELECT CASE ( nmsh )
-            CASE ( 1 )                
-               CALL iom_close( inum0 )
-            CASE ( 2 )
-               CALL iom_close( inum1 )
-               CALL iom_close( inum2 )
-            CASE ( 3 )
-               CALL iom_close( inum2 )
-               CALL iom_close( inum3 )
-               CALL iom_close( inum4 )
-         END SELECT
-
+      !                                     ! ============================
+      !                                     !        close the files 
+      !                                     ! ============================
+      SELECT CASE ( nmsh )
+         CASE ( 1 )                
+            CALL iom_close( inum0 )
+         CASE ( 2 )
+            CALL iom_close( inum1 )
+            CALL iom_close( inum2 )
+         CASE ( 3 )
+            CALL iom_close( inum2 )
+            CALL iom_close( inum3 )
+            CALL iom_close( inum4 )
+      END SELECT
+      !
    END SUBROUTINE dom_rea