source: NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/ISF/isfparmlt.F90 @ 11987

Last change on this file since 11987 was 11987, checked in by mathiot, 11 months ago

ENHANCE-02_ISF_nemo: changes needed after Dave's review

File size: 11.2 KB
RevLine 
[11395]1MODULE isfparmlt
2   !!======================================================================
[11403]3   !!                       ***  MODULE  isfparmlt  ***
[11541]4   !! Ice shelf parametrisation module :  update surface ocean boundary condition under ice
5   !!                   shelf using an ice shelf melt parametrisation
[11395]6   !!======================================================================
7   !! History :  4.0  ! original code
8   !!----------------------------------------------------------------------
9
[11987]10   USE isf_oce                  ! ice shelf
[11852]11   USE isftbl , ONLY: isf_tbl   ! ice shelf depth average
[11395]12
[11852]13   USE dom_oce                  ! ocean space and time domain
14   USE oce    , ONLY: tsn       ! ocean dynamics and tracers
15   USE phycst , ONLY: rcp, rau0 ! physical constants
16   USE eosbn2 , ONLY: eos_fzp   ! equation of state
[11395]17
[11852]18   USE in_out_manager              ! I/O manager
19   USE iom        , ONLY: iom_put  ! I/O library
[11931]20   USE fldread    , ONLY: fld_read, FLD, FLD_N !
[11852]21   USE lib_fortran, ONLY: glob_sum !
22   USE lib_mpp    , ONLY: ctl_stop !
23
[11395]24   IMPLICIT NONE
25
26   PRIVATE
27
28   PUBLIC  isfpar_mlt 
29   
30   !!----------------------------------------------------------------------
31   !! NEMO/OCE 4.0 , NEMO Consortium (2018)
32   !! $Id: sbcisf.F90 10536 2019-01-16 19:21:09Z mathiot $
33   !! Software governed by the CeCILL license (see ./LICENSE)
34   !!----------------------------------------------------------------------
35CONTAINS
36
37! -------------------------------------------------------------------------------------------------------
38! -------------------------------- PUBLIC SUBROUTINE ----------------------------------------------------
39! -------------------------------------------------------------------------------------------------------
40
[11494]41  SUBROUTINE isfpar_mlt( kt, pqhc, pqoce, pqfwf )
[11395]42      !!---------------------------------------------------------------------
[11494]43      !!                  ***  ROUTINE isfpar_mlt  ***
[11395]44      !!
45      !! ** Purpose : Compute Salt and Heat fluxes related to ice_shelf
46      !!              melting and freezing
47      !!
[11403]48      !! ** Method  :  2 parameterizations are available according
49      !!                        1 : Specified melt flux
[11395]50      !!                        2 : Beckmann & Goose parameterization
[11931]51      !!
[11395]52      !!-------------------------- OUT -------------------------------------
[11494]53      REAL(wp), DIMENSION(jpi,jpj), INTENT(inout) :: pqfwf, pqoce, pqhc  ! fresh water, ice-ocean heat and heat content fluxes
[11395]54      !!-------------------------- IN  -------------------------------------
55      INTEGER, INTENT(in) ::   kt   ! ocean time step
56      !!---------------------------------------------------------------------
57      !
58      ! Choose among the available ice shelf parametrisation
59      SELECT CASE ( cn_isfpar_mlt )
60      CASE ( 'spe' )    ! specified runoff in depth (Mathiot et al., 2017 in preparation)
61         CALL isfpar_mlt_spe(kt, pqhc, pqoce, pqfwf)
62      CASE ( 'bg03' )    ! Beckmann and Goosse parametrisation
63         CALL isfpar_mlt_bg03(kt, pqhc, pqoce, pqfwf)
64      CASE ( 'oasis' )
[11403]65         CALL isfpar_mlt_oasis( kt, pqhc, pqoce, pqfwf)
[11395]66      CASE DEFAULT
67         CALL ctl_stop('STOP', 'unknown isf melt formulation : cn_isfpar (should not see this)')
68      END SELECT
69      !
70   END SUBROUTINE isfpar_mlt
71
72! -------------------------------------------------------------------------------------------------------
73! -------------------------------- PRIVATE SUBROUTINE ---------------------------------------------------
74! -------------------------------------------------------------------------------------------------------
75
[11494]76   SUBROUTINE isfpar_mlt_spe(kt, pqhc, pqoce, pqfwf)
[11395]77      !!---------------------------------------------------------------------
[11403]78      !!                  ***  ROUTINE isfpar_mlt_spe  ***
[11395]79      !!
80      !! ** Purpose : prescribed ice shelf melting in case ice shelf cavities are closed.
81      !!              data read into a forcing files.
82      !!
83      !!-------------------------- OUT -------------------------------------
[11494]84      REAL(wp), DIMENSION(jpi,jpj), INTENT(inout) :: pqhc, pqfwf, pqoce  ! fresh water and ice-ocean heat fluxes
[11395]85      !!-------------------------- IN  -------------------------------------
86      INTEGER,  INTENT(in) :: kt
87      !!--------------------------------------------------------------------
88      INTEGER :: jk
89      REAL(wp), DIMENSION(jpi,jpj,jpk)  :: ztfrz3d
90      REAL(wp), DIMENSION(jpi,jpj)      :: ztfrz
91      !!--------------------------------------------------------------------
92      !
[11403]93      ! 0. ------------Read specified runoff
[11852]94      CALL fld_read ( kt, 1, sf_isfpar_fwf   )
[11395]95      !
96      ! compute ptfrz
[11403]97      ! 1. ------------Mean freezing point
[11395]98      DO jk = 1,jpk
99         CALL eos_fzp(tsn(:,:,jk,jp_sal), ztfrz3d(:,:,jk), gdept_n(:,:,jk))
100      END DO
[11495]101      CALL isf_tbl(ztfrz3d, ztfrz, 'T', misfkt_par, rhisf_tbl_par, misfkb_par, rfrac_tbl_par )
[11395]102      !
103      pqfwf(:,:) = - sf_isfpar_fwf(1)%fnow(:,:,1)      ! fresh water flux from the isf (fwfisf <0 mean melting)
104      pqoce(:,:) =   pqfwf(:,:) * rLfusisf             ! ocean/ice shelf flux assume to be equal to latent heat flux
105      pqhc (:,:) =   pqfwf(:,:) * ztfrz(:,:) * rcp     ! heat content flux
[11931]106      !
[11876]107      CALL iom_put('isftfrz_par', ztfrz )
[11395]108      !
109   END SUBROUTINE isfpar_mlt_spe
110
111   SUBROUTINE isfpar_mlt_bg03(kt, pqhc, pqoce, pqfwf)
112      !!---------------------------------------------------------------------
[11403]113      !!                  ***  ROUTINE isfpar_mlt_bg03  ***
[11395]114      !!
[11403]115      !! ** Purpose : compute an estimate of ice shelf melting and
116      !!              latent, ocean-ice and heat content heat fluxes
117      !!              in case cavities are closed based on the far fields T and S properties.
[11395]118      !!
[11403]119      !! ** Method  : The ice shelf melt is computed as proportional to the differences between the
120      !!              mean temperature and mean freezing point in front of the ice shelf averaged
121      !!              over the ice shelf min ice shelf draft and max ice shelf draft and the freezing point
[11395]122      !!
123      !! ** Reference : Beckmann and Goosse (2003), "A parameterization of ice shelf-ocean
[11403]124      !!                interaction for climate models", Ocean Modelling 5(2003) 157-170.
[11395]125      !!----------------------------------------------------------------------
126      !!-------------------------- OUT -------------------------------------
[11494]127      REAL(wp), DIMENSION(jpi,jpj), INTENT(inout) :: pqhc, pqfwf, pqoce  ! fresh water and ice-ocean heat fluxes
[11395]128      !!-------------------------- IN  -------------------------------------
129      INTEGER,  INTENT(in) :: kt
130      !!--------------------------------------------------------------------
131      INTEGER :: jk
132      REAL(wp), DIMENSION(jpi,jpj,jpk) :: ztfrz3d        ! freezing point
133      REAL(wp), DIMENSION(jpi,jpj)     :: ztfrz          ! freezing point
134      REAL(wp), DIMENSION(jpi,jpj)     :: ztavg          ! temperature avg
135      !!----------------------------------------------------------------------
136      !
137      ! 0. ------------Mean freezing point
138      DO jk = 1,jpk
139         CALL eos_fzp(tsn(:,:,jk,jp_sal), ztfrz3d(:,:,jk), gdept_n(:,:,jk))
140      END DO
[11495]141      CALL isf_tbl(ztfrz3d, ztfrz, 'T', misfkt_par, rhisf_tbl_par, misfkb_par, rfrac_tbl_par )
[11395]142      !
143      ! 1. ------------Mean temperature
[11495]144      CALL isf_tbl(tsn(:,:,jk,jp_tem), ztavg, 'T', misfkt_par, rhisf_tbl_par, misfkb_par, rfrac_tbl_par )
[11395]145      !
146      ! 2. ------------Net heat flux and fresh water flux due to the ice shelf
[11403]147      pqoce(:,:) =   rau0 * rcp * rn_gammat0 * risfLeff(:,:) * e1t(:,:) * ( ztavg(:,:) - ztfrz(:,:) ) * r1_e1e2t(:,:)
[11541]148      pqfwf(:,:) = - pqoce(:,:) / rLfusisf             ! derived from the latent heat flux
[11395]149      pqhc (:,:) =   pqfwf(:,:) * ztfrz(:,:) * rcp     ! heat content flux
150      !
[11403]151      ! 3. ------------BG03 output
[11395]152      ! output ttbl
[11521]153      CALL iom_put('ttbl_par', ztavg(:,:) * mskisf_par(:,:) )
[11395]154      !
155      ! output thermal driving
[11521]156      CALL iom_put('isfthermald_par',( ztfrz(:,:) - ztavg(:,:) ) * mskisf_par(:,:))
[11931]157      !
158      ! output freezing point used to define the thermal driving and heat content fluxes
[11876]159      CALL iom_put('isftfrz_par', ztfrz )
[11395]160      !
161   END SUBROUTINE isfpar_mlt_bg03
162
[11403]163   SUBROUTINE isfpar_mlt_oasis(kt, pqhc , pqoce, pqfwf )
164      !!----------------------------------------------------------------------
[11541]165      !!                  ***  ROUTINE isfpar_mlt_oasis  ***
[11403]166      !!
167      !! ** Purpose    : scale the fwf read from input file by the total amount received by the sbccpl interface
168      !!
[11494]169      !! ** Purpose    : - read ice shelf melt from forcing file and scale it by the input file total amount => pattern
170      !!                 - compute total amount of fwf given by sbccpl (fwfisf_oasis)
[11403]171      !!                 - scale fwf and compute heat fluxes
172      !!
173      !!---------------------------------------------------------------------
174      !!-------------------------- OUT -------------------------------------
175      REAL(wp), DIMENSION(jpi,jpj), INTENT(  out) :: pqhc, pqoce, pqfwf  ! heat content, latent heat and fwf fluxes
176      !!-------------------------- IN  -------------------------------------
177      INTEGER                     , INTENT(in   ) :: kt                  ! current time step
178      !!--------------------------------------------------------------------
179      INTEGER                           :: jk                            ! loop index
[11541]180      REAL(wp)                          :: zfwf_fld, zfwf_oasis          ! total fwf in the forcing fields (pattern) and from the cpl interface (amount)
[11403]181      REAL(wp), DIMENSION(jpi,jpj)      :: ztfrz                         ! tbl freezing temperature
182      REAL(wp), DIMENSION(jpi,jpj)      :: zfwf                          ! 2d fwf map after scaling
183      REAL(wp), DIMENSION(jpi,jpj,jpk)  :: ztfrz3d
184      !!--------------------------------------------------------------------
185      !
186      ! 0. ------------Read specified runoff
[11852]187      CALL fld_read ( kt, 1, sf_isfpar_fwf   )
[11403]188      !
[11541]189      ! 1. ------------Mean freezing point (needed for heat content flux)
[11403]190      DO jk = 1,jpk
191         CALL eos_fzp(tsn(:,:,jk,jp_sal), ztfrz3d(:,:,jk), gdept_n(:,:,jk))
192      END DO
[11495]193      CALL isf_tbl(ztfrz3d, ztfrz, 'T', misfkt_par, rhisf_tbl_par, misfkb_par, rfrac_tbl_par )
[11403]194      !
[11541]195      ! 2. ------------Scale isf melt pattern with total amount from oasis
[11403]196      ! ice shelf 2d map
197      zfwf(:,:) = - sf_isfpar_fwf(1)%fnow(:,:,1)
198      !
199      ! compute glob sum from input file
[11425]200      ! (PM) should we consider delay sum as in fwb ? (it will offset by 1 time step if I understood well)
201      zfwf_fld = glob_sum('isfcav_mlt', e1e2t(:,:) * zfwf(:,:))
[11403]202      !
203      ! compute glob sum from atm->oce ice shelf fwf
[11425]204      ! (PM) should we consider delay sum as in fwb ?
205      zfwf_oasis = glob_sum('isfcav_mlt', e1e2t(:,:) * fwfisf_oasis(:,:))
[11403]206      !
207      ! scale fwf
[11423]208      zfwf(:,:) = zfwf(:,:) * zfwf_oasis / zfwf_fld
[11403]209      !
[11931]210      ! 3. -----------Define fwf and qoce
[11403]211      ! ocean heat flux is assume to be equal to the latent heat
212      pqfwf(:,:) =   zfwf(:,:)                         ! fwf                ( >0 out )
213      pqoce(:,:) = - pqfwf(:,:) * rLfusisf             ! ocean heat flux    ( >0 out ) (assumed to be the latent heat flux)
214      pqhc (:,:) =   pqfwf(:,:) * ztfrz(:,:) * rcp     ! heat content flux  ( >0 out )
215      !
[11876]216      CALL iom_put('isftfrz_par', ztfrz )
217      !
[11395]218   END SUBROUTINE isfpar_mlt_oasis
219
220END MODULE isfparmlt
Note: See TracBrowser for help on using the repository browser.