Changeset 291 for trunk/NEMO/OPA_SRC/DOM
- Timestamp:
- 2005-09-22T12:47:37+02:00 (19 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/DOM/dommsk.F90
r255 r291 173 173 END DO 174 174 175 176 175 #if defined key_zdfkpp 177 ! land point on Bab el Mandeb zonal section 178 ij0 = 87 ; ij1 = 88 179 ii0 = 160 ; ii1 = 161 180 tmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , 1:jpk ) = 0.e0 176 IF( cp_cfg == 'orca' ) THEN 177 IF( jp_cfg = 2 ) THEN 178 ! land point on Bab el Mandeb zonal section 179 ij0 = 87 ; ij1 = 88 180 ii0 = 160 ; ii1 = 161 181 tmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , 1:jpk ) = 0.e0 182 ELSE 183 IF(lwp) WRITE(numout,*) 184 IF(lwp) WRITE(numout,cform_war) 185 IF(lwp) WRITE(numout,*) 186 IF(lwp) WRITE(numout,*)' A mask must be applied on Bab el Mandeb strait' 187 IF(lwp) WRITE(numout,*)' in case of ORCAs configurations' 188 IF(lwp) WRITE(numout,*)' This is a problem which is not yet solved' 189 IF(lwp) WRITE(numout,*) 190 ENDIF 191 ENDIF 181 192 #endif 193 182 194 ! Interior domain mask (used for global sum) 183 195 ! -------------------- … … 200 212 tpol(jpiglo/2+1:jpiglo) = 0.e0 201 213 fpol( 1 :jpiglo) = 0.e0 214 ! T-point pivot: only half of the nlcj-1 row 215 IF( mjg(nlej) == jpjglo ) THEN 216 DO ji = iif+1, iil-1 217 tmask_i(ji,nlej-1) = tmask_i(ji,nlej-1) * tpol(mig(ji)) 218 END DO 219 ENDIF 202 220 ENDIF 203 221 IF( jperio == 5 .OR. jperio == 6 ) THEN ! F-point pivot … … 205 223 fpol(jpiglo/2+1:jpiglo) = 0.e0 206 224 ENDIF 207 208 IF( nperio == 3 .OR. nperio == 4 ) THEN ! T-point pivot: only half of the nlcj-1 row209 DO ji = iif+1, iil-1210 tmask_i(ji,ijl-1) = tmask_i(ji,ijl-1) * tpol(mig(ji))211 END DO212 ENDIF213 214 225 215 226 ! 2. Ocean/land mask at u-, v-, and z-points (computed from tmask)
Note: See TracChangeset
for help on using the changeset viewer.