- Timestamp:
- 2018-08-08T17:12:35+02:00 (6 years ago)
- Location:
- branches/NERC/dev_r5518_GO6_under_ice_relax/NEMOGCM/NEMO/OPA_SRC/DOM
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/NERC/dev_r5518_GO6_under_ice_relax/NEMOGCM/NEMO/OPA_SRC/DOM/dom_oce.F90
r8427 r10047 254 254 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: bathy !: ocean depth (meters) 255 255 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: tmask_i, umask_i, vmask_i, fmask_i !: interior domain T-point mask 256 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: tmask_i_diag !: partial mask for use in T diagnostic mask calc. 257 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: umask_i_diag !: partial mask for use in U diagnostic mask calc. 258 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: vmask_i_diag !: partial mask for use in V diagnostic mask calc. 256 259 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: bmask !: land/ocean mask of barotropic stream function 257 260 … … 406 409 407 410 ALLOCATE( tmask(jpi,jpj,jpk) , umask(jpi,jpj,jpk), & 408 & vmask(jpi,jpj,jpk) , fmask(jpi,jpj,jpk), STAT=ierr(11) ) 411 & vmask(jpi,jpj,jpk) , fmask(jpi,jpj,jpk), & 412 tmask_i_diag(jpi,jpj,jpk), & 413 umask_i_diag(jpi,jpj,jpk), & 414 vmask_i_diag(jpi,jpj,jpk), & 415 STAT=ierr(11) ) 409 416 410 417 ALLOCATE( wmask(jpi,jpj,jpk) , wumask(jpi,jpj,jpk), wvmask(jpi,jpj,jpk) , STAT=ierr(12) ) -
branches/NERC/dev_r5518_GO6_under_ice_relax/NEMOGCM/NEMO/OPA_SRC/DOM/dommsk.F90
r8280 r10047 1 1 2 MODULE dommsk 2 3 !!====================================================================== … … 30 31 USE dynspg_oce ! choice/control of key cpp for surface pressure gradient 31 32 USE wrk_nemo ! Memory allocation 33 USE domwri 32 34 USE timing ! Timing 33 35 … … 138 140 REAL(wp) :: zphi_drake_passage, zshlat_antarc 139 141 REAL(wp), POINTER, DIMENSION(:,:) :: zwf 142 REAL(wp) :: uvt(jpi,jpj) ! dummy array for masking purposes. 140 143 !! 141 144 NAMELIST/namlbc/ rn_shlat, ln_vorlat … … 223 226 ! -------------------- 224 227 tmask_i(:,:) = ssmask(:,:) ! (ISH) tmask_i = 1 even on the ice shelf 228 225 229 iif = jpreci ! ??? 226 230 iil = nlci - jpreci + 1 … … 246 250 ENDIF 247 251 ENDIF 252 253 248 254 IF( jperio == 5 .OR. jperio == 6 ) THEN ! F-point pivot 249 255 tpol( 1 :jpiglo) = 0._wp … … 263 269 & * tmask(ji,jj+1,jk) * tmask(ji+1,jj+1,jk) 264 270 END DO 265 END DO 271 END DO 266 272 END DO 267 273 ! (ISF) MIN(1,SUM(umask)) is here to check if you have effectively at least 1 wet u point … … 282 288 CALL lbc_lnk( vmask_i, 'V', 1._wp ) 283 289 CALL lbc_lnk( fmask_i, 'F', 1._wp ) 290 291 292 ! Set up mask for diagnostics on T points, to exclude duplicate 293 ! data points in wrap and N-fold regions. 294 CALL dom_uniq( uvt, 'T' ) 295 DO jk = 1, jpk 296 tmask_i_diag(:,:,jk) = tmask(:,:,jk) * uvt(:,:) 297 END DO 298 299 ! Set up mask for diagnostics on U points, to exclude duplicate 300 ! data points in wrap and N-fold regions. 301 umask_i_diag(:,:,:) = 1.0 302 umask_i_diag(2:jpim1,:,:) = tmask(2:jpim1,:,:) + tmask(3:jpi,:,:) 303 CALL lbc_lnk( umask_i_diag, 'U', 1. ) 304 305 ! Now mask out any duplicate points 306 CALL dom_uniq( uvt, 'U' ) 307 DO jk = 1, jpk 308 umask_i_diag(:,:,jk) = umask_i_diag(:,:,jk) * uvt(:,:) 309 END DO 310 311 312 ! Set up mask for diagnostics on V points, to exclude duplicate 313 ! data points in wrap and N-fold regions. 314 vmask_i_diag(:,:,:) = 1.0 315 vmask_i_diag(:,2:jpjm1,:) = tmask(:,2:jpjm1,:) + tmask(:,3:jpj,:) 316 CALL lbc_lnk( vmask_i_diag, 'V', 1. ) 317 318 ! Now mask out any duplicate points 319 CALL dom_uniq( uvt, 'V' ) 320 DO jk = 1, jpk 321 vmask_i_diag(:,:,jk) = vmask_i_diag(:,:,jk) * uvt(:,:) 322 END DO 323 324 284 325 285 326 ! 3. Ocean/land mask at wu-, wv- and w points -
branches/NERC/dev_r5518_GO6_under_ice_relax/NEMOGCM/NEMO/OPA_SRC/DOM/domwri.F90
r9321 r10047 26 26 PRIVATE 27 27 28 PUBLIC dom_wri ! routine called by inidom.F9028 PUBLIC dom_wri, dom_uniq ! routines called by inidom.F90 and iom.F90 29 29 30 30 !! * Substitutions
Note: See TracChangeset
for help on using the changeset viewer.