- Timestamp:
- 2018-10-03T16:49:50+02:00 (5 years ago)
- Location:
- NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/tests
- Files:
-
- 16 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/tests/BENCH/MY_SRC/diawri.F90
r10166 r10170 181 181 END DO 182 182 END DO 183 CALL lbc_lnk( z2d, 'T', 1. )183 CALL lbc_lnk( 'diawri', z2d, 'T', 1. ) 184 184 CALL iom_put( "taubot", z2d ) 185 185 ENDIF … … 237 237 END DO 238 238 END DO 239 CALL lbc_lnk( z2d, 'T', 1. )239 CALL lbc_lnk( 'diawri', z2d, 'T', 1. ) 240 240 CALL iom_put( "sstgrad2", z2d ) ! square of module of sst gradient 241 241 z2d(:,:) = SQRT( z2d(:,:) ) … … 281 281 END DO 282 282 END DO 283 CALL lbc_lnk( z3d, 'T', 1. )283 CALL lbc_lnk( 'diawri', z3d, 'T', 1. ) 284 284 CALL iom_put( "eken", z3d ) ! kinetic energy 285 285 ENDIF … … 307 307 END DO 308 308 END DO 309 CALL lbc_lnk( z2d, 'U', -1. )309 CALL lbc_lnk( 'diawri', z2d, 'U', -1. ) 310 310 CALL iom_put( "u_heattr", 0.5*rcp * z2d ) ! heat transport in i-direction 311 311 ENDIF … … 320 320 END DO 321 321 END DO 322 CALL lbc_lnk( z2d, 'U', -1. )322 CALL lbc_lnk( 'diawri', z2d, 'U', -1. ) 323 323 CALL iom_put( "u_salttr", 0.5 * z2d ) ! heat transport in i-direction 324 324 ENDIF … … 342 342 END DO 343 343 END DO 344 CALL lbc_lnk( z2d, 'V', -1. )344 CALL lbc_lnk( 'diawri', z2d, 'V', -1. ) 345 345 CALL iom_put( "v_heattr", 0.5*rcp * z2d ) ! heat transport in j-direction 346 346 ENDIF … … 355 355 END DO 356 356 END DO 357 CALL lbc_lnk( z2d, 'V', -1. )357 CALL lbc_lnk( 'diawri', z2d, 'V', -1. ) 358 358 CALL iom_put( "v_salttr", 0.5 * z2d ) ! heat transport in j-direction 359 359 ENDIF … … 368 368 END DO 369 369 END DO 370 CALL lbc_lnk( z2d, 'T', -1. )370 CALL lbc_lnk( 'diawri', z2d, 'T', -1. ) 371 371 CALL iom_put( "tosmint", rau0 * z2d ) ! Vertical integral of temperature 372 372 ENDIF … … 380 380 END DO 381 381 END DO 382 CALL lbc_lnk( z2d, 'T', -1. )382 CALL lbc_lnk( 'diawri', z2d, 'T', -1. ) 383 383 CALL iom_put( "somint", rau0 * z2d ) ! Vertical integral of salinity 384 384 ENDIF -
NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/tests/BENCH/MY_SRC/usrdef_istate.F90
r10166 r10170 80 80 ENDDO 81 81 ! 82 CALL lbc_lnk( pssh, 'T', 1. ) ! apply boundary conditions83 CALL lbc_lnk( pts, 'T', 1. ) ! apply boundary conditions84 CALL lbc_lnk( pu, 'U', -1. ) ! apply boundary conditions85 CALL lbc_lnk( pv, 'V', -1. ) ! apply boundary conditions82 CALL lbc_lnk('usrdef_istate', pssh, 'T', 1. ) ! apply boundary conditions 83 CALL lbc_lnk( 'usrdef_istate', pts, 'T', 1. ) ! apply boundary conditions 84 CALL lbc_lnk( 'usrdef_istate', pu, 'U', -1. ) ! apply boundary conditions 85 CALL lbc_lnk( 'usrdef_istate', pv, 'V', -1. ) ! apply boundary conditions 86 86 87 87 END SUBROUTINE usr_def_istate -
NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/tests/BENCH/MY_SRC/usrdef_zgr.F90
r10166 r10170 194 194 z2d(:,:) = REAL( jpkm1 , wp ) ! flat bottom 195 195 ! 196 CALL lbc_lnk( z2d, 'T', 1. ) ! set surrounding land to zero (here jperio=0 ==>> closed)196 CALL lbc_lnk( 'usrdef_zgr', z2d, 'T', 1. ) ! set surrounding land to zero (here jperio=0 ==>> closed) 197 197 ! 198 198 k_bot(:,:) = INT( z2d(:,:) ) ! =jpkm1 over the ocean point, =0 elsewhere -
NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/tests/CANAL/MY_SRC/diawri.F90
r10074 r10170 185 185 END DO 186 186 END DO 187 CALL lbc_lnk( z2d, 'T', 1. )187 CALL lbc_lnk( 'diawri', z2d, 'T', 1. ) 188 188 CALL iom_put( "taubot", z2d ) 189 189 ENDIF … … 244 244 END DO 245 245 END DO 246 CALL lbc_lnk( z3d, 'T', 1. )246 CALL lbc_lnk( 'diawri', z3d, 'T', 1. ) 247 247 CALL iom_put( "salgrad2", z3d ) ! square of module of sal gradient 248 248 z3d(:,:,:) = SQRT( z3d(:,:,:) ) … … 260 260 END DO 261 261 END DO 262 CALL lbc_lnk( z2d, 'T', 1. )262 CALL lbc_lnk( 'diawri', z2d, 'T', 1. ) 263 263 CALL iom_put( "sstgrad2", z2d ) ! square of module of sst gradient 264 264 z2d(:,:) = SQRT( z2d(:,:) ) … … 314 314 END DO 315 315 END DO 316 CALL lbc_lnk( z3d, 'T', 1. )316 CALL lbc_lnk( 'diawri', z3d, 'T', 1. ) 317 317 CALL iom_put( "eken", z3d ) ! kinetic energy 318 318 ENDIF … … 335 335 END DO 336 336 337 CALL lbc_lnk( z3d, 'T', 1. )337 CALL lbc_lnk( 'diawri', z3d, 'T', 1. ) 338 338 CALL iom_put( "ke", z3d ) ! kinetic energy 339 339 … … 363 363 END DO 364 364 END DO 365 CALL lbc_lnk( z3d, 'F', 1. )365 CALL lbc_lnk( 'diawri', z3d, 'F', 1. ) 366 366 CALL iom_put( "relvor", z3d ) ! relative vorticity 367 367 … … 387 387 END DO 388 388 END DO 389 CALL lbc_lnk( z3d, 'F', 1. )389 CALL lbc_lnk( 'diawri', z3d, 'F', 1. ) 390 390 CALL iom_put( "potvor", z3d ) ! potential vorticity 391 391 … … 413 413 END DO 414 414 END DO 415 CALL lbc_lnk( z2d, 'U', -1. )415 CALL lbc_lnk( 'diawri', z2d, 'U', -1. ) 416 416 CALL iom_put( "u_heattr", 0.5*rcp * z2d ) ! heat transport in i-direction 417 417 ENDIF … … 426 426 END DO 427 427 END DO 428 CALL lbc_lnk( z2d, 'U', -1. )428 CALL lbc_lnk( 'diawri', z2d, 'U', -1. ) 429 429 CALL iom_put( "u_salttr", 0.5 * z2d ) ! heat transport in i-direction 430 430 ENDIF … … 448 448 END DO 449 449 END DO 450 CALL lbc_lnk( z2d, 'V', -1. )450 CALL lbc_lnk( 'diawri', z2d, 'V', -1. ) 451 451 CALL iom_put( "v_heattr", 0.5*rcp * z2d ) ! heat transport in j-direction 452 452 ENDIF … … 461 461 END DO 462 462 END DO 463 CALL lbc_lnk( z2d, 'V', -1. )463 CALL lbc_lnk( 'diawri', z2d, 'V', -1. ) 464 464 CALL iom_put( "v_salttr", 0.5 * z2d ) ! heat transport in j-direction 465 465 ENDIF … … 474 474 END DO 475 475 END DO 476 CALL lbc_lnk( z2d, 'T', -1. )476 CALL lbc_lnk( 'diawri', z2d, 'T', -1. ) 477 477 CALL iom_put( "tosmint", rau0 * z2d ) ! Vertical integral of temperature 478 478 ENDIF … … 486 486 END DO 487 487 END DO 488 CALL lbc_lnk( z2d, 'T', -1. )488 CALL lbc_lnk( 'diawri', z2d, 'T', -1. ) 489 489 CALL iom_put( "somint", rau0 * z2d ) ! Vertical integral of salinity 490 490 ENDIF -
NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/tests/CANAL/MY_SRC/domvvl.F90
r10074 r10170 385 385 ! ! d - thickness diffusion transport: boundary conditions 386 386 ! (stored for tracer advction and continuity equation) 387 CALL lbc_lnk_multi( un_td , 'U' , -1._wp, vn_td , 'V' , -1._wp)387 CALL lbc_lnk_multi( 'domvvl', un_td , 'U' , -1._wp, vn_td , 'V' , -1._wp) 388 388 389 389 ! 4 - Time stepping of baroclinic scale factors … … 396 396 z2dt = 2.0_wp * rdt 397 397 ENDIF 398 CALL lbc_lnk( tilde_e3t_a(:,:,:), 'T', 1._wp )398 CALL lbc_lnk( 'domvvl', tilde_e3t_a(:,:,:), 'T', 1._wp ) 399 399 tilde_e3t_a(:,:,:) = tilde_e3t_b(:,:,:) + z2dt * tmask(:,:,:) * tilde_e3t_a(:,:,:) 400 400 … … 688 688 END DO 689 689 END DO 690 CALL lbc_lnk( pe3_out(:,:,:), 'U', 1._wp )690 CALL lbc_lnk( 'domvvl', pe3_out(:,:,:), 'U', 1._wp ) 691 691 pe3_out(:,:,:) = pe3_out(:,:,:) + e3u_0(:,:,:) 692 692 ! … … 701 701 END DO 702 702 END DO 703 CALL lbc_lnk( pe3_out(:,:,:), 'V', 1._wp )703 CALL lbc_lnk( 'domvvl', pe3_out(:,:,:), 'V', 1._wp ) 704 704 pe3_out(:,:,:) = pe3_out(:,:,:) + e3v_0(:,:,:) 705 705 ! … … 715 715 END DO 716 716 END DO 717 CALL lbc_lnk( pe3_out(:,:,:), 'F', 1._wp )717 CALL lbc_lnk( 'domvvl', pe3_out(:,:,:), 'F', 1._wp ) 718 718 pe3_out(:,:,:) = pe3_out(:,:,:) + e3f_0(:,:,:) 719 719 ! -
NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/tests/CANAL/MY_SRC/trazdf.F90
r10074 r10170 92 92 END DO 93 93 !!gm this should be moved in trdtra.F90 and done on all trends 94 CALL lbc_lnk_multi( ztrdt, 'T', 1. , ztrds, 'T', 1. )94 CALL lbc_lnk_multi( 'trazdf', ztrdt, 'T', 1. , ztrds, 'T', 1. ) 95 95 !!gm 96 96 CALL trd_tra( kt, 'TRA', jp_tem, jptra_zdf, ztrdt ) -
NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/tests/CANAL/MY_SRC/usrdef_istate.F90
r10074 r10170 169 169 END DO 170 170 END DO 171 CALL lbc_lnk( pssh, 'T', 1. )171 CALL lbc_lnk( 'usrdef_istate', pssh, 'T', 1. ) 172 172 END DO 173 173 … … 293 293 pssh(:,:) = pssh(:,:) + ( 0.1 * zrandom(:,:) - 0.05 ) 294 294 END IF 295 CALL lbc_lnk( pssh, 'T', 1. )296 CALL lbc_lnk( pts, 'T', 1. )297 CALL lbc_lnk( pu, 'U', -1. )298 CALL lbc_lnk( pv, 'V', -1. )295 CALL lbc_lnk( 'usrdef_istate', pssh, 'T', 1. ) 296 CALL lbc_lnk( 'usrdef_istate', pts, 'T', 1. ) 297 CALL lbc_lnk( 'usrdef_istate', pu, 'U', -1. ) 298 CALL lbc_lnk( 'usrdef_istate', pv, 'V', -1. ) 299 299 300 300 END SUBROUTINE usr_def_istate -
NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/tests/CANAL/MY_SRC/usrdef_zgr.F90
r10074 r10170 204 204 END SELECT 205 205 ! 206 CALL lbc_lnk( z2d, 'T', 1. ) ! set surrounding land to zero (here jperio=0 ==>> closed)206 CALL lbc_lnk( 'usrdef_zgr', z2d, 'T', 1. ) ! set surrounding land to zero (here jperio=0 ==>> closed) 207 207 ! 208 208 k_bot(:,:) = INT( z2d(:,:) ) ! =jpkm1 over the ocean point, =0 elsewhere -
NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/tests/ISOMIP/MY_SRC/usrdef_zgr.F90
r10074 r10170 89 89 ! the ocean basin surrounded by land (1 grid-point) set through lbc_lnk call as jperio=0 90 90 z2d(:,:) = 1._wp ! surface ocean is the 1st level 91 CALL lbc_lnk( z2d, 'T', 1. ) ! closed basin since jperio = 0 (see userdef_nam.F90)91 CALL lbc_lnk( 'usrdef_zgr', z2d, 'T', 1. ) ! closed basin since jperio = 0 (see userdef_nam.F90) 92 92 zmsk(:,:) = NINT( z2d(:,:) ) 93 93 ! … … 177 177 END DO 178 178 END DO 179 CALL lbc_lnk( pe3v , 'V', 1._wp ) ; CALL lbc_lnk(pe3vw, 'V', 1._wp )180 CALL lbc_lnk( pe3f , 'F', 1._wp )179 CALL lbc_lnk( 'usrdef_zgr', pe3v , 'V', 1._wp ) ; CALL lbc_lnk( 'usrdef_zgr', pe3vw, 'V', 1._wp ) 180 CALL lbc_lnk( 'usrdef_zgr', pe3f , 'F', 1._wp ) 181 181 DO jk = 1, jpk 182 182 ! set to z-scale factor if zero (i.e. along closed boundaries) because of lbclnk -
NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/tests/LOCK_EXCHANGE/MY_SRC/usrdef_zgr.F90
r10074 r10170 88 88 ! the ocean basin surrounded by land (1 grid-point) set through lbc_lnk call as jperio=0 89 89 z2d(:,:) = 1._wp ! surface ocean is the 1st level 90 CALL lbc_lnk( z2d, 'T', 1. ) ! closed basin since jperio = 0 (see userdef_nam.F90)90 CALL lbc_lnk( 'usrdef_zgr', z2d, 'T', 1. ) ! closed basin since jperio = 0 (see userdef_nam.F90) 91 91 k_top(:,:) = NINT( z2d(:,:) ) 92 92 ! -
NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/tests/OVERFLOW/MY_SRC/usrdef_zgr.F90
r10074 r10170 93 93 zhu(ji,:) = 0.5_wp * ( zht(ji,:) + zht(ji+1,:) ) 94 94 END DO 95 CALL lbc_lnk( zhu, 'U', 1. ) ! boundary condition: this mask the surrouding grid-points95 CALL lbc_lnk( 'usrdef_zgr', zhu, 'U', 1. ) ! boundary condition: this mask the surrouding grid-points 96 96 ! ! ==>>> set by hand non-zero value on first/last columns & rows 97 97 DO ji = mi0(1), mi1(1) ! first row of global domain only … … 112 112 ! the ocean basin surrounded by land (1 grid-point) set through lbc_lnk call as jperio=0 113 113 z2d(:,:) = 1._wp ! surface ocean is the 1st level 114 CALL lbc_lnk( z2d, 'T', 1. ) ! closed basin since jperio = 0 (see userdef_nam.F90)114 CALL lbc_lnk( 'usrdef_zgr', z2d, 'T', 1. ) ! closed basin since jperio = 0 (see userdef_nam.F90) 115 115 k_top(:,:) = NINT( z2d(:,:) ) 116 116 ! -
NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/tests/VORTEX/MY_SRC/domvvl.F90
r10074 r10170 408 408 ! ! d - thickness diffusion transport: boundary conditions 409 409 ! (stored for tracer advction and continuity equation) 410 CALL lbc_lnk_multi( un_td , 'U' , -1._wp, vn_td , 'V' , -1._wp)410 CALL lbc_lnk_multi( 'domvvl', un_td , 'U' , -1._wp, vn_td , 'V' , -1._wp) 411 411 412 412 ! 4 - Time stepping of baroclinic scale factors … … 419 419 z2dt = 2.0_wp * rdt 420 420 ENDIF 421 CALL lbc_lnk( tilde_e3t_a(:,:,:), 'T', 1._wp )421 CALL lbc_lnk( 'domvvl', tilde_e3t_a(:,:,:), 'T', 1._wp ) 422 422 tilde_e3t_a(:,:,:) = tilde_e3t_b(:,:,:) + z2dt * tmask(:,:,:) * tilde_e3t_a(:,:,:) 423 423 … … 711 711 END DO 712 712 END DO 713 CALL lbc_lnk( pe3_out(:,:,:), 'U', 1._wp )713 CALL lbc_lnk( 'domvvl', pe3_out(:,:,:), 'U', 1._wp ) 714 714 pe3_out(:,:,:) = pe3_out(:,:,:) + e3u_0(:,:,:) 715 715 ! … … 724 724 END DO 725 725 END DO 726 CALL lbc_lnk( pe3_out(:,:,:), 'V', 1._wp )726 CALL lbc_lnk( 'domvvl', pe3_out(:,:,:), 'V', 1._wp ) 727 727 pe3_out(:,:,:) = pe3_out(:,:,:) + e3v_0(:,:,:) 728 728 ! … … 738 738 END DO 739 739 END DO 740 CALL lbc_lnk( pe3_out(:,:,:), 'F', 1._wp )740 CALL lbc_lnk( 'domvvl', pe3_out(:,:,:), 'F', 1._wp ) 741 741 pe3_out(:,:,:) = pe3_out(:,:,:) + e3f_0(:,:,:) 742 742 ! -
NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/tests/VORTEX/MY_SRC/usrdef_istate.F90
r10074 r10170 136 136 END DO 137 137 138 CALL lbc_lnk( pu, 'U', -1. )139 CALL lbc_lnk( pv, 'V', -1. )138 CALL lbc_lnk( 'usrdef_istate', pu, 'U', -1. ) 139 CALL lbc_lnk( 'usrdef_istate', pv, 'V', -1. ) 140 140 ! 141 141 END SUBROUTINE usr_def_istate -
NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/tests/VORTEX/MY_SRC/usrdef_zgr.F90
r10074 r10170 192 192 z2d(:,:) = REAL( jpkm1 , wp ) ! flat bottom 193 193 ! 194 CALL lbc_lnk( z2d, 'T', 1. ) ! set surrounding land to zero (here jperio=0 ==>> closed)194 CALL lbc_lnk( 'usrdef_zgr', z2d, 'T', 1. ) ! set surrounding land to zero (here jperio=0 ==>> closed) 195 195 ! 196 196 k_bot(:,:) = INT( z2d(:,:) ) ! =jpkm1 over the ocean point, =0 elsewhere -
NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/tests/WAD/MY_SRC/bdyini.F90
r10074 r10170 1133 1133 END DO 1134 1134 END DO 1135 CALL lbc_lnk_multi( bdyumask, 'U', 1. , bdyvmask, 'V', 1. ) ! Lateral boundary cond.1135 CALL lbc_lnk_multi( 'bdyini', bdyumask, 'U', 1. , bdyvmask, 'V', 1. ) ! Lateral boundary cond. 1136 1136 1137 1137 ! bdy masks are now set to zero on boundary points: … … 1169 1169 1170 1170 ! Lateral boundary conditions 1171 CALL lbc_lnk( zfmask, 'F', 1. )1172 CALL lbc_lnk_multi( bdyumask, 'U', 1. , bdyvmask, 'V', 1., bdytmask, 'T', 1. )1171 CALL lbc_lnk( 'bdyini', zfmask, 'F', 1. ) 1172 CALL lbc_lnk_multi( 'bdyini', bdyumask, 'U', 1. , bdyvmask, 'V', 1., bdytmask, 'T', 1. ) 1173 1173 DO ib_bdy = 1, nb_bdy ! Indices and directions of rim velocity components 1174 1174 -
NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/tests/WAD/MY_SRC/usrdef_zgr.F90
r10074 r10170 234 234 zhu(ji,:) = 0.5_wp * ( zht(ji,:) + zht(ji+1,:) ) 235 235 END DO 236 CALL lbc_lnk( zhu, 'U', 1. ) ! boundary condition: this mask the surrounding grid-points236 CALL lbc_lnk( 'usrdef_zgr', zhu, 'U', 1. ) ! boundary condition: this mask the surrounding grid-points 237 237 ! ! ==>>> set by hand non-zero value on first/last columns & rows 238 238 DO ji = mi0(1), mi1(1) ! first row of global domain only … … 247 247 zhv(:,jj) = 0.5_wp * ( zht(:,jj) + zht(:,jj+1) ) 248 248 END DO 249 CALL lbc_lnk( zhv, 'V', 1. ) ! boundary condition: this mask the surrounding grid-points249 CALL lbc_lnk( 'usrdef_zgr', zhv, 'V', 1. ) ! boundary condition: this mask the surrounding grid-points 250 250 DO jj = mj0(1), mj1(1) ! first row of global domain only 251 251 zhv(:,jj) = zht(:,jj) … … 272 272 273 273 274 CALL lbc_lnk( z2d, 'T', 1. ) ! closed basin since jperio = 0 (see userdef_nam.F90)274 CALL lbc_lnk( 'usrdef_zgr', z2d, 'T', 1. ) ! closed basin since jperio = 0 (see userdef_nam.F90) 275 275 k_top(:,:) = NINT( z2d(:,:) ) 276 276 ! … … 311 311 END DO 312 312 END DO 313 CALL lbc_lnk( pdept, 'T', 1. )314 CALL lbc_lnk( pdepw, 'T', 1. )315 CALL lbc_lnk( pe3t , 'T', 1. )316 CALL lbc_lnk( pe3w , 'T', 1. )317 CALL lbc_lnk( pe3u , 'U', 1. )318 CALL lbc_lnk( pe3uw, 'U', 1. )319 CALL lbc_lnk( pe3f , 'F', 1. )320 CALL lbc_lnk( pe3v , 'V', 1. )321 CALL lbc_lnk( pe3vw, 'V', 1. )313 CALL lbc_lnk( 'usrdef_zgr', pdept, 'T', 1. ) 314 CALL lbc_lnk( 'usrdef_zgr', pdepw, 'T', 1. ) 315 CALL lbc_lnk( 'usrdef_zgr', pe3t , 'T', 1. ) 316 CALL lbc_lnk( 'usrdef_zgr', pe3w , 'T', 1. ) 317 CALL lbc_lnk( 'usrdef_zgr', pe3u , 'U', 1. ) 318 CALL lbc_lnk( 'usrdef_zgr', pe3uw, 'U', 1. ) 319 CALL lbc_lnk( 'usrdef_zgr', pe3f , 'F', 1. ) 320 CALL lbc_lnk( 'usrdef_zgr', pe3v , 'V', 1. ) 321 CALL lbc_lnk( 'usrdef_zgr', pe3vw, 'V', 1. ) 322 322 WHERE( pe3t (:,:,:) == 0._wp ) pe3t (:,:,:) = 1._wp 323 323 WHERE( pe3u (:,:,:) == 0._wp ) pe3u (:,:,:) = 1._wp
Note: See TracChangeset
for help on using the changeset viewer.