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 13540 for NEMO/branches/2020/r12377_ticket2386/src/OCE/ASM/asminc.F90 – NEMO

Ignore:
Timestamp:
2020-09-29T12:41:06+02:00 (4 years ago)
Author:
andmirek
Message:

Ticket #2386: update to latest trunk

Location:
NEMO/branches/2020/r12377_ticket2386
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • NEMO/branches/2020/r12377_ticket2386

    • Property svn:externals

      •  

        old new  
        33^/utils/build/mk@HEAD         mk 
        44^/utils/tools@HEAD            tools 
        5 ^/vendors/AGRIF/dev@HEAD      ext/AGRIF 
         5^/vendors/AGRIF/dev_r12970_AGRIF_CMEMS      ext/AGRIF 
        66^/vendors/FCM@HEAD            ext/FCM 
        77^/vendors/IOIPSL@HEAD         ext/IOIPSL 
        88 
        99# SETTE 
        10 ^/utils/CI/sette@HEAD         sette 
         10^/utils/CI/sette@13507        sette 

  • NEMO/branches/2020/r12377_ticket2386/src/OCE/ASM/asminc.F90

    r12511 r13540  
    9595   !! * Substitutions 
    9696#  include "do_loop_substitute.h90" 
     97#  include "domzgr_substitute.h90" 
    9798   !!---------------------------------------------------------------------- 
    9899   !! NEMO/OCE 4.0 , NEMO Consortium (2018) 
    … …  
    359360 
    360361         IF ( ln_trainc ) THEN    
    361             CALL iom_get( inum, jpdom_autoglo, 'bckint', t_bkginc, 1 ) 
    362             CALL iom_get( inum, jpdom_autoglo, 'bckins', s_bkginc, 1 ) 
     362            CALL iom_get( inum, jpdom_auto, 'bckint', t_bkginc, 1 ) 
     363            CALL iom_get( inum, jpdom_auto, 'bckins', s_bkginc, 1 ) 
    363364            ! Apply the masks 
    364365            t_bkginc(:,:,:) = t_bkginc(:,:,:) * tmask(:,:,:) 
    … …  
    371372 
    372373         IF ( ln_dyninc ) THEN    
    373             CALL iom_get( inum, jpdom_autoglo, 'bckinu', u_bkginc, 1 )               
    374             CALL iom_get( inum, jpdom_autoglo, 'bckinv', v_bkginc, 1 )               
     374            CALL iom_get( inum, jpdom_auto, 'bckinu', u_bkginc, 1 )               
     375            CALL iom_get( inum, jpdom_auto, 'bckinv', v_bkginc, 1 )               
    375376            ! Apply the masks 
    376377            u_bkginc(:,:,:) = u_bkginc(:,:,:) * umask(:,:,:) 
    … …  
    383384         
    384385         IF ( ln_sshinc ) THEN 
    385             CALL iom_get( inum, jpdom_autoglo, 'bckineta', ssh_bkginc, 1 ) 
     386            CALL iom_get( inum, jpdom_auto, 'bckineta', ssh_bkginc, 1 ) 
    386387            ! Apply the masks 
    387388            ssh_bkginc(:,:) = ssh_bkginc(:,:) * tmask(:,:,1) 
    … …  
    392393 
    393394         IF ( ln_seaiceinc ) THEN 
    394             CALL iom_get( inum, jpdom_autoglo, 'bckinseaice', seaice_bkginc, 1 ) 
     395            CALL iom_get( inum, jpdom_auto, 'bckinseaice', seaice_bkginc, 1 ) 
    395396            ! Apply the masks 
    396397            seaice_bkginc(:,:) = seaice_bkginc(:,:) * tmask(:,:,1) 
    … …  
    413414            DO jk = 1, jpkm1           ! zhdiv = e1e1 * div 
    414415               zhdiv(:,:) = 0._wp 
    415                DO_2D_00_00 
     416               DO_2D( 0, 0, 0, 0 ) 
    416417                  zhdiv(ji,jj) = (  e2u(ji  ,jj) * e3u(ji  ,jj,jk,Kmm) * u_bkginc(ji  ,jj,jk)    & 
    417418                     &            - e2u(ji-1,jj) * e3u(ji-1,jj,jk,Kmm) * u_bkginc(ji-1,jj,jk)    & 
    418419                     &            + e1v(ji,jj  ) * e3v(ji,jj  ,jk,Kmm) * v_bkginc(ji,jj  ,jk)    & 
    419                      &            - e1v(ji,jj-1) * e3v(ji,jj-1,jk,Kmm) * v_bkginc(ji,jj-1,jk)  ) / e3t(ji,jj,jk,Kmm) 
     420                     &            - e1v(ji,jj-1) * e3v(ji,jj-1,jk,Kmm) * v_bkginc(ji,jj-1,jk)  ) & 
     421                     &            / e3t(ji,jj,jk,Kmm) 
    420422               END_2D 
    421                CALL lbc_lnk( 'asminc', zhdiv, 'T', 1. )   ! lateral boundary cond. (no sign change) 
     423               CALL lbc_lnk( 'asminc', zhdiv, 'T', 1.0_wp )   ! lateral boundary cond. (no sign change) 
    422424               ! 
    423                DO_2D_00_00 
     425               DO_2D( 0, 0, 0, 0 ) 
    424426                  u_bkginc(ji,jj,jk) = u_bkginc(ji,jj,jk)                         & 
    425427                     &               + 0.2_wp * ( zhdiv(ji+1,jj) - zhdiv(ji  ,jj) ) * r1_e1u(ji,jj) * umask(ji,jj,jk) 
    … …  
    465467         ! 
    466468         IF ( ln_trainc ) THEN    
    467             CALL iom_get( inum, jpdom_autoglo, 'tn', t_bkg ) 
    468             CALL iom_get( inum, jpdom_autoglo, 'sn', s_bkg ) 
     469            CALL iom_get( inum, jpdom_auto, 'tn', t_bkg ) 
     470            CALL iom_get( inum, jpdom_auto, 'sn', s_bkg ) 
    469471            t_bkg(:,:,:) = t_bkg(:,:,:) * tmask(:,:,:) 
    470472            s_bkg(:,:,:) = s_bkg(:,:,:) * tmask(:,:,:) 
    … …  
    472474         ! 
    473475         IF ( ln_dyninc ) THEN    
    474             CALL iom_get( inum, jpdom_autoglo, 'un', u_bkg ) 
    475             CALL iom_get( inum, jpdom_autoglo, 'vn', v_bkg ) 
     476            CALL iom_get( inum, jpdom_auto, 'un', u_bkg, cd_type = 'U', psgn = 1._wp ) 
     477            CALL iom_get( inum, jpdom_auto, 'vn', v_bkg, cd_type = 'V', psgn = 1._wp ) 
    476478            u_bkg(:,:,:) = u_bkg(:,:,:) * umask(:,:,:) 
    477479            v_bkg(:,:,:) = v_bkg(:,:,:) * vmask(:,:,:) 
    … …  
    479481         ! 
    480482         IF ( ln_sshinc ) THEN 
    481             CALL iom_get( inum, jpdom_autoglo, 'sshn', ssh_bkg ) 
     483            CALL iom_get( inum, jpdom_auto, 'sshn', ssh_bkg ) 
    482484            ssh_bkg(:,:) = ssh_bkg(:,:) * tmask(:,:,1) 
    483485         ENDIF 
    … …  
    758760            ! 
    759761            ssh(:,:,Kbb) = ssh(:,:,Kmm)                        ! Update before fields 
     762#if ! defined key_qco 
    760763            e3t(:,:,:,Kbb) = e3t(:,:,:,Kmm) 
     764#endif 
    761765!!gm why not e3u(:,:,:,Kbb), e3v(:,:,:,Kbb), gdept(:,:,:,Kbb) ???? 
    762766            ! 
    … …  
    896900         IF ( kt == nitdin_r ) THEN 
    897901            ! 
    898             l_1st_euler = 0              ! Force Euler forward step 
     902            l_1st_euler = .TRUE.              ! Force Euler forward step 
    899903            ! 
    900904            ! Sea-ice : SI3 case 
    … …  
    970974!           ! set to bottom of a level  
    971975!                 DO jk = jpk-1, 2, -1 
    972 !                   IF ((mld > gdepw(ji,jj,jk)) .and. (mld < gdepw(ji,jj,jk+1))) THEN  
    973 !                     mld=gdepw(ji,jj,jk+1) 
     976!                   IF ((mld > gdepw(ji,jj,jk,Kmm)) .and. (mld < gdepw(ji,jj,jk+1,Kmm))) THEN 
     977!                     mld=gdepw(ji,jj,jk+1,Kmm) 
    974978!                     jkmax=jk 
    975979!                   ENDIF 
Note: See TracChangeset for help on using the changeset viewer.