- Timestamp:
- 2020-06-22T11:13:03+02:00 (4 years ago)
- Location:
- NEMO/branches/2020/dev_r12558_HPC-08_epico_Extra_Halo/src/OCE
- Files:
-
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r12558_HPC-08_epico_Extra_Halo/src/OCE/DOM/dommsk.F90
r13065 r13138 138 138 END_2D 139 139 ! 140 ! the following is mandatory 141 ! it masks boundaries (bathy=0) where needed depending on the configuration (closed, periodic...) 142 !!$ IF( .NOT. (nbondj == 1 .OR. nbondj == 0 .OR. l_Jperio) ) THEN 143 !!$ tmask(mi0( 1 ):mi1(jpiglo),mj0(Njs0):mj1(Njs0 ),:) = 0._wp ! line number Njs0 at 0 144 !!$ ENDIF 145 !!$ IF( .NOT. (nbondi == 1 .OR. nbondi == 0 .OR. l_Iperio) ) THEN 146 !!$ tmask(mi0(Nis0):mi1( Nis0),mj0( 1 ):mj1(jpjglo),:) = 0._wp ! column number Nis0 at 0 147 !!$ ENDIF 148 !!$ CALL lbc_lnk( 'dommsk', tmask, 'T', 1._wp ) ! Lateral boundary conditions 149 150 ! Mask corrections for bdy (read in mppini2) 140 ! Mask corrections for bdy (read in mppini2) 151 141 READ ( numnam_ref, nambdy, IOSTAT = ios, ERR = 903) 152 142 903 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nambdy in reference namelist' ) -
NEMO/branches/2020/dev_r12558_HPC-08_epico_Extra_Halo/src/OCE/DOM/domzgr.F90
r13065 r13138 111 111 ENDIF 112 112 ! 113 ! the following is mandatory 113 114 ! make sure that closed boundaries are correctly defined in k_top that will be used to compute all mask arrays 114 115 ! -
NEMO/branches/2020/dev_r12558_HPC-08_epico_Extra_Halo/src/OCE/DYN/sshwzv.F90
r12980 r13138 200 200 ENDIF 201 201 ! 202 #if defined key_agrif203 202 IF( .NOT. AGRIF_Root() ) THEN 204 203 ! 205 204 ! Mask vertical velocity at first/last columns/row 206 205 ! inside computational domain (cosmetic) 207 ! --- West --- ! 208 DO ji = mi0(2), mi1(2) 209 DO jj = 1, jpj 210 pww(ji,jj,:) = 0._wp 211 ENDDO 212 ENDDO 213 ! 214 ! --- East --- ! 215 DO ji = mi0(jpiglo-1), mi1(jpiglo-1) 216 DO jj = 1, jpj 217 pww(ji,jj,:) = 0._wp 218 ENDDO 219 ENDDO 220 ! 221 ! --- South --- ! 222 DO jj = mj0(2), mj1(2) 223 DO ji = 1, jpi 224 pww(ji,jj,:) = 0._wp 225 ENDDO 226 ENDDO 227 ! 228 ! --- North --- ! 229 DO jj = mj0(jpjglo-1), mj1(jpjglo-1) 230 DO ji = 1, jpi 231 pww(ji,jj,:) = 0._wp 232 ENDDO 233 ENDDO 206 DO jk = 1, jpkm1 207 ! --- West --- ! 208 DO ji = mi0(2+nn_hls), mi1(2+nn_hls) 209 DO jj = 1, jpj 210 pww(ji,jj,jk) = 0._wp 211 END DO 212 END DO 213 ! 214 ! --- East --- ! 215 DO ji = mi0(jpiglo-1-nn_hls), mi1(jpiglo-1-nn_hls) 216 DO jj = 1, jpj 217 pww(ji,jj,jk) = 0._wp 218 END DO 219 END DO 220 ! 221 ! --- South --- ! 222 DO jj = mj0(2+nn_hls), mj1(2+nn_hls) 223 DO ji = 1, jpi 224 pww(ji,jj,jk) = 0._wp 225 END DO 226 END DO 227 ! 228 ! --- North --- ! 229 DO jj = mj0(jpjglo-1-nn_hls), mj1(jpjglo-1-nn_hls) 230 DO ji = 1, jpi 231 pww(ji,jj,jk) = 0._wp 232 END DO 233 END DO 234 END DO 234 235 ! 235 236 ENDIF 236 #endif237 237 ! 238 238 IF( ln_timing ) CALL timing_stop('wzv') -
NEMO/branches/2020/dev_r12558_HPC-08_epico_Extra_Halo/src/OCE/USR/usrdef_fmask.F90
r12960 r13138 76 76 ! 77 77 IF(lwp) WRITE(numout,*) ' Bab el Mandeb ' 78 ij0 = 87 + nn_hls ; ij1 = 79 ii0 = 160 + nn_hls - 1; ii1 = 160 + nn_hls - 178 ij0 = 87 + nn_hls ; ij1 = 88 + nn_hls ! Bab el Mandeb : partial slip (pfmsk=1) 79 ii0 = 160 + nn_hls - 1 ; ii1 = 160 + nn_hls - 1 80 80 pfmsk( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , 1:jpk ) = 1._wp 81 81 ij0 = 88 + nn_hls ; ij1 = 88 + nn_hls … … 98 98 !!gm ! Currently these hard-wired indices relate to configuration with extend grid (jpjglo=332) 99 99 ! 100 isrow = 332 - Nj0glo - 2100 isrow = 332 - (Nj0glo + 1) ! was 332 - jpjglo -> jpjglo_old_version = Nj0glo + 1 101 101 ! 102 102 IF(lwp) WRITE(numout,*) 103 103 IF(lwp) WRITE(numout,*) ' orca_r1: increase friction near the following straits : ' 104 104 IF(lwp) WRITE(numout,*) ' Gibraltar ' 105 ii0 = 282 ; ii1 = 283 ! Gibraltar Strait 106 ij0 = 241 - isrow ; ij1 = 241 - isrow ; pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp 105 ii0 = 282 + nn_hls - 1 ; ii1 = 283 + nn_hls - 1 ! Gibraltar Strait 106 ij0 = 241 + nn_hls - isrow ; ij1 = 241 + nn_hls - isrow 107 pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp 107 108 ! 108 109 IF(lwp) WRITE(numout,*) ' Bhosporus ' 109 ii0 = 314 ; ii1 = 315 ! Bhosporus Strait 110 ij0 = 248 - isrow ; ij1 = 248 - isrow ; pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp 110 ii0 = 314 + nn_hls - 1 ; ii1 = 315 + nn_hls - 1 ! Bhosporus Strait 111 ij0 = 248 + nn_hls - isrow ; ij1 = 248 + nn_hls - isrow 112 pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp 111 113 ! 112 114 IF(lwp) WRITE(numout,*) ' Makassar (Top) ' 113 ii0 = 48 ; ii1 = 48 ! Makassar Strait (Top) 114 ij0 = 189 - isrow ; ij1 = 190 - isrow ; pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 3._wp 115 ii0 = 48 + nn_hls - 1 ; ii1 = 48 + nn_hls - 1 ! Makassar Strait (Top) 116 ij0 = 189 + nn_hls - isrow ; ij1 = 190 + nn_hls - isrow 117 pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 3._wp 115 118 ! 116 119 IF(lwp) WRITE(numout,*) ' Lombok ' 117 ii0 = 44 ; ii1 = 44 ! Lombok Strait 118 ij0 = 164 - isrow ; ij1 = 165 - isrow ; pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp 120 ii0 = 44 + nn_hls - 1 ; ii1 = 44 + nn_hls - 1 ! Lombok Strait 121 ij0 = 164 + nn_hls - isrow ; ij1 = 165 + nn_hls - isrow 122 pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp 119 123 ! 120 124 IF(lwp) WRITE(numout,*) ' Ombai ' 121 ii0 = 53 ; ii1 = 53 ! Ombai Strait 122 ij0 = 164 - isrow ; ij1 = 165 - isrow ; pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp 125 ii0 = 53 + nn_hls - 1 ; ii1 = 53 + nn_hls - 1 ! Ombai Strait 126 ij0 = 164 + nn_hls - isrow ; ij1 = 165 + nn_hls - isrow 127 pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp 123 128 ! 124 129 IF(lwp) WRITE(numout,*) ' Timor Passage ' 125 ii0 = 56 ; ii1 = 56 ! Timor Passage 126 ij0 = 164 - isrow ; ij1 = 165 - isrow ; pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp 130 ii0 = 56 + nn_hls - 1 ; ii1 = 56 + nn_hls - 1 ! Timor Passage 131 ij0 = 164 + nn_hls - isrow ; ij1 = 165 + nn_hls - isrow 132 pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 2._wp 127 133 ! 128 134 IF(lwp) WRITE(numout,*) ' West Halmahera ' 129 ii0 = 58 ; ii1 = 58 ! West Halmahera Strait 130 ij0 = 181 - isrow ; ij1 = 182 - isrow ; pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 3._wp 135 ii0 = 58 + nn_hls - 1 ; ii1 = 58 + nn_hls - 1 ! West Halmahera Strait 136 ij0 = 181 + nn_hls - isrow ; ij1 = 182 + nn_hls - isrow 137 pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 3._wp 131 138 ! 132 139 IF(lwp) WRITE(numout,*) ' East Halmahera ' 133 ii0 = 55 ; ii1 = 55 ! East Halmahera Strait 134 ij0 = 181 - isrow ; ij1 = 182 - isrow ; pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 3._wp 140 ii0 = 55 + nn_hls - 1 ; ii1 = 55 + nn_hls - 1 ! East Halmahera Strait 141 ij0 = 181 + nn_hls - isrow ; ij1 = 182 + nn_hls - isrow 142 pfmsk( mi0(ii0):mi1(ii1),mj0(ij0):mj1(ij1),1:jpk ) = 3._wp 135 143 ! 136 144 CASE DEFAULT -
NEMO/branches/2020/dev_r12558_HPC-08_epico_Extra_Halo/src/OCE/USR/usrdef_zgr.F90
r12939 r13138 199 199 ! 200 200 z2d(:,:) = REAL( jpkm1 , wp ) ! flat bottom 201 z2d(mi0( 1 ):mi1(jpiglo),mj0(Njs0):mj1( Njs0 )) = 0._wp ! line number Njs0 at 0202 z2d(mi0(Nis0):mi1( Nis0 ),mj0( 1 ):mj1(jpjglo)) = 0._wp ! column number Nis0 at 0203 !204 CALL lbc_lnk( 'usrdef_zgr', z2d, 'T', 1. ) ! set surrounding land to zero (here jperio=0 ==>> closed)205 201 ! 206 202 k_bot(:,:) = NINT( z2d(:,:) ) ! =jpkm1 over the ocean point, =0 elsewhere
Note: See TracChangeset
for help on using the changeset viewer.