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 9734 for branches – NEMO

Changeset 9734 for branches


Ignore:
Timestamp:
2018-06-04T15:58:26+02:00 (6 years ago)
Author:
dancopsey
Message:

Changed tn_ice to t1_ice in sbccpl

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_merge_2017_CICE_interface/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90

    r9499 r9734  
    18951895 
    18961896      ! Heat content per unit mass of snow (J/kg) 
    1897       WHERE( SUM( a_i, dim=3 ) > 1.e-10 )   ;   zcptsnw(:,:) = cpic * SUM( (tn_ice - rt0) * a_i, dim=3 ) / SUM( a_i, dim=3 ) 
     1897      WHERE( SUM( a_i, dim=3 ) > 1.e-10 )   ;   zcptsnw(:,:) = cpic * SUM( (t1_ice - rt0) * a_i, dim=3 ) / SUM( a_i, dim=3 ) 
    18981898      ELSEWHERE                             ;   zcptsnw(:,:) = zcptn(:,:) 
    18991899      ENDWHERE 
    19001900      ! Heat content per unit mass of rain (J/kg) 
    1901       zcptrain(:,:) = rcp * ( SUM( (tn_ice(:,:,:) - rt0) * a_i(:,:,:), dim=3 ) + sst_m(:,:) * ziceld(:,:) )  
     1901      zcptrain(:,:) = rcp * ( SUM( (t1_ice(:,:,:) - rt0) * a_i(:,:,:), dim=3 ) + sst_m(:,:) * ziceld(:,:) )  
    19021902 
    19031903      ! --- enthalpy of snow precip over ice in J/m3 (to be used in 1D-thermo) --- ! 
     
    21402140               SELECT CASE( sn_snd_temp%clcat ) 
    21412141               CASE( 'yes' )    
    2142                   ztmp3(:,:,1:jpl) = tn_ice(:,:,1:jpl) 
     2142                  ztmp3(:,:,1:jpl) = t1_ice(:,:,1:jpl) 
    21432143               CASE( 'no' ) 
    21442144                  WHERE( SUM( a_i, dim=3 ) /= 0. ) 
    2145                      ztmp3(:,:,1) = SUM( tn_ice * a_i, dim=3 ) / SUM( a_i, dim=3 ) 
     2145                     ztmp3(:,:,1) = SUM( t1_ice * a_i, dim=3 ) / SUM( a_i, dim=3 ) 
    21462146                  ELSEWHERE 
    21472147                     ztmp3(:,:,1) = rt0 
     
    21522152               SELECT CASE( sn_snd_temp%clcat ) 
    21532153               CASE( 'yes' )    
    2154                   ztmp3(:,:,1:jpl) = tn_ice(:,:,1:jpl) * a_i(:,:,1:jpl) 
     2154                  ztmp3(:,:,1:jpl) = t1_ice(:,:,1:jpl) * a_i(:,:,1:jpl) 
    21552155               CASE( 'no' ) 
    21562156                  ztmp3(:,:,:) = 0.0 
    21572157                  DO jl=1,jpl 
    2158                      ztmp3(:,:,1) = ztmp3(:,:,1) + tn_ice(:,:,jl) * a_i(:,:,jl) 
     2158                     ztmp3(:,:,1) = ztmp3(:,:,1) + t1_ice(:,:,jl) * a_i(:,:,jl) 
    21592159                  ENDDO 
    21602160               CASE default                  ;   CALL ctl_stop( 'sbc_cpl_snd: wrong definition of sn_snd_temp%clcat' ) 
     
    21632163               SELECT CASE( sn_snd_temp%clcat )  
    21642164               CASE( 'yes' )     
    2165                   ztmp3(:,:,1:jpl) = tn_ice(:,:,1:jpl) * a_i(:,:,1:jpl)  
     2165                  ztmp3(:,:,1:jpl) = t1_ice(:,:,1:jpl) * a_i(:,:,1:jpl)  
    21662166               CASE( 'no' )  
    21672167                  ztmp3(:,:,:) = 0.0  
    21682168                  DO jl=1,jpl  
    2169                      ztmp3(:,:,1) = ztmp3(:,:,1) + tn_ice(:,:,jl) * a_i(:,:,jl)  
     2169                     ztmp3(:,:,1) = ztmp3(:,:,1) + t1_ice(:,:,jl) * a_i(:,:,jl)  
    21702170                  ENDDO  
    21712171               CASE default                  ;   CALL ctl_stop( 'sbc_cpl_snd: wrong definition of sn_snd_temp%clcat' )  
     
    21742174               ztmp1(:,:) = ( ztmp1(:,:) + rt0 ) * zfr_l(:,:)  
    21752175               DO jl=1,jpl 
    2176                   ztmp1(:,:) = ztmp1(:,:) + tn_ice(:,:,jl) * a_i(:,:,jl) 
     2176                  ztmp1(:,:) = ztmp1(:,:) + t1_ice(:,:,jl) * a_i(:,:,jl) 
    21772177               ENDDO 
    21782178            CASE default                     ;   CALL ctl_stop( 'sbc_cpl_snd: wrong definition of sn_snd_temp%cldes' ) 
Note: See TracChangeset for help on using the changeset viewer.