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 9885 – NEMO

Changeset 9885


Ignore:
Timestamp:
2018-07-05T18:45:20+02:00 (6 years ago)
Author:
clem
Message:

add melt ponds in BDY (very roughly)

Location:
NEMO/trunk/src
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/src/ICE/icedyn_adv.F90

    r9880 r9885  
    9999      !---------------------------- 
    100100      ! clem: The 2 advection schemes above are not strictly positive. 
    101       !       In Prather, advected fields are bounded by 0 in the routine with a MAX(0,field) ==> likely conservation issues 
     101      !       In Prather, advected fields are bounded by 0 (not anymore?) in the routine with a MAX(0,field) ==> likely conservation issues 
    102102      !       In UMx    , advected fields are not bounded and negative values can appear. 
    103103      !                   These values are usually very small but in some occasions they can also be non-negligible 
  • NEMO/trunk/src/OCE/BDY/bdyice.F90

    r9880 r9885  
    7474      ! 
    7575      CALL ice_cor( kt , 0 )      ! -- In case categories are out of bounds, do a remapping 
    76       !                           !    i.e. inputs have not the same ice thickness distribution  
    77       !                           !    (set by rn_himean) than the regional simulation 
     76      !                           !    i.e. inputs have not the same ice thickness distribution (set by rn_himean) 
     77      !                           !        than the regional simulation 
    7878      CALL ice_var_agg(1) 
    7979      ! 
     
    189189               END DO 
    190190               ! 
     191               ! Ice ponds 
     192               a_ip(ji,jj,jl) = 0._wp 
     193               v_ip(ji,jj,jl) = 0._wp 
     194               ! 
    191195            CASE( 1 )   ! velocity is outward 
    192196               ! 
     
    202206                  sz_i(ji,jj,jk,jl) = rswitch * sz_i(ii,ij,jk,jl) + ( 1.0 - rswitch ) * rn_simin 
    203207               END DO 
     208               ! 
     209               ! Ice ponds 
     210               a_ip(ji,jj,jl) = rswitch * a_ip(ii,ij,jl) 
     211               v_ip(ji,jj,jl) = rswitch * v_ip(ii,ij,jl) 
    204212               ! 
    205213            END SELECT 
     
    250258         END DO 
    251259         ! 
     260         CALL lbc_bdy_lnk( a_ip(:,:,jl), 'T', 1., ib_bdy ) 
     261         CALL lbc_bdy_lnk( v_ip(:,:,jl), 'T', 1., ib_bdy ) 
     262         ! 
    252263      END DO !jl 
    253       ! 
    254 !!      ! --- In case categories are out of bounds, do a remapping --- ! 
    255 !!      !     i.e. inputs have not the same ice thickness distribution  
    256 !!      !          (set by rn_himean) than the regional simulation 
    257 !!      IF( jpl > 1 )   CALL ice_itd_reb( kt ) 
    258264      !       
    259265   END SUBROUTINE bdy_ice_frs 
Note: See TracChangeset for help on using the changeset viewer.