Changeset 14627 for branches/UKMO/jonnywilliams-u-br854/NEMOGCM/NEMO/OPA_SRC/SBC/sbc_oce.F90

Timestamp:

2021-03-22T05:00:21+01:00 (3 years ago)

Author:

jonnywilliams

Message:

runs successfully on maui

File:

• branches/UKMO/jonnywilliams-u-br854/NEMOGCM/NEMO/OPA_SRC/SBC/sbc_oce.F90 (modified) (2 diffs)

branches/UKMO/jonnywilliams-u-br854/NEMOGCM/NEMO/OPA_SRC/SBC/sbc_oce.F90

@@ -153 +153 @@
    ! sbccpl namelist parameters associated with icesheet freshwater input code. Included here rather than in sbccpl.F90 to 
    ! avoid circular dependencies.
+   ! SM 24.05.2018 - added defn for rn_antarctica_total_fw_flux_roc, and comment for option for nn_coupled_iceshelf_fluxes to =3
    INTEGER, PUBLIC     ::   nn_coupled_iceshelf_fluxes     ! =0 : total freshwater input from iceberg calving and ice shelf basal melting 
                                                            ! taken from climatologies used (no action in coupling routines).
@@ -159 +160 @@
                                                            ! =2 :  specify constant freshwater inputs in this namelist to set the combined
                                                            ! magnitude of iceberg calving and iceshelf melting freshwater fluxes.
-   LOGICAL, PUBLIC     ::   ln_iceshelf_init_atmos         ! If true force ocean to initialise iceshelf masses from atmospheric values rather
-                                                           ! than values in ocean restart (applicable if nn_coupled_iceshelf_fluxes=1).
-   REAL(wp), PUBLIC    ::   rn_greenland_total_fw_flux    ! Constant total rate of freshwater input (kg/s) for Greenland (if nn_coupled_iceshelf_fluxes=2) 
-   REAL(wp), PUBLIC    ::   rn_greenland_calving_fraction  ! Set fraction of total freshwater flux for iceberg calving - remainder goes to iceshelf melting.
-   REAL(wp), PUBLIC    ::   rn_antarctica_total_fw_flux   ! Constant total rate of freshwater input (kg/s) for Antarctica (if nn_coupled_iceshelf_fluxes=2) 
-   REAL(wp), PUBLIC    ::   rn_antarctica_calving_fraction ! Set fraction of total freshwater flux for iceberg calving - remainder goes to iceshelf melting.
-   REAL(wp), PUBLIC    ::   rn_iceshelf_fluxes_tolerance   ! Absolute tolerance for detecting differences in icesheet masses. 
+                                             ! =3 : as for =2, but with a rate of change also specified in the namelist
+   LOGICAL, PUBLIC     ::   ln_iceshelf_init_atmos          ! If true force ocean to initialise iceshelf masses from atmospheric values rather
+                                                            ! than values in ocean restart (applicable if nn_coupled_iceshelf_fluxes=1).
+   REAL(wp), PUBLIC    ::   rn_greenland_total_fw_flux      ! Constant total rate of freshwater input (kg/s) for Greenland 
+                                                            ! (if nn_coupled_iceshelf_fluxes=2), or starting rate if nn_coupled_iceshelf_fluxes=3
+   REAL(wp), PUBLIC    ::   rn_greenland_calving_fraction   ! Set fraction of total freshwater flux for iceberg calving - remainder goes to iceshelf melting.
+   REAL(wp), PUBLIC    ::   rn_greenland_total_fw_flux_roc  ! As rn_antarctica_total_fw_flux_roc but for Greenland Ice Sheet
+   REAL(wp), PUBLIC    ::   rn_antarctica_total_fw_flux     ! Constant total rate of freshwater input (kg/s) for Antarctica if 
+                                                            ! nn_coupled_iceshelf_fluxes=2, or starting rate if nn_coupled_iceshelf_fluxes=3
+   REAL(wp), PUBLIC    ::   rn_antarctica_calving_fraction  ! Set fraction of total freshwater flux for iceberg calving - remainder goes to iceshelf melting.
+   REAL(wp), PUBLIC    ::   rn_iceshelf_fluxes_tolerance    ! Absolute tolerance for detecting differences in icesheet masses. 
+   REAL(wp), PUBLIC    ::   rn_antarctica_total_fw_flux_roc ! Rate of change for ice sheet->ocean freshwater flux (% of orig flux after 1 yr, 
+                                                            ! so 1.1 is 10% inc/yr, 1.0 is no change, -1.1 is a 10% decr/yr) if nn_coupled_iceshelf_fluxes=3 
 
    !! * Substitutions