Changeset 12576 for NEMO/branches/UKMO/NEMO_4.0.1_icesheet_and_river_coupling/src/OCE/ICB/icbclv.F90

Timestamp:

2020-03-20T12:39:37+01:00 (4 years ago)

Author:

dancopsey

Message:

Merge in iceberg calving stuff from dev_r5518_coupling_GSI7_GSI8_landice from the start of the branch to revision 6023.

File:

• NEMO/branches/UKMO/NEMO_4.0.1_icesheet_and_river_coupling/src/OCE/ICB/icbclv.F90 (modified) (3 diffs)

NEMO/branches/UKMO/NEMO_4.0.1_icesheet_and_river_coupling/src/OCE/ICB/icbclv.F90

@@ -26 +26 @@
    USE icb_oce        ! iceberg parameters 
 
+   USE sbc_oce        ! for icesheet freshwater input variables 
+   USE in_out_manager 
+   USE iom 
+
    IMPLICIT NONE
    PRIVATE
@@ -49 +53 @@
       !
       REAL(wp)      ::   zcalving_used, zdist, zfact
+      REAL(wp)      ::   zgreenland_calving_sum, zantarctica_calving_sum 
       INTEGER       ::   jn, ji, jj                    ! loop counters
       INTEGER       ::   imx                           ! temporary integer for max berg class
@@ -63 +68 @@
       ! Heat in units of W/m2, and mask (just in case)
       berg_grid%calving_hflx(:,:) = src_calving_hflx(:,:) * tmask_i(:,:) * tmask(:,:,1)
+
+      IF( lk_oasis) THEN
+        ! nn_coupled_iceshelf_fluxes uninitialised unless lk_oasis=true
+        IF( nn_coupled_iceshelf_fluxes .gt. 0 ) THEN
+          ll_write = ((MOD( kt, sn_cfctl%ptimincr ) == 0) .OR. ( kt == nitend )) .AND. lwp
+          ! Adjust total calving rates so that sum of iceberg calving and iceshelf melting in the northern
+          ! and southern hemispheres equals rate of increase of mass of greenland and antarctic ice sheets
+          ! to preserve total freshwater conservation in coupled models without an active ice sheet model.
+
+           zgreenland_calving_sum = SUM( berg_grid%calving(:,:) * greenland_icesheet_mask(:,:) )
+           IF( lk_mpp ) CALL mpp_sum( zgreenland_calving_sum )
+           WHERE( greenland_icesheet_mask(:,:) == 1.0 )                                                                                 &
+          &    berg_grid%calving(:,:) = berg_grid%calving(:,:) * greenland_icesheet_mass_rate_of_change * rn_greenland_calving_fraction &
+          &                                     / ( zgreenland_calving_sum + 1.0e-10_wp )
+
+           ! check
+           IF(ll_write) WRITE(numout, *) 'Greenland iceberg calving climatology (kg/s) : ',zgreenland_calving_sum
+           zgreenland_calving_sum = SUM( berg_grid%calving(:,:) * greenland_icesheet_mask(:,:) )
+           IF( lk_mpp ) CALL mpp_sum( zgreenland_calving_sum )
+           IF(ll_write) WRITE(numout, *) 'Greenland iceberg calving adjusted value (kg/s) : ',zgreenland_calving_sum
+
+           zantarctica_calving_sum = SUM( berg_grid%calving(:,:) * antarctica_icesheet_mask(:,:) )
+           IF( lk_mpp ) CALL mpp_sum( zantarctica_calving_sum )
+           WHERE( antarctica_icesheet_mask(:,:) == 1.0 )                                                                              &
+           berg_grid%calving(:,:) = berg_grid%calving(:,:) * antarctica_icesheet_mass_rate_of_change * rn_antarctica_calving_fraction &
+          &                           / ( zantarctica_calving_sum + 1.0e-10_wp )
+
+           ! check
+           IF(ll_write) WRITE(numout, *) 'Antarctica iceberg calving climatology (kg/s) : ',zantarctica_calving_sum
+           zantarctica_calving_sum = SUM( berg_grid%calving(:,:) * antarctica_icesheet_mask(:,:) )
+           IF( lk_mpp ) CALL mpp_sum( zantarctica_calving_sum )
+           IF(ll_write) WRITE(numout, *) 'Antarctica iceberg calving adjusted value (kg/s) : ',zantarctica_calving_sum
+
+        ENDIF
+      ENDIF
+   
+      CALL iom_put( 'berg_calve', berg_grid%calving(:,:) )
+
 
       IF( ll_first_call .AND. .NOT. l_restarted_bergs ) THEN      ! This is a hack to simplify initialization