Ignore:
Timestamp:
2014-12-15T17:42:49+01:00 (6 years ago)
Author:
timgraham
Message:

Merged branches/2014/dev_MERGE_2014 back onto the trunk as follows:

In the working copy of branch ran:
svn merge svn+ssh://forge.ipsl.jussieu.fr/ipsl/forge/projets/nemo/svn/trunk@HEAD
1 conflict in LIM_SRC_3/limdiahsb.F90
Resolved by keeping the version from dev_MERGE_2014 branch
and commited at r4989

In working copy run:
svn switch svn+ssh://forge.ipsl.jussieu.fr/ipsl/forge/projets/nemo/svn/trunk
to switch working copy

Run:
svn merge —reintegrate svn+ssh://forge.ipsl.jussieu.fr/ipsl/forge/projets/nemo/svn/branches/2014/dev_MERGE_2014
to merge the branch into the trunk - no conflicts at this stage.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/OPA_SRC/BDY/bdyice_lim.F90

    r4767 r4990  
    4545   PUBLIC   bdy_ice_lim_dyn ! routine called in limrhg 
    4646 
    47    REAL(wp) ::   epsi20 = 1.e-20_wp  ! module constants 
    48    REAL(wp) ::   epsi10 = 1.e-10_wp  ! min area allowed by ice model 
    4947   !!---------------------------------------------------------------------- 
    5048   !! NEMO/OPA 3.3 , NEMO Consortium (2010) 
     
    9896      INTEGER  ::   ji, jj, ii, ij     ! local scalar 
    9997      REAL(wp) ::   zwgt, zwgt1        ! local scalar 
    100       REAL(wp) ::   zinda, ztmelts, zdh 
     98      REAL(wp) ::   ztmelts, zdh 
    10199 
    102100      !!------------------------------------------------------------------------------ 
     
    115113         hicif(ji,jj) = ( hicif(ji,jj) * zwgt1 + dta%hicif(jb) * zwgt ) * tmask(ji,jj,1)     ! Ice depth  
    116114         hsnif(ji,jj) = ( hsnif(ji,jj) * zwgt1 + dta%hsnif(jb) * zwgt ) * tmask(ji,jj,1)     ! Snow depth 
    117  
    118 !         zinda = 1.0 - MAX( 0.0_wp , SIGN ( 1.0_wp , - frld(ji,jj) ) ) ! 0 if no ice 
    119 !         !------------------------------ 
    120 !         ! Sea ice surface temperature 
    121 !         !------------------------------ 
    122 !         sist(ji,jj)   = zinda * 270.0 + ( 1.0 - zinda ) * tfu(ji,jj) 
    123 !         !----------------------------------------------- 
    124 !         ! Ice/snow temperatures and energy stored in brines  
    125 !         !----------------------------------------------- 
    126 !         !!! TO BE CONTIUNED (as LIM3 below) !!! 
    127 !            zindhe = MAX( 0.e0, SIGN( 1.e0, fcor(1,jj) ) )              ! = 0 for SH, =1 for NH 
    128 !  
    129 !               ! Recover in situ values. 
    130 !               zindb         = MAX( rzero, SIGN( rone, zs0a(ji,jj) - epsi06 ) ) 
    131 !               zacrith       = 1.0 - ( zindhe * acrit(1) + ( 1.0 - zindhe ) * acrit(2) ) 
    132 !               zs0a (ji,jj)  = zindb * MIN( zs0a(ji,jj), zacrith ) 
    133 !               hsnif(ji,jj)  = zindb * ( zs0sn(ji,jj) /MAX( zs0a(ji,jj), epsi16 ) ) 
    134 !               hicif(ji,jj)  = zindb * ( zs0ice(ji,jj)/MAX( zs0a(ji,jj), epsi16 ) ) 
    135 !               zindsn        = MAX( rzero, SIGN( rone, hsnif(ji,jj) - epsi06 ) ) 
    136 !               zindic        = MAX( rzero, SIGN( rone, hicif(ji,jj) - epsi03 ) ) 
    137 !               zindb         = MAX( zindsn, zindic ) 
    138 !               zs0a (ji,jj)  = zindb * zs0a(ji,jj) 
    139 !               frld (ji,jj)  = 1.0 - zs0a(ji,jj) 
    140 !               hsnif(ji,jj)  = zindsn * hsnif(ji,jj) 
    141 !               hicif(ji,jj)  = zindic * hicif(ji,jj) 
    142 !               zusvosn       = 1.0/MAX( hsnif(ji,jj) * zs0a(ji,jj), epsi16 ) 
    143 !               zusvoic       = 1.0/MAX( hicif(ji,jj) * zs0a(ji,jj), epsi16 ) 
    144 !               zignm         = MAX( rzero,  SIGN( rone, hsndif - hsnif(ji,jj) ) ) 
    145 !               zrtt          = 173.15 * rone  
    146 !               ztsn          =          zignm   * tbif(ji,jj,1)  & 
    147 !                              + ( 1.0 - zignm ) * MIN( MAX( zrtt, rt0_snow * zusvosn * zs0c0(ji,jj)) , tfu(ji,jj) )  
    148 !               ztic1          = MIN( MAX( zrtt, rt0_ice * zusvoic * zs0c1(ji,jj) ) , tfu(ji,jj) ) 
    149 !               ztic2          = MIN( MAX( zrtt, rt0_ice * zusvoic * zs0c2(ji,jj) ) , tfu(ji,jj) ) 
    150 !  
    151 !               tbif(ji,jj,1) = zindsn * ztsn  + ( 1.0 - zindsn ) * tfu(ji,jj)                
    152 !               tbif(ji,jj,2) = zindic * ztic1 + ( 1.0 - zindic ) * tfu(ji,jj) 
    153 !               tbif(ji,jj,3) = zindic * ztic2 + ( 1.0 - zindic ) * tfu(ji,jj) 
    154 !               qstoif(ji,jj) = zindb  * xlic * zs0st(ji,jj) /  MAX( zs0a(ji,jj), epsi16 ) 
    155  
    156115      END DO  
    157116 
     
    215174            IF ( v_ice(ji  ,jj-1) > 0. .AND. vmask(ji  ,jj+1,1) == 0. ) jpbound = 1; ii = ji  ; ij = jj-1 
    216175 
    217             zinda = 1.0 - MAX( 0.0_wp , SIGN ( 1.0_wp , - at_i(ii,ij) + 0.01 ) ) ! 0 if no ice 
     176            rswitch = 1.0 - MAX( 0.0_wp , SIGN ( 1.0_wp , - at_i(ii,ij) + 0.01 ) ) ! 0 if no ice 
    218177 
    219178            ! concentration and thickness 
    220             a_i (ji,jj,jl) = a_i (ii,ij,jl) * zinda 
    221             ht_i(ji,jj,jl) = ht_i(ii,ij,jl) * zinda 
    222             ht_s(ji,jj,jl) = ht_s(ii,ij,jl) * zinda 
     179            a_i (ji,jj,jl) = a_i (ii,ij,jl) * rswitch 
     180            ht_i(ji,jj,jl) = ht_i(ii,ij,jl) * rswitch 
     181            ht_s(ji,jj,jl) = ht_s(ii,ij,jl) * rswitch 
    223182 
    224183            ! Ice and snow volumes 
     
    231190 
    232191               ! Ice salinity, age, temperature 
    233                sm_i(ji,jj,jl)   = zinda * rn_ice_sal(ib_bdy)  + ( 1.0 - zinda ) * s_i_min 
    234                o_i(ji,jj,jl)    = zinda * rn_ice_age(ib_bdy)  + ( 1.0 - zinda ) 
    235                t_su(ji,jj,jl)   = zinda * rn_ice_tem(ib_bdy)  + ( 1.0 - zinda ) * rn_ice_tem(ib_bdy) 
     192               sm_i(ji,jj,jl)   = rswitch * rn_ice_sal(ib_bdy)  + ( 1.0 - rswitch ) * s_i_min 
     193               o_i(ji,jj,jl)    = rswitch * rn_ice_age(ib_bdy)  + ( 1.0 - rswitch ) 
     194               t_su(ji,jj,jl)   = rswitch * rn_ice_tem(ib_bdy)  + ( 1.0 - rswitch ) * rn_ice_tem(ib_bdy) 
    236195               DO jk = 1, nlay_s 
    237                   t_s(ji,jj,jk,jl) = zinda * rn_ice_tem(ib_bdy) + ( 1.0 - zinda ) * rtt 
     196                  t_s(ji,jj,jk,jl) = rswitch * rn_ice_tem(ib_bdy) + ( 1.0 - rswitch ) * rtt 
    238197               END DO 
    239198               DO jk = 1, nlay_i 
    240                   t_i(ji,jj,jk,jl) = zinda * rn_ice_tem(ib_bdy) + ( 1.0 - zinda ) * rtt  
    241                   s_i(ji,jj,jk,jl) = zinda * rn_ice_sal(ib_bdy) + ( 1.0 - zinda ) * s_i_min 
     199                  t_i(ji,jj,jk,jl) = rswitch * rn_ice_tem(ib_bdy) + ( 1.0 - rswitch ) * rtt  
     200                  s_i(ji,jj,jk,jl) = rswitch * rn_ice_sal(ib_bdy) + ( 1.0 - rswitch ) * s_i_min 
    242201               END DO 
    243202                
     
    245204  
    246205               ! Ice salinity, age, temperature 
    247                sm_i(ji,jj,jl)   = zinda * sm_i(ii,ij,jl)  + ( 1.0 - zinda ) * s_i_min 
    248                o_i(ji,jj,jl)    = zinda * o_i(ii,ij,jl)   + ( 1.0 - zinda ) 
    249                t_su(ji,jj,jl)   = zinda * t_su(ii,ij,jl)  + ( 1.0 - zinda ) * rtt 
     206               sm_i(ji,jj,jl)   = rswitch * sm_i(ii,ij,jl)  + ( 1.0 - rswitch ) * s_i_min 
     207               o_i(ji,jj,jl)    = rswitch * o_i(ii,ij,jl)   + ( 1.0 - rswitch ) 
     208               t_su(ji,jj,jl)   = rswitch * t_su(ii,ij,jl)  + ( 1.0 - rswitch ) * rtt 
    250209               DO jk = 1, nlay_s 
    251                   t_s(ji,jj,jk,jl) = zinda * t_s(ii,ij,jk,jl) + ( 1.0 - zinda ) * rtt 
     210                  t_s(ji,jj,jk,jl) = rswitch * t_s(ii,ij,jk,jl) + ( 1.0 - rswitch ) * rtt 
    252211               END DO 
    253212               DO jk = 1, nlay_i 
    254                   t_i(ji,jj,jk,jl) = zinda * t_i(ii,ij,jk,jl) + ( 1.0 - zinda ) * rtt 
    255                   s_i(ji,jj,jk,jl) = zinda * s_i(ii,ij,jk,jl) + ( 1.0 - zinda ) * s_i_min 
     213                  t_i(ji,jj,jk,jl) = rswitch * t_i(ii,ij,jk,jl) + ( 1.0 - rswitch ) * rtt 
     214                  s_i(ji,jj,jk,jl) = rswitch * s_i(ii,ij,jk,jl) + ( 1.0 - rswitch ) * s_i_min 
    256215               END DO 
    257216 
     
    269228            DO jk = 1, nlay_s 
    270229               ! Snow energy of melting 
    271                e_s(ji,jj,jk,jl) = zinda * rhosn * ( cpic * ( rtt - t_s(ji,jj,jk,jl) ) + lfus ) 
     230               e_s(ji,jj,jk,jl) = rswitch * rhosn * ( cpic * ( rtt - t_s(ji,jj,jk,jl) ) + lfus ) 
    272231               ! Change dimensions 
    273232               e_s(ji,jj,jk,jl) = e_s(ji,jj,jk,jl) / unit_fac 
     
    278237               ztmelts          = - tmut * s_i(ji,jj,jk,jl) + rtt !Melting temperature in K                   
    279238               ! heat content per unit volume 
    280                e_i(ji,jj,jk,jl) = zinda * rhoic * & 
     239               e_i(ji,jj,jk,jl) = rswitch * rhoic * & 
    281240                  (   cpic    * ( ztmelts - t_i(ji,jj,jk,jl) ) & 
    282241                  +   lfus    * ( 1.0 - (ztmelts-rtt) / MIN((t_i(ji,jj,jk,jl)-rtt),-epsi20) ) & 
     
    335294      INTEGER  ::   ji, jj             ! local scalar 
    336295      INTEGER  ::   ib_bdy ! Loop index 
    337       REAL(wp) ::   zmsk1, zmsk2, zflag, zinda  
     296      REAL(wp) ::   zmsk1, zmsk2, zflag 
    338297     !!------------------------------------------------------------------------------ 
    339298      ! 
     
    375334                  ENDIF 
    376335                  ! mask ice velocities 
    377                   zinda = 1.0 - MAX( 0.0_wp , SIGN ( 1.0_wp , - at_i(ji,jj) + 0.01 ) ) ! 0 if no ice 
    378                   u_ice(ji,jj) = zinda * u_ice(ji,jj) 
     336                  rswitch = 1.0 - MAX( 0.0_wp , SIGN ( 1.0_wp , - at_i(ji,jj) + 0.01 ) ) ! 0 if no ice 
     337                  u_ice(ji,jj) = rswitch * u_ice(ji,jj) 
    379338                   
    380339               ENDDO 
     
    404363                  ENDIF 
    405364                  ! mask ice velocities 
    406                   zinda = 1.0 - MAX( 0.0_wp , SIGN ( 1.0_wp , - at_i(ji,jj) + 0.01 ) ) ! 0 if no ice 
    407                   v_ice(ji,jj) = zinda * v_ice(ji,jj) 
     365                  rswitch = 1.0 - MAX( 0.0_wp , SIGN ( 1.0_wp , - at_i(ji,jj) + 0.01 ) ) ! 0 if no ice 
     366                  v_ice(ji,jj) = rswitch * v_ice(ji,jj) 
    408367                   
    409368               ENDDO 
Note: See TracChangeset for help on using the changeset viewer.