Changeset 13909

Timestamp:

2020-11-30T10:28:01+01:00 (4 years ago)

Author:

vancop

Message:

pond routine easy infrastructure part

File:

• NEMO/branches/2020/SI3_martin_ponds/src/ICE/icethd_pnd.F90 (modified) (6 diffs)

NEMO/branches/2020/SI3_martin_ponds/src/ICE/icethd_pnd.F90

@@ -20 +20 @@
    USE ice1D          ! sea-ice: thermodynamics variables
    USE icetab         ! sea-ice: 1D <==> 2D transformation
+   USE sbc_ice        ! surface energy budget
    !
    USE in_out_manager ! I/O manager
+   USE iom            ! I/O manager library
    USE lib_mpp        ! MPP library
    USE lib_fortran    ! fortran utilities (glob_sum + no signed zero)
@@ -37 +39 @@
    INTEGER, PARAMETER ::   np_pndCST = 1   ! Constant ice pond scheme
    INTEGER, PARAMETER ::   np_pndLEV = 2   ! Level ice pond scheme
+   INTEGER, PARAMETER ::   np_pndTOPO = 3   ! Level ice pond scheme
+
+   REAL(wp), PARAMETER :: &   ! shared parameters for topographic melt ponds
+      zhi_min = 0.1_wp        , & ! minimum ice thickness with ponds (m) 
+      zTd     = 0.15_wp       , & ! temperature difference for freeze-up (C)
+      zvp_min = 1.e-4_wp          ! minimum pond volume (m)
+
+   !--------------------------------------------------------------------------
+   ! 
+   ! Pond volume per area budget diags
+   !  
+   ! The idea of diags is the volume of ponds per grid cell area is
+   !
+   ! dV/dt = mlt + drn + lid + rnf
+   ! mlt   = input from surface melting
+   ! drn   = drainage through brine network
+   ! lid   = lid growth & melt
+   ! rnf   = runoff (water directly removed out of surface melting + overflow)
+   !
+   ! In topo mode, the pond water lost because it is in the snow is not included in the budget
+   !
+   ! In level mode, all terms are incorporated
+
+   REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: & ! pond volume budget diagnostics
+      diag_dvpn_mlt,  &     !! meltwater pond volume input      [m/day]
+      diag_dvpn_drn,  &     !! pond volume lost by drainage     [m/day]
+      diag_dvpn_lid,  &     !! exchange with lid / refreezing   [m/day]
+      diag_dvpn_rnf         !! meltwater pond lost to runoff    [m/day]
+      
+   REAL(wp), ALLOCATABLE, DIMENSION(:) ::   & ! pond volume budget diagnostics (1d)
+      diag_dvpn_mlt_1d,  &  !! meltwater pond volume input      [m/day]
+      diag_dvpn_drn_1d,  &  !! pond volume lost by drainage     [m/day]
+      diag_dvpn_lid_1d,  &  !! exchange with lid / refreezing   [m/day]
+      diag_dvpn_rnf_1d      !! meltwater pond lost to runoff    [m/day]
+   !
 
    !!----------------------------------------------------------------------
@@ -46 +83 @@
 
    SUBROUTINE ice_thd_pnd
+
       !!-------------------------------------------------------------------
       !!               ***  ROUTINE ice_thd_pnd   ***
       !!               
       !! ** Purpose :   change melt pond fraction and thickness
-      !!                
-      !!-------------------------------------------------------------------
-      !
+      !!
+      !! Note: Melt ponds affect only radiative transfer for now
+      !!
+      !!       No heat, no salt.
+      !!       The melt water they carry is collected but 
+      !!       not removed from fw budget or released to the ocean
+      !!
+      !!       A wfx_pnd has been coded for diagnostic purposes
+      !!       It is not fully consistent yet.
+      !!
+      !!       The current diagnostic lacks a contribution from drainage
+      !!
+      !!-------------------------------------------------------------------
+      !!
+      
+      ALLOCATE( diag_dvpn_mlt(jpi,jpj), diag_dvpn_lid(jpi,jpj), diag_dvpn_drn(jpi,jpj), diag_dvpn_rnf(jpi,jpj) )
+      ALLOCATE( diag_dvpn_mlt_1d(jpij), diag_dvpn_lid_1d(jpij), diag_dvpn_drn_1d(jpij), diag_dvpn_rnf_1d(jpij) )
+
+      diag_dvpn_mlt(:,:) = 0._wp    ;   diag_dvpn_drn(:,:)  = 0._wp
+      diag_dvpn_lid(:,:) = 0._wp    ;   diag_dvpn_rnf(:,:)  = 0._wp
+      diag_dvpn_mlt_1d(:) = 0._wp   ;   diag_dvpn_drn_1d(:) = 0._wp
+      diag_dvpn_lid_1d(:) = 0._wp   ;   diag_dvpn_rnf_1d(:) = 0._wp
+
       SELECT CASE ( nice_pnd )
       !
-      CASE (np_pndCST)   ;   CALL pnd_CST    !==  Constant melt ponds  ==!
+      CASE (np_pndCST)   ;   CALL pnd_CST                              !==  Constant melt ponds  ==!
          !
-      CASE (np_pndLEV)   ;   CALL pnd_LEV    !==  Level ice melt ponds  ==!
+      CASE (np_pndLEV)   ;   CALL pnd_LEV                              !==  Level ice melt ponds  ==!
+         !
+      CASE (np_pndTOPO)  ;   CALL pnd_TOPO                             !==  Topographic melt ponds  ==!
          !
       END SELECT
       !
+
+      IF( iom_use('dvpn_mlt'  ) )   CALL iom_put( 'dvpn_mlt', diag_dvpn_mlt * 100._wp * 86400._wp ) ! input from melting
+      IF( iom_use('dvpn_lid'  ) )   CALL iom_put( 'dvpn_lid', diag_dvpn_lid * 100._wp * 86400._wp ) ! exchanges with lid
+      IF( iom_use('dvpn_drn'  ) )   CALL iom_put( 'dvpn_drn', diag_dvpn_drn * 100._wp * 86400._wp ) ! vertical drainage
+      IF( iom_use('dvpn_rnf'  ) )   CALL iom_put( 'dvpn_rnf', diag_dvpn_rnf * 100._wp * 86400._wp ) ! runoff + overflow
+
+      DEALLOCATE( diag_dvpn_mlt, diag_dvpn_lid, diag_dvpn_drn, diag_dvpn_rnf )
+      DEALLOCATE( diag_dvpn_mlt_1d, diag_dvpn_lid_1d, diag_dvpn_drn_1d, diag_dvpn_rnf_1d )
+      
    END SUBROUTINE ice_thd_pnd 
+
 
 
@@ -315 +385 @@
       INTEGER  ::   ios, ioptio   ! Local integer
       !!
-      NAMELIST/namthd_pnd/  ln_pnd, ln_pnd_LEV , rn_apnd_min, rn_apnd_max, &
+      NAMELIST/namthd_pnd/  ln_pnd, ln_pnd_LEV , rn_apnd_min, rn_apnd_max, rn_pnd_flush, &
          &                          ln_pnd_CST , rn_apnd, rn_hpnd,         &
+         &                          ln_pnd_TOPO,                           &
          &                          ln_pnd_lids, ln_pnd_alb
       !!-------------------------------------------------------------------
@@ -332 +403 @@
          WRITE(numout,*) '   Namelist namicethd_pnd:'
          WRITE(numout,*) '      Melt ponds activated or not                                 ln_pnd       = ', ln_pnd
+         WRITE(numout,*) '         Topographic melt pond scheme                             ln_pnd_TOPO  = ', ln_pnd_TOPO
          WRITE(numout,*) '         Level ice melt pond scheme                               ln_pnd_LEV   = ', ln_pnd_LEV
          WRITE(numout,*) '            Minimum ice fraction that contributes to melt ponds   rn_apnd_min  = ', rn_apnd_min
          WRITE(numout,*) '            Maximum ice fraction that contributes to melt ponds   rn_apnd_max  = ', rn_apnd_max
+         WRITE(numout,*) '            Pond flushing efficiency                              rn_pnd_flush = ', rn_pnd_flus
          WRITE(numout,*) '         Constant ice melt pond scheme                            ln_pnd_CST   = ', ln_pnd_CST
          WRITE(numout,*) '            Prescribed pond fraction                              rn_apnd      = ', rn_apnd
@@ -347 +420 @@
       IF( ln_pnd_CST  ) THEN   ;   ioptio = ioptio + 1   ;   nice_pnd = np_pndCST    ;   ENDIF
       IF( ln_pnd_LEV  ) THEN   ;   ioptio = ioptio + 1   ;   nice_pnd = np_pndLEV    ;   ENDIF
+      IF( ln_pnd_TOPO ) THEN   ;   ioptio = ioptio + 1   ;   nice_pnd = np_pndTOPO   ;   ENDIF
       IF( ioptio /= 1 )   &
-         & CALL ctl_stop( 'ice_thd_pnd_init: choose either none (ln_pnd=F) or only one pond scheme (ln_pnd_LEV or ln_pnd_CST)' )
+         & CALL ctl_stop( 'ice_thd_pnd_init: choose either none (ln_pnd=F) or only one pond scheme (ln_pnd_LEV, ln_pnd_CST or ln_pnd_TOPO)' )
       !
       SELECT CASE( nice_pnd )