Changeset 5260 for branches/2014/dev_r4650_UKMO10_Tidally_Meaned_Diagnostics/NEMOGCM/NEMO/OFF_SRC/dommsk.F90
- Timestamp:
- 2015-05-12T12:37:15+02:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_r4650_UKMO10_Tidally_Meaned_Diagnostics/NEMOGCM/NEMO/OFF_SRC/dommsk.F90
r3294 r5260 15 15 USE in_out_manager ! I/O manager 16 16 USE wrk_nemo 17 USE lbclnk 17 18 18 19 IMPLICIT NONE … … 48 49 !!---------------------------------------------------------------------- 49 50 ! 50 INTEGER :: ji, j k ! dummy loop indices51 INTEGER :: ji, jj, jk ! dummy loop indices 51 52 INTEGER :: iif, iil, ijf, ijl ! local integers 52 53 INTEGER, POINTER, DIMENSION(:,:) :: imsk … … 58 59 ! Interior domain mask (used for global sum) 59 60 ! -------------------- 61 ssmask(:,:) = tmask(:,:,1) 60 62 tmask_i(:,:) = tmask(:,:,1) 61 63 iif = jpreci ! thickness of exchange halos in i-axis … … 82 84 ENDIF 83 85 ! 86 ! (ISF) MIN(1,SUM(umask)) is here to check if you have effectively at 87 ! least 1 wet u point 88 DO jj = 1, jpjm1 89 DO ji = 1, fs_jpim1 ! vector loop 90 umask_i(ji,jj) = ssmask(ji,jj) * ssmask(ji+1,jj ) * MIN(1._wp,SUM(umask(ji,jj,:))) 91 vmask_i(ji,jj) = ssmask(ji,jj) * ssmask(ji ,jj+1) * MIN(1._wp,SUM(vmask(ji,jj,:))) 92 END DO 93 DO ji = 1, jpim1 ! NO vector opt. 94 fmask_i(ji,jj) = ssmask(ji,jj ) * ssmask(ji+1,jj ) & 95 & * ssmask(ji,jj+1) * ssmask(ji+1,jj+1) * MIN(1._wp,SUM(fmask(ji,jj,:))) 96 END DO 97 END DO 98 CALL lbc_lnk( umask_i, 'U', 1._wp ) ! Lateral boundary conditions 99 CALL lbc_lnk( vmask_i, 'V', 1._wp ) 100 CALL lbc_lnk( fmask_i, 'F', 1._wp ) 101 102 ! 3. Ocean/land mask at wu-, wv- and w points 103 !---------------------------------------------- 104 wmask (:,:,1) = tmask(:,:,1) ! ???????? 105 wumask(:,:,1) = umask(:,:,1) ! ???????? 106 wvmask(:,:,1) = vmask(:,:,1) ! ???????? 107 DO jk=2,jpk 108 wmask (:,:,jk)=tmask(:,:,jk) * tmask(:,:,jk-1) 109 wumask(:,:,jk)=umask(:,:,jk) * umask(:,:,jk-1) 110 wvmask(:,:,jk)=vmask(:,:,jk) * vmask(:,:,jk-1) 111 END DO 112 ! 84 113 IF( nprint == 1 .AND. lwp ) THEN ! Control print 85 114 imsk(:,:) = INT( tmask_i(:,:) )
Note: See TracChangeset
for help on using the changeset viewer.