Ignore:
Timestamp:
2019-01-25T13:45:45+01:00 (19 months ago)
Author:
clem
Message:

add a comment

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/src/ICE/icectl.F90

    r10576 r10581  
    186186      !!------------------------------------------------------------------- 
    187187      CHARACTER(len=*), INTENT(in)    :: cd_routine    ! name of the routine 
    188       REAL(wp)                        :: zhfx, zsfx, zvfx 
     188      REAL(wp)                        :: zqmass, zhfx, zsfx, zvfx 
    189189      REAL(wp)                        :: zarea, zv_sill, zs_sill, zt_sill 
    190190      REAL(wp), PARAMETER             :: zconv = 1.e-9 ! convert W to GW and kg to Mt 
     
    199199      ! heat flux 
    200200      ! clem: not the good formulation 
    201       !!zhfx  = glob_sum( 'icectl', ( qt_atm_oi - qt_oce_ai - diag_heat ) * e1e2t ) * zconv 
     201      !!zhfx  = glob_sum( 'icectl', ( qt_oce_ai - qt_atm_oi + diag_heat + hfx_thd + hfx_dyn + hfx_res + hfx_sub + hfx_spr  & 
     202      !!   &                        ) * e1e2t ) * zconv 
    202203 
    203204      ! set threshold values and calculate the ice area (+epsi10 to set a threshold > 0 when there is no ice)  
Note: See TracChangeset for help on using the changeset viewer.