Changeset 10727 for utils/tools_AGRIF_CMEMS_2020/DOMAINcfg/src/domhgr.F90
- Timestamp:
- 2019-02-27T17:02:02+01:00 (5 years ago)
- File:
-
- 1 moved
Legend:
- Unmodified
- Added
- Removed
-
utils/tools_AGRIF_CMEMS_2020/DOMAINcfg/src/domhgr.F90
r10725 r10727 112 112 !!---------------------------------------------------------------------- 113 113 ! 114 IF( nn_timing == 1 ) CALL timing_start('dom_hgr')114 ! IF( nn_timing == 1 ) CALL timing_start('dom_hgr') 115 115 ! 116 116 IF(lwp) THEN … … 131 131 CASE ( 0 ) !== read in coordinate.nc file ==! 132 132 ! 133 #if defined key_agrif 134 IF (agrif_root()) THEN 135 #endif 133 136 IF(lwp) WRITE(numout,*) 134 137 IF(lwp) WRITE(numout,*) ' curvilinear coordinate on the sphere read in "coordinate" file' … … 143 146 e1e2v (:,:) = e1v(:,:) * e2v(:,:) 144 147 ENDIF 148 #if defined key_agrif 149 ELSE 150 CALL Agrif_InitValues_cont() 151 ENDIF 152 #endif 145 153 ! 146 154 CASE ( 1 ) !== geographical mesh on the sphere with regular (in degree) grid-spacing ==! … … 272 280 ze1 = 106000. / REAL( jp_cfg , wp ) 273 281 ! benchmark: forced the resolution to be about 100 km 274 IF( nbench /= 0 ) ze1 = 106000._wp282 ! IF( nbench /= 0 ) ze1 = 106000._wp 275 283 zsin_alpha = - SQRT( 2._wp ) * 0.5_wp 276 284 zcos_alpha = SQRT( 2._wp ) * 0.5_wp 277 285 ze1deg = ze1 / (ra * rad) 278 IF( nbench /= 0 ) ze1deg = ze1deg / REAL( jp_cfg , wp ) ! benchmark: keep the lat/+lon286 ! IF( nbench /= 0 ) ze1deg = ze1deg / REAL( jp_cfg , wp ) ! benchmark: keep the lat/+lon 279 287 ! ! at the right jp_cfg resolution 280 288 glam0 = zlam1 + zcos_alpha * ze1deg * REAL( jpjglo-2 , wp ) … … 395 403 zminff=ff_f(nldi,nldj) 396 404 zmaxff=ff_f(nldi,nlej) 397 CALL mpp_min( zminff ) ! min over the global domain398 CALL mpp_max( zmaxff ) ! max over the global domain405 CALL mpp_min( 'toto',zminff ) ! min over the global domain 406 CALL mpp_max( 'toto',zmaxff ) ! max over the global domain 399 407 IF(lwp) WRITE(numout,*) ' Coriolis parameter varies globally from ', zminff,' to ', zmaxff 400 408 END IF … … 418 426 zminff=ff_f(nldi,nldj) 419 427 zmaxff=ff_f(nldi,nlej) 420 CALL mpp_min( zminff ) ! min over the global domain421 CALL mpp_max( zmaxff ) ! max over the global domain428 CALL mpp_min('toto', zminff ) ! min over the global domain 429 CALL mpp_max( 'toto',zmaxff ) ! max over the global domain 422 430 IF(lwp) WRITE(numout,*) ' Coriolis parameter varies globally from ', zminff,' to ', zmaxff 423 431 END IF … … 430 438 ! The equator line must be the latitude coordinate axe 431 439 432 IF( nperio == 2 ) THEN433 znorme = SQRT( SUM( gphiu(:,2) * gphiu(:,2) ) ) / REAL( jpi )434 IF( znorme > 1.e-13 ) CALL ctl_stop( ' ===>>>> : symmetrical condition: rerun with good equator line' )435 ENDIF436 ! 437 IF( nn_timing == 1 ) CALL timing_stop('dom_hgr')440 ! IF( nperio == 2 ) THEN 441 ! znorme = SQRT( SUM( gphiu(:,2) * gphiu(:,2) ) ) / REAL( jpi ) 442 ! IF( znorme > 1.e-13 ) CALL ctl_stop( ' ===>>>> : symmetrical condition: rerun with good equator line' ) 443 ! ENDIF 444 ! 445 ! IF( nn_timing == 1 ) CALL timing_stop('dom_hgr') 438 446 ! 439 447 END SUBROUTINE dom_hgr
Note: See TracChangeset
for help on using the changeset viewer.