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 14072 for NEMO/trunk/tests/STATION_ASF/MY_SRC/nemogcm.F90 – NEMO

Ignore:
Timestamp:
2020-12-04T08:48:38+01:00 (3 years ago)
Author:
laurent
Message:

Merging branch "2020/dev_r13648_ASINTER-04_laurent_bulk_ice", ticket #2369

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/tests/STATION_ASF/MY_SRC/nemogcm.F90

    r13286 r14072  
    3030   USE step_c1d       ! Time stepping loop for the 1D configuration 
    3131   ! 
    32    USE prtctl         ! Print control 
    3332   USE in_out_manager ! I/O manager 
    3433   USE lib_mpp        ! distributed memory computing 
     
    4746   CHARACTER(lc) ::   cform_aaa="( /, 'AAAAAAAA', / ) "     ! flag for output listing 
    4847 
     48#if defined key_mpp_mpi 
     49   ! need MPI_Wtime 
     50   INCLUDE 'mpif.h' 
     51#endif 
     52 
    4953   !!---------------------------------------------------------------------- 
    5054   !! NEMO/OCE 4.0 , NEMO Consortium (2018) 
    51    !! $Id: nemogcm.F90 12489 2020-02-28 15:55:11Z davestorkey $ 
     55   !! $Id: nemogcm.F90 13286 2020-07-09 15:48:29Z smasson $ 
    5256   !! Software governed by the CeCILL license (see ./LICENSE) 
    5357   !!---------------------------------------------------------------------- 
     
    110114      ! 
    111115#if defined key_iomput 
    112                                     CALL xios_finalize  ! end mpp communications with xios 
     116      CALL xios_finalize  ! end mpp communications with xios 
    113117#else 
    114118      IF( lk_mpp )                  CALL mppstop      ! end mpp communications 
     
    146150#if defined key_iomput 
    147151      IF( Agrif_Root() ) THEN 
    148             CALL xios_initialize( "for_xios_mpi_id", return_comm=ilocal_comm )   ! nemo local communicator given by xios 
     152         CALL xios_initialize( "for_xios_mpi_id", return_comm=ilocal_comm )   ! nemo local communicator given by xios 
    149153      ENDIF 
    150154      CALL mpp_start( ilocal_comm ) 
    151155#else 
    152          CALL mpp_start( ) 
     156      CALL mpp_start( ) 
    153157#endif 
    154158      ! 
     
    163167      IF( lwm )   CALL ctl_opn(     numout,        'ocean.output', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, -1, .FALSE. ) 
    164168      ! open reference and configuration namelist files 
    165                   CALL load_nml( numnam_ref,        'namelist_ref',                                           -1, lwm ) 
    166                   CALL load_nml( numnam_cfg,        'namelist_cfg',                                           -1, lwm ) 
     169      CALL load_nml( numnam_ref,        'namelist_ref',                                           -1, lwm ) 
     170      CALL load_nml( numnam_cfg,        'namelist_cfg',                                           -1, lwm ) 
    167171      IF( lwm )   CALL ctl_opn(     numond, 'output.namelist.dyn', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, -1, .FALSE. ) 
    168172      ! open /dev/null file to be able to supress output write easily 
    169173      IF( Agrif_Root() ) THEN 
    170                   CALL ctl_opn(     numnul,           '/dev/null', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, -1, .FALSE. ) 
     174         CALL ctl_opn(     numnul,           '/dev/null', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, -1, .FALSE. ) 
    171175#ifdef key_agrif 
    172176      ELSE 
    173                   numnul = Agrif_Parent(numnul)    
    174 #endif 
    175       ENDIF 
    176       ! 
     177         numnul = Agrif_Parent(numnul) 
     178#endif 
     179      ENDIF 
    177180      !                             !--------------------! 
    178181      !                             ! Open listing units !  -> need sn_cfctl from namctl to define lwp 
     
    215218         WRITE(numout,*) "     ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ " 
    216219         WRITE(numout,*) 
     220 
     221         ! Print the working precision to ocean.output 
     222         IF (wp == dp) THEN 
     223            WRITE(numout,*) "Working precision = double-precision" 
     224         ELSE 
     225            WRITE(numout,*) "Working precision = single-precision" 
     226         ENDIF 
     227         WRITE(numout,*) 
    217228         ! 
    218229         WRITE(numout,cform_aaa)                                        ! Flag AAAAAAA 
     
    229240903   IF( ios /= 0 )   CALL ctl_nam ( ios , 'namcfg in reference namelist' ) 
    230241      READ  ( numnam_cfg, namcfg, IOSTAT = ios, ERR = 904 ) 
    231 904   IF( ios >  0 )   CALL ctl_nam ( ios , 'namcfg in configuration namelist' )    
     242904   IF( ios >  0 )   CALL ctl_nam ( ios , 'namcfg in configuration namelist' ) 
    232243      ! 
    233244      IF( ln_read_cfg ) THEN            ! Read sizes in domain configuration file 
     
    263274                           CALL     eos_init        ! Equation of state 
    264275      IF( lk_c1d       )   CALL     c1d_init        ! 1D column configuration 
    265                            CALL     dom_init( Nbb, Nnn, Naa, "OPA") ! Domain 
     276                           CALL     dom_init( Nbb, Nnn, Naa ) ! Domain 
    266277      IF( sn_cfctl%l_prtctl )   & 
    267278         &                 CALL prt_ctl_init        ! Print control 
    268       ! 
    269        
     279 
    270280                           CALL  istate_init( Nbb, Nnn, Naa )    ! ocean initial state (Dynamics and tracers) 
    271281 
    272       !                                      ! external forcing  
     282      !                                      ! external forcing 
    273283                           CALL     sbc_init( Nbb, Nnn, Naa )    ! surface boundary conditions (including sea-ice) 
     284 
     285      !#LB: 
     286#if defined key_si3 
     287      IF(lwp) WRITE(numout,*) 'LOLO: nemo_init@nemogcm.F90: shape of fr_i ==>', SIZE(fr_i,1), SIZE(fr_i,2) 
     288      fr_i(:,:) = 0._wp 
     289#endif 
     290      !#LB. 
    274291 
    275292      ! 
     
    302319         WRITE(numout,*) '                              sn_cfctl%l_prttrc  = ', sn_cfctl%l_prttrc 
    303320         WRITE(numout,*) '                              sn_cfctl%l_oasout  = ', sn_cfctl%l_oasout 
    304          WRITE(numout,*) '                              sn_cfctl%procmin   = ', sn_cfctl%procmin   
    305          WRITE(numout,*) '                              sn_cfctl%procmax   = ', sn_cfctl%procmax   
    306          WRITE(numout,*) '                              sn_cfctl%procincr  = ', sn_cfctl%procincr  
    307          WRITE(numout,*) '                              sn_cfctl%ptimincr  = ', sn_cfctl%ptimincr  
     321         WRITE(numout,*) '                              sn_cfctl%procmin   = ', sn_cfctl%procmin 
     322         WRITE(numout,*) '                              sn_cfctl%procmax   = ', sn_cfctl%procmax 
     323         WRITE(numout,*) '                              sn_cfctl%procincr  = ', sn_cfctl%procincr 
     324         WRITE(numout,*) '                              sn_cfctl%ptimincr  = ', sn_cfctl%ptimincr 
    308325         WRITE(numout,*) '      timing by routine               ln_timing  = ', ln_timing 
    309326         WRITE(numout,*) '      CFL diagnostics                 ln_diacfl  = ', ln_diacfl 
     
    366383      !!---------------------------------------------------------------------- 
    367384      ! 
    368       ierr =        oce_alloc    ()    ! ocean  
     385      ierr =        oce_alloc    ()    ! ocean 
    369386      ierr = ierr + dia_wri_alloc() 
    370387      ierr = ierr + dom_oce_alloc()    ! ocean domain 
     
    375392   END SUBROUTINE nemo_alloc 
    376393 
    377     
     394 
    378395   SUBROUTINE nemo_set_cfctl(sn_cfctl, setto ) 
    379396      !!---------------------------------------------------------------------- 
     
    399416   !!====================================================================== 
    400417END MODULE nemogcm 
    401  
Note: See TracChangeset for help on using the changeset viewer.