- Timestamp:
- 2013-07-24T07:36:14+02:00 (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/NEMO/OPA_SRC/DIA/diadct.F90
r3452 r3988 42 42 #endif 43 43 #if defined key_lim3 44 USE ice_3 44 USE par_ice 45 USE ice 45 46 #endif 46 47 USE domvvl … … 566 567 567 568 !! * Local variables 568 INTEGER :: jk,jseg,jclass, &!loop on level/segment/classes569 INTEGER :: jk,jseg,jclass,jl, &!loop on level/segment/classes/ice type 569 570 isgnu , isgnv ! 570 571 INTEGER :: ii, ij ! local integer … … 816 817 zTnorm=zumid_ice*e2u(k%I,k%J)+zvmid_ice*e1v(k%I,k%J) 817 818 819 #if defined key_lim2 818 820 IF( zTnorm .GE. 0)THEN 819 821 zice_vol_pos = (zTnorm)* & … … 835 837 +zice_surf_neg 836 838 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 838 863 zsum(11,1) = zsum(11,1)+zice_vol_pos 839 864 zsum(12,1) = zsum(12,1)+zice_vol_neg
Note: See TracChangeset
for help on using the changeset viewer.