Changeset 5654 for branches/UKMO/2015_V36_STABLE_CO6_CO5_zenv_pomsdwl/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90

Timestamp:

2015-07-30T18:45:28+02:00 (9 years ago)

Author:

deazer

Message:

Added in changes to allow CO5 like run to verify CO6 with NEMO STABLE VN3.6

File:

• branches/UKMO/2015_V36_STABLE_CO6_CO5_zenv_pomsdwl/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90 (modified) (8 diffs)

branches/UKMO/2015_V36_STABLE_CO6_CO5_zenv_pomsdwl/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90

@@ -102 +102 @@
       INTEGER ::   ios
       !
-      NAMELIST/namzgr/ ln_zco, ln_zps, ln_sco, ln_isfcav
+      NAMELIST/namzgr/ ln_zco, ln_zps, ln_sco, ln_isfcav, ln_read_zenv
       !!----------------------------------------------------------------------
       !
@@ -125 +125 @@
          WRITE(numout,*) '             s- or hybrid z-s-coordinate    ln_sco    = ', ln_sco
          WRITE(numout,*) '             ice shelf cavities             ln_isfcav = ', ln_isfcav
+         WRITE(numout,*) '             Read zenv from Bathy T/F ln_read_zenv    = ', ln_read_zenv
       ENDIF
 
@@ -518 +519 @@
             ELSE
                CALL iom_get  ( inum, jpdom_data, 'Bathymetry'    , bathy, lrowattr=ln_use_jattr  )
+               IF ( ln_read_zenv ) THEN                  ! Whether we should read zenv or not
+                   CALL iom_get  ( inum, jpdom_data, 'zenv', zenv )
+               ENDIF
             END IF
             CALL iom_close( inum )
@@ -1798 +1802 @@
       !
       REAL(wp), POINTER, DIMENSION(:,:  ) :: ztmpi1, ztmpi2, ztmpj1, ztmpj2
-      REAL(wp), POINTER, DIMENSION(:,:  ) :: zenv, ztmp, zmsk, zri, zrj, zhbat
+      REAL(wp), POINTER, DIMENSION(:,:  ) :: ztmp, zmsk, zri, zrj, zhbat
 
       NAMELIST/namzgr_sco/ln_s_sh94, ln_s_sf12, ln_sigcrit, rn_sbot_min, rn_sbot_max, rn_hc, rn_rmax,rn_theta, &
@@ -1860 +1864 @@
       !                                        ! =============================
       ! use r-value to create hybrid coordinates
+      scosrf(:,:) = 0._wp             ! ocean surface depth (here zero: no under ice-shelf sea)
+      scobot(:,:) = bathy(:,:)        ! ocean bottom  depth
+      IF( ln_read_zenv) THEN
+           WRITE(numout,*) '      Zenv is not calculated but read from Bathy File ln_read_zenv        = ', ln_read_zenv
+      ELSE
       zenv(:,:) = bathy(:,:)
       !
@@ -1881 +1890 @@
       ! 
       ! smooth the bathymetry (if required)
-      scosrf(:,:) = 0._wp             ! ocean surface depth (here zero: no under ice-shelf sea)
-      scobot(:,:) = bathy(:,:)        ! ocean bottom  depth
       !
       jl = 0
@@ -1946 +1953 @@
       END DO                                                !     End loop     !
       !                                                     ! ================ !
+      ENDIF ! End of IF Block for reading in zenv
+      CALL lbc_lnk( zenv, 'T', 1._wp, 'no0' )
       DO jj = 1, jpj
          DO ji = 1, jpi
@@ -2221 +2230 @@
       END DO
       !
-      CALL wrk_dealloc( jpi, jpj, zenv, ztmp, zmsk, zri, zrj, zhbat , ztmpi1, ztmpi2, ztmpj1, ztmpj2 )
+      CALL wrk_dealloc( jpi, jpj, ztmp, zmsk, zri, zrj, zhbat , ztmpi1, ztmpi2, ztmpj1, ztmpj2 )
       !
       IF( nn_timing == 1 )  CALL timing_stop('zgr_sco')