Changeset 3988


Ignore:
Timestamp:
2013-07-24T07:36:14+02:00 (7 years ago)
Author:
cbricaud
Message:

bugfix diadct for ice transport ( LIM3 case ) ; see ticket 1080

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/NEMO/OPA_SRC/DIA/diadct.F90

    r3452 r3988  
    4242#endif 
    4343#if defined key_lim3 
    44   USE ice_3 
     44  USE par_ice 
     45  USE ice 
    4546#endif 
    4647  USE domvvl 
     
    566567     
    567568     !! * Local variables 
    568      INTEGER             :: jk,jseg,jclass,   &!loop on level/segment/classes  
     569     INTEGER             :: jk,jseg,jclass,jl,   &!loop on level/segment/classes/ice type  
    569570                            isgnu  , isgnv     ! 
    570571     INTEGER :: ii, ij ! local integer 
     
    816817              zTnorm=zumid_ice*e2u(k%I,k%J)+zvmid_ice*e1v(k%I,k%J) 
    817818    
     819#if defined key_lim2    
    818820              IF( zTnorm .GE. 0)THEN 
    819821                 zice_vol_pos = (zTnorm)*   & 
     
    835837                                     +zice_surf_neg 
    836838              ENDIF 
    837     
     839#endif 
     840#if defined key_lim3 
     841              DO jl=1,jpl 
     842 
     843                 IF( zTnorm .GE. 0)THEN 
     844                    zice_vol_pos = zice_vol_pos + zTnorm*   & 
     845                                      a_i(sec%listPoint(jseg)%I,sec%listPoint(jseg)%J,jl)     * & 
     846                                     ( ht_i(sec%listPoint(jseg)%I,sec%listPoint(jseg)%J,jl) +   & 
     847                                       ht_s(sec%listPoint(jseg)%I,sec%listPoint(jseg)%J,jl) ) 
     848 
     849                    zice_surf_pos = zice_surf_pos + zTnorm * a_i(sec%listPoint(jseg)%I,sec%listPoint(jseg)%J,jl)  
     850                 ELSE 
     851 
     852                    zice_vol_neg = zice_vol_neg + zTnorm*   & 
     853                                      a_i(sec%listPoint(jseg)%I,sec%listPoint(jseg)%J,jl)     * & 
     854                                     ( ht_i(sec%listPoint(jseg)%I,sec%listPoint(jseg)%J,jl) +   & 
     855                                       ht_s(sec%listPoint(jseg)%I,sec%listPoint(jseg)%J,jl) ) 
     856                                  
     857                    zice_surf_neg = zice_surf_neg + zTnorm * a_i(sec%listPoint(jseg)%I,sec%listPoint(jseg)%J,jl) 
     858                 ENDIF 
     859 
     860              ENDDO 
     861#endif 
     862  
    838863              zsum(11,1) = zsum(11,1)+zice_vol_pos 
    839864              zsum(12,1) = zsum(12,1)+zice_vol_neg 
Note: See TracChangeset for help on using the changeset viewer.