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/OCE/ICB/icbini.F90 – 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:
2 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/OCE/ICB/icbini.F90

    r13786 r13899  
    8080      ssu_e(:,:) = 0._wp   ;   ssv_e(:,:) = 0._wp   ; 
    8181      ua_e(:,:)  = 0._wp   ;   va_e(:,:)  = 0._wp   ; 
    82       ff_e(:,:)  = 0._wp   ;   sst_e(:,:)  = 0._wp   ; 
    83       fr_e(:,:)  = 0._wp   ; 
     82      ff_e(:,:)  = 0._wp   ;   sst_e(:,:) = 0._wp   ; 
     83      fr_e(:,:)  = 0._wp   ;   sss_e(:,:) = 0._wp   ; 
    8484      ! 
    8585      IF ( ln_M2016 ) THEN 
    … …  
    137137      nicbfldproc(:) = -1 
    138138 
    139       DO_2D_11_11 
     139      DO_2D( 1, 1, 1, 1 ) 
    140140         src_calving_hflx(ji,jj) = narea 
    141141         src_calving     (ji,jj) = nicbpack * mjg(jj) + mig(ji) 
    … …  
    147147      ! first entry with narea for this processor is left hand interior index 
    148148      ! last  entry                               is right hand interior index 
    149       jj = nlcj/2 
     149      jj = jpj/2 
    150150      nicbdi = -1 
    151151      nicbei = -1 
    … …  
    163163      ! 
    164164      ! repeat for j direction 
    165       ji = nlci/2 
     165      ji = jpi/2 
    166166      nicbdj = -1 
    167167      nicbej = -1 
    … …  
    180180      ! special for east-west boundary exchange we save the destination index 
    181181      i1 = MAX( nicbdi-1, 1) 
    182       i3 = INT( src_calving(i1,nlcj/2) ) 
     182      i3 = INT( src_calving(i1,jpj/2) ) 
    183183      jj = INT( i3/nicbpack ) 
    184184      ricb_left = REAL( i3 - nicbpack*jj, wp ) 
    185185      i1 = MIN( nicbei+1, jpi ) 
    186       i3 = INT( src_calving(i1,nlcj/2) ) 
     186      i3 = INT( src_calving(i1,jpj/2) ) 
    187187      jj = INT( i3/nicbpack ) 
    188188      ricb_right = REAL( i3 - nicbpack*jj, wp ) 
    … …  
    217217         WRITE(numicb,*) 'processor ', narea 
    218218         WRITE(numicb,*) 'jpi, jpj   ', jpi, jpj 
    219          WRITE(numicb,*) 'nldi, nlei ', nldi, nlei 
    220          WRITE(numicb,*) 'nldj, nlej ', nldj, nlej 
     219         WRITE(numicb,*) 'Nis0, Nie0 ', Nis0, Nie0 
     220         WRITE(numicb,*) 'Njs0, Nje0 ', Njs0, Nje0 
    221221         WRITE(numicb,*) 'berg i interior ', nicbdi, nicbei 
    222222         WRITE(numicb,*) 'berg j interior ', nicbdj, nicbej 
    223223         WRITE(numicb,*) 'berg left       ', ricb_left 
    224224         WRITE(numicb,*) 'berg right      ', ricb_right 
    225          jj = nlcj/2 
     225         jj = jpj/2 
    226226         WRITE(numicb,*) "central j line:" 
    227227         WRITE(numicb,*) "i processor" 
    … …  
    229229         WRITE(numicb,*) "i point" 
    230230         WRITE(numicb,*) (INT(src_calving(ji,jj)), ji=1,jpi) 
    231          ji = nlci/2 
     231         ji = jpi/2 
    232232         WRITE(numicb,*) "central i line:" 
    233233         WRITE(numicb,*) "j processor" 
    … …  
    280280         ivar = iom_varid( inum, 'maxclass', ldstop=.FALSE. ) 
    281281         IF( ivar > 0 ) THEN 
    282             CALL iom_get  ( inum, jpdom_data, 'maxclass', src_calving )   ! read the max distribution array 
     282            CALL iom_get  ( inum, jpdom_global, 'maxclass', src_calving )   ! read the max distribution array 
    283283            berg_grid%maxclass(:,:) = INT( src_calving ) 
    284284            src_calving(:,:) = 0._wp 
Note: See TracChangeset for help on using the changeset viewer.