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

changes in style - part5 - almost done

File:
1 edited

Legend:

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

    r8500 r8514  
    287287         !------------------------------------------------------------------------------| 
    288288         ! 
    289          IF( nn_ice_thcon == 0 ) THEN      ! Untersteiner (1964) formula 
     289         IF( ln_cndi_U64 ) THEN         !-- Untersteiner (1964) formula: k = k0 + beta.S/T 
    290290            DO ji = 1 , nidx 
    291291               ztcond_i(ji,0) = rcdic + zbeta * s_i_1d(ji,1) / MIN( -epsi10, t_i_1d(ji,1) - rt0 ) 
     
    299299               END DO 
    300300            END DO 
    301          ENDIF 
    302  
    303          IF( nn_ice_thcon == 1 ) THEN      ! Pringle et al formula included: 2.11 + 0.09 S/T - 0.011.T 
     301 
     302         ELSEIF( ln_cndi_P07 ) THEN     !-- Pringle et al formula: k = k0 + beta1.S/T - beta2.T 
    304303            DO ji = 1 , nidx 
    305304               ztcond_i(ji,0) = rcdic + 0.090_wp * s_i_1d(ji,1) / MIN( -epsi10, t_i_1d(ji,1) - rt0 )   & 
     
    347346 
    348347               ! Effective thickness he (zhe) 
    349                zfac     = 1._wp / ( rn_cdsn + zkimean ) 
    350                zratio_s = rn_cdsn   * zfac 
     348               zfac     = 1._wp / ( rn_cnd_s + zkimean ) 
     349               zratio_s = rn_cnd_s   * zfac 
    351350               zratio_i = zkimean * zfac 
    352351               zhe      = zratio_s * ht_i_1d(ji) + zratio_i * ht_s_1d(ji) 
     
    371370         DO ji = 1, nidx 
    372371            zfac                  =  1. / MAX( epsi10 , zh_s(ji) ) 
    373             zkappa_s(ji,0)        = zghe(ji) * rn_cdsn * zfac 
    374             zkappa_s(ji,nlay_s)   = zghe(ji) * rn_cdsn * zfac 
     372            zkappa_s(ji,0)        = zghe(ji) * rn_cnd_s * zfac 
     373            zkappa_s(ji,nlay_s)   = zghe(ji) * rn_cnd_s * zfac 
    375374         END DO 
    376375 
    377376         DO jk = 1, nlay_s-1 
    378377            DO ji = 1 , nidx 
    379                zkappa_s(ji,jk)    = zghe(ji) * 2.0 * rn_cdsn / MAX( epsi10, 2.0 * zh_s(ji) ) 
     378               zkappa_s(ji,jk)    = zghe(ji) * 2.0 * rn_cnd_s / MAX( epsi10, 2.0 * zh_s(ji) ) 
    380379            END DO 
    381380         END DO 
     
    393392            zkappa_i(ji,0)        = zghe(ji) * ztcond_i(ji,0) * zfac 
    394393            zkappa_i(ji,nlay_i)   = zghe(ji) * ztcond_i(ji,nlay_i) * zfac 
    395             zkappa_s(ji,nlay_s)   = zghe(ji) * zghe(ji) * 2.0 * rn_cdsn * ztcond_i(ji,0) / &  
    396            &                        MAX( epsi10, ( zghe(ji) * ztcond_i(ji,0) * zh_s(ji) + zghe(ji) * rn_cdsn * zh_i(ji) ) ) 
     394            zkappa_s(ji,nlay_s)   = zghe(ji) * zghe(ji) * 2.0 * rn_cnd_s * ztcond_i(ji,0) / &  
     395           &                        MAX( epsi10, ( zghe(ji) * ztcond_i(ji,0) * zh_s(ji) + zghe(ji) * rn_cnd_s * zh_i(ji) ) ) 
    397396            zkappa_i(ji,0)        = zkappa_s(ji,nlay_s) * isnow(ji) + zkappa_i(ji,0) * ( 1._wp - isnow(ji) ) 
    398397         END DO 
     
    423422         !------------------------------------------------------------------------------| 
    424423         ! 
    425          IF ( ln_dqnsice ) THEN  
     424         IF ( ln_dqns_i ) THEN  
    426425            DO ji = 1 , nidx 
    427426               ! update of the non solar flux according to the update in T_su 
     
    708707      !--- Snow-ice interfacial temperature (diagnostic SIMIP) 
    709708      DO ji = 1, nidx 
    710          zfac        = 1. / MAX( epsi10 , rn_cdsn * zh_i(ji) + ztcond_i(ji,1) * zh_s(ji) ) 
     709         zfac        = 1. / MAX( epsi10 , rn_cnd_s * zh_i(ji) + ztcond_i(ji,1) * zh_s(ji) ) 
    711710         IF( zh_s(ji) >= 1.e-3 ) THEN 
    712             t_si_1d(ji) = ( rn_cdsn        * zh_i(ji) * t_s_1d(ji,1) + & 
     711            t_si_1d(ji) = ( rn_cnd_s       * zh_i(ji) * t_s_1d(ji,1) + & 
    713712               &            ztcond_i(ji,1) * zh_s(ji) * t_i_1d(ji,1) ) * zfac 
    714713         ELSE 
     
    718717      ! END MV SIMIP 2016 
    719718 
    720       IF( ln_limctl .AND. lwp ) THEN 
     719      IF( ln_icectl .AND. lwp ) THEN 
    721720         WRITE(numout,*) ' zdti_max : ', zdti_max 
    722721         WRITE(numout,*) ' iconv    : ', iconv 
     
    748747 
    749748      ! --- diagnose the change in non-solar flux due to surface temperature change --- ! 
    750       IF ( ln_dqnsice ) THEN 
     749      IF ( ln_dqns_i ) THEN 
    751750         DO ji = 1, nidx 
    752751            hfx_err_dif_1d(ji) = hfx_err_dif_1d(ji) - ( qns_ice_1d(ji)  - zqns_ice_b(ji) ) * a_i_1d(ji)  
Note: See TracChangeset for help on using the changeset viewer.