- Timestamp:
- 2015-10-22T10:17:51+02:00 (9 years ago)
- Location:
- branches/UKMO/icebergs_latent_heat/NEMOGCM
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/icebergs_latent_heat/NEMOGCM/CONFIG/SHARED/field_def.xml
r5517 r5821 563 563 <!-- variables available with iceberg trajectories --> 564 564 565 <field_group id="icbvar" domain_ref="grid_T" > 565 <field_group id="icbvar" grid_ref="grid_T_2D" > 566 <field id="berg_total_melt" long_name="icb melt rate 2 of icebergs" unit="kg/m2/s" /> 567 <field id="berg_total_heat_flux" long_name="icb latent heat of melting of icebergs" unit="kg/m2/s" /> 566 568 <field id="berg_melt" long_name="icb melt rate of icebergs" unit="kg/m2/s" /> 567 569 <field id="berg_buoy_melt" long_name="icb buoyancy component of iceberg melt rate" unit="kg/m2/s" /> -
branches/UKMO/icebergs_latent_heat/NEMOGCM/NEMO/OPA_SRC/ICB/icbdia.F90
r5808 r5821 371 371 IF( .NOT. ln_bergdia ) RETURN !!gm useless iom will control whether it is output or not 372 372 ! 373 CALL iom_put( "berg_total_melt" , berg_grid%floating_melt(:,:) ) ! Total melt flux to ocean [kg/m2/s] 374 CALL iom_put( "berg_total_heat_flux" , berg_grid%calving_hflx(:,:) ) ! Total iceberg-ocean heat flux [W/m2] 373 375 CALL iom_put( "berg_melt" , berg_melt (:,:) ) ! Melt rate of icebergs [kg/m2/s] 374 376 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
r5808 r5821 18 18 USE dom_oce ! NEMO domain 19 19 USE in_out_manager ! NEMO IO routines, numout in particular 20 USE iom 20 21 USE lib_mpp ! NEMO MPI routines, ctl_stop in particular 21 22 USE phycst ! NEMO physical constants … … 160 161 zmelt = ( zdM - ( zdMbitsE - zdMbitsM ) ) * z1_dt ! kg/s 161 162 berg_grid%floating_melt(ii,ij) = berg_grid%floating_melt(ii,ij) + zmelt * z1_e1e2 ! kg/m2/s 162 zheat = zmelt * pt%heat_density ! kg/s x J/kg = J/s 163 ! zheat = zmelt * pt%heat_density ! kg/s x J/kg = J/s 164 zheat = zmelt * lfus !rma kg/s x J/kg (latent heat of fusion) = J/s 163 165 berg_grid%calving_hflx (ii,ij) = berg_grid%calving_hflx (ii,ij) + zheat * z1_e1e2 ! W/m2 164 166 CALL icb_dia_melt( ii, ij, zMnew, zheat, this%mass_scaling, & … … 208 210 IF(.NOT. ln_passive_mode ) THEN 209 211 emp (:,:) = emp (:,:) - berg_grid%floating_melt(:,:) 210 !! qns (:,:) = qns (:,:) + berg_grid%calving_hflx (:,:) !!gm heat flux not yet properly coded ==>> need it, SOLVE that! 212 qns (:,:) = qns (:,:) - berg_grid%calving_hflx (:,:) 211 213 ENDIF 212 214 !
Note: See TracChangeset
for help on using the changeset viewer.