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/SBC/sbcice_cice.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/SBC/sbcice_cice.F90

    r13237 r13899  
    217217! T point to U point 
    218218! T point to V point 
    219       DO_2D_10_10 
     219      DO_2D( 1, 0, 1, 0 ) 
    220220         fr_iu(ji,jj)=0.5*(fr_i(ji,jj)+fr_i(ji+1,jj))*umask(ji,jj,1) 
    221221         fr_iv(ji,jj)=0.5*(fr_i(ji,jj)+fr_i(ji,jj+1))*vmask(ji,jj,1) 
    … …  
    312312! x comp of wind stress (CI_1) 
    313313! U point to F point 
    314          DO_2D_10_11 
     314         DO_2D( 1, 0, 1, 1 ) 
    315315            ztmp(ji,jj) = 0.5 * (  fr_iu(ji,jj) * utau(ji,jj)      & 
    316316                                 + fr_iu(ji,jj+1) * utau(ji,jj+1) ) * fmask(ji,jj,1) 
    … …  
    320320! y comp of wind stress (CI_2) 
    321321! V point to F point 
    322          DO_2D_11_10 
     322         DO_2D( 1, 1, 1, 0 ) 
    323323            ztmp(ji,jj) = 0.5 * (  fr_iv(ji,jj) * vtau(ji,jj)      & 
    324324                                 + fr_iv(ji+1,jj) * vtau(ji+1,jj) ) * fmask(ji,jj,1) 
    … …  
    335335            qla_ice(:,:,1)= - ( emp_ice(:,:)+sprecip(:,:) ) * rLsub 
    336336! End of temporary code 
    337             DO_2D_11_11 
     337            DO_2D( 1, 1, 1, 1 ) 
    338338               IF(fr_i(ji,jj).eq.0.0) THEN 
    339339                  DO jl=1,ncat 
    … …  
    437437! x comp and y comp of surface ocean current 
    438438! U point to F point 
    439       DO_2D_10_11 
     439      DO_2D( 1, 0, 1, 1 ) 
    440440         ztmp(ji,jj)=0.5*(ssu_m(ji,jj)+ssu_m(ji,jj+1))*fmask(ji,jj,1) 
    441441      END_2D 
    … …  
    443443 
    444444! V point to F point 
    445       DO_2D_11_10 
     445      DO_2D( 1, 1, 1, 0 ) 
    446446         ztmp(ji,jj)=0.5*(ssv_m(ji,jj)+ssv_m(ji+1,jj))*fmask(ji,jj,1) 
    447447      END_2D 
    … …  
    467467! x comp and y comp of sea surface slope (on F points) 
    468468! T point to F point 
    469       DO_2D_10_10 
     469      DO_2D( 1, 0, 1, 0 ) 
    470470         ztmp(ji,jj)=0.5 * (  (zpice(ji+1,jj  )-zpice(ji,jj  )) * r1_e1u(ji,jj  )    & 
    471471            &               + (zpice(ji+1,jj+1)-zpice(ji,jj+1)) * r1_e1u(ji,jj+1)  ) * fmask(ji,jj,1) 
    … …  
    474474 
    475475! T point to F point 
    476       DO_2D_10_10 
     476      DO_2D( 1, 0, 1, 0 ) 
    477477         ztmp(ji,jj)=0.5 * (  (zpice(ji  ,jj+1)-zpice(ji  ,jj)) * r1_e2v(ji  ,jj)    & 
    478478            &               + (zpice(ji+1,jj+1)-zpice(ji+1,jj)) * r1_e2v(ji+1,jj)  ) *  fmask(ji,jj,1) 
    … …  
    503503      ss_iou(:,:)=0.0 
    504504! F point to U point 
    505       DO_2D_00_00 
     505      DO_2D( 0, 0, 0, 0 ) 
    506506         ss_iou(ji,jj) = 0.5 * ( ztmp1(ji,jj-1) + ztmp1(ji,jj) ) * umask(ji,jj,1) 
    507507      END_2D 
    … …  
    513513! F point to V point 
    514514 
    515       DO_2D_10_00 
     515      DO_2D( 1, 0, 0, 0 ) 
    516516         ss_iov(ji,jj) = 0.5 * ( ztmp1(ji-1,jj) + ztmp1(ji,jj) ) * vmask(ji,jj,1) 
    517517      END_2D 
    … …  
    597597      CALL lbc_lnk( 'sbcice_cice', qsr , 'T', 1.0_wp ) 
    598598 
    599       DO_2D_11_11 
     599      DO_2D( 1, 1, 1, 1 ) 
    600600         nfrzmlt(ji,jj)=MAX(nfrzmlt(ji,jj),0.0) 
    601601      END_2D 
    … …  
    621621! T point to U point 
    622622! T point to V point 
    623       DO_2D_10_10 
     623      DO_2D( 1, 0, 1, 0 ) 
    624624         fr_iu(ji,jj)=0.5*(fr_i(ji,jj)+fr_i(ji+1,jj))*umask(ji,jj,1) 
    625625         fr_iv(ji,jj)=0.5*(fr_i(ji,jj)+fr_i(ji,jj+1))*vmask(ji,jj,1) 
    … …  
    880880!        pcg(:,:)=0.0 
    881881         DO jn=1,jpnij 
    882             DO jj=nldjt(jn),nlejt(jn) 
    883                DO ji=nldit(jn),nleit(jn) 
     882            DO jj=njs0all(jn),nje0all(jn) 
     883               DO ji=nis0all(jn),nie0all(jn) 
    884884                  png2(ji+nimppt(jn)-1,jj+njmppt(jn)-1)=png(ji,jj,jn) 
    885885               ENDDO 
    … …  
    981981 
    982982      pn(:,:)=0.0 
    983       DO_2D_10_10 
     983      DO_2D( 1, 0, 1, 0 ) 
    984984         pn(ji,jj)=pc(ji+1-ji_off,jj+1-jj_off,1) 
    985985      END_2D 
    … …  
    10011001         png(:,:,:)=0.0 
    10021002         DO jn=1,jpnij 
    1003             DO jj=nldjt(jn),nlejt(jn) 
    1004                DO ji=nldit(jn),nleit(jn) 
     1003            DO jj=njs0all(jn),nje0all(jn) 
     1004               DO ji=nis0all(jn),nie0all(jn) 
    10051005                  png(ji,jj,jn)=pcg(ji+nimppt(jn)-1-ji_off,jj+njmppt(jn)-1-jj_off) 
    10061006               ENDDO 
Note: See TracChangeset for help on using the changeset viewer.