- Timestamp:
- 2017-04-26T09:02:56+02:00 (7 years ago)
- Location:
- trunk/NEMOGCM
- Files:
-
- 4 edited
-
CONFIG/SHARED/field_def_nemo-opa.xml (modified) (4 diffs)
-
NEMO/OPA_SRC/SBC/sbccpl.F90 (modified) (12 diffs)
-
NEMO/OPA_SRC/SBC/sbcisf.F90 (modified) (1 diff)
-
NEMO/OPA_SRC/SBC/sbcrnf.F90 (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/CONFIG/SHARED/field_def_nemo-opa.xml
r7828 r7968 207 207 208 208 <!-- * variable related to ice shelf forcing * --> 209 <field id="fwfisf" long_name="Ice shelf melting" unit="Kg/m2/s" />209 <field id="fwfisf" long_name="Ice shelf melting" unit="kg/m2/s" /> 210 210 <field id="fwfisf3d" long_name="Ice shelf melting" unit="kg/m2/s" grid_ref="grid_T_3D" /> 211 211 <field id="qlatisf" long_name="Ice shelf latent heat flux" unit="W/m2" /> … … 215 215 <field id="isfgammat" long_name="transfert coefficient for isf (temperature) " unit="m/s" /> 216 216 <field id="isfgammas" long_name="transfert coefficient for isf (salinity) " unit="m/s" /> 217 217 <field id="stbl" long_name="salinity in the Losh tbl " unit="PSU" /> 218 218 <field id="ttbl" long_name="temperature in the Losh tbl " unit="C" /> 219 219 <field id="utbl" long_name="zonal current in the Losh tbl at T point " unit="m/s" /> … … 238 238 <field id="calving_cea" long_name="Calving" standard_name="water_flux_into_sea_water_from_icebergs" unit="kg/m2/s" /> 239 239 <field id="iceberg_cea" long_name="Iceberg" standard_name="water_flux_into_sea_water_from_icebergs" unit="kg/m2/s" /> 240 <field id="iceshelf_cea" long_name="Iceshelf" standard_name="water_flux_into_sea_water_from_iceshelf" unit="kg/m2/s" /> 240 <field id="iceshelf_cea" long_name="Iceshelf" standard_name="water_flux_into_sea_water_from_iceshelf" unit="kg/m2/s" /> 241 241 242 242 243 <!-- available if key_oasis3 + conservative method --> … … 253 254 <field id="hflx_ice_cea" long_name="heat flux due to ice thermodynamics" standard_name="heat_flux_into_sea_water_due_to_sea_ice_thermodynamics" unit="W/m2" /> 254 255 <field id="hflx_rnf_cea" long_name="heat flux due to runoffs" standard_name="temperature_flux_due_to_runoff_expressed_as_heat_flux_into_sea_water" unit="W/m2" /> 255 <field id="hflx_cal_cea" long_name="heat flux due to calving" standard_name="heat_flux_into_sea_water_due_to_iceberg_thermodynamics" unit="W/m2" /> 256 <field id="hflx_cal_cea" long_name="heat flux due to calving" standard_name="heat_flux_into_sea_water_due_to_calving" unit="W/m2" /> 257 <field id="hflx_icb_cea" long_name="heat flux due to iceberg" standard_name="heat_flux_into_sea_water_due_to_icebergs" unit="W/m2" /> 258 <field id="hflx_isf_cea" long_name="heat flux due to iceshelf" standard_name="heat_flux_into_sea_water_due_to_iceshelf" unit="W/m2" /> 256 259 <field id="bicemel_cea" long_name="Rate of Melt at Sea Ice Base (cell average)" standard_name="tendency_of_sea_ice_amount_due_to_basal_melting" unit="kg/m2/s" /> 257 260 <field id="licepro_cea" long_name="Lateral Sea Ice Growth Rate (cell average)" standard_name="tendency_of_sea_ice_amount_due_to_lateral_growth_of_ice_floes" unit="kg/m2/s" /> -
trunk/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90
r7952 r7968 1566 1566 !! qsr_tot = pfrld * qsr_oce + ( 1 - pfrld ) * qsr_ice => provided 1567 1567 !! 1568 !! emp_tot = emp_oce + emp_ice => calving is provided and added to emp_tot (and emp_oce) 1569 !! river runoff (rnf) is provided but not included here 1570 !! 1568 !! emp_tot = emp_oce + emp_ice => calving is provided and added to emp_tot (and emp_oce). 1569 !! runoff (which includes rivers+icebergs) and iceshelf 1570 !! are provided but not included in emp here. Only runoff will 1571 !! be included in emp in other parts of NEMO code 1571 1572 !! ** Action : update at each nf_ice time step: 1572 1573 !! qns_tot, qsr_tot non-solar and solar total heat fluxes … … 1584 1585 ! 1585 1586 INTEGER :: jl ! dummy loop index 1586 REAL(wp), POINTER, DIMENSION(:,: ) :: zcptn, z tmp, zcptrain, zcptsnw, zicefr, zmsk, zsnw1587 REAL(wp), POINTER, DIMENSION(:,: ) :: zcptn, zcptrain, zcptsnw, zicefr, zmsk, zsnw 1587 1588 REAL(wp), POINTER, DIMENSION(:,: ) :: zemp_tot, zemp_ice, zemp_oce, ztprecip, zsprecip, zevap_oce, zevap_ice, zdevap_ice 1588 1589 REAL(wp), POINTER, DIMENSION(:,: ) :: zqns_tot, zqns_oce, zqsr_tot, zqsr_oce, zqprec_ice, zqemp_oce, zqemp_ice … … 1592 1593 IF( nn_timing == 1 ) CALL timing_start('sbc_cpl_ice_flx') 1593 1594 ! 1594 CALL wrk_alloc( jpi,jpj, zcptn, z tmp, zcptrain, zcptsnw, zicefr, zmsk, zsnw )1595 CALL wrk_alloc( jpi,jpj, zcptn, zcptrain, zcptsnw, zicefr, zmsk, zsnw ) 1595 1596 CALL wrk_alloc( jpi,jpj, zemp_tot, zemp_ice, zemp_oce, ztprecip, zsprecip, zevap_oce, zevap_ice, zdevap_ice ) 1596 1597 CALL wrk_alloc( jpi,jpj, zqns_tot, zqns_oce, zqsr_tot, zqsr_oce, zqprec_ice, zqemp_oce, zqemp_ice ) … … 1615 1616 zemp_tot(:,:) = frcv(jpr_tevp)%z3(:,:,1) - ztprecip(:,:) 1616 1617 zemp_ice(:,:) = ( frcv(jpr_ievp)%z3(:,:,1) - frcv(jpr_snow)%z3(:,:,1) ) * zicefr(:,:) 1617 IF( iom_use('precip') ) &1618 & CALL iom_put( 'precip' , frcv(jpr_rain)%z3(:,:,1) + frcv(jpr_snow)%z3(:,:,1) ) ! total precipitation1619 IF( iom_use('rain') ) &1620 & CALL iom_put( 'rain' , frcv(jpr_rain)%z3(:,:,1) ) ! liquid precipitation1621 IF( iom_use('rain_ao_cea') ) &1622 & CALL iom_put( 'rain_ao_cea' , frcv(jpr_rain)%z3(:,:,1)* p_frld(:,:) * tmask(:,:,1) ) ! liquid precipitation1623 IF( iom_use('hflx_rain_cea') ) &1624 CALL iom_put( 'hflx_rain_cea', frcv(jpr_rain)%z3(:,:,1) * zcptn(:,:) * tmask(:,:,1)) ! heat flux from liq. precip.1625 IF( iom_use('hflx_prec_cea') ) &1626 CALL iom_put( 'hflx_prec_cea', ztprecip * zcptn(:,:) * tmask(:,:,1) * p_frld(:,:) ) ! heat content flux from all precip (cell avg)1627 IF( iom_use('evap_ao_cea') .OR. iom_use('hflx_evap_cea') ) &1628 ztmp(:,:) = frcv(jpr_tevp)%z3(:,:,1) - frcv(jpr_ievp)%z3(:,:,1) * zicefr(:,:)1629 IF( iom_use('evap_ao_cea' ) ) &1630 CALL iom_put( 'evap_ao_cea' , ztmp * tmask(:,:,1) ) ! ice-free oce evap (cell average)1631 IF( iom_use('hflx_evap_cea') ) &1632 CALL iom_put( 'hflx_evap_cea', ztmp(:,:) * zcptn(:,:) * tmask(:,:,1) ) ! heat flux from from evap (cell average)1633 1618 CASE( 'oce and ice' ) ! received fields: jpr_sbpr, jpr_semp, jpr_oemp, jpr_ievp 1634 1619 zemp_tot(:,:) = p_frld(:,:) * frcv(jpr_oemp)%z3(:,:,1) + zicefr(:,:) * frcv(jpr_sbpr)%z3(:,:,1) … … 1639 1624 1640 1625 #if defined key_lim3 1641 ! zsnw = snow fraction over ice after wind blowing 1626 ! zsnw = snow fraction over ice after wind blowing (=zicefr if no blowing) 1642 1627 zsnw(:,:) = 0._wp ; CALL lim_thd_snwblow( p_frld, zsnw ) 1643 1628 … … 1652 1637 zevap_ice(:,:) = frcv(jpr_ievp)%z3(:,:,1) 1653 1638 ! since the sensitivity of evap to temperature (devap/dT) is not prescribed by the atmosphere, we set it to 0 1654 ! therefore, sublimation is not redistributed over the ice categories in caseno subgrid scale fluxes are provided by atm.1639 ! therefore, sublimation is not redistributed over the ice categories when no subgrid scale fluxes are provided by atm. 1655 1640 zdevap_ice(:,:) = 0._wp 1656 1641 1657 ! --- runoffs (included in emp later on)--- !1658 IF( srcv(jpr_rnf)%laction ) rnf(:,:) = frcv(jpr_rnf)%z3(:,:,1)1659 1660 ! --- calving (put in emp_tot and emp_oce) --- !1661 IF( srcv(jpr_cal)%laction ) THEN 1642 ! --- Continental fluxes --- ! 1643 IF( srcv(jpr_rnf)%laction ) THEN ! runoffs (included in emp later on) 1644 rnf(:,:) = frcv(jpr_rnf)%z3(:,:,1) 1645 ENDIF 1646 IF( srcv(jpr_cal)%laction ) THEN ! calving (put in emp_tot and emp_oce) 1662 1647 zemp_tot(:,:) = zemp_tot(:,:) - frcv(jpr_cal)%z3(:,:,1) 1663 1648 zemp_oce(:,:) = zemp_oce(:,:) - frcv(jpr_cal)%z3(:,:,1) 1664 CALL iom_put( 'calving_cea', frcv(jpr_cal)%z3(:,:,1) ) 1665 ENDIF 1666 1667 IF( srcv(jpr_icb)%laction ) THEN 1649 ENDIF 1650 IF( srcv(jpr_icb)%laction ) THEN ! iceberg added to runoffs 1668 1651 fwficb(:,:) = frcv(jpr_icb)%z3(:,:,1) 1669 rnf(:,:) = rnf(:,:) + fwficb(:,:) ! iceberg added to runoffs 1670 CALL iom_put( 'iceberg_cea', frcv(jpr_icb)%z3(:,:,1) ) 1671 ENDIF 1672 IF( srcv(jpr_isf)%laction ) THEN 1673 fwfisf(:,:) = - frcv(jpr_isf)%z3(:,:,1) ! fresh water flux from the isf (fwfisf <0 mean melting) 1674 CALL iom_put( 'iceshelf_cea', frcv(jpr_isf)%z3(:,:,1) ) 1675 ENDIF 1676 1652 rnf(:,:) = rnf(:,:) + fwficb(:,:) 1653 ENDIF 1654 IF( srcv(jpr_isf)%laction ) THEN ! iceshelf (fwfisf <0 mean melting) 1655 fwfisf(:,:) = - frcv(jpr_isf)%z3(:,:,1) 1656 ENDIF 1677 1657 1678 1658 IF( ln_mixcpl ) THEN … … 1698 1678 ENDIF 1699 1679 1700 IF( iom_use('subl_ai_cea') ) CALL iom_put( 'subl_ai_cea', zevap_ice(:,:) * zicefr(:,:) ) ! Sublimation over sea-ice (cell average)1701 CALL iom_put( 'snowpre' , sprecip(:,:) ) ! Snow1702 IF( iom_use('snow_ao_cea') ) CALL iom_put( 'snow_ao_cea', sprecip(:,:) * ( 1._wp - zsnw(:,:) ) ) ! Snow over ice-free ocean (cell average)1703 IF( iom_use('snow_ai_cea') ) CALL iom_put( 'snow_ai_cea', sprecip(:,:) * zsnw(:,:) ) ! Snow over sea-ice (cell average)1704 1680 #else 1705 ! runoffs and calving (put in emp_tot) 1706 IF( srcv(jpr_rnf)%laction ) rnf(:,:) = frcv(jpr_rnf)%z3(:,:,1) 1707 IF( iom_use('hflx_rnf_cea') ) & 1708 CALL iom_put( 'hflx_rnf_cea' , rnf(:,:) * zcptn(:,:) ) 1709 IF( srcv(jpr_cal)%laction ) THEN 1681 zsnw(:,:) = zicefr(:,:) 1682 ! --- Continental fluxes --- ! 1683 IF( srcv(jpr_rnf)%laction ) THEN ! runoffs (included in emp later on) 1684 rnf(:,:) = frcv(jpr_rnf)%z3(:,:,1) 1685 ENDIF 1686 IF( srcv(jpr_cal)%laction ) THEN ! calving (put in emp_tot) 1710 1687 zemp_tot(:,:) = zemp_tot(:,:) - frcv(jpr_cal)%z3(:,:,1) 1711 CALL iom_put( 'calving_cea', frcv(jpr_cal)%z3(:,:,1) ) 1712 ENDIF 1713 1714 1715 IF( srcv(jpr_icb)%laction ) THEN 1688 ENDIF 1689 IF( srcv(jpr_icb)%laction ) THEN ! iceberg added to runoffs 1716 1690 fwficb(:,:) = frcv(jpr_icb)%z3(:,:,1) 1717 rnf(:,:) = rnf(:,:) + fwficb(:,:) ! iceberg added to runoffs 1718 CALL iom_put( 'iceberg_cea', frcv(jpr_icb)%z3(:,:,1) ) 1719 ENDIF 1720 IF( srcv(jpr_isf)%laction ) THEN 1721 fwfisf(:,:) = - frcv(jpr_isf)%z3(:,:,1) ! fresh water flux from the isf (fwfisf <0 mean melting) 1722 CALL iom_put( 'iceshelf_cea', frcv(jpr_isf)%z3(:,:,1) ) 1723 ENDIF 1724 1691 rnf(:,:) = rnf(:,:) + fwficb(:,:) 1692 ENDIF 1693 IF( srcv(jpr_isf)%laction ) THEN ! iceshelf (fwfisf <0 mean melting) 1694 fwfisf(:,:) = - frcv(jpr_isf)%z3(:,:,1) 1695 ENDIF 1725 1696 1726 1697 IF( ln_mixcpl ) THEN … … 1736 1707 ENDIF 1737 1708 1738 IF( iom_use('subl_ai_cea') ) CALL iom_put( 'subl_ai_cea', frcv(jpr_ievp)%z3(:,:,1) * zicefr(:,:) ) ! Sublimation over sea-ice (cell average)1739 CALL iom_put( 'snowpre' , sprecip(:,:) ) ! Snow1740 IF( iom_use('snow_ao_cea') ) CALL iom_put( 'snow_ao_cea', sprecip(:,:) * p_frld(:,:) ) ! Snow over ice-free ocean (cell average)1741 IF( iom_use('snow_ai_cea') ) CALL iom_put( 'snow_ai_cea', sprecip(:,:) * zicefr(:,:) ) ! Snow over sea-ice (cell average)1742 1709 #endif 1743 1710 ! outputs 1711 !! IF( srcv(jpr_rnf)%laction ) CALL iom_put( 'runoffs' , rnf(:,:) * tmask(:,:,1) ) ! runoff 1712 !! IF( srcv(jpr_isf)%laction ) CALL iom_put( 'iceshelf_cea', -fwfisf(:,:) * tmask(:,:,1) ) ! iceshelf 1713 IF( srcv(jpr_cal)%laction ) CALL iom_put( 'calving_cea' , frcv(jpr_cal)%z3(:,:,1) * tmask(:,:,1) ) ! calving 1714 IF( srcv(jpr_icb)%laction ) CALL iom_put( 'iceberg_cea' , frcv(jpr_icb)%z3(:,:,1) * tmask(:,:,1) ) ! icebergs 1715 IF( iom_use('snowpre') ) CALL iom_put( 'snowpre' , sprecip(:,:) ) ! Snow 1716 IF( iom_use('precip') ) CALL iom_put( 'precip' , tprecip(:,:) ) ! total precipitation 1717 IF( iom_use('rain') ) CALL iom_put( 'rain' , tprecip(:,:) - sprecip(:,:) ) ! liquid precipitation 1718 IF( iom_use('snow_ao_cea') ) CALL iom_put( 'snow_ao_cea' , sprecip(:,:) * ( 1._wp - zsnw(:,:) ) ) ! Snow over ice-free ocean (cell average) 1719 IF( iom_use('snow_ai_cea') ) CALL iom_put( 'snow_ai_cea' , sprecip(:,:) * zsnw(:,:) ) ! Snow over sea-ice (cell average) 1720 IF( iom_use('subl_ai_cea') ) CALL iom_put( 'subl_ai_cea' , frcv(jpr_ievp)%z3(:,:,1) * zicefr(:,:) * tmask(:,:,1) ) ! Sublimation over sea-ice (cell average) 1721 IF( iom_use('evap_ao_cea') ) CALL iom_put( 'evap_ao_cea' , ( frcv(jpr_tevp)%z3(:,:,1) & 1722 & - frcv(jpr_ievp)%z3(:,:,1) * zicefr(:,:) ) * tmask(:,:,1) ) ! ice-free oce evap (cell average) 1723 ! note: runoff output is done in sbcrnf (which includes icebergs too) and iceshelf output is done in sbcisf 1724 ! 1744 1725 ! ! ========================= ! 1745 1726 SELECT CASE( TRIM( sn_rcv_qns%cldes ) ) ! non solar heat fluxes ! (qns) … … 1777 1758 & + pist(:,:,1) * zicefr(:,:) ) ) 1778 1759 END SELECT 1779 !!gm1780 !! currently it is taken into account in leads budget but not in the zqns_tot, and thus not in1781 !! the flux that enter the ocean....1782 !! moreover 1 - it is not diagnose anywhere....1783 !! 2 - it is unclear for me whether this heat lost is taken into account in the atmosphere or not...1784 !!1785 !! similar job should be done for snow and precipitation temperature1786 1760 ! 1787 IF( srcv(jpr_cal)%laction ) THEN ! Iceberg melting 1788 zqns_tot(:,:) = zqns_tot(:,:) - frcv(jpr_cal)%z3(:,:,1) * lfus ! add the latent heat of iceberg melting 1789 ! we suppose it melts at 0deg, though it should be temp. of surrounding ocean 1790 IF( iom_use('hflx_cal_cea') ) CALL iom_put( 'hflx_cal_cea', - frcv(jpr_cal)%z3(:,:,1) * lfus ) ! heat flux from calving 1791 ENDIF 1792 1793 !!chris 1794 !! The heat content associated to the ice shelf in removed in the routine sbcisf.F90 1795 ! 1796 IF( srcv(jpr_icb)%laction ) zqns_tot(:,:) = zqns_tot(:,:) - frcv(jpr_icb)%z3(:,:,1) * lfus ! remove heat content associated to iceberg melting 1797 1761 ! --- calving (removed from qns_tot) --- ! 1762 IF( srcv(jpr_cal)%laction ) zqns_tot(:,:) = zqns_tot(:,:) - frcv(jpr_cal)%z3(:,:,1) * lfus ! remove latent heat of calving 1763 ! we suppose it melts at 0deg, though it should be temp. of surrounding ocean 1764 ! --- iceberg (removed from qns_tot) --- ! 1765 IF( srcv(jpr_icb)%laction ) zqns_tot(:,:) = zqns_tot(:,:) - frcv(jpr_icb)%z3(:,:,1) * lfus ! remove latent heat of iceberg melting 1798 1766 1799 1767 #if defined key_lim3 … … 1850 1818 ENDIF 1851 1819 1852 ! some more outputs1853 IF( iom_use('hflx_snow_cea') ) CALL iom_put('hflx_snow_cea', sprecip(:,:) * ( zcptsnw(:,:) - Lfus ) ) ! heat flux from snow (cell average)1854 IF( iom_use('hflx_rain_cea') ) CALL iom_put('hflx_rain_cea', ( tprecip(:,:) - sprecip(:,:) ) * zcptrain(:,:) ) ! heat flux from rain (cell average)1855 IF( iom_use('hflx_snow_ao_cea') ) CALL iom_put('hflx_snow_ao_cea',sprecip(:,:) * ( zcptsnw(:,:) - Lfus ) * (1._wp - zsnw(:,:)) ) ! heat flux from snow (over ocean)1856 IF( iom_use('hflx_snow_ai_cea') ) CALL iom_put('hflx_snow_ai_cea',sprecip(:,:) * ( zcptsnw(:,:) - Lfus ) * zsnw(:,:) ) ! heat flux from snow (over ice)1857 1858 1820 #else 1821 zcptsnw (:,:) = zcptn(:,:) 1822 zcptrain(:,:) = zcptn(:,:) 1823 1859 1824 ! clem: this formulation is certainly wrong... but better than it was... 1860 zqns_tot(:,:) = zqns_tot(:,:) &! zqns_tot update over free ocean with:1861 & - ztmp(:,:) &! remove the latent heat flux of solid precip. melting1862 & - ( zemp_tot(:,:) &! remove the heat content of mass flux (assumed to be at SST)1825 zqns_tot(:,:) = zqns_tot(:,:) & ! zqns_tot update over free ocean with: 1826 & - ( p_frld(:,:) * zsprecip(:,:) * lfus ) & ! remove the latent heat flux of solid precip. melting 1827 & - ( zemp_tot(:,:) & ! remove the heat content of mass flux (assumed to be at SST) 1863 1828 & - zemp_ice(:,:) ) * zcptn(:,:) 1864 1829 … … 1873 1838 qns_ice(:,:,:) = zqns_ice(:,:,:) 1874 1839 ENDIF 1840 1875 1841 #endif 1876 1842 ! outputs 1843 IF( srcv(jpr_cal)%laction ) CALL iom_put('hflx_cal_cea' , - frcv(jpr_cal)%z3(:,:,1) * lfus ) ! latent heat from calving 1844 IF( srcv(jpr_icb)%laction ) CALL iom_put('hflx_icb_cea' , - frcv(jpr_icb)%z3(:,:,1) * lfus ) ! latent heat from icebergs melting 1845 IF( iom_use('hflx_snow_cea') ) CALL iom_put('hflx_snow_cea', sprecip(:,:) * ( zcptsnw(:,:) - Lfus ) ) ! heat flux from snow (cell average) 1846 IF( iom_use('hflx_rain_cea') ) CALL iom_put('hflx_rain_cea',( tprecip(:,:) - sprecip(:,:) ) * zcptrain(:,:) ) ! heat flux from rain (cell average) 1847 IF( iom_use('hflx_evap_cea') ) CALL iom_put('hflx_evap_cea',(frcv(jpr_tevp)%z3(:,:,1) - frcv(jpr_ievp)%z3(:,:,1) * zicefr(:,:) & ! heat flux from from evap (cell average) 1848 & ) * zcptn(:,:) * tmask(:,:,1) ) 1849 IF( iom_use('hflx_snow_ao_cea') ) CALL iom_put('hflx_snow_ao_cea',sprecip(:,:) * (zcptsnw(:,:) - Lfus) * (1._wp - zsnw(:,:)) ) ! heat flux from snow (over ocean) 1850 IF( iom_use('hflx_snow_ai_cea') ) CALL iom_put('hflx_snow_ai_cea',sprecip(:,:) * (zcptsnw(:,:) - Lfus) * zsnw(:,:) ) ! heat flux from snow (over ice) 1851 ! note: hflx for runoff and iceshelf are done in sbcrnf and sbcisf resp. 1852 ! 1877 1853 ! ! ========================= ! 1878 1854 SELECT CASE( TRIM( sn_rcv_qsr%cldes ) ) ! solar heat fluxes ! (qsr) … … 1980 1956 fr2_i0(:,:) = ( 0.82 * ( 1.0 - cldf_ice ) + 0.65 * cldf_ice ) 1981 1957 1982 CALL wrk_dealloc( jpi,jpj, zcptn, z tmp, zcptrain, zcptsnw, zicefr, zmsk, zsnw )1958 CALL wrk_dealloc( jpi,jpj, zcptn, zcptrain, zcptsnw, zicefr, zmsk, zsnw ) 1983 1959 CALL wrk_dealloc( jpi,jpj, zemp_tot, zemp_ice, zemp_oce, ztprecip, zsprecip, zevap_oce, zevap_ice, zdevap_ice ) 1984 1960 CALL wrk_dealloc( jpi,jpj, zqns_tot, zqns_oce, zqsr_tot, zqsr_oce, zqprec_ice, zqemp_oce, zqemp_ice ) -
trunk/NEMOGCM/NEMO/OPA_SRC/SBC/sbcisf.F90
r7816 r7968 168 168 169 169 ! output 170 CALL iom_put('qlatisf' , qisf) 171 CALL iom_put('fwfisf' , fwfisf) 170 IF( iom_use('iceshelf_cea') ) CALL iom_put( 'iceshelf_cea', -fwfisf(:,:) ) ! isf mass flux 171 IF( iom_use('hflx_isf_cea') ) CALL iom_put( 'hflx_isf_cea', risf_tsc(:,:,jp_tem) * rau0 * rcp ) ! isf sensible+latent heat (W/m2) 172 IF( iom_use('qlatisf' ) ) CALL iom_put( 'qlatisf' , qisf(:,:) ) ! isf latent heat 173 IF( iom_use('fwfisf' ) ) CALL iom_put( 'fwfisf' , fwfisf(:,:) ) ! isf mass flux (opposite sign) 172 174 173 175 ! Diagnostics -
trunk/NEMOGCM/NEMO/OPA_SRC/SBC/sbcrnf.F90
r7753 r7968 138 138 IF( ln_rnf_sal ) rnf_tsc(:,:,jp_sal) = ( sf_s_rnf(1)%fnow(:,:,1) ) * rnf(:,:) * r1_rau0 139 139 ! ! else use S=0 for runoffs (done one for all in the init) 140 CALL iom_put( "runoffs", rnf ) ! output runoffs arrays 140 IF( iom_use('runoffs') ) CALL iom_put( 'runoffs' , rnf(:,:) ) ! output runoff mass flux 141 IF( iom_use('hflx_rnf_cea') ) CALL iom_put( 'hflx_rnf_cea', rnf_tsc(:,:,jp_tem) * rau0 * rcp ) ! output runoff sensible heat (W/m2) 141 142 ENDIF 142 143 !
Note: See TracChangeset
for help on using the changeset viewer.
