Ignore:
Timestamp:
2015-07-31T10:55:56+02:00 (5 years ago)
Author:
timgraham
Message:

Merge of AGRIF branch (branches/2014/dev_r4765_CNRS_agrif) onto the trunk

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/OPA_SRC/nemogcm.F90

    r5536 r5656  
    161161          ENDIF 
    162162 
     163#if defined key_agrif 
     164          CALL Agrif_Regrid() 
     165#endif 
     166 
    163167         DO WHILE ( istp <= nitend .AND. nstop == 0 ) 
    164168#if defined key_agrif 
    165             CALL Agrif_Step( stp )           ! AGRIF: time stepping 
     169            CALL stp                         ! AGRIF: time stepping 
    166170#else 
    167171            CALL stp( istp )                 ! standard time stepping 
     
    187191      ! 
    188192#if defined key_agrif 
    189       IF ( Agrif_Level() < Agrif_MaxLevel() ) THEN  
     193      IF( .NOT. Agrif_Root() ) THEN 
    190194         CALL Agrif_ParentGrid_To_ChildGrid() 
    191195         IF( lk_diaobs ) CALL dia_obs_wri 
     
    336340         jpj = ( jpjglo-2*jprecj + (jpnj-1) ) / jpnj + 2*jprecj   ! second dim. 
    337341#endif 
    338       ENDIF 
     342      ENDIF          
    339343         jpk = jpkdta                                             ! third dim 
     344#if defined key_agrif 
     345         ! simple trick to use same vertical grid as parent 
     346         ! but different number of levels:  
     347         ! Save maximum number of levels in jpkdta, then define all vertical grids 
     348         ! with this number. 
     349         ! Suppress once vertical online interpolation is ok 
     350         IF(.NOT.Agrif_Root()) jpkdta = Agrif_Parent(jpkdta) 
     351#endif 
    340352         jpim1 = jpi-1                                            ! inner domain indices 
    341353         jpjm1 = jpj-1                                            !   "           " 
     
    712724      INTEGER :: ifac, jl, inu 
    713725      INTEGER, PARAMETER :: ntest = 14 
    714       INTEGER :: ilfax(ntest) 
    715       ! 
     726      INTEGER, DIMENSION(ntest) :: ilfax 
     727      !!---------------------------------------------------------------------- 
    716728      ! lfax contains the set of allowed factors. 
    717       data (ilfax(jl),jl=1,ntest) / 16384, 8192, 4096, 2048, 1024, 512, 256,  & 
    718          &                            128,   64,   32,   16,    8,   4,   2  / 
    719       !!---------------------------------------------------------------------- 
     729      ilfax(:) = (/(2**jl,jl=ntest,1,-1)/) 
    720730 
    721731      ! Clear the error flag and initialise output vars 
Note: See TracChangeset for help on using the changeset viewer.