Changeset 7320 for branches/2015/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/OPA_SRC/CRS/crsini.F90

Timestamp:

2016-11-23T14:24:07+01:00 (8 years ago)

Author:

cbricaud

Message:

create lib_fortran routine for CRS grid and interface it in TOP

File:

• branches/2015/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/OPA_SRC/CRS/crsini.F90 (modified) (2 diffs)

branches/2015/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/OPA_SRC/CRS/crsini.F90

@@ -67 +67 @@
       !!-------------------------------------------------------------------
       !! Local variables
-      INTEGER  :: ji,jj,jk      ! dummy indices
-      INTEGER  :: ierr                                ! allocation error status
-      INTEGER  ::   ios                 ! Local integer output status for namelist read
+      INTEGER  :: ji,jj,jk           ! dummy indices
+      INTEGER  :: ierr               ! allocation error status
+      INTEGER  :: ios                ! Local integer output status for namelist read
+      INTEGER  :: iif, iil, ijf, ijl ! indices for tmask_i_crs construction
       REAL(wp) :: zmin,zmax
       REAL(wp), DIMENSION(:,:,:), POINTER :: zfse3t, zfse3u, zfse3v, zfse3w
@@ -309 +310 @@
 
      !---------------------------------------------------------
-     ! 5.  Coarse grid ocean volume and averaging weights
-     !---------------------------------------------------------
-     !CALL dom_grid_crs   ! Save the parent grid information  & Switch to coarse grid domain
-     !CALL ldf_tra_crs_init
-     !CALL dom_grid_glo   ! Return to parent grid domain
-
+     ! 5.  Build tmask_i_crs ( for crsdomwri and lib_fortran_crs (glob_sum)
+     !---------------------------------------------------------
+     tmask_i_crs(:,:) = tmask_crs(:,:,1)
+     iif = jpreci
+     iil = nlci_crs - jpreci + 1
+     ijf = jpreci
+     ijl = nlcj_crs - jprecj + 1
+
+     tmask_i_crs( 1:iif ,    :  ) = 0._wp
+     tmask_i_crs(iil:jpi_crs,    :  ) = 0._wp
+     tmask_i_crs(   :   , 1:ijf ) = 0._wp
+     tmask_i_crs(   :   ,ijl:jpj_crs) = 0._wp
+
+     tpol_crs(1:jpiglo_crs,:) = 1._wp
+     fpol_crs(1:jpiglo_crs,:) = 1._wp
+     IF( jperio == 3 .OR. jperio == 4 ) THEN
+        tpol_crs(jpiglo_crs/2+1:jpiglo_crs,:) = 0._wp
+        fpol_crs(       1      :jpiglo_crs,:) = 0._wp
+        IF( mjg_crs(nlej_crs) == jpiglo_crs ) THEN
+           DO ji = iif+1, iil-1
+              tmask_i_crs(ji,nlej_crs-1) = tmask_i_crs(ji,nlej_crs-1) &
+              & * tpol_crs(mig_crs(ji),1)
+           ENDDO
+        ENDIF
+     ENDIF
+     IF( jperio == 5 .OR. jperio == 6 ) THEN
+        tpol_crs(      1       :jpiglo_crs,:)=0._wp
+        fpol_crs(jpiglo_crs/2+1:jpiglo_crs,:)=0._wp
+     ENDIF
 
      !