Changeset 12449 for NEMO/branches/UKMO/NEMO_4.0_add_pond_lids_prints/src/ICE/icectl.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/icectl.F90 (modified) (2 diffs)

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

@@ -87 +87 @@
          pdiag_fs = glob_sum( 'icectl',                                                                     &
             &                  ( sfx_bri(:,:) + sfx_bog(:,:) + sfx_bom(:,:) + sfx_sum(:,:) + sfx_sni(:,:) +  &
-            &                    sfx_opw(:,:) + sfx_res(:,:) + sfx_dyn(:,:) + sfx_sub(:,:) + sfx_lam(:,:)    &
-            &                  ) *  e1e2t(:,:) ) * zconv 
+            &                    sfx_opw(:,:) + sfx_res(:,:) + sfx_dyn(:,:) + sfx_sub(:,:) + sfx_lam(:,:) +  &
+            &                    sfx_pnd(:,:) ) *  e1e2t(:,:) ) * zconv 
          !
          !                          ! heat flux
@@ -116 +116 @@
          zfs = glob_sum( 'icectl',                                                                       &
             &              ( sfx_bri(:,:) + sfx_bog(:,:) + sfx_bom(:,:) + sfx_sum(:,:) + sfx_sni(:,:) +  &
-            &                sfx_opw(:,:) + sfx_res(:,:) + sfx_dyn(:,:) + sfx_sub(:,:) + sfx_lam(:,:)    & 
-            &              ) * e1e2t(:,:) ) * zconv - pdiag_fs
+            &                sfx_opw(:,:) + sfx_res(:,:) + sfx_dyn(:,:) + sfx_sub(:,:) + sfx_lam(:,:) +  & 
+            &                sfx_pnd(:,:) ) * e1e2t(:,:) ) * zconv - pdiag_fs
 
          ! heat flux