Changeset 5821 for branches/UKMO/icebergs_latent_heat/NEMOGCM/NEMO/OPA_SRC/ICB

Timestamp:

2015-10-22T10:17:51+02:00 (9 years ago)

Author:

davestorkey

Message:

UKMO/icebergs_latent_heat branch:

1. Include latent heat term.
2. New diagnostics.

Location:

branches/UKMO/icebergs_latent_heat/NEMOGCM/NEMO/OPA_SRC/ICB

Files:

• icbdia.F90 (modified) (1 diff)
• icbthm.F90 (modified) (3 diffs)

branches/UKMO/icebergs_latent_heat/NEMOGCM/NEMO/OPA_SRC/ICB/icbdia.F90

@@ -371 +371 @@
       IF( .NOT. ln_bergdia )   RETURN            !!gm useless iom will control whether it is output or not
       !
+      CALL iom_put( "berg_total_melt"      , berg_grid%floating_melt(:,:)   )   ! Total melt flux to ocean      [kg/m2/s]
+      CALL iom_put( "berg_total_heat_flux" , berg_grid%calving_hflx(:,:)    )   ! Total iceberg-ocean heat flux [W/m2]
       CALL iom_put( "berg_melt"        , berg_melt   (:,:)   )   ! Melt rate of icebergs                     [kg/m2/s]
       CALL iom_put( "berg_buoy_melt"   , buoy_melt   (:,:)   )   ! Buoyancy component of iceberg melt rate   [kg/m2/s]

branches/UKMO/icebergs_latent_heat/NEMOGCM/NEMO/OPA_SRC/ICB/icbthm.F90

@@ -18 +18 @@
    USE dom_oce        ! NEMO domain
    USE in_out_manager ! NEMO IO routines, numout in particular
+   USE iom
    USE lib_mpp        ! NEMO MPI routines, ctl_stop in particular
    USE phycst         ! NEMO physical constants
@@ -160 +161 @@
             zmelt    = ( zdM - ( zdMbitsE - zdMbitsM ) ) * z1_dt   ! kg/s
             berg_grid%floating_melt(ii,ij) = berg_grid%floating_melt(ii,ij) + zmelt    * z1_e1e2    ! kg/m2/s
-            zheat = zmelt * pt%heat_density              ! kg/s x J/kg = J/s
+!            zheat = zmelt * pt%heat_density              ! kg/s x J/kg = J/s
+            zheat = zmelt * lfus                           !rma kg/s x J/kg (latent heat of fusion) = J/s
             berg_grid%calving_hflx (ii,ij) = berg_grid%calving_hflx (ii,ij) + zheat    * z1_e1e2    ! W/m2
             CALL icb_dia_melt( ii, ij, zMnew, zheat, this%mass_scaling,       &
@@ -208 +210 @@
       IF(.NOT. ln_passive_mode ) THEN
          emp (:,:) = emp (:,:) - berg_grid%floating_melt(:,:)
-!!       qns (:,:) = qns (:,:) + berg_grid%calving_hflx (:,:)  !!gm heat flux not yet properly coded ==>> need it, SOLVE that!
+         qns (:,:) = qns (:,:) - berg_grid%calving_hflx (:,:)  
       ENDIF
       !