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 8565 for branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icethd_dh.F90 – NEMO

Ignore:
Timestamp:
2017-09-27T12:09:10+02:00 (7 years ago)
Author:
clem
Message:

trying to respect naming convention

File:
1 edited

Legend:

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

    r8564 r8565  
    117117 
    118118      ! initialize layer thicknesses and enthalpies 
    119       h_i_old (1:nidx,0:nlay_i+1) = 0._wp 
    120       eh_i_old(1:nidx,0:nlay_i+1) = 0._wp 
     119      h_i_old (1:npti,0:nlay_i+1) = 0._wp 
     120      eh_i_old(1:npti,0:nlay_i+1) = 0._wp 
    121121      DO jk = 1, nlay_i 
    122          DO ji = 1, nidx 
     122         DO ji = 1, npti 
    123123            h_i_old (ji,jk) = h_i_1d(ji) * r1_nlay_i 
    124124            eh_i_old(ji,jk) = e_i_1d(ji,jk) * h_i_old(ji,jk) 
     
    130130      !------------------------------------------------------------------------------! 
    131131      ! 
    132       DO ji = 1, nidx 
     132      DO ji = 1, npti 
    133133         zdum       = qns_ice_1d(ji) + ( 1._wp - i0(ji) ) * qsr_ice_1d(ji) - fc_su(ji)  
    134134         zf_tt(ji)  = fc_bo_i(ji) + fhtur_1d(ji) + fhld_1d(ji)  
     
    142142      ! (should not happen but sometimes it does) 
    143143      !------------------------------------------------------------------------------! 
    144       DO ji = 1, nidx 
     144      DO ji = 1, npti 
    145145         IF( t_s_1d(ji,1) > rt0 ) THEN !!! Internal melting 
    146146            ! Contribution to heat flux to the ocean [W.m-2], < 0   
     
    159159      !------------------------------------------------------------! 
    160160      DO jk = 1, nlay_i 
    161          DO ji = 1, nidx 
     161         DO ji = 1, npti 
    162162            zh_i(ji,jk) = h_i_1d(ji) * r1_nlay_i 
    163163            zeh_i(ji)   = zeh_i(ji) + e_i_1d(ji,jk) * zh_i(ji,jk) 
     
    183183      ! Martin Vancoppenolle, December 2006 
    184184 
    185       CALL ice_thd_snwblow( 1. - at_i_1d(1:nidx), zsnw(1:nidx) ) ! snow distribution over ice after wind blowing 
    186  
    187       zdeltah(1:nidx,:) = 0._wp 
    188       DO ji = 1, nidx 
     185      CALL ice_thd_snwblow( 1. - at_i_1d(1:npti), zsnw(1:npti) ) ! snow distribution over ice after wind blowing 
     186 
     187      zdeltah(1:npti,:) = 0._wp 
     188      DO ji = 1, npti 
    189189         !----------- 
    190190         ! Snow fall 
     
    220220 
    221221      ! If heat still available (zq_su > 0), then melt more snow 
    222       zdeltah(1:nidx,:) = 0._wp 
    223       zdh_s_mel(1:nidx) = 0._wp 
     222      zdeltah(1:npti,:) = 0._wp 
     223      zdh_s_mel(1:npti) = 0._wp 
    224224      DO jk = 1, nlay_s 
    225          DO ji = 1, nidx 
     225         DO ji = 1, npti 
    226226            ! thickness change 
    227227            rswitch          = 1._wp - MAX( 0._wp, SIGN( 1._wp, - h_s_1d(ji) ) )  
     
    245245      ! qla_ice is always >=0 (upwards), heat goes to the atmosphere, therefore snow sublimates 
    246246      ! clem comment: not counted in mass/heat exchange in iceupdate.F90 since this is an exchange with atm. (not ocean) 
    247       zdeltah(1:nidx,:) = 0._wp 
    248       DO ji = 1, nidx 
     247      zdeltah(1:npti,:) = 0._wp 
     248      DO ji = 1, npti 
    249249         zdh_s_sub(ji)  = MAX( - h_s_1d(ji) , - evap_ice_1d(ji) * r1_rhosn * rdt_ice ) 
    250250         ! remaining evap in kg.m-2 (used for ice melting later on) 
     
    265265       
    266266      ! --- Update snow diags --- ! 
    267       DO ji = 1, nidx 
     267      DO ji = 1, npti 
    268268         dh_s_tot(ji) = zdh_s_mel(ji) + zdh_s_pre(ji) + zdh_s_sub(ji) 
    269269      END DO 
     
    274274      ! new temp and enthalpy of the snow (remaining snow precip + remaining pre-existing snow) 
    275275      DO jk = 1, nlay_s 
    276          DO ji = 1,nidx 
     276         DO ji = 1,npti 
    277277            rswitch       = MAX( 0._wp , SIGN( 1._wp, h_s_1d(ji) - epsi20 ) ) 
    278278            e_s_1d(ji,jk) = rswitch / MAX( h_s_1d(ji), epsi20 ) *           & 
     
    285285      ! 3.4 Surface ice ablation  
    286286      !-------------------------- 
    287       zdeltah(1:nidx,:) = 0._wp ! important 
     287      zdeltah(1:npti,:) = 0._wp ! important 
    288288      DO jk = 1, nlay_i 
    289          DO ji = 1, nidx 
     289         DO ji = 1, npti 
    290290            ztmelts           = - tmut * sz_i_1d(ji,jk)          ! Melting point of layer k [C] 
    291291             
     
    375375      END DO 
    376376      ! update ice thickness 
    377       DO ji = 1, nidx 
     377      DO ji = 1, npti 
    378378         h_i_1d(ji) =  MAX( 0._wp , h_i_1d(ji) + dh_i_surf(ji) + dh_i_sub(ji) ) 
    379379      END DO 
    380380 
    381381      ! remaining "potential" evap is sent to ocean 
    382       DO ji = 1, nidx 
     382      DO ji = 1, npti 
    383383         wfx_err_sub_1d(ji) = wfx_err_sub_1d(ji) - zevap_rema(ji) * a_i_1d(ji) * r1_rdtice  ! <=0 (net evap for the ocean in kg.m-2.s-1) 
    384384      END DO 
     
    407407       
    408408      ! Iterative procedure 
    409       DO ji = 1, nidx 
     409      DO ji = 1, npti 
    410410         IF(  zf_tt(ji) < 0._wp  ) THEN 
    411411            DO iter = 1, num_iter_max 
     
    478478      ! 4.2 Basal melt 
    479479      !---------------- 
    480       zdeltah(1:nidx,:) = 0._wp ! important 
     480      zdeltah(1:npti,:) = 0._wp ! important 
    481481      DO jk = nlay_i, 1, -1 
    482          DO ji = 1, nidx 
     482         DO ji = 1, npti 
    483483            IF(  zf_tt(ji)  >  0._wp  .AND. jk > icount(ji,jk) ) THEN   ! do not calculate where layer has already disappeared by surface melting  
    484484 
     
    554554      ! Update temperature, energy 
    555555      !------------------------------------------- 
    556       DO ji = 1, nidx 
     556      DO ji = 1, npti 
    557557         h_i_1d(ji) =  MAX( 0._wp , h_i_1d(ji) + dh_i_bott(ji) ) 
    558558      END DO   
     
    563563      ! If heat still available for melting and snow remains, then melt more snow 
    564564      !------------------------------------------- 
    565       zdeltah(1:nidx,:) = 0._wp ! important 
    566       DO ji = 1, nidx 
     565      zdeltah(1:npti,:) = 0._wp ! important 
     566      DO ji = 1, npti 
    567567         zq_rema(ji)     = zq_su(ji) + zq_bo(ji)  
    568568         rswitch         = 1._wp - MAX( 0._wp, SIGN( 1._wp, - h_s_1d(ji) ) )   ! =1 if snow 
     
    592592      ! flooding of seawater transforms snow into ice dh_snowice is positive for the ice 
    593593      z1_rho = 1._wp / ( rhosn+rau0-rhoic ) 
    594       DO ji = 1, nidx 
     594      DO ji = 1, npti 
    595595         ! 
    596596         dh_snowice(ji) = MAX(  0._wp , ( rhosn * h_s_1d(ji) + (rhoic-rau0) * h_i_1d(ji) ) * z1_rho ) 
     
    631631      ! Update temperature, energy 
    632632      !------------------------------------------- 
    633       DO ji = 1, nidx 
     633      DO ji = 1, npti 
    634634         rswitch     =  1.0 - MAX( 0._wp , SIGN( 1._wp , - h_i_1d(ji) ) )  
    635635         t_su_1d(ji) =  rswitch * t_su_1d(ji) + ( 1.0 - rswitch ) * rt0 
     
    637637 
    638638      DO jk = 1, nlay_s 
    639          DO ji = 1,nidx 
     639         DO ji = 1,npti 
    640640            ! mask enthalpy 
    641641            rswitch       = 1._wp - MAX(  0._wp , SIGN( 1._wp, - h_s_1d(ji) )  ) 
     
    647647 
    648648      ! --- ensure that a_i = 0 where h_i = 0 --- 
    649       WHERE( h_i_1d(1:nidx) == 0._wp )   a_i_1d(1:nidx) = 0._wp 
     649      WHERE( h_i_1d(1:npti) == 0._wp )   a_i_1d(1:npti) = 0._wp 
    650650      ! 
    651651   END SUBROUTINE ice_thd_dh 
Note: See TracChangeset for help on using the changeset viewer.