- Timestamp:
- 2016-04-20T11:33:10+02:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/nemogcm.F90
r6486 r6487 161 161 ENDIF 162 162 163 #if defined key_agrif 164 CALL Agrif_Regrid() 165 #endif 166 163 167 DO WHILE ( istp <= nitend .AND. nstop == 0 ) 164 168 #if defined key_agrif 165 CALL Agrif_Step( stp )! AGRIF: time stepping169 CALL stp ! AGRIF: time stepping 166 170 #else 167 171 CALL stp( istp ) ! standard time stepping … … 187 191 ! 188 192 #if defined key_agrif 189 CALL Agrif_ParentGrid_To_ChildGrid() 190 IF( lk_diaobs ) CALL dia_obs_wri 191 IF( nn_timing == 1 ) CALL timing_finalize 192 CALL Agrif_ChildGrid_To_ParentGrid() 193 IF( .NOT. Agrif_Root() ) THEN 194 CALL Agrif_ParentGrid_To_ChildGrid() 195 IF( lk_diaobs ) CALL dia_obs_wri 196 IF( nn_timing == 1 ) CALL timing_finalize 197 CALL Agrif_ChildGrid_To_ParentGrid() 198 ENDIF 193 199 #endif 194 200 IF( nn_timing == 1 ) CALL timing_finalize … … 334 340 jpj = ( jpjglo-2*jprecj + (jpnj-1) ) / jpnj + 2*jprecj ! second dim. 335 341 #endif 336 ENDIF 342 ENDIF 337 343 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 338 352 jpim1 = jpi-1 ! inner domain indices 339 353 jpjm1 = jpj-1 ! " " … … 710 724 INTEGER :: ifac, jl, inu 711 725 INTEGER, PARAMETER :: ntest = 14 712 INTEGER :: ilfax(ntest) 713 ! 714 ! lfax contains the set of allowed factors. 715 data (ilfax(jl),jl=1,ntest) / 16384, 8192, 4096, 2048, 1024, 512, 256, & 716 & 128, 64, 32, 16, 8, 4, 2 / 717 !!---------------------------------------------------------------------- 726 INTEGER, DIMENSION(ntest) :: ilfax 727 ! 728 ! ilfax contains the set of allowed factors. 729 ilfax(:) = (/(2**jl,jl=ntest,1,-1)/) 730 !!---------------------------------------------------------------------- 731 ! ilfax contains the set of allowed factors. 732 ilfax(:) = (/(2**jl,jl=ntest,1,-1)/) 718 733 719 734 ! Clear the error flag and initialise output vars
Note: See TracChangeset
for help on using the changeset viewer.