Changeset 6140 for trunk/NEMOGCM/NEMO/LIM_SRC_2/limthd_2.F90
- Timestamp:
- 2015-12-21T12:35:23+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/LIM_SRC_2/limthd_2.F90
r5836 r6140 18 18 USE phycst ! physical constants 19 19 USE dom_oce ! ocean space and time domain variables 20 USE domvvl 21 USE lbclnk 20 USE domvvl ! ocean domain 21 USE ice_2 ! LIM sea-ice variables 22 USE sbc_oce ! surface boundary condition: ocean 23 USE sbc_ice ! surface boundary condition: sea-ice 24 USE thd_ice_2 ! LIM thermodynamic sea-ice variables 25 USE dom_ice_2 ! LIM sea-ice domain 26 USE limthd_zdf_2 ! 27 USE limthd_lac_2 ! 28 USE limtab_2 ! 29 ! 22 30 USE in_out_manager ! I/O manager 23 USE lib_mpp 31 USE lbclnk ! 32 USE lib_mpp ! 24 33 USE wrk_nemo ! work arrays 25 34 USE iom ! IOM library 26 USE ice_2 ! LIM sea-ice variables27 USE sbc_oce !28 USE sbc_ice !29 USE thd_ice_2 ! LIM thermodynamic sea-ice variables30 USE dom_ice_2 ! LIM sea-ice domain31 USE limthd_zdf_232 USE limthd_lac_233 USE limtab_234 35 USE prtctl ! Print control 35 36 USE lib_fortran ! Fortran utilities (allows no signed zero when 'key_nosignedzero' defined) … … 43 44 REAL(wp) :: epsi16 = 1.e-16 ! 44 45 REAL(wp) :: epsi04 = 1.e-04 ! 45 REAL(wp) :: rzero = 0. e0!46 REAL(wp) :: rone = 1. e0!46 REAL(wp) :: rzero = 0._wp ! 47 REAL(wp) :: rone = 1._wp ! 47 48 48 49 !! * Substitutions 49 # include "domzgr_substitute.h90"50 50 # include "vectopt_loop_substitute.h90" 51 51 !!-------- ------------------------------------------------------------- … … 75 75 !!--------------------------------------------------------------------- 76 76 INTEGER, INTENT(in) :: kt ! number of iteration 77 ! !77 ! 78 78 INTEGER :: ji, jj ! dummy loop indices 79 79 INTEGER :: nbpb ! nb of icy pts for thermo. cal. … … 233 233 234 234 ! energy needed to bring ocean surface layer until its freezing 235 qcmif (ji,jj) = rau0 * rcp * fse3t_m(ji,jj) * ( tfu(ji,jj) - sst_m(ji,jj) - rt0 ) * ( 1 - zinda )235 qcmif (ji,jj) = rau0 * rcp * e3t_m(ji,jj) * ( tfu(ji,jj) - sst_m(ji,jj) - rt0 ) * ( 1 - zinda ) 236 236 237 237 ! calculate oceanic heat flux.
Note: See TracChangeset
for help on using the changeset viewer.