Ignore:
Timestamp:
2017-09-22T16:55:24+02:00 (3 years ago)
Author:
clem
Message:

almost useless commits

File:
1 edited

Legend:

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

    r8534 r8559  
    293293      DO jk = 1, nlay_i 
    294294         DO ji = 1, nidx 
    295             ztmelts           = - tmut * s_i_1d(ji,jk) + rt0          ! Melting point of layer k [K] 
    296              
    297             IF( t_i_1d(ji,jk) >= ztmelts ) THEN !!! Internal melting 
     295            ztmelts           = - tmut * s_i_1d(ji,jk)          ! Melting point of layer k [C] 
     296             
     297            IF( t_i_1d(ji,jk) >= (ztmelts+rt0) ) THEN !!! Internal melting 
    298298 
    299299               zEi            = - e_i_1d(ji,jk) * r1_rhoic            ! Specific enthalpy of layer k [J/kg, <0]        
     
    320320                
    321321               zEi            = - e_i_1d(ji,jk) * r1_rhoic            ! Specific enthalpy of layer k [J/kg, <0] 
    322                zEw            =    rcp * ( ztmelts - rt0 )            ! Specific enthalpy of resulting meltwater [J/kg, <0] 
     322               zEw            =    rcp * ztmelts                      ! Specific enthalpy of resulting meltwater [J/kg, <0] 
    323323               zdE            =    zEi - zEw                          ! Specific enthalpy difference < 0 
    324324                
     
    430430                                  + ( 1. - zswitch_sal ) * sm_i_1d(ji)  
    431431               ! New ice growth 
    432                ztmelts            = - tmut * s_i_new(ji) + rt0          ! New ice melting point (K) 
     432               ztmelts            = - tmut * s_i_new(ji)                 ! New ice melting point (C) 
    433433 
    434434               zt_i_new           = zswitch_sal * t_bo_1d(ji) + ( 1. - zswitch_sal) * t_i_1d(ji, nlay_i) 
    435435                
    436                zEi                = cpic * ( zt_i_new - ztmelts ) &     ! Specific enthalpy of forming ice (J/kg, <0)       
    437                   &               - lfus * ( 1.0 - ( ztmelts - rt0 ) / ( zt_i_new - rt0 ) )   & 
    438                   &               + rcp  * ( ztmelts-rt0 )           
     436               zEi                = cpic * ( zt_i_new - (ztmelts+rt0) ) &     ! Specific enthalpy of forming ice (J/kg, <0)       
     437                  &               - lfus * ( 1.0 - ztmelts / ( zt_i_new - rt0 ) )   & 
     438                  &               + rcp  * ztmelts           
    439439 
    440440               zEw                = rcp  * ( t_bo_1d(ji) - rt0 )         ! Specific enthalpy of seawater (J/kg, < 0) 
     
    450450            zfmdt          = - rhoic * dh_i_bott(ji)             ! Mass flux x time step (kg/m2, < 0) 
    451451 
    452             ztmelts        = - tmut * s_i_new(ji) + rt0          ! New ice melting point (K) 
     452            ztmelts        = - tmut * s_i_new(ji)                ! New ice melting point (C) 
    453453             
    454454            zt_i_new       = zswitch_sal * t_bo_1d(ji) + ( 1. - zswitch_sal) * t_i_1d(ji, nlay_i) 
    455455             
    456             zEi            = cpic * ( zt_i_new - ztmelts ) &     ! Specific enthalpy of forming ice (J/kg, <0)       
    457                &               - lfus * ( 1.0 - ( ztmelts - rt0 ) / ( zt_i_new - rt0 ) )   & 
    458                &               + rcp  * ( ztmelts-rt0 )           
     456            zEi            = cpic * ( zt_i_new - (ztmelts+rt0) ) &     ! Specific enthalpy of forming ice (J/kg, <0)       
     457               &               - lfus * ( 1.0 - ztmelts / ( zt_i_new - rt0 ) )   & 
     458               &               + rcp  * ztmelts           
    459459             
    460460            zEw            = rcp  * ( t_bo_1d(ji) - rt0 )         ! Specific enthalpy of seawater (J/kg, < 0) 
     
    490490            IF(  zf_tt(ji)  >  0._wp  .AND. jk > icount(ji,jk) ) THEN   ! do not calculate where layer has already disappeared by surface melting  
    491491 
    492                ztmelts = - tmut * s_i_1d(ji,jk) + rt0  ! Melting point of layer jk (K) 
    493  
    494                IF( t_i_1d(ji,jk) >= ztmelts ) THEN !!! Internal melting 
     492               ztmelts = - tmut * s_i_1d(ji,jk)  ! Melting point of layer jk (C) 
     493 
     494               IF( t_i_1d(ji,jk) >= (ztmelts+rt0) ) THEN !!! Internal melting 
    495495 
    496496                  zEi               = - e_i_1d(ji,jk) * r1_rhoic    ! Specific enthalpy of melting ice (J/kg, <0) 
     
    520520 
    521521                  zEi             = - e_i_1d(ji,jk) * r1_rhoic ! Specific enthalpy of melting ice (J/kg, <0) 
    522                   zEw             = rcp * ( ztmelts - rt0 )    ! Specific enthalpy of meltwater (J/kg, <0) 
     522                  zEw             = rcp * ztmelts              ! Specific enthalpy of meltwater (J/kg, <0) 
    523523                  zdE             = zEi - zEw                  ! Specific enthalpy difference   (J/kg, <0) 
    524524 
Note: See TracChangeset for help on using the changeset viewer.