- Timestamp:
- 2020-04-03T18:54:55+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/DIA/diahsb.F90
r12625 r12680 57 57 !! Software governed by the CeCILL license (see ./LICENSE) 58 58 !!---------------------------------------------------------------------- 59 60 59 CONTAINS 61 60 … … 159 158 ! 160 159 DO jk = 1, jpkm1 ! volume variation (calculated with scale factors) 161 zwrk(:,:,jk) = surf (:,:) * e3t(:,:,jk,Kmm)*tmask(:,:,jk) &160 zwrk(:,:,jk) = surf (:,:) * e3t (:,:,jk,Kmm)*tmask (:,:,jk) & 162 161 & - surf_ini(:,:) * e3t_ini(:,:,jk )*tmask_ini(:,:,jk) 163 162 END DO 164 163 zdiff_v2 = glob_sum_full( 'diahsb', zwrk(:,:,:) ) ! glob_sum_full needed as tmask and tmask_ini could be different 165 164 DO jk = 1, jpkm1 ! heat content variation 166 zwrk(:,:,jk) = ( surf (:,:)*e3t(:,:,jk,Kmm)*ts(:,:,jk,jp_tem,Kmm) &167 & - surf_ini(:,:) *hc_loc_ini(:,:,jk) )165 zwrk(:,:,jk) = ( surf (:,:) * e3t(:,:,jk,Kmm)*ts(:,:,jk,jp_tem,Kmm) & 166 & - surf_ini(:,:) * hc_loc_ini(:,:,jk) ) 168 167 END DO 169 168 zdiff_hc = glob_sum_full( 'diahsb', zwrk(:,:,:) ) 170 169 DO jk = 1, jpkm1 ! salt content variation 171 zwrk(:,:,jk) = ( surf (:,:)*e3t(:,:,jk,Kmm)*ts(:,:,jk,jp_sal,Kmm) &172 & - surf_ini(:,:) *sc_loc_ini(:,:,jk) )170 zwrk(:,:,jk) = ( surf (:,:) * e3t(:,:,jk,Kmm)*ts(:,:,jk,jp_sal,Kmm) & 171 & - surf_ini(:,:) * sc_loc_ini(:,:,jk) ) 173 172 END DO 174 173 zdiff_sc = glob_sum_full( 'diahsb', zwrk(:,:,:) )
Note: See TracChangeset
for help on using the changeset viewer.