- Timestamp:
- 2015-02-17T10:06:39+01:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_r4650_UKMO3_masked_damping/NEMOGCM/NEMO/OPA_SRC/nemogcm.F90
r4645 r5086 42 42 !!---------------------------------------------------------------------- 43 43 USE step_oce ! module used in the ocean time stepping module 44 USE sbc_oce ! surface boundary condition: ocean45 44 USE cla ! cross land advection (tra_cla routine) 46 45 USE domcfg ! domain configuration (dom_cfg routine) … … 51 50 #endif 52 51 USE tideini ! tidal components initialization (tide_ini routine) 53 USE bdyini ! open boundary cond. initialization(bdy_init routine)54 USE bdydta ! open boundary cond. initialization(bdy_dta_init routine)55 USE bdytides ! open boundary cond. initialization(bdytide_init routine)52 USE bdyini ! open boundary cond. setting (bdy_init routine) 53 USE bdydta ! open boundary cond. setting (bdy_dta_init routine) 54 USE bdytides ! open boundary cond. setting (bdytide_init routine) 56 55 USE istate ! initial state setting (istate_init routine) 57 56 USE ldfdyn ! lateral viscosity setting (ldfdyn_init routine) … … 59 58 USE zdfini ! vertical physics setting (zdf_init routine) 60 59 USE phycst ! physical constant (par_cst routine) 61 USE trd mod ! momentum/tracers trends (trd_mod_init routine)60 USE trdini ! dyn/tra trends initialization (trd_init routine) 62 61 USE asminc ! assimilation increments 63 62 USE asmbkg ! writing out state trajectory … … 69 68 USE icbini ! handle bergs, initialisation 70 69 USE icbstp ! handle bergs, calving, themodynamics and transport 71 #if defined key_oasis372 70 USE cpl_oasis3 ! OASIS3 coupling 73 #elif defined key_oasis474 USE cpl_oasis4 ! OASIS4 coupling (not working)75 #endif76 71 USE c1d ! 1D configuration 77 72 USE step_c1d ! Time stepping loop for the 1D configuration … … 86 81 USE sbctide, ONLY: lk_tide 87 82 USE crsini ! initialise grid coarsening utility 88 USE lbcnfd, ONLY: isendto, nsndto ! Setup of north fold exchanges83 USE lbcnfd, ONLY: isendto, nsndto, nfsloop, nfeloop ! Setup of north fold exchanges 89 84 90 85 IMPLICIT NONE … … 121 116 !!---------------------------------------------------------------------- 122 117 ! 123 124 118 #if defined key_agrif 125 119 CALL Agrif_Init_Grids() ! AGRIF: set the meshes … … 139 133 # endif 140 134 #endif 141 142 135 ! check that all process are still there... If some process have an error, 143 136 ! they will never enter in step and other processes will wait until the end of the cpu time! … … 166 159 167 160 DO WHILE ( istp <= nitend .AND. nstop == 0 ) 168 169 161 #if defined key_agrif 170 162 CALL Agrif_Step( stp ) ! AGRIF: time stepping … … 172 164 CALL stp( istp ) ! standard time stepping 173 165 #endif 174 175 166 istp = istp + 1 176 167 IF( lk_mpp ) CALL mpp_max( nstop ) … … 201 192 ! 202 193 CALL nemo_closefile 194 ! 203 195 #if defined key_iomput 204 196 CALL xios_finalize ! end mpp communications with xios 205 # if defined key_oasis3 || defined key_oasis4 206 CALL cpl_prism_finalize ! end coupling and mpp communications with OASIS 207 # endif 197 IF( lk_cpl ) CALL cpl_finalize ! end coupling and mpp communications with OASIS 208 198 #else 209 # if defined key_oasis3 || defined key_oasis4 210 CALL cpl_prism_finalize! end coupling and mpp communications with OASIS211 # else 212 IF( lk_mpp ) CALL mppstop! end mpp communications213 # endif 199 IF( lk_cpl ) THEN 200 CALL cpl_finalize ! end coupling and mpp communications with OASIS 201 ELSE 202 IF( lk_mpp ) CALL mppstop ! end mpp communications 203 ENDIF 214 204 #endif 215 205 ! … … 227 217 INTEGER :: ios 228 218 CHARACTER(len=80), DIMENSION(16) :: cltxt 229 ! !230 NAMELIST/namctl/ ln_ctl , nn_print, nn_ictls, nn_ictle, &219 ! 220 NAMELIST/namctl/ ln_ctl , nn_print, nn_ictls, nn_ictle, & 231 221 & nn_isplt, nn_jsplt, nn_jctls, nn_jctle, & 232 222 & nn_bench, nn_timing … … 281 271 #if defined key_iomput 282 272 IF( Agrif_Root() ) THEN 283 # if defined key_oasis3 || defined key_oasis4 284 CALL cpl_prism_init( ilocal_comm )! nemo local communicator given by oasis285 CALL xios_initialize( "oceanx",local_comm=ilocal_comm )286 # else 287 CALL xios_initialize( "nemo",return_comm=ilocal_comm )288 # endif 273 IF( lk_cpl ) THEN 274 CALL cpl_init( ilocal_comm ) ! nemo local communicator given by oasis 275 CALL xios_initialize( "oceanx",local_comm=ilocal_comm ) ! send nemo communicator to xios 276 ELSE 277 CALL xios_initialize( "nemo",return_comm=ilocal_comm ) ! nemo local communicator given by xios 278 ENDIF 289 279 ENDIF 290 280 narea = mynode( cltxt, numnam_ref, numnam_cfg, numond , nstop, ilocal_comm ) ! Nodes selection 291 281 #else 292 # if defined key_oasis3 || defined key_oasis4 293 IF( Agrif_Root() ) THEN294 CALL cpl_prism_init( ilocal_comm )! nemo local communicator given by oasis295 ENDIF296 narea = mynode( cltxt, numnam_ref, numnam_cfg, numond , nstop, ilocal_comm ) ! Nodes selection (control print return in cltxt)297 # else 298 ilocal_comm = 0299 narea = mynode( cltxt, numnam_ref, numnam_cfg, numond , nstop )! Nodes selection (control print return in cltxt)300 # endif 282 IF( lk_cpl ) THEN 283 IF( Agrif_Root() ) THEN 284 CALL cpl_init( ilocal_comm ) ! nemo local communicator given by oasis 285 ENDIF 286 narea = mynode( cltxt, numnam_ref, numnam_cfg, numond , nstop, ilocal_comm ) ! Nodes selection (control print return in cltxt) 287 ELSE 288 ilocal_comm = 0 289 narea = mynode( cltxt, numnam_ref, numnam_cfg, numond , nstop ) ! Nodes selection (control print return in cltxt) 290 ENDIF 301 291 #endif 302 292 narea = narea + 1 ! mynode return the rank of proc (0 --> jpnij -1 ) … … 393 383 IF( lk_tide ) CALL tide_init( nit000 ) ! Initialisation of the tidal harmonics 394 384 395 IF( lk_bdy ) CALL bdy_init! Open boundaries initialisation396 IF( lk_bdy ) CALL bdy_dta_init! Open boundaries initialisation of external data arrays385 IF( lk_bdy ) CALL bdy_init ! Open boundaries initialisation 386 IF( lk_bdy ) CALL bdy_dta_init ! Open boundaries initialisation of external data arrays 397 387 IF( lk_bdy .AND. lk_tide ) & 398 & CALL bdytide_init! Open boundaries initialisation of tidal harmonic forcing388 & CALL bdytide_init ! Open boundaries initialisation of tidal harmonic forcing 399 389 400 390 CALL dyn_nept_init ! simplified form of Neptune effect … … 406 396 CALL sbc_init ! Forcings : surface module 407 397 ! ! Vertical physics 408 409 398 CALL zdf_init ! namelist read 410 411 399 CALL zdf_bfr_init ! bottom friction 412 413 400 IF( lk_zdfric ) CALL zdf_ric_init ! Richardson number dependent Kz 414 401 IF( lk_zdftke ) CALL zdf_tke_init ! TKE closure scheme … … 449 436 CALL trc_init 450 437 #endif 451 ! 452 453 ! Diagnostics 438 ! ! Diagnostics 454 439 IF( lk_floats ) CALL flo_init ! drifting Floats 455 440 IF( lk_diaar5 ) CALL dia_ar5_init ! ar5 diag … … 457 442 IF( lk_diadct ) CALL dia_dct_init ! Sections tranports 458 443 CALL dia_hsb_init ! heat content, salt content and volume budgets 459 CALL trd_mod_init ! Mixed-layer/Vorticity/Integral constraints trends444 CALL trd_init ! Mixed-layer/Vorticity/Integral constraints trends 460 445 IF( lk_diaobs ) THEN ! Observation & model comparison 461 446 CALL dia_obs_init ! Initialize observational data 462 447 CALL dia_obs( nit000 - 1 ) ! Observation operator for restart 463 448 ENDIF 449 464 450 ! ! Assimilation increments 465 451 IF( lk_asminc ) CALL asm_inc_init ! Initialize assimilation increments … … 568 554 ENDIF 569 555 ! 570 IF( lk_c1d .AND. .NOT.lk_iomput ) CALL ctl_stop( 'nemo_ctl: The 1D configuration must be used ', &571 & 'with the IOM Input/Output manager. ' , &572 & 'Compile with key_iomput enabled' )573 !574 556 IF( 1_wp /= SIGN(1._wp,-0._wp) ) CALL ctl_stop( 'nemo_ctl: The intrinsec SIGN function follows ', & 575 557 & 'f2003 standard. ' , & … … 663 645 !! ** Method : 664 646 !!---------------------------------------------------------------------- 665 INTEGER, INTENT(in) :: num_pes! The number of MPI processes we have647 INTEGER, INTENT(in) :: num_pes ! The number of MPI processes we have 666 648 ! 667 649 INTEGER, PARAMETER :: nfactmax = 20 … … 672 654 INTEGER, DIMENSION(nfactmax) :: ifact ! Array of factors 673 655 !!---------------------------------------------------------------------- 674 656 ! 675 657 ierr = 0 676 658 ! 677 659 CALL factorise( ifact, nfactmax, nfact, num_pes, ierr ) 678 660 ! 679 661 IF( nfact <= 1 ) THEN 680 662 WRITE (numout, *) 'WARNING: factorisation of number of PEs failed' … … 718 700 INTEGER, PARAMETER :: ntest = 14 719 701 INTEGER :: ilfax(ntest) 720 702 ! 721 703 ! lfax contains the set of allowed factors. 722 704 data (ilfax(jl),jl=1,ntest) / 16384, 8192, 4096, 2048, 1024, 512, 256, & … … 769 751 770 752 #if defined key_mpp_mpi 753 771 754 SUBROUTINE nemo_northcomms 772 755 !!====================================================================== … … 803 786 !loop over the other north-fold processes to find the processes 804 787 !managing the points belonging to the sxT-dxT range 805 DO jn = jpnij - jpni +1, jpnij806 IF ( njmppt(jn) == njmppmax ) THEN788 789 DO jn = 1, jpni 807 790 !sxT is the first point (in the global domain) of the jn 808 791 !process 809 sxT = n imppt(jn)792 sxT = nfiimpp(jn, jpnj) 810 793 !dxT is the last point (in the global domain) of the jn 811 794 !process 812 dxT = n imppt(jn) + nlcit(jn) - 1795 dxT = nfiimpp(jn, jpnj) + nfilcit(jn, jpnj) - 1 813 796 IF ((sxM .gt. sxT) .AND. (sxM .lt. dxT)) THEN 814 797 nsndto = nsndto + 1 815 isendto(nsndto) = jn798 isendto(nsndto) = jn 816 799 ELSEIF ((sxM .le. sxT) .AND. (dxM .ge. dxT)) THEN 817 800 nsndto = nsndto + 1 818 isendto(nsndto) = jn801 isendto(nsndto) = jn 819 802 ELSEIF ((dxM .lt. dxT) .AND. (sxT .lt. dxM)) THEN 820 803 nsndto = nsndto + 1 821 isendto(nsndto) = jn804 isendto(nsndto) = jn 822 805 END IF 823 END IF824 806 END DO 807 nfsloop = 1 808 nfeloop = nlci 809 DO jn = 2,jpni-1 810 IF(nfipproc(jn,jpnj) .eq. (narea - 1)) THEN 811 IF (nfipproc(jn - 1 ,jpnj) .eq. -1) THEN 812 nfsloop = nldi 813 ENDIF 814 IF (nfipproc(jn + 1,jpnj) .eq. -1) THEN 815 nfeloop = nlei 816 ENDIF 817 ENDIF 818 END DO 819 825 820 ENDIF 826 821 l_north_nogather = .TRUE. … … 831 826 END SUBROUTINE nemo_northcomms 832 827 #endif 828 833 829 !!====================================================================== 834 830 END MODULE nemogcm
Note: See TracChangeset
for help on using the changeset viewer.
