Changeset 12449 for NEMO/branches/UKMO/NEMO_4.0_add_pond_lids_prints/src/ICE/iceupdate.F90

Timestamp:

2020-02-24T18:04:40+01:00 (4 years ago)

Author:

dancopsey

Message:

• Snow into meltponds contributes to snow to ice diagnostics affecting thickness and salinity of the ice.
• Avoid divide by zero errors by:
  • Removing meltponds from small ice areas.
  • Removing divide by pond fractions apart from lid shrinking there a minimum pond fraction is imnposed
• Ponds leaking into the ocean have a salinity flux to the ocean (also a new diagnostic).
• Add vertical fluxing of pond water.

File:

• NEMO/branches/UKMO/NEMO_4.0_add_pond_lids_prints/src/ICE/iceupdate.F90 (modified) (3 diffs)

NEMO/branches/UKMO/NEMO_4.0_add_pond_lids_prints/src/ICE/iceupdate.F90

@@ -165 +165 @@
             !------------------------------------------
             sfx(ji,jj) = sfx_bog(ji,jj) + sfx_bom(ji,jj) + sfx_sum(ji,jj) + sfx_sni(ji,jj) + sfx_opw(ji,jj)   &
-               &       + sfx_res(ji,jj) + sfx_dyn(ji,jj) + sfx_bri(ji,jj) + sfx_sub(ji,jj) + sfx_lam(ji,jj)
+               &       + sfx_res(ji,jj) + sfx_dyn(ji,jj) + sfx_bri(ji,jj) + sfx_sub(ji,jj) + sfx_lam(ji,jj) + sfx_pnd(ji,jj)
             
             ! Mass of snow and ice per unit area   
@@ -197 +197 @@
       !
       ! --- salt fluxes [kg/m2/s] --- !
-      !                           ! sfxice =  sfxbog + sfxbom + sfxsum + sfxsni + sfxopw + sfxres + sfxdyn + sfxbri + sfxsub + sfxlam
+      !                           ! sfxice =  sfxbog + sfxbom + sfxsum + sfxsni + sfxopw + sfxres + sfxdyn + sfxbri + sfxsub + sfxlam + sfxpnd
       IF( iom_use('sfxice'  ) )   CALL iom_put( "sfxice", sfx     * 1.e-03 )   ! salt flux from total ice growth/melt
       IF( iom_use('sfxbog'  ) )   CALL iom_put( "sfxbog", sfx_bog * 1.e-03 )   ! salt flux from bottom growth
@@ -209 +209 @@
       IF( iom_use('sfxres'  ) )   CALL iom_put( "sfxres", sfx_res * 1.e-03 )   ! salt flux from undiagnosed processes
       IF( iom_use('sfxsub'  ) )   CALL iom_put( "sfxsub", sfx_sub * 1.e-03 )   ! salt flux from sublimation
+      IF( iom_use('sfxpnd'  ) )   CALL iom_put( "sfxpnd", sfx_pnd * 1.e-03 )   ! salt flux from melt ponds leaking into the ocean
 
       ! --- mass fluxes [kg/m2/s] --- !