Changeset 9885

Timestamp:

2018-07-05T18:45:20+02:00 (5 years ago)

Author:

clem

Message:

add melt ponds in BDY (very roughly)

Location:

NEMO/trunk/src

Files:

• ICE/icedyn_adv.F90 (modified) (1 diff)
• OCE/BDY/bdyice.F90 (modified) (4 diffs)

NEMO/trunk/src/ICE/icedyn_adv.F90

@@ -99 +99 @@
       !----------------------------
       ! clem: The 2 advection schemes above are not strictly positive.
-      !       In Prather, advected fields are bounded by 0 in the routine with a MAX(0,field) ==> likely conservation issues
+      !       In Prather, advected fields are bounded by 0 (not anymore?) in the routine with a MAX(0,field) ==> likely conservation issues
       !       In UMx    , advected fields are not bounded and negative values can appear.
       !                   These values are usually very small but in some occasions they can also be non-negligible

NEMO/trunk/src/OCE/BDY/bdyice.F90

@@ -74 +74 @@
       !
       CALL ice_cor( kt , 0 )      ! -- In case categories are out of bounds, do a remapping
-      !                           !    i.e. inputs have not the same ice thickness distribution 
-      !                           !    (set by rn_himean) than the regional simulation
+      !                           !    i.e. inputs have not the same ice thickness distribution (set by rn_himean)
+      !                           !         than the regional simulation
       CALL ice_var_agg(1)
       !
@@ -189 +189 @@
                END DO
                !
+               ! Ice ponds
+               a_ip(ji,jj,jl) = 0._wp
+               v_ip(ji,jj,jl) = 0._wp
+               !
             CASE( 1 )   ! velocity is outward
                !
@@ -202 +206 @@
                   sz_i(ji,jj,jk,jl) = rswitch * sz_i(ii,ij,jk,jl) + ( 1.0 - rswitch ) * rn_simin
                END DO
+               !
+               ! Ice ponds
+               a_ip(ji,jj,jl) = rswitch * a_ip(ii,ij,jl)
+               v_ip(ji,jj,jl) = rswitch * v_ip(ii,ij,jl)
                !
             END SELECT
@@ -250 +258 @@
          END DO
          !
+         CALL lbc_bdy_lnk( a_ip(:,:,jl), 'T', 1., ib_bdy )
+         CALL lbc_bdy_lnk( v_ip(:,:,jl), 'T', 1., ib_bdy )
+         !
       END DO !jl
-      !
-!!      ! --- In case categories are out of bounds, do a remapping --- !
-!!      !     i.e. inputs have not the same ice thickness distribution 
-!!      !          (set by rn_himean) than the regional simulation
-!!      IF( jpl > 1 )   CALL ice_itd_reb( kt )
       !      
    END SUBROUTINE bdy_ice_frs