New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 10727 for utils/tools_AGRIF_CMEMS_2020/DOMAINcfg/src/domhgr.F90 – NEMO

Ignore:
Timestamp:
2019-02-27T17:02:02+01:00 (5 years ago)
Author:
rblod
Message:

new nesting tools (attempt) and brutal cleaning of DOMAINcfg, see ticket #2129

File:
1 moved

Legend:

Unmodified
Added
Removed
  • utils/tools_AGRIF_CMEMS_2020/DOMAINcfg/src/domhgr.F90

    r10725 r10727  
    112112      !!---------------------------------------------------------------------- 
    113113      ! 
    114       IF( nn_timing == 1 )  CALL timing_start('dom_hgr') 
     114  !    IF( nn_timing == 1 )  CALL timing_start('dom_hgr') 
    115115      ! 
    116116      IF(lwp) THEN 
     
    131131      CASE ( 0 )                     !==  read in coordinate.nc file  ==! 
    132132         ! 
     133#if defined key_agrif 
     134         IF (agrif_root()) THEN 
     135#endif 
    133136         IF(lwp) WRITE(numout,*) 
    134137         IF(lwp) WRITE(numout,*) '          curvilinear coordinate on the sphere read in "coordinate" file' 
     
    143146            e1e2v (:,:) = e1v(:,:) * e2v(:,:)  
    144147         ENDIF 
     148#if defined key_agrif 
     149         ELSE 
     150           CALL Agrif_InitValues_cont() 
     151         ENDIF 
     152#endif 
    145153         ! 
    146154      CASE ( 1 )                     !==  geographical mesh on the sphere with regular (in degree) grid-spacing  ==! 
     
    272280         ze1 = 106000. / REAL( jp_cfg , wp )             
    273281         ! benchmark: forced the resolution to be about 100 km 
    274          IF( nbench /= 0 )   ze1 = 106000._wp      
     282       !  IF( nbench /= 0 )   ze1 = 106000._wp      
    275283         zsin_alpha = - SQRT( 2._wp ) * 0.5_wp 
    276284         zcos_alpha =   SQRT( 2._wp ) * 0.5_wp 
    277285         ze1deg = ze1 / (ra * rad) 
    278          IF( nbench /= 0 )   ze1deg = ze1deg / REAL( jp_cfg , wp )   ! benchmark: keep the lat/+lon 
     286       !  IF( nbench /= 0 )   ze1deg = ze1deg / REAL( jp_cfg , wp )   ! benchmark: keep the lat/+lon 
    279287         !                                                           ! at the right jp_cfg resolution 
    280288         glam0 = zlam1 + zcos_alpha * ze1deg * REAL( jpjglo-2 , wp ) 
     
    395403            zminff=ff_f(nldi,nldj) 
    396404            zmaxff=ff_f(nldi,nlej) 
    397             CALL mpp_min( zminff )   ! min over the global domain 
    398             CALL mpp_max( zmaxff )   ! max over the global domain 
     405            CALL mpp_min( 'toto',zminff )   ! min over the global domain 
     406            CALL mpp_max( 'toto',zmaxff )   ! max over the global domain 
    399407            IF(lwp) WRITE(numout,*) '          Coriolis parameter varies globally from ', zminff,' to ', zmaxff 
    400408         END IF 
     
    418426            zminff=ff_f(nldi,nldj) 
    419427            zmaxff=ff_f(nldi,nlej) 
    420             CALL mpp_min( zminff )   ! min over the global domain 
    421             CALL mpp_max( zmaxff )   ! max over the global domain 
     428            CALL mpp_min('toto', zminff )   ! min over the global domain 
     429            CALL mpp_max( 'toto',zmaxff )   ! max over the global domain 
    422430            IF(lwp) WRITE(numout,*) '          Coriolis parameter varies globally from ', zminff,' to ', zmaxff 
    423431         END IF 
     
    430438      ! The equator line must be the latitude coordinate axe 
    431439 
    432       IF( nperio == 2 ) THEN 
    433          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       ENDIF 
    436       ! 
    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') 
    438446      ! 
    439447   END SUBROUTINE dom_hgr 
Note: See TracChangeset for help on using the changeset viewer.