New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 1855 for branches/DEV_r1837_mass_heat_salt_fluxes/NEMO/LIM_SRC_2/limdmp_2.F90 – NEMO

Ignore:
Timestamp:
2010-04-30T17:49:04+02:00 (14 years ago)
Author:
gm
Message:

ticket:#665 style change only, with the suppression of thd_ice_2 (merged in ice_2)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/DEV_r1837_mass_heat_salt_fluxes/NEMO/LIM_SRC_2/limdmp_2.F90

    r1715 r1855  
    44   !!  Ice model : restoring Ice thickness and Fraction leads 
    55   !!====================================================================== 
    6    !! History :   2.0  !  04-04 (S. Theetten) Original code 
     6   !! History :   2.0  ! 2004-04 (S. Theetten) Original code 
    77   !!---------------------------------------------------------------------- 
    88#if defined key_lim2   &&   defined key_tradmp 
     
    1111   !!   'key_tradmp'                                                Damping 
    1212   !!---------------------------------------------------------------------- 
     13   !!   lim_dmp_2       : ice model damping 
    1314   !!---------------------------------------------------------------------- 
    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              ! 
    2322    
    2423   IMPLICIT NONE 
     
    2726   PUBLIC   lim_dmp_2     ! called by ice_step_2 
    2827    
    29    INTEGER                        ::   nice1, nice2,  &  ! first and second record used 
    30       &                                inumice_dmp       ! logical unit for ice variables (damping) 
    31    REAL(wp), DIMENSION(jpi,jpj)   ::   hicif_dta  ,   &  ! ice thickness at a given time 
    32       &                                frld_dta          ! fraction lead at a given time 
    33    REAL(wp), DIMENSION(jpi,jpj,2) ::   hicif_data ,   &  ! ice thickness data at two consecutive times 
    34       &                                frld_data         ! fraction lead data at two consecutive times 
     28   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 
    3534 
    3635   !! * Substitution 
    3736#  include "vectopt_loop_substitute.h90" 
    3837   !!---------------------------------------------------------------------- 
    39    !!   LIM 2.0 , UCL-LOCEAN-IPSL  (2006) 
     38   !! NEMO/LIM 3.3,  UCL-LOCEAN-IPSL (2010)  
    4039   !! $Id$ 
    4140   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) 
     
    4443CONTAINS 
    4544 
    46    SUBROUTINE lim_dmp_2(kt) 
     45   SUBROUTINE lim_dmp_2( kt ) 
    4746      !!------------------------------------------------------------------- 
    4847      !!                   ***  ROUTINE lim_dmp_2 *** 
     
    5352      !! ** method  : the key_tradmp must be used to compute resto(:,:) coef. 
    5453      !!--------------------------------------------------------------------- 
    55       INTEGER, INTENT(in) ::   kt     ! ocean time-step 
    56       ! 
    57       INTEGER             ::   ji, jj         ! dummy loop indices 
     54      INTEGER, INTENT(in) ::   kt   ! ocean time-step 
     55      !! 
     56      INTEGER ::   ji, jj   ! dummy loop indices 
    5857      !!--------------------------------------------------------------------- 
    5958      ! 
    6059      CALL dta_lim_2( kt ) 
    61  
     60      ! 
    6261      DO jj = 2, jpjm1 
    6362         DO ji = fs_2, fs_jpim1   ! vector opt. 
    6463            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) )   
    6665         END DO 
    6766      END DO 
Note: See TracChangeset for help on using the changeset viewer.