Changeset 3294 for trunk/NEMOGCM/NEMO/OPA_SRC/DIA/diahsb.F90
- Timestamp:
- 2012-01-28T17:44:18+01:00 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/OPA_SRC/DIA/diahsb.F90
r2528 r3294 18 18 USE lib_mpp ! distributed memory computing library 19 19 USE trabbc ! bottom boundary condition 20 USE obc_par ! (for lk_obc) 20 21 USE bdy_par ! (for lk_bdy) 21 USE obc_par ! (for lk_obc)22 USE timing ! preformance summary 22 23 23 24 IMPLICIT NONE … … 72 73 REAL(dp) :: z_frc_trd_v ! - - 73 74 !!--------------------------------------------------------------------------- 75 IF( nn_timing == 1 ) CALL timing_start('dia_hsb') 74 76 75 77 ! ------------------------- ! … … 107 109 ! heat content variation 108 110 zdiff_hc = zdiff_hc + SUM( surf(:,:) * tmask(:,:,jk) & 109 & * ( fse3t_n(:,:,jk) * t n(:,:,jk) &111 & * ( fse3t_n(:,:,jk) * tsn(:,:,jk,jp_tem) & 110 112 & - hc_loc_ini(:,:,jk) ) ) 111 113 ! salt content variation 112 114 zdiff_sc = zdiff_sc + SUM( surf(:,:) * tmask(:,:,jk) & 113 & * ( fse3t_n(:,:,jk) * sn(:,:,jk) &115 & * ( fse3t_n(:,:,jk) * tsn(:,:,jk,jp_sal) & 114 116 & - sc_loc_ini(:,:,jk) ) ) 115 117 ENDDO … … 139 141 IF ( kt == nitend ) CLOSE( numhsb ) 140 142 143 IF( nn_timing == 1 ) CALL timing_stop('dia_hsb') 144 141 145 9020 FORMAT(I5,11D15.7) 142 146 ! … … 205 209 WRITE(numout,*) "dia_hsb: heat salt volume budgets activated" 206 210 WRITE(numout,*) "~~~~~~~ output written in the 'heat_salt_volume_budgets.txt' ASCII file" 207 IF( lk_obc . OR. lk_bdy) THEN211 IF( lk_obc .or. lk_bdy ) THEN 208 212 CALL ctl_warn( 'dia_hsb does not take open boundary fluxes into account' ) 209 213 ENDIF … … 248 252 ! 4 - initial conservation variables ! 249 253 ! ---------------------------------- ! 250 ssh_ini(:,:) = sshn(:,:) ! initial ssh254 ssh_ini(:,:) = sshn(:,:) ! initial ssh 251 255 DO jk = 1, jpk 252 e3t_ini (:,:,jk) = fse3t_n(:,:,jk) ! initial vertical scale factors253 hc_loc_ini(:,:,jk) = t n(:,:,jk) * fse3t_n(:,:,jk) ! initial heat content254 sc_loc_ini(:,:,jk) = sn(:,:,jk) * fse3t_n(:,:,jk) ! initial salt content256 e3t_ini (:,:,jk) = fse3t_n(:,:,jk) ! initial vertical scale factors 257 hc_loc_ini(:,:,jk) = tsn(:,:,jk,jp_tem) * fse3t_n(:,:,jk) ! initial heat content 258 sc_loc_ini(:,:,jk) = tsn(:,:,jk,jp_sal) * fse3t_n(:,:,jk) ! initial salt content 255 259 END DO 256 260 frc_v = 0.d0 ! volume trend due to forcing
Note: See TracChangeset
for help on using the changeset viewer.