Ignore:
Timestamp:
2008-04-11T19:05:03+02:00 (13 years ago)
Author:
ctlod
Message:

merge dev_001_SBC branche with the trunk to include the New Surface Module package, see ticket: #113

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/LIM_SRC_3/limthd_sal.F90

    r869 r888  
    11MODULE limthd_sal 
    2 #if defined key_lim3 
    32   !!---------------------------------------------------------------------- 
    43   !!   'key_lim3'                                      LIM3 sea-ice model 
     
    98   !!                               the ice 
    109   !!====================================================================== 
    11  
     10#if defined key_lim3 
    1211   !!---------------------------------------------------------------------- 
    1312   !!   lim_thd_sal : salinity variations in the ice 
     
    1615   USE phycst           ! physical constants (ocean directory) 
    1716   USE ice_oce          ! ice variables 
     17   USE sbc_oce          ! Surface boundary condition: ocean fields 
    1818   USE thd_ice 
    1919   USE iceini 
     
    4040   CONTAINS 
    4141 
    42    SUBROUTINE lim_thd_sal(kideb,kiut,jl) 
     42   SUBROUTINE lim_thd_sal(kideb,kiut) 
    4343      !!------------------------------------------------------------------- 
    4444      !!                ***  ROUTINE lim_thd_sal  ***        
     
    7676      !! * Local variables 
    7777      INTEGER, INTENT(in) :: & 
    78          kideb, kiut, jl         !: thickness category index 
     78         kideb, kiut             !: thickness category index 
    7979 
    8080      INTEGER ::             & 
     
    318318            zjj                 = ( npb(ji) - 1 ) / jpi + 1 
    319319            fseqv_1d(ji) = fseqv_1d(ji)              + &  
    320                            ( sss_io(zji,zjj) - bulk_sal    ) * &  
     320                           ( sss_m(zji,zjj) - bulk_sal    ) * &  
    321321                           rhoic * a_i_b(ji) * & 
    322322                           MAX( dh_i_bott(ji) , 0.0 ) / rdt_ice 
     
    327327            zjj                 = ( npb(ji) - 1 ) / jpi + 1 
    328328            fseqv_1d(ji) = fseqv_1d(ji)              + &  
    329                            ( sss_io(zji,zjj) - s_i_new(ji) ) * &  
     329                           ( sss_m(zji,zjj) - s_i_new(ji) ) * &  
    330330                             rhoic * a_i_b(ji) * & 
    331331                             MAX( dh_i_bott(ji) , 0.0 ) / rdt_ice 
Note: See TracChangeset for help on using the changeset viewer.