- Timestamp:
- 2010-04-30T17:49:04+02:00 (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/DEV_r1837_mass_heat_salt_fluxes/NEMO/LIM_SRC_2/limdmp_2.F90
r1715 r1855 4 4 !! Ice model : restoring Ice thickness and Fraction leads 5 5 !!====================================================================== 6 !! History : 2.0 ! 6 !! History : 2.0 ! 2004-04 (S. Theetten) Original code 7 7 !!---------------------------------------------------------------------- 8 8 #if defined key_lim2 && defined key_tradmp … … 11 11 !! 'key_tradmp' Damping 12 12 !!---------------------------------------------------------------------- 13 !! lim_dmp_2 : ice model damping 13 14 !!---------------------------------------------------------------------- 14 !! lim_dmp_2 : ice model damping 15 !!---------------------------------------------------------------------- 16 USE in_out_manager ! I/O manager 17 USE phycst ! physical constants 18 USE ice_2 19 USE tradmp 20 USE dom_oce 21 USE oce 22 USE iom 15 USE oce ! ocean variables 16 USE dom_oce ! ocean domain 17 USE phycst ! physical constants 18 USE ice_2 ! LIM-2 variables 19 USE tradmp ! traceur damping 20 USE in_out_manager ! I/O manager 21 USE iom ! 23 22 24 23 IMPLICIT NONE … … 27 26 PUBLIC lim_dmp_2 ! called by ice_step_2 28 27 29 INTEGER :: nice1, nice2 , &! first and second record used30 & inumice_dmp! logical unit for ice variables (damping)31 REAL(wp), DIMENSION(jpi,jpj) :: hicif_dta , &! ice thickness at a given time32 & frld_dta! fraction lead at a given time33 REAL(wp), DIMENSION(jpi,jpj,2) :: hicif_data , &! ice thickness data at two consecutive times34 & frld_data! fraction lead data at two consecutive times28 INTEGER :: nice1, nice2 ! first and second record used 29 INTEGER :: inumice_dmp ! logical unit for ice variables (damping) 30 REAL(wp), DIMENSION(jpi,jpj) :: hicif_dta ! ice thickness at a given time 31 REAL(wp), DIMENSION(jpi,jpj) :: frld_dta ! fraction lead at a given time 32 REAL(wp), DIMENSION(jpi,jpj,2) :: hicif_data ! ice thickness data at two consecutive times 33 REAL(wp), DIMENSION(jpi,jpj,2) :: frld_data ! fraction lead data at two consecutive times 35 34 36 35 !! * Substitution 37 36 # include "vectopt_loop_substitute.h90" 38 37 !!---------------------------------------------------------------------- 39 !! LIM 2.0 , UCL-LOCEAN-IPSL (2006)38 !! NEMO/LIM 3.3, UCL-LOCEAN-IPSL (2010) 40 39 !! $Id$ 41 40 !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) … … 44 43 CONTAINS 45 44 46 SUBROUTINE lim_dmp_2( kt)45 SUBROUTINE lim_dmp_2( kt ) 47 46 !!------------------------------------------------------------------- 48 47 !! *** ROUTINE lim_dmp_2 *** … … 53 52 !! ** method : the key_tradmp must be used to compute resto(:,:) coef. 54 53 !!--------------------------------------------------------------------- 55 INTEGER, INTENT(in) :: kt 56 ! 57 INTEGER :: ji, jj! dummy loop indices54 INTEGER, INTENT(in) :: kt ! ocean time-step 55 !! 56 INTEGER :: ji, jj ! dummy loop indices 58 57 !!--------------------------------------------------------------------- 59 58 ! 60 59 CALL dta_lim_2( kt ) 61 60 ! 62 61 DO jj = 2, jpjm1 63 62 DO ji = fs_2, fs_jpim1 ! vector opt. 64 63 hicif(ji,jj) = hicif(ji,jj) - rdt_ice * resto(ji,jj,1) * ( hicif(ji,jj) - hicif_dta(ji,jj) ) 65 frld(ji,jj) = frld (ji,jj) - rdt_ice * resto(ji,jj,1) * ( frld (ji,jj)- frld_dta (ji,jj) )64 frld(ji,jj) = frld (ji,jj) - rdt_ice * resto(ji,jj,1) * ( frld (ji,jj) - frld_dta (ji,jj) ) 66 65 END DO 67 66 END DO
Note: See TracChangeset
for help on using the changeset viewer.