Changeset 4234 for branches/2013/dev_LOCEAN_CMCC_INGV_2013
- Timestamp:
- 2013-11-18T14:30:46+01:00 (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2013/dev_LOCEAN_CMCC_INGV_2013/NEMOGCM/NEMO/OPA_SRC/DIA/diahsb.F90
r4205 r4234 112 112 DO jk = 1, jpkm1 113 113 ! volume variation (calculated with scale factors) 114 zdiff_v2 = zdiff_v2 + glob_sum( zsurf(:,:) * tmask(:,:,jk) * ( fse3t_n(:,:,jk) - e3t_ini(:,:,jk) ) ) 114 zdiff_v2 = zdiff_v2 & 115 & + glob_sum( zsurf(:,:) * tmask(:,:,jk) * ( fse3t_n(:,:,jk) - e3t_ini(:,:,jk) ) ) 115 116 ! heat content variation 116 zdiff_hc = zdiff_hc + glob_sum( zsurf(:,:) * tmask(:,:,jk) * ( fse3t_n(:,:,jk) * tsn(:,:,jk,jp_tem) & 117 zdiff_hc = zdiff_hc & 118 & + glob_sum( zsurf(:,:) * tmask(:,:,jk) * ( fse3t_n(:,:,jk) * tsn(:,:,jk,jp_tem) & 117 119 & - hc_loc_ini(:,:,jk) ) ) 118 120 ! salt content variation 119 zdiff_sc = zdiff_sc + glob_sum( zsurf(:,:) * tmask(:,:,jk) * ( fse3t_n(:,:,jk) * tsn(:,:,jk,jp_sal) & 121 zdiff_sc = zdiff_sc & 122 & + glob_sum( zsurf(:,:) * tmask(:,:,jk) * ( fse3t_n(:,:,jk) * tsn(:,:,jk,jp_sal) & 120 123 & - sc_loc_ini(:,:,jk) ) ) 121 124 ENDDO … … 145 148 DO jk = 1, jpkm1 146 149 ! volume (calculated with scale factors) 147 z_v2 = z_v2 + glob_sum( zsurf(:,:) * tmask(:,:,jk) * fse3t_n(:,:,jk) ) 150 z_v2 = z_v2 & 151 & + glob_sum( zsurf(:,:) * tmask(:,:,jk) * fse3t_n(:,:,jk) ) 148 152 ! heat content 149 z_hc = z_hc + glob_sum( zsurf(:,:) * tmask(:,:,jk) * fse3t_n(:,:,jk) * tsn(:,:,jk,jp_tem) ) 153 z_hc = z_hc & 154 & + glob_sum( zsurf(:,:) * tmask(:,:,jk) * fse3t_n(:,:,jk) * tsn(:,:,jk,jp_tem) ) 150 155 ! salt content 151 z_sc = z_sc + glob_sum( zsurf(:,:) * tmask(:,:,jk) * fse3t_n(:,:,jk) * tsn(:,:,jk,jp_sal) ) 156 z_sc = z_sc & 157 & + glob_sum( zsurf(:,:) * tmask(:,:,jk) * fse3t_n(:,:,jk) * tsn(:,:,jk,jp_sal) ) 152 158 ENDDO 153 159 ! add ssh if not vvl
Note: See TracChangeset
for help on using the changeset viewer.