New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 8563 for branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icethd_da.F90 – NEMO

Ignore:
Timestamp:
2017-09-26T15:24:17+02:00 (7 years ago)
Author:
clem
Message:

change variable names (ht_s => h_s & ht_i => h_i)

File:
1 edited

Legend:

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

    r8534 r8563  
    135135             
    136136            ! Contribution to salt flux 
    137             sfx_lam_1d(ji) = sfx_lam_1d(ji) + rhoic *  ht_i_1d(ji) * zda * sm_i_1d(ji) * r1_rdtice 
     137            sfx_lam_1d(ji) = sfx_lam_1d(ji) + rhoic *  h_i_1d(ji) * zda * sm_i_1d(ji) * r1_rdtice 
    138138             
    139139            ! Contribution to heat flux into the ocean [W.m-2], (<0)   
    140             hfx_thd_1d(ji) = hfx_thd_1d(ji) - zda * r1_rdtice * ( ht_i_1d(ji) * r1_nlay_i * SUM( e_i_1d(ji,1:nlay_i) )  & 
    141                                                                 + ht_s_1d(ji) * r1_nlay_s * SUM( e_s_1d(ji,1:nlay_s) ) )  
     140            hfx_thd_1d(ji) = hfx_thd_1d(ji) - zda * r1_rdtice * ( h_i_1d(ji) * r1_nlay_i * SUM( e_i_1d(ji,1:nlay_i) )  & 
     141                                                                + h_s_1d(ji) * r1_nlay_s * SUM( e_s_1d(ji,1:nlay_s) ) )  
    142142             
    143143            ! Contribution to mass flux 
    144             wfx_lam_1d(ji) =  wfx_lam_1d(ji) + zda * r1_rdtice * ( rhoic * ht_i_1d(ji) + rhosn * ht_s_1d(ji) ) 
     144            wfx_lam_1d(ji) =  wfx_lam_1d(ji) + zda * r1_rdtice * ( rhoic * h_i_1d(ji) + rhosn * h_s_1d(ji) ) 
    145145             
    146146            ! new concentration 
    147147            a_i_1d(ji) = a_i_1d(ji) - zda 
    148148 
    149             ! ensure that ht_i = 0 where a_i = 0 
     149            ! ensure that h_i = 0 where a_i = 0 
    150150            IF( a_i_1d(ji) == 0._wp ) THEN 
    151                ht_i_1d(ji) = 0._wp 
    152                ht_s_1d(ji) = 0._wp 
     151               h_i_1d(ji) = 0._wp 
     152               h_s_1d(ji) = 0._wp 
    153153            ENDIF 
    154154         ENDIF 
Note: See TracChangeset for help on using the changeset viewer.