Changeset 9448


Ignore:
Timestamp:
2018-03-28T18:42:49+02:00 (2 years ago)
Author:
clem
Message:

correct minor bugs mostly for outputs

Location:
branches/2017/dev_merge_2017/NEMOGCM/NEMO/LIM_SRC_3
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/LIM_SRC_3/iceistate.F90

    r9416 r9448  
    114114      DO jl = 1, jpl 
    115115         t_su   (:,:,jl) = rt0 * tmask(:,:,1)  ! temp at the surface 
    116          t_si   (:,:,jl) = rt0 * tmask(:,:,1)  ! temp at the ice-snow interface 
    117116         cnd_ice(:,:,jl) = 0._wp               ! initialisation of the effective conductivity at the top of ice/snow (Jules coupling) 
    118117      END DO 
  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/LIM_SRC_3/icestp.F90

    r9409 r9448  
    409409      ! SIMIP diagnostics 
    410410      diag_fc_bo(:,:) = 0._wp ; diag_fc_su(:,:) = 0._wp 
     411      t_si(:,:,:) = rt0       ! temp at the ice-snow interface 
    411412 
    412413      tau_icebfr(:,:)   = 0._wp   ! landfast ice param only (clem: important to keep the init here) 
  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/LIM_SRC_3/icevar.F90

    r9433 r9448  
    464464         ! Zap ice energy and use ocean heat to melt ice 
    465465         !----------------------------------------------------------------- 
    466          WHERE( a_i(:,:,jl) > epsi20 )   ;   h_i(:,:,jl) = v_i(:,:,jl) / a_i(:,:,jl) 
     466         WHERE( a_i(:,:,jl) > epsi10 )   ;   h_i(:,:,jl) = v_i(:,:,jl) / a_i(:,:,jl) 
    467467         ELSEWHERE                       ;   h_i(:,:,jl) = 0._wp 
    468468         END WHERE 
     
    504504               ! zap ice and snow volume, add water and salt to ocean 
    505505               !----------------------------------------------------------------- 
    506                ato_i(ji,jj)    = a_i (ji,jj,jl) * ( 1._wp - zswitch(ji,jj) ) + ato_i(ji,jj) 
    507506               a_i  (ji,jj,jl) = a_i (ji,jj,jl) * zswitch(ji,jj) 
    508507               v_i  (ji,jj,jl) = v_i (ji,jj,jl) * zswitch(ji,jj) 
Note: See TracChangeset for help on using the changeset viewer.