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 13899 for NEMO/branches/2020/tickets_icb_1900/src/SAS – NEMO

Ignore:
Timestamp:
2020-11-27T17:26:33+01:00 (4 years ago)
Author:
mathiot
Message:

ticket #1900: update branch to trunk and add ICB test case

Location:
NEMO/branches/2020/tickets_icb_1900
Files:
4 edited

Legend:

Unmodified
Added
Removed

  • NEMO/branches/2020/tickets_icb_1900

    • Property svn:externals

      •  

        old new  
        22^/utils/build/makenemo@HEAD   makenemo 
        33^/utils/build/mk@HEAD         mk 
        4 ^/utils/tools/@HEAD           tools 
         4^/utils/tools@HEAD            tools 
        55^/vendors/AGRIF/dev_r12970_AGRIF_CMEMS      ext/AGRIF 
        66^/vendors/FCM@HEAD            ext/FCM 
        … …  
        88 
        99# SETTE 
        10 ^/utils/CI/sette@12931        sette 
         10^/utils/CI/sette@13559        sette 

  • NEMO/branches/2020/tickets_icb_1900/src/SAS/nemogcm.F90

    r13365 r13899  
    22   !!====================================================================== 
    33   !!                       ***  MODULE nemogcm   *** 
    4    !! StandAlone Surface module : surface fluxes + sea-ice + iceberg floats 
     4   !! StandAlone Surface module : surface fluxes + sea-ice + iceberg floats + ABL 
    55   !!====================================================================== 
    66   !! History :  3.6  ! 2011-11  (S. Alderson, G. Madec) original code 
    … …  
    3636   USE icb_oce        ! icebergs 
    3737   ! 
     38   USE prtctl         ! Print control 
    3839   USE in_out_manager ! I/O manager 
    3940   USE lib_mpp        ! distributed memory computing 
    4041   USE mppini         ! shared/distributed memory setting (mpp_init routine) 
    41    USE lbcnfd  , ONLY : isendto, nsndto, nfsloop, nfeloop  ! Setup of north fold exchanges 
     42   USE lbcnfd  , ONLY : isendto, nsndto ! Setup of north fold exchanges 
    4243   USE lib_fortran    ! Fortran utilities (allows no signed zero when 'key_nosignedzero' defined) 
    4344#if defined key_iomput 
    … …  
    4748   USE agrif_ice_update ! ice update 
    4849#endif 
     50   USE halo_mng 
    4951 
    5052   IMPLICIT NONE 
    … …  
    5759 
    5860#if defined key_mpp_mpi 
     61   ! need MPI_Wtime 
    5962   INCLUDE 'mpif.h' 
    6063#endif 
    … …  
    8285      !!---------------------------------------------------------------------- 
    8386      INTEGER ::   istp   ! time step index 
     87      REAL(wp)::   zstptiming   ! elapsed time for 1 time step 
    8488      !!---------------------------------------------------------------------- 
    8589      ! 
    … …  
    9296#if defined key_agrif 
    9397      Kbb_a = Nbb; Kmm_a = Nnn; Krhs_a = Nrhs   ! agrif_oce module copies of time level indices 
    94       CALL Agrif_Declare_Var       !  "      "   "   "      "  DYN/TRA  
     98      CALL Agrif_Declare_Var       !  "      "   "   "      "  DYN/TRA 
    9599# if defined key_top 
    96100      CALL Agrif_Declare_Var_top   !  "      "   "   "      "  TOP 
    … …  
    106110      !                            !==   time stepping   ==! 
    107111      !                            !-----------------------! 
     112      ! 
     113      !                                               !== set the model time-step  ==! 
     114      ! 
    108115      istp = nit000 
    109116      ! 
    … …  
    123130      END DO 
    124131      ! 
    125 #else 
     132# else 
    126133      ! 
    127134      IF( .NOT.ln_diurnal_only ) THEN                 !==  Standard time-stepping  ==! 
    128135         ! 
    129136         DO WHILE( istp <= nitend .AND. nstop == 0 ) 
    130 #if defined key_mpp_mpi 
     137 
    131138            ncom_stp = istp 
    132             IF ( istp == ( nit000 + 1 ) ) elapsed_time = MPI_Wtime() 
    133             IF ( istp ==         nitend ) elapsed_time = MPI_Wtime() - elapsed_time 
    134 #endif 
     139            IF( ln_timing ) THEN 
     140               zstptiming = MPI_Wtime() 
     141               IF ( istp == ( nit000 + 1 ) ) elapsed_time = zstptiming 
     142               IF ( istp ==         nitend ) elapsed_time = zstptiming - elapsed_time 
     143            ENDIF 
     144             
    135145            CALL stp        ( istp )  
    136146            istp = istp + 1 
     147 
     148            IF( lwp .AND. ln_timing )   WRITE(numtime,*) 'timing step ', istp-1, ' : ', MPI_Wtime() - zstptiming 
     149 
    137150         END DO 
    138151         ! 
    … …  
    198211      INTEGER ::   ios, ilocal_comm   ! local integers 
    199212      !! 
    200       NAMELIST/namctl/ sn_cfctl, nn_print, nn_ictls, nn_ictle,              & 
    201          &             nn_isplt , nn_jsplt, nn_jctls, nn_jctle,             & 
    202          &             ln_timing, ln_diacfl 
     213      NAMELIST/namctl/ sn_cfctl, ln_timing, ln_diacfl,                                & 
     214         &             nn_isplt,  nn_jsplt,  nn_ictls, nn_ictle, nn_jctls, nn_jctle             
    203215      NAMELIST/namcfg/ ln_read_cfg, cn_domcfg, ln_closea, ln_write_cfg, cn_domcfg_out, ln_use_jattr 
    204216      !!---------------------------------------------------------------------- 
    … …  
    207219      ELSE                  ;   cxios_context = 'nemo' 
    208220      ENDIF 
     221      nn_hls = 1 
    209222      ! 
    210223      !                             !-------------------------------------------------! 
    … …  
    304317         WRITE(numout,*) "       )  )       \) |`\ \)  '.   \      (   (   " 
    305318         WRITE(numout,*) "      (  (           \_/       '-._\      )   )  " 
    306          WRITE(numout,*) "       )  ) jgs                    `     (   (   " 
     319         WRITE(numout,*) "       )  ) jgs                     `    (   (   " 
    307320         WRITE(numout,*) "     ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ " 
    308321         WRITE(numout,*) 
    … …  
    325338      ! 
    326339      IF( ln_read_cfg ) THEN            ! Read sizes in domain configuration file 
    327          CALL domain_cfg ( cn_cfg, nn_cfg, jpiglo, jpjglo, jpkglo, jperio ) 
     340         CALL domain_cfg ( cn_cfg, nn_cfg, Ni0glo, Nj0glo, jpkglo, jperio ) 
    328341      ELSE                              ! user-defined namelist 
    329          CALL usr_def_nam( cn_cfg, nn_cfg, jpiglo, jpjglo, jpkglo, jperio ) 
     342         CALL usr_def_nam( cn_cfg, nn_cfg, Ni0glo, Nj0glo, jpkglo, jperio ) 
    330343      ENDIF 
    331344      ! 
    … …  
    337350      CALL mpp_init 
    338351 
     352      CALL halo_mng_init() 
    339353      ! Now we know the dimensions of the grid and numout has been set: we can allocate arrays 
    340354      CALL nemo_alloc() 
    … …  
    353367      ! 
    354368      !                                      ! General initialization 
    355       IF( ln_timing    )   CALL timing_init     ! timing 
     369      IF( ln_timing    )   CALL timing_init ( 'timing_sas.output' ) 
    356370      IF( ln_timing    )   CALL timing_start( 'nemo_init') 
    357371 
    … …  
    365379         &                 CALL prt_ctl_init        ! Print control 
    366380       
     381      IF( ln_rstart )      CALL rst_read_open 
    367382                           CALL day_init        ! model calendar (using both namelist and restart infos) 
    368       IF( ln_rstart )      CALL rst_read_open 
    369  
     383 
     384#if defined key_agrif 
     385      uu(:,:,:,:) = 0.0_wp   ;   vv(:,:,:,:) = 0.0_wp   ;   ts(:,:,:,:,:) = 0.0_wp   ! needed for interp done at initialization phase 
     386#endif  
    370387      !                                      ! external forcing  
    371388                           CALL sbc_init( Nbb, Nnn, Naa )  ! Forcings : surface module  
    … …  
    416433         WRITE(numout,*) '                              sn_cfctl%procincr  = ', sn_cfctl%procincr  
    417434         WRITE(numout,*) '                              sn_cfctl%ptimincr  = ', sn_cfctl%ptimincr  
    418          WRITE(numout,*) '      level of print                  nn_print   = ', nn_print 
    419          WRITE(numout,*) '      Start i indice for SUM control  nn_ictls   = ', nn_ictls 
    420          WRITE(numout,*) '      End i indice for SUM control    nn_ictle   = ', nn_ictle 
    421          WRITE(numout,*) '      Start j indice for SUM control  nn_jctls   = ', nn_jctls 
    422          WRITE(numout,*) '      End j indice for SUM control    nn_jctle   = ', nn_jctle 
    423          WRITE(numout,*) '      number of proc. following i     nn_isplt   = ', nn_isplt 
    424          WRITE(numout,*) '      number of proc. following j     nn_jsplt   = ', nn_jsplt 
    425435         WRITE(numout,*) '      timing by routine               ln_timing  = ', ln_timing 
    426436         WRITE(numout,*) '      CFL diagnostics                 ln_diacfl  = ', ln_diacfl 
    427437      ENDIF 
    428438      ! 
    429       nprint    = nn_print          ! convert DOCTOR namelist names into OLD names 
    430       nictls    = nn_ictls 
    431       nictle    = nn_ictle 
    432       njctls    = nn_jctls 
    433       njctle    = nn_jctle 
    434       isplt     = nn_isplt 
    435       jsplt     = nn_jsplt 
    436  
     439      IF( .NOT.ln_read_cfg )   ln_closea = .FALSE.   ! dealing possible only with a domcfg file 
    437440      IF(lwp) THEN                  ! control print 
    438441         WRITE(numout,*) 
    … …  
    445448         WRITE(numout,*) '      use file attribute if exists as i/p j-start   ln_use_jattr     = ', ln_use_jattr 
    446449      ENDIF 
    447       IF( .NOT.ln_read_cfg )   ln_closea = .false.   ! dealing possible only with a domcfg file 
    448       ! 
    449       !                             ! Parameter control 
    450       ! 
    451       IF( sn_cfctl%l_prtctl .OR. sn_cfctl%l_prttrc ) THEN              ! sub-domain area indices for the control prints 
    452          IF( lk_mpp .AND. jpnij > 1 ) THEN 
    453             isplt = jpni   ;   jsplt = jpnj   ;   ijsplt = jpni*jpnj   ! the domain is forced to the real split domain 
    454          ELSE 
    455             IF( isplt == 1 .AND. jsplt == 1  ) THEN 
    456                CALL ctl_warn( ' - isplt & jsplt are equal to 1',   & 
    457                   &           ' - the print control will be done over the whole domain' ) 
    458             ENDIF 
    459             ijsplt = isplt * jsplt            ! total number of processors ijsplt 
    460          ENDIF 
    461          IF(lwp) WRITE(numout,*)'          - The total number of processors over which the' 
    462          IF(lwp) WRITE(numout,*)'            print control will be done is ijsplt : ', ijsplt 
    463          ! 
    464          !                              ! indices used for the SUM control 
    465          IF( nictls+nictle+njctls+njctle == 0 )   THEN    ! print control done over the default area 
    466             lsp_area = .FALSE. 
    467          ELSE                                             ! print control done over a specific  area 
    468             lsp_area = .TRUE. 
    469             IF( nictls < 1 .OR. nictls > jpiglo )   THEN 
    470                CALL ctl_warn( '          - nictls must be 1<=nictls>=jpiglo, it is forced to 1' ) 
    471                nictls = 1 
    472             ENDIF 
    473             IF( nictle < 1 .OR. nictle > jpiglo )   THEN 
    474                CALL ctl_warn( '          - nictle must be 1<=nictle>=jpiglo, it is forced to jpiglo' ) 
    475                nictle = jpiglo 
    476             ENDIF 
    477             IF( njctls < 1 .OR. njctls > jpjglo )   THEN 
    478                CALL ctl_warn( '          - njctls must be 1<=njctls>=jpjglo, it is forced to 1' ) 
    479                njctls = 1 
    480             ENDIF 
    481             IF( njctle < 1 .OR. njctle > jpjglo )   THEN 
    482                CALL ctl_warn( '          - njctle must be 1<=njctle>=jpjglo, it is forced to jpjglo' ) 
    483                njctle = jpjglo 
    484             ENDIF 
    485          ENDIF 
    486       ENDIF 
    487450      ! 
    488451      IF( 1._wp /= SIGN(1._wp,-0._wp)  )   CALL ctl_stop( 'nemo_ctl: The intrinsec SIGN function follows f2003 standard.',  & 
    … …  
    538501      ierr =        dia_wri_alloc() 
    539502      ierr = ierr + dom_oce_alloc()          ! ocean domain 
    540       ierr = ierr + oce_alloc    ()          ! (tsn...) needed for agrif and/or SI3 and bdy 
     503      ierr = ierr + oce_alloc    ()          ! (ts...) needed for agrif and/or SI3 and bdy 
    541504      ierr = ierr + bdy_oce_alloc()          ! bdy masks (incl. initialization) 
    542505      ! 

  • NEMO/branches/2020/tickets_icb_1900/src/SAS/sbcssm.F90

    r13365 r13899  
    294294            !                                         ! fill sf with slf_i and control print 
    295295            CALL fld_fill( sf_ssm_3d, slf_3d, cn_dir, 'sbc_ssm_init', '3D Data in file', 'namsbc_ssm' ) 
     296            sf_ssm_3d(jf_usp)%cltype = 'U'   ;   sf_ssm_3d(jf_usp)%zsgn = -1._wp 
     297            sf_ssm_3d(jf_vsp)%cltype = 'V'   ;   sf_ssm_3d(jf_vsp)%zsgn = -1._wp 
    296298         ENDIF 
    297299         ! 
    … …  
    310312            ! 
    311313            CALL fld_fill( sf_ssm_2d, slf_2d, cn_dir, 'sbc_ssm_init', '2D Data in file', 'namsbc_ssm' ) 
     314            IF( .NOT. ln_3d_uve ) THEN 
     315               sf_ssm_2d(jf_usp)%cltype = 'U'   ;   sf_ssm_2d(jf_usp)%zsgn = -1._wp 
     316               sf_ssm_2d(jf_vsp)%cltype = 'V'   ;   sf_ssm_2d(jf_vsp)%zsgn = -1._wp 
     317            ENDIF 
    312318         ENDIF 
    313319         ! 

  • NEMO/branches/2020/tickets_icb_1900/src/SAS/stpctl.F90

    r13136 r13899  
    2020   USE dom_oce         ! ocean space and time domain variables  
    2121   USE ice      , ONLY : vt_i, u_ice, tm_i 
     22   USE phycst   , ONLY : rt0 
     23   USE sbc_oce  , ONLY : lk_oasis 
    2224   ! 
    2325   USE diawri          ! Standard run outputs       (dia_wri_state routine) 
    … …  
    4850      !! 
    4951      !! ** Method  : - Save the time step in numstp 
    50       !!              - Print it each 50 time steps 
    5152      !!              - Stop the run IF problem encountered by setting nstop > 0 
    5253      !!                Problems checked: ice thickness maximum > 100 m 
    … …  
    6768      REAL(wp)                        ::   zzz                                   ! local real  
    6869      REAL(wp), DIMENSION(4)          ::   zmax, zmaxlocal 
    69       LOGICAL                         ::   ll_wrtstp, ll_colruns, ll_wrtruns 
     70      LOGICAL                         ::   ll_wrtstp, ll_colruns, ll_wrtruns, ll_0oce 
    7071      LOGICAL, DIMENSION(jpi,jpj)     ::   llmsk 
    7172      CHARACTER(len=20)               ::   clname 
    … …  
    8586         ENDIF 
    8687         !                                ! open time.step    ascii file, done only by 1st subdomain 
    87          IF( lwm )   CALL ctl_opn( numstp, 'time.step', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp, narea ) 
     88         IF( lk_oasis ) THEN   ;   clname = 'time_sas.step' 
     89         ELSE                  ;   clname = 'time.step' 
     90         ENDIF 
     91         IF( lwm )   CALL ctl_opn( numstp, clname, 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp, narea ) 
    8892         ! 
    8993         IF( ll_wrtruns ) THEN 
     94            IF( lk_oasis ) THEN   ;   clname = 'run_sas.stat' 
     95            ELSE                  ;   clname = 'run.stat' 
     96            ENDIF 
    9097            !                             ! open run.stat     ascii file, done only by 1st subdomain 
    91             CALL ctl_opn( numrun, 'run.stat', 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp, narea ) 
     98            CALL ctl_opn( numrun, clname, 'REPLACE', 'FORMATTED', 'SEQUENTIAL', -1, numout, lwp, narea ) 
    9299            !                             ! open run.stat.nc netcdf file, done only by 1st subdomain 
    93             clname = 'run.stat.nc' 
     100            clname = TRIM(clname)//'.nc' 
    94101            IF( .NOT. Agrif_Root() )   clname = TRIM(Agrif_CFixed())//"_"//TRIM(clname) 
    95102            istatus = NF90_CREATE( TRIM(clname), NF90_CLOBBER, nrunid ) 
    … …  
    111118      !                                   !==            test of local extrema           ==! 
    112119      !                                   !==  done by all processes at every time step  ==! 
    113       llmsk(:,:) = tmask(:,:,1) == 1._wp 
    114       IF( COUNT( llmsk(:,:) ) > 0 ) THEN   ! avoid huge values sent back for land processors... 
    115          zmax(1) = MAXVAL(      vt_i (:,:)            , mask = llmsk )   ! max ice thickness 
    116          zmax(2) = MAXVAL( ABS( u_ice(:,:) )          , mask = llmsk )   ! max ice velocity (zonal only) 
    117          zmax(3) = MAXVAL(     -tm_i (:,:) + 273.15_wp, mask = llmsk )   ! min ice temperature 
    118       ELSE 
    119          IF( ll_colruns ) THEN    ! default value: must not be kept when calling mpp_max -> must be as small as possible 
    120             zmax(1:3) = -HUGE(1._wp) 
    121          ELSE                     ! default value: must not give true for any of the tests bellow (-> avoid manipulating HUGE...) 
    122             zmax(1:3) = 0._wp 
    123          ENDIF 
    124       ENDIF 
    125       zmax(4) = REAL( nstop, wp )                                     ! stop indicator 
     120      ! 
     121      llmsk(   1:Nis1,:) = .FALSE.                                              ! exclude halos from the checked region 
     122      llmsk(Nie1: jpi,:) = .FALSE. 
     123      llmsk(:,   1:Njs1) = .FALSE. 
     124      llmsk(:,Nje1: jpj) = .FALSE. 
     125      ! 
     126      llmsk(Nis0:Nie0,Njs0:Nje0) = tmask(Nis0:Nie0,Njs0:Nje0,1) == 1._wp        ! test only the inner domain 
     127      ! 
     128      ll_0oce = .NOT. ANY( llmsk(:,:) )                                         ! no ocean point in the inner domain? 
     129      ! 
     130      zmax(1) = MAXVAL(      vt_i (:,:)      , mask = llmsk )                   ! max ice thickness 
     131      zmax(2) = MAXVAL( ABS( u_ice(:,:) )    , mask = llmsk )                   ! max ice velocity (zonal only) 
     132      zmax(3) = MAXVAL(     -tm_i (:,:) + rt0, mask = llmsk )                   ! min ice temperature (in degC) 
     133      zmax(4) = REAL( nstop, wp )                                               ! stop indicator 
     134      ! 
    126135      !                                   !==               get global extrema             ==! 
    127136      !                                   !==  done by all processes if writting run.stat  ==! 
    128137      IF( ll_colruns ) THEN 
    129138         zmaxlocal(:) = zmax(:) 
    130          CALL mpp_max( "stpctl", zmax )          ! max over the global domain 
     139         CALL mpp_max( "stpctl", zmax )          ! max over the global domain: ok even of ll_0oce = .true. 
    131140         nstop = NINT( zmax(4) )                 ! update nstop indicator (now sheared among all local domains) 
    132       ENDIF 
     141      ELSE 
     142         ! if no ocean point: MAXVAL returns -HUGE => we must overwrite this value to avoid error handling bellow. 
     143         IF( ll_0oce )   zmax(1:3) = 0._wp       ! default "valid" values... 
     144      ENDIF 
     145      ! 
     146      zmax(3) = -zmax(3)                              ! move back from max(-zz) to min(zz) : easier to manage! 
     147      IF( ll_colruns ) zmaxlocal(3) = -zmaxlocal(3)   ! move back from max(-zz) to min(zz) : easier to manage! 
     148      ! 
    133149      !                                   !==              write "run.stat" files              ==! 
    134150      !                                   !==  done only by 1st subdomain at writting timestep  ==! 
    135151      IF( ll_wrtruns ) THEN 
    136          WRITE(numrun,9500) kt, zmax(1), zmax(2), -zmax(3) 
    137          istatus = NF90_PUT_VAR( nrunid, nvarid(1), (/ zmax(1)/), (/kt/), (/1/) ) 
    138          istatus = NF90_PUT_VAR( nrunid, nvarid(2), (/ zmax(2)/), (/kt/), (/1/) ) 
    139          istatus = NF90_PUT_VAR( nrunid, nvarid(3), (/-zmax(3)/), (/kt/), (/1/) ) 
     152         WRITE(numrun,9500) kt, zmax(1), zmax(2), zmax(3) 
     153         DO ji = 1, 3 
     154            istatus = NF90_PUT_VAR( nrunid, nvarid(ji), (/zmax(ji)/), (/kt/), (/1/) ) 
     155         END DO 
    140156         IF( kt == nitend )   istatus = NF90_CLOSE(nrunid) 
    141157      END IF 
    … …  
    145161      IF(   zmax(1) >  100._wp .OR.   &                   ! too large ice thickness maximum ( > 100 m) 
    146162         &  zmax(2) >   10._wp .OR.   &                   ! too large ice velocity ( > 10 m/s) 
    147          &  zmax(3) >  101._wp .OR.   &                   ! too cold ice temperature ( < -100 degC) 
     163         &  zmax(3) < -101._wp .OR.   &                   ! too cold ice temperature ( < -100 degC) 
    148164         &  ISNAN( zmax(1) + zmax(2) + zmax(3) ) .OR.   &               ! NaN encounter in the tests 
    149165         &  ABS(   zmax(1) + zmax(2) + zmax(3) ) > HUGE(1._wp) ) THEN   ! Infinity encounter in the tests 
    … …  
    154170            IF( lwm .AND. kt /= nitend )   istatus = NF90_CLOSE(nrunid) 
    155171            ! get global loc on the min/max 
    156             CALL mpp_maxloc( 'stpctl',      vt_i(:,:)            , tmask(:,:,1), zzz, iloc(1:2,1) )   ! mpp_maxloc ok if mask = F  
    157             CALL mpp_maxloc( 'stpctl',ABS( u_ice(:,:) )          , tmask(:,:,1), zzz, iloc(1:2,2) ) 
    158             CALL mpp_minloc( 'stpctl',      tm_i(:,:) - 273.15_wp, tmask(:,:,1), zzz, iloc(1:2,3) ) 
     172            CALL mpp_maxloc( 'stpctl',      vt_i(:,:)      , llmsk, zzz, iloc(1:2,1) )   ! mpp_maxloc ok if mask = F  
     173            CALL mpp_maxloc( 'stpctl',ABS( u_ice(:,:) )    , llmsk, zzz, iloc(1:2,2) ) 
     174            CALL mpp_minloc( 'stpctl',      tm_i(:,:) - rt0, llmsk, zzz, iloc(1:2,3) ) 
    159175            ! find which subdomain has the max. 
    160176            iareamin(:) = jpnij+1   ;   iareamax(:) = 0   ;   iareasum(:) = 0 
    … …  
    169185         ELSE                    ! find local min and max locations: 
    170186            ! if we are here, this means that the subdomain contains some oce points -> no need to test the mask used in maxloc 
    171             iloc(1:2,1) = MAXLOC(       vt_i(:,:)            , mask = llmsk ) + (/ nimpp - 1, njmpp - 1/) 
    172             iloc(1:2,2) = MAXLOC( ABS( u_ice(:,:) )          , mask = llmsk ) + (/ nimpp - 1, njmpp - 1/) 
    173             iloc(1:2,3) = MINLOC(       tm_i(:,:) - 273.15_wp, mask = llmsk ) + (/ nimpp - 1, njmpp - 1/) 
     187            iloc(1:2,1) = MAXLOC(       vt_i(:,:)      , mask = llmsk ) 
     188            iloc(1:2,2) = MAXLOC( ABS( u_ice(:,:) )    , mask = llmsk ) 
     189            iloc(1:2,3) = MINLOC(       tm_i(:,:) - rt0, mask = llmsk ) 
     190            DO ji = 1, 3   ! local domain indices ==> global domain indices, excluding halos 
     191               iloc(1:2,ji) = (/ mig0(iloc(1,ji)), mjg0(iloc(2,ji)) /) 
     192            END DO 
    174193            iareamin(:) = narea   ;   iareamax(:) = narea   ;   iareasum(:) = 1         ! this is local information 
    175194         ENDIF 
    176195         ! 
    177196         WRITE(ctmp1,*) ' stp_ctl: ice_thick > 100 m or |ice_vel| > 10 m/s or ice_temp < -100 degC or NaN encounter in the tests' 
    178          CALL wrt_line( ctmp2, kt, 'ice_thick max',  zmax(1), iloc(:,1), iareasum(1), iareamin(1), iareamax(1) ) 
    179          CALL wrt_line( ctmp3, kt, '|ice_vel| max',  zmax(2), iloc(:,2), iareasum(2), iareamin(2), iareamax(2) ) 
    180          CALL wrt_line( ctmp4, kt, 'ice_temp  min', -zmax(3), iloc(:,3), iareasum(3), iareamin(3), iareamax(3) ) 
     197         CALL wrt_line( ctmp2, kt, 'ice_thick max', zmax(1), iloc(:,1), iareasum(1), iareamin(1), iareamax(1) ) 
     198         CALL wrt_line( ctmp3, kt, '|ice_vel| max', zmax(2), iloc(:,2), iareasum(2), iareamin(2), iareamax(2) ) 
     199         CALL wrt_line( ctmp4, kt, 'ice_temp  min', zmax(3), iloc(:,3), iareasum(3), iareamin(3), iareamax(3) ) 
    181200         IF( Agrif_Root() ) THEN 
    182201            WRITE(ctmp6,*) '      ===> output of last computed fields in output.abort* files' 
Note: See TracChangeset for help on using the changeset viewer.