Ignore:
Timestamp:
2017-09-14T17:52:02+02:00 (3 years ago)
Author:
clem
Message:

changes in style - part6 - ice diffusion (hope the scheme still converges)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icewri.F90

    r8518 r8522  
    179179      ENDIF 
    180180 
    181       ztmp = rday / rhosn 
     181      ztmp = rday * r1_rhosn 
    182182      CALL iom_put( "vfxspr"     , wfx_spr * ztmp       )        ! precip (snow) 
    183183      CALL iom_put( "vfxsnw"     , wfx_snw * ztmp       )        ! total snw growth/melt  
     
    237237      ! Add-ons for SIMIP 
    238238      !-------------------------------- 
    239       zrho1 = ( rau0 - rhoic ) / rau0; zrho2 = rhosn / rau0 
     239      zrho1 = ( rau0 - rhoic ) * r1_rau0; zrho2 = rhosn * r1_rau0 
    240240 
    241241      IF ( iom_use( "icepres"  ) ) CALL iom_put( "icepres"     , zswi(:,:)                     )                                ! Ice presence (1 or 0)  
     
    278278 
    279279      IF ( iom_use( "dmsspr"   ) ) CALL iom_put( "dmsspr"      , - wfx_spr                  )                            ! Snow mass change through snow fall 
    280       IF ( iom_use( "dmsssi"   ) ) CALL iom_put( "dmsssi"      ,   wfx_sni*rhosn/rhoic      )                            ! Snow mass change through snow-to-ice conversion 
     280      IF ( iom_use( "dmsssi"   ) ) CALL iom_put( "dmsssi"      ,   wfx_sni*rhosn*r1_rhoic   )                            ! Snow mass change through snow-to-ice conversion 
    281281 
    282282      IF ( iom_use( "dmsmel"   ) ) CALL iom_put( "dmsmel"      , - wfx_snw_sum              )                            ! Snow mass change through melt 
Note: See TracChangeset for help on using the changeset viewer.