Changeset 7277 for branches/2016/dev_CNRS_2016/NEMOGCM/NEMO/OPA_SRC/C1D/domc1d.F90

Timestamp:

2016-11-21T09:55:07+01:00 (8 years ago)

Author:

flavoni

Message:

update 2016 branch with simplif-2

File:

• branches/2016/dev_CNRS_2016/NEMOGCM/NEMO/OPA_SRC/C1D/domc1d.F90 (modified) (4 diffs)

branches/2016/dev_CNRS_2016/NEMOGCM/NEMO/OPA_SRC/C1D/domc1d.F90

@@ -24 +24 @@
 
    PUBLIC   dom_c1d   ! called in domcfg.F90
+
+   INTEGER ::   jpizoom = 1      !: left bottom (i,j) indices of the zoom
+   INTEGER ::   jpjzoom = 1      !: in data domain indices
 
    !!----------------------------------------------------------------------
@@ -82 +85 @@
       !  mesh, only glamt and gphit   !
       ! ============================= !
-
+      !
       SELECT CASE( jphgr_msh )   ! type of horizontal mesh
-
+      !
       CASE ( 0 )                 !  curvilinear coordinate on the sphere read in coordinate.nc file
-
+         !
          CALL iom_open( 'coordinates', inum )
          CALL iom_get( inum, jpdom_unknown, 'glamt', glamdta ) ! mig, mjg undefined at this point
          CALL iom_get( inum, jpdom_unknown, 'gphit', gphidta ) ! so use jpdom_unknown not jpdom_data
          CALL iom_close ( inum )
-
+         !
       CASE ( 1 )                 ! geographical mesh on the sphere with regular grid-spacing
-
+         !
          DO jj = 1, jpjdta
             DO ji = 1, jpidta
                zti = FLOAT( ji - 1 + nimpp - 1 )
                ztj = FLOAT( jj - 1 + njmpp - 1 )
-
+               !
                glamdta(ji,jj) = ppglam0 + ppe1_deg * zti
                gphidta(ji,jj) = ppgphi0 + ppe2_deg * ztj
             END DO
          END DO
-
+         !
       CASE ( 2:3 )               ! f- or beta-plane with regular grid-spacing
-         
+         !
          glam0 = 0.e0
          gphi0 = - ppe2_m * 1.e-3
-
+         !
          DO jj = 1, jpjdta
             DO ji = 1, jpidta
@@ -115 +118 @@
             END DO
          END DO
-
+         !
       CASE ( 4 )                 ! geographical mesh on the sphere, isotropic MERCATOR type
-
+         !
          IF( ppgphi0 == -90 )   CALL ctl_stop( ' Mercator grid cannot start at south pole !!!! ' )
-
+         !
          zarg = rpi / 4. - rpi / 180. * ppgphi0 / 2.
          ijeq = ABS( 180. / rpi * LOG( COS( zarg ) / SIN( zarg ) ) / ppe1_deg )
          IF( ppgphi0 > 0 )   ijeq = -ijeq
-
+         !
          DO jj = 1, jpjdta
             DO ji = 1, jpidta
                zti = FLOAT( ji - 1    + nimpp - 1 )
                ztj = FLOAT( jj - ijeq + njmpp - 1 )
-
+               !
                glamdta(ji,jj) = ppglam0 + ppe1_deg * zti
                gphidta(ji,jj) = 1. / rad * ASIN ( TANH( ppe1_deg * rad * ztj ) )
             END DO
          END DO
-
+         !
       CASE ( 5 )                 ! beta-plane with regular grid-spacing and rotated domain (GYRE configuration)
-   
+         !
          zlam1 = -85
          zphi1 = 29
-         ze1 = 106000. / FLOAT(jp_cfg)
- 
+         ze1 = 106000. / REAL( nn_cfg , wp )
+         !
          zsin_alpha = - SQRT( 2. ) / 2.
          zcos_alpha =   SQRT( 2. ) / 2.
          ze1deg = ze1 / (ra * rad)
-
+         !
          glam0 = zlam1 + zcos_alpha * ze1deg * FLOAT( jpjdta-2 ) ! Force global
          gphi0 = zphi1 + zsin_alpha * ze1deg * FLOAT( jpjdta-2 )
-
+         !
          DO jj = 1, jpjdta
             DO ji = 1, jpidta
@@ -156 +159 @@
             END DO
          END DO
-
+         !
       CASE DEFAULT
-
+         !
          WRITE(ctmp1,*) '          bad flag value for jphgr_msh = ', jphgr_msh
          CALL ctl_stop( ctmp1 )
-
+         !
       END SELECT