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 4153 for branches/2013/dev_LOCEAN_2013/NEMOGCM/NEMO/OPA_SRC/DIA/diadct.F90 – NEMO

Ignore:
Timestamp:
2013-11-05T13:25:45+01:00 (10 years ago)
Author:
cetlod
Message:

dev_LOCEAN_2013: merge in trunk changes between r3940 and r4028, see ticket #1169

File:
1 edited

Legend:

Unmodified
Added
Removed

  • branches/2013/dev_LOCEAN_2013/NEMOGCM/NEMO/OPA_SRC/DIA/diadct.F90

    r4147 r4153  
    4242#endif 
    4343#if defined key_lim3 
    44   USE ice_3 
     44  USE par_ice 
     45  USE ice 
    4546#endif 
    4647  USE domvvl 
    … …  
    490491                 ijglo = secs(jsec)%listPoint(jpt)%J + jpjzoom - 1 + njmpp - 1 
    491492                 WRITE(numout,*)'         # I J : ',iiglo,ijglo 
     493                 CALL FLUSH(numout) 
    492494              ENDDO 
    493495           ENDIF 
    … …  
    612614     
    613615     !! * Local variables 
    614      INTEGER             :: jk, jseg, jclass,                    &!loop on level/segment/classes   
     616     INTEGER             :: jk, jseg, jclass,jl,                 &!loop on level/segment/classes/ice categories 
    615617                            isgnu, isgnv                          !  
    616618     REAL(wp)            :: zumid, zvmid,                        &!U/V velocity on a cell segment  
    … …  
    777779    
    778780              zTnorm=zumid_ice*e2u(k%I,k%J)+zvmid_ice*e1v(k%I,k%J) 
    779     
     781 
     782#if defined key_lim2    
    780783              transports_2d(1,jsec,jseg) = transports_2d(1,jsec,jseg) + (zTnorm)*   &  
    781784                                   (1.0 - frld(sec%listPoint(jseg)%I,sec%listPoint(jseg)%J))  &  
    … …  
    784787              transports_2d(2,jsec,jseg) = transports_2d(2,jsec,jseg) + (zTnorm)*   &  
    785788                                    (1.0 -  frld(sec%listPoint(jseg)%I,sec%listPoint(jseg)%J)) 
     789#endif 
     790#if defined key_lim3 
     791              DO jl=1,jpl 
     792                 transports_2d(1,jsec,jseg) = transports_2d(1,jsec,jseg) + (zTnorm)*     & 
     793                                   a_i(sec%listPoint(jseg)%I,sec%listPoint(jseg)%J,jl) * & 
     794                                  ( ht_i(sec%listPoint(jseg)%I,sec%listPoint(jseg)%J,jl) +  & 
     795                                    ht_s(sec%listPoint(jseg)%I,sec%listPoint(jseg)%J,jl) ) 
     796                                    
     797                 transports_2d(2,jsec,jseg) = transports_2d(2,jsec,jseg) + (zTnorm)*   & 
     798                                   a_i(sec%listPoint(jseg)%I,sec%listPoint(jseg)%J,jl) 
     799              ENDDO 
     800#endif 
    786801    
    787802           ENDIF !end of ice case 
Note: See TracChangeset for help on using the changeset viewer.