Changeset 6596 for branches/2016/dev_r6409_SIMPLIF_2_usrdef/NEMOGCM/NEMO/OPA_SRC/CRS

Timestamp:

2016-05-22T20:11:06+02:00 (8 years ago)

Author:

gm

Message:

#1692 - branch SIMPLIF_2_usrdef: remove from namcfg and namdom many obsolete variables ; remove izoom/jzoom option

Location:

branches/2016/dev_r6409_SIMPLIF_2_usrdef/NEMOGCM/NEMO/OPA_SRC/CRS

Files:

• crsdom.F90 (modified) (1 diff)
• crsini.F90 (modified) (2 diffs)

branches/2016/dev_r6409_SIMPLIF_2_usrdef/NEMOGCM/NEMO/OPA_SRC/CRS/crsdom.F90

@@ -2023 +2023 @@
          nimpp_crs = nimppt_crs(nproc + 1)
 
-         ! No coarsening with zoom
-         IF( jpizoom /= 1 .OR. jpjzoom /= 1)    STOP 
-
          DO ji = 1, jpi_crs
             mig_crs(ji) = ji + nimpp_crs - 1

branches/2016/dev_r6409_SIMPLIF_2_usrdef/NEMOGCM/NEMO/OPA_SRC/CRS/crsini.F90

@@ -11 +11 @@
    !!----------------------------------------------------------------------
    USE par_kind, ONLY: wp
-   USE par_oce                  ! For parameter jpi,jpj,jphgr_msh
+   USE par_oce                  ! For parameter jpi,jpj
    USE dom_oce                  ! For parameters in par_oce
    USE crs                      ! Coarse grid domain
@@ -170 +170 @@
      !      3.c.2 Coriolis factor  
 
-      SELECT CASE( jphgr_msh )   ! type of horizontal mesh
-
-      CASE ( 0, 1, 4 )           ! mesh on the sphere
-
-         ff_crs(:,:) = 2. * omega * SIN( rad * gphif_crs(:,:) )
-
-      CASE DEFAULT 
-
-       IF(lwp)    WRITE(numout,*) 'crsini.F90. crs_init. Only jphgr_msh = 0, 1 or 4 supported' 
+!!gm  Not sure CRS needs Coriolis parameter....
+!!gm  If needed, then update this to have Coriolis at both f- and t-points
+
+      ff_crs(:,:) = 2. * omega * SIN( rad * gphif_crs(:,:) )
+
+      CALL ctl_warn( 'crsini: CAUTION, CRS only designed for Coriolis defined on the sphere' ) 
  
-      END SELECT 
 
      !    3.d.1 mbathy ( vertical k-levels of bathymetry )