Changeset 2528 for trunk/NEMOGCM/NEMO/OPA_SRC/DOM/dommsk.F90
- Timestamp:
- 2010-12-27T18:33:53+01:00 (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/OPA_SRC/DOM/dommsk.F90
- Property svn:eol-style deleted
r1707 r2528 5 5 !!====================================================================== 6 6 !! History : OPA ! 1987-07 (G. Madec) Original code 7 !! -! 1993-03 (M. Guyon) symetrical conditions (M. Guyon)8 !! -! 1996-01 (G. Madec) suppression of common work arrays7 !! 6.0 ! 1993-03 (M. Guyon) symetrical conditions (M. Guyon) 8 !! 7.0 ! 1996-01 (G. Madec) suppression of common work arrays 9 9 !! - ! 1996-05 (G. Madec) mask computed from tmask and sup- 10 10 !! ! pression of the double computation of bmask 11 !! -! 1997-02 (G. Madec) mesh information put in domhgr.F12 !! -! 1997-07 (G. Madec) modification of mbathy and fmask11 !! 8.0 ! 1997-02 (G. Madec) mesh information put in domhgr.F 12 !! 8.1 ! 1997-07 (G. Madec) modification of mbathy and fmask 13 13 !! - ! 1998-05 (G. Roullet) free surface 14 !! -! 2000-03 (G. Madec) no slip accurate14 !! 8.2 ! 2000-03 (G. Madec) no slip accurate 15 15 !! - ! 2001-09 (J.-M. Molines) Open boundaries 16 16 !! NEMO 1.0 ! 2002-08 (G. Madec) F90: Free form and module … … 44 44 !! NEMO/OPA 3.2 , LODYC-IPSL (2009) 45 45 !! $Id$ 46 !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) 47 !!---------------------------------------------------------------------- 48 46 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 47 !!---------------------------------------------------------------------- 49 48 CONTAINS 50 49 … … 56 55 !! zontal velocity points (u & v), vorticity points (f) and baro- 57 56 !! tropic stream function points (b). 58 !! Set mbathy to the number of non-zero w-levels of a water column59 57 !! 60 58 !! ** Method : The ocean/land mask is computed from the basin bathy- … … 73 71 !! or mbathy(ji+1,jj) or mbathy(ji+1,jj+1) =< 0 74 72 !! 1. IF mbathy( ji ,jj) and mbathy( ji ,jj+1) 75 !! and mbathy(ji+1,jj) and mbathy(ji+1,jj+1) >= jk.73 !! and mbathy(ji+1,jj) and mbathy(ji+1,jj+1) >= jk. 76 74 !! b-point : the same definition as for f-point of the first ocean 77 75 !! level (surface level) but with 0 along coastlines. 76 !! tmask_i : interior ocean mask at t-point, i.e. excluding duplicated 77 !! rows/lines due to cyclic or North Fold boundaries as well 78 !! as MPP halos. 78 79 !! 79 80 !! The lateral friction is set through the value of fmask along … … 99 100 !! - bmask is set to 0 on the open boundaries. 100 101 !! 101 !! Set mbathy to the number of non-zero w-levels of a water column102 !! mbathy = min( mbathy, 1 ) + 1103 !! (note that the minimum value of mbathy is 2).104 !!105 102 !! ** Action : tmask : land/ocean mask at t-point (=0. or 1.) 106 103 !! umask : land/ocean mask at u-point (=0. or 1.) … … 110 107 !! bmask : land/ocean mask at barotropic stream 111 108 !! function point (=0. or 1.) and set to 0 along lateral boundaries 112 !! mbathy : number of non-zero w-levels109 !! tmask_i : interior ocean mask 113 110 !!---------------------------------------------------------------------- 114 111 INTEGER :: ji, jj, jk ! dummy loop indices … … 132 129 ENDIF 133 130 134 IF ( rn_shlat == 0. ) THEN ; IF(lwp) WRITE(numout,*) ' ocean lateral free-slip '131 IF ( rn_shlat == 0. ) THEN ; IF(lwp) WRITE(numout,*) ' ocean lateral free-slip ' 135 132 ELSEIF ( rn_shlat == 2. ) THEN ; IF(lwp) WRITE(numout,*) ' ocean lateral no-slip ' 136 133 ELSEIF ( 0. < rn_shlat .AND. rn_shlat < 2. ) THEN ; IF(lwp) WRITE(numout,*) ' ocean lateral partial-slip ' … … 145 142 ! N.B. tmask has already the right boundary conditions since mbathy is ok 146 143 ! 147 tmask(:,:,:) = 0. e0144 tmask(:,:,:) = 0._wp 148 145 DO jk = 1, jpk 149 146 DO jj = 1, jpj 150 147 DO ji = 1, jpi 151 IF( REAL( mbathy(ji,jj) - jk ) +.1 >= 0.e0 ) tmask(ji,jj,jk) = 1.e0148 IF( REAL( mbathy(ji,jj) - jk, wp ) + 0.1_wp >= 0._wp ) tmask(ji,jj,jk) = 1._wp 152 149 END DO 153 150 END DO … … 160 157 ij0 = 87 ; ij1 = 88 161 158 ii0 = 160 ; ii1 = 161 162 tmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , 1:jpk ) = 0. e0159 tmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , 1:jpk ) = 0._wp 163 160 ELSE 164 161 IF(lwp) WRITE(numout,*) … … 182 179 ijl = nlcj - jprecj + 1 183 180 184 tmask_i( 1 :iif, : ) = 0. e0! first columns185 tmask_i(iil:jpi, : ) = 0. e0! last columns (including mpp extra columns)186 tmask_i( : , 1 :ijf) = 0. e0! first rows187 tmask_i( : ,ijl:jpj) = 0. e0! last rows (including mpp extra rows)181 tmask_i( 1 :iif, : ) = 0._wp ! first columns 182 tmask_i(iil:jpi, : ) = 0._wp ! last columns (including mpp extra columns) 183 tmask_i( : , 1 :ijf) = 0._wp ! first rows 184 tmask_i( : ,ijl:jpj) = 0._wp ! last rows (including mpp extra rows) 188 185 189 186 ! north fold mask 190 187 ! --------------- 191 tpol(1:jpiglo) = 1. e0192 fpol(1:jpiglo) = 1. e0188 tpol(1:jpiglo) = 1._wp 189 fpol(1:jpiglo) = 1._wp 193 190 IF( jperio == 3 .OR. jperio == 4 ) THEN ! T-point pivot 194 tpol(jpiglo/2+1:jpiglo) = 0. e0195 fpol( 1 :jpiglo) = 0. e0191 tpol(jpiglo/2+1:jpiglo) = 0._wp 192 fpol( 1 :jpiglo) = 0._wp 196 193 IF( mjg(nlej) == jpjglo ) THEN ! only half of the nlcj-1 row 197 194 DO ji = iif+1, iil-1 … … 201 198 ENDIF 202 199 IF( jperio == 5 .OR. jperio == 6 ) THEN ! F-point pivot 203 tpol( 1 :jpiglo) = 0. e0204 fpol(jpiglo/2+1:jpiglo) = 0. e0200 tpol( 1 :jpiglo) = 0._wp 201 fpol(jpiglo/2+1:jpiglo) = 0._wp 205 202 ENDIF 206 203 … … 219 216 END DO 220 217 END DO 221 CALL lbc_lnk( umask, 'U', 1. ) ! Lateral boundary conditions222 CALL lbc_lnk( vmask, 'V', 1. )223 CALL lbc_lnk( fmask, 'F', 1. )218 CALL lbc_lnk( umask, 'U', 1._wp ) ! Lateral boundary conditions 219 CALL lbc_lnk( vmask, 'V', 1._wp ) 220 CALL lbc_lnk( fmask, 'F', 1._wp ) 224 221 225 222 … … 231 228 ! ! cyclic east-west : bmask must be set to 0. on rows 1 and jpi 232 229 IF( nperio == 1 .OR. nperio == 4 .OR. nperio == 6 ) THEN 233 bmask( 1 ,:) = 0. e0234 bmask(jpi,:) = 0. e0230 bmask( 1 ,:) = 0._wp 231 bmask(jpi,:) = 0._wp 235 232 ENDIF 236 233 IF( nperio == 2 ) THEN ! south symmetric : bmask must be set to 0. on row 1 237 bmask(:, 1 ) = 0. e0234 bmask(:, 1 ) = 0._wp 238 235 ENDIF 239 236 ! ! north fold : … … 242 239 ii = ji + nimpp - 1 243 240 bmask(ji,jpj-1) = bmask(ji,jpj-1) * tpol(ii) 244 bmask(ji,jpj ) = 0. e0241 bmask(ji,jpj ) = 0._wp 245 242 END DO 246 243 ENDIF 247 244 IF( nperio == 5 .OR. nperio == 6 ) THEN ! F-pt pivot and T-pt elliptic eq. : bmask set to 0. on row jpj 248 bmask(:,jpj) = 0. e0245 bmask(:,jpj) = 0._wp 249 246 ENDIF 250 247 ! 251 248 IF( lk_mpp ) THEN ! mpp specificities 252 249 ! ! bmask is set to zero on the overlap region 253 IF( nbondi /= -1 .AND. nbondi /= 2 ) bmask( 1 :jpreci,:) = 0. e0254 IF( nbondi /= 1 .AND. nbondi /= 2 ) bmask(nlci:jpi ,:) = 0. e0255 IF( nbondj /= -1 .AND. nbondj /= 2 ) bmask(:, 1 :jprecj) = 0. e0256 IF( nbondj /= 1 .AND. nbondj /= 2 ) bmask(:,nlcj:jpj ) = 0. e0250 IF( nbondi /= -1 .AND. nbondi /= 2 ) bmask( 1 :jpreci,:) = 0._wp 251 IF( nbondi /= 1 .AND. nbondi /= 2 ) bmask(nlci:jpi ,:) = 0._wp 252 IF( nbondj /= -1 .AND. nbondj /= 2 ) bmask(:, 1 :jprecj) = 0._wp 253 IF( nbondj /= 1 .AND. nbondj /= 2 ) bmask(:,nlcj:jpj ) = 0._wp 257 254 ! 258 255 IF( npolj == 3 .OR. npolj == 4 ) THEN ! north fold : bmask must be set to 0. on rows jpj-1 and jpj … … 260 257 ii = ji + nimpp - 1 261 258 bmask(ji,nlcj-1) = bmask(ji,nlcj-1) * tpol(ii) 262 bmask(ji,nlcj ) = 0. e0259 bmask(ji,nlcj ) = 0._wp 263 260 END DO 264 261 ENDIF 265 262 IF( npolj == 5 .OR. npolj == 6 ) THEN ! F-pt pivot and T-pt elliptic eq. : bmask set to 0. on row jpj 266 263 DO ji = 1, nlci 267 bmask(ji,nlcj ) = 0. e0264 bmask(ji,nlcj ) = 0._wp 268 265 END DO 269 266 ENDIF … … 282 279 DO jj = 2, jpjm1 283 280 DO ji = fs_2, fs_jpim1 ! vector opt. 284 IF( fmask(ji,jj,jk) == 0. ) THEN285 fmask(ji,jj,jk) = rn_shlat * MIN( 1. , MAX( zwf(ji+1,jj), zwf(ji,jj+1), &286 & zwf(ji-1,jj), zwf(ji,jj-1) ) )281 IF( fmask(ji,jj,jk) == 0._wp ) THEN 282 fmask(ji,jj,jk) = rn_shlat * MIN( 1._wp , MAX( zwf(ji+1,jj), zwf(ji,jj+1), & 283 & zwf(ji-1,jj), zwf(ji,jj-1) ) ) 287 284 ENDIF 288 285 END DO 289 286 END DO 290 287 DO jj = 2, jpjm1 291 IF( fmask(1,jj,jk) == 0. ) THEN292 fmask(1 ,jj,jk) = rn_shlat * MIN( 1. , MAX( zwf(2,jj), zwf(1,jj+1), zwf(1,jj-1) ) )293 ENDIF 294 IF( fmask(jpi,jj,jk) == 0. ) THEN295 fmask(jpi,jj,jk) = rn_shlat * MIN( 1. , MAX( zwf(jpi,jj+1), zwf(jpim1,jj), zwf(jpi,jj-1) ) )288 IF( fmask(1,jj,jk) == 0._wp ) THEN 289 fmask(1 ,jj,jk) = rn_shlat * MIN( 1._wp , MAX( zwf(2,jj), zwf(1,jj+1), zwf(1,jj-1) ) ) 290 ENDIF 291 IF( fmask(jpi,jj,jk) == 0._wp ) THEN 292 fmask(jpi,jj,jk) = rn_shlat * MIN( 1._wp , MAX( zwf(jpi,jj+1), zwf(jpim1,jj), zwf(jpi,jj-1) ) ) 296 293 ENDIF 297 294 END DO 298 295 DO ji = 2, jpim1 299 IF( fmask(ji,1,jk) == 0. ) THEN300 fmask(ji, 1 ,jk) = rn_shlat * MIN( 1. , MAX( zwf(ji+1,1), zwf(ji,2), zwf(ji-1,1) ) )301 ENDIF 302 IF( fmask(ji,jpj,jk) == 0. ) THEN303 fmask(ji,jpj,jk) = rn_shlat * MIN( 1. , MAX( zwf(ji+1,jpj), zwf(ji-1,jpj), zwf(ji,jpjm1) ) )296 IF( fmask(ji,1,jk) == 0._wp ) THEN 297 fmask(ji, 1 ,jk) = rn_shlat * MIN( 1._wp , MAX( zwf(ji+1,1), zwf(ji,2), zwf(ji-1,1) ) ) 298 ENDIF 299 IF( fmask(ji,jpj,jk) == 0._wp ) THEN 300 fmask(ji,jpj,jk) = rn_shlat * MIN( 1._wp , MAX( zwf(ji+1,jpj), zwf(ji-1,jpj), zwf(ji,jpjm1) ) ) 304 301 ENDIF 305 302 END DO … … 308 305 IF( cp_cfg == "orca" .AND. jp_cfg == 2 ) THEN ! ORCA_R2 configuration 309 306 ! ! Increased lateral friction near of some straits 310 IF( n _cla == 0 ) THEN307 IF( nn_cla == 0 ) THEN 311 308 ! ! Gibraltar strait : partial slip (fmask=0.5) 312 309 ij0 = 101 ; ij1 = 101 313 ii0 = 139 ; ii1 = 140 ; fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , 1:jpk ) = 0.5 e0310 ii0 = 139 ; ii1 = 140 ; fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , 1:jpk ) = 0.5_wp 314 311 ij0 = 102 ; ij1 = 102 315 ii0 = 139 ; ii1 = 140 ; fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , 1:jpk ) = 0.5 e0312 ii0 = 139 ; ii1 = 140 ; fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , 1:jpk ) = 0.5_wp 316 313 ! 317 314 ! ! Bab el Mandeb : partial slip (fmask=1) 318 315 ij0 = 87 ; ij1 = 88 319 ii0 = 160 ; ii1 = 160 ; fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , 1:jpk ) = 1. e0316 ii0 = 160 ; ii1 = 160 ; fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , 1:jpk ) = 1._wp 320 317 ij0 = 88 ; ij1 = 88 321 ii0 = 159 ; ii1 = 159 ; fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , 1:jpk ) = 1. e0318 ii0 = 159 ; ii1 = 159 ; fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , 1:jpk ) = 1._wp 322 319 ! 323 320 ENDIF 324 325 321 ! ! Danish straits : strong slip (fmask > 2) 326 322 ! We keep this as an example but it is instable in this case 327 323 ! ij0 = 115 ; ij1 = 115 328 ! ii0 = 145 ; ii1 = 146 ; fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , 1:jpk ) = 4. 0e0324 ! ii0 = 145 ; ii1 = 146 ; fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , 1:jpk ) = 4._wp 329 325 ! ij0 = 116 ; ij1 = 116 330 ! ii0 = 145 ; ii1 = 146 ; fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , 1:jpk ) = 4. 0e0326 ! ii0 = 145 ; ii1 = 146 ; fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , 1:jpk ) = 4._wp 331 327 ! 332 328 ENDIF 333 329 ! 334 CALL lbc_lnk( fmask, 'F', 1. ) ! Lateral boundary conditions on fmask 335 336 337 ! Mbathy set to the number of w-level (minimum value 2) 338 ! ----------------------------------- 339 DO jj = 1, jpj 340 DO ji = 1, jpi 341 mbathy(ji,jj) = MAX( 1, mbathy(ji,jj) ) + 1 342 END DO 343 END DO 344 330 IF( cp_cfg == "orca" .AND. jp_cfg == 1 ) THEN ! ORCA R1 configuration 331 ! ! Increased lateral friction near of some straits 332 IF(lwp) WRITE(numout,*) 333 IF(lwp) WRITE(numout,*) ' orca_r1: increase friction near the following straits : ' 334 IF(lwp) WRITE(numout,*) ' Gibraltar ' 335 ii0 = 283 ; ii1 = 284 ! Gibraltar Strait 336 ij0 = 200 ; ij1 = 200 ; fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) = 2._wp 337 338 IF(lwp) WRITE(numout,*) ' Bhosporus ' 339 ii0 = 314 ; ii1 = 315 ! Bhosporus Strait 340 ij0 = 208 ; ij1 = 208 ; fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) = 2._wp 341 342 IF(lwp) WRITE(numout,*) ' Makassar (Top) ' 343 ii0 = 48 ; ii1 = 48 ! Makassar Strait (Top) 344 ij0 = 149 ; ij1 = 150 ; fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) = 3._wp 345 346 IF(lwp) WRITE(numout,*) ' Lombok ' 347 ii0 = 44 ; ii1 = 44 ! Lombok Strait 348 ij0 = 124 ; ij1 = 125 ; fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) = 2._wp 349 350 IF(lwp) WRITE(numout,*) ' Ombai ' 351 ii0 = 53 ; ii1 = 53 ! Ombai Strait 352 ij0 = 124 ; ij1 = 125 ; fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) = 2._wp 353 354 IF(lwp) WRITE(numout,*) ' Timor Passage ' 355 ii0 = 56 ; ii1 = 56 ! Timor Passage 356 ij0 = 124 ; ij1 = 125 ; fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) = 2._wp 357 358 IF(lwp) WRITE(numout,*) ' West Halmahera ' 359 ii0 = 58 ; ii1 = 58 ! West Halmahera Strait 360 ij0 = 141 ; ij1 = 142 ; fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) = 3._wp 361 362 IF(lwp) WRITE(numout,*) ' East Halmahera ' 363 ii0 = 55 ; ii1 = 55 ! East Halmahera Strait 364 ij0 = 141 ; ij1 = 142 ; fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1), 1:jpk ) = 3._wp 365 ! 366 ENDIF 367 ! 368 CALL lbc_lnk( fmask, 'F', 1._wp ) ! Lateral boundary conditions on fmask 369 370 345 371 IF( nprint == 1 .AND. lwp ) THEN ! Control print 346 372 imsk(:,:) = INT( tmask_i(:,:) ) … … 385 411 imsk(:,:) = INT( bmask(:,:) ) 386 412 CALL prihin( imsk(:,:), jpi, jpj, 1, jpi, 1, & 387 &1, jpj, 1, 1, numout )413 & 1, jpj, 1, 1, numout ) 388 414 ENDIF 389 415 ! … … 404 430 !! 405 431 !! ** Action : 406 !!407 432 !!---------------------------------------------------------------------- 408 433 INTEGER :: ji, jj, jk, jl ! dummy loop indices … … 448 473 zaa = tmask(ji ,jj,jk) + tmask(ji ,jj+1,jk) & 449 474 &+ tmask(ji+1,jj,jk) + tmask(ji+1,jj+1,jk) 450 IF( ABS(zaa-3. ) <= 0.1 ) fmask(ji,jj,jk) = 1.475 IF( ABS(zaa-3._wp) <= 0.1_wp ) fmask(ji,jj,jk) = 1._wp 451 476 END DO 452 477 END DO … … 461 486 DO ji = 2, jpim1 462 487 zaa = tmask(ji+1,jj,jk) + tmask(ji+1,jj+1,jk) 463 IF( ABS(zaa-2. ) <= 0.1 .AND. fmask(ji,jj,jk) == 0) THEN488 IF( ABS(zaa-2._wp) <= 0.1_wp .AND. fmask(ji,jj,jk) == 0._wp ) THEN 464 489 inw = inw + 1 465 490 nicoa(inw,1,jk) = ji … … 468 493 ENDIF 469 494 zaa = tmask(ji,jj,jk) + tmask(ji,jj+1,jk) 470 IF( ABS(zaa-2. ) <= 0.1 .AND. fmask(ji,jj,jk) == 0) THEN495 IF( ABS(zaa-2._wp) <= 0.1_wp .AND. fmask(ji,jj,jk) == 0._wp ) THEN 471 496 ine = ine + 1 472 497 nicoa(ine,2,jk) = ji … … 488 513 DO ji =2, jpim1 489 514 zaa = tmask(ji,jj+1,jk) + tmask(ji+1,jj+1,jk) 490 IF( ABS(zaa-2. ) <= 0.1 .AND. fmask(ji,jj,jk) == 0) THEN515 IF( ABS(zaa-2._wp) <= 0.1_wp .AND. fmask(ji,jj,jk) == 0._wp ) THEN 491 516 ins = ins + 1 492 517 nicoa(ins,3,jk) = ji … … 495 520 ENDIF 496 521 zaa = tmask(ji+1,jj,jk) + tmask(ji,jj,jk) 497 IF( ABS(zaa-2. ) <= 0.1 .AND. fmask(ji,jj,jk) == 0) THEN522 IF( ABS(zaa-2._wp) <= 0.1_wp .AND. fmask(ji,jj,jk) == 0._wp ) THEN 498 523 inn = inn + 1 499 524 nicoa(inn,4,jk) = ji … … 524 549 iind = 0 525 550 ijnd = 0 526 IF( nperio == 1 .OR. nperio == 4 .OR. nperio == 6 ) iind = 2527 IF( nperio == 3 .OR. nperio == 4 .OR. nperio == 5 .OR. nperio == 6 ) ijnd = 2551 IF( nperio == 1 .OR. nperio == 4 .OR. nperio == 6 ) iind = 2 552 IF( nperio == 3 .OR. nperio == 4 .OR. nperio == 5 .OR. nperio == 6 ) ijnd = 2 528 553 DO jk = 1, jpk 529 554 DO jl = 1, npcoa(1,jk) … … 551 576 ENDIF 552 577 END DO 553 DO jl =1,npcoa(4,jk)578 DO jl = 1, npcoa(4,jk) 554 579 IF( njcoa(jl,4,jk)-2 < 1) THEN 555 ierror=ierror +1556 icoord(ierror,1) =nicoa(jl,4,jk)557 icoord(ierror,2) =njcoa(jl,4,jk)558 icoord(ierror,3) =jk580 ierror=ierror + 1 581 icoord(ierror,1) = nicoa(jl,4,jk) 582 icoord(ierror,2) = njcoa(jl,4,jk) 583 icoord(ierror,3) = jk 559 584 ENDIF 560 585 END DO
Note: See TracChangeset
for help on using the changeset viewer.