- Timestamp:
- 2017-09-26T15:24:17+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90
r8422 r8563 1998 1998 ! we must send the surface potential temperature 1999 1999 IF( l_useCT ) THEN ; ztmp1(:,:) = eos_pt_from_ct( tsn(:,:,1,jp_tem), tsn(:,:,1,jp_sal) ) 2000 ELSE 2000 ELSE ; ztmp1(:,:) = tsn(:,:,1,jp_tem) 2001 2001 ENDIF 2002 2002 ! … … 2111 2111 SELECT CASE( sn_snd_thick%clcat ) 2112 2112 CASE( 'yes' ) 2113 ztmp3(:,:,1:jpl) = h t_i(:,:,1:jpl) * a_i(:,:,1:jpl)2114 ztmp4(:,:,1:jpl) = h t_s(:,:,1:jpl) * a_i(:,:,1:jpl)2113 ztmp3(:,:,1:jpl) = h_i(:,:,1:jpl) * a_i(:,:,1:jpl) 2114 ztmp4(:,:,1:jpl) = h_s(:,:,1:jpl) * a_i(:,:,1:jpl) 2115 2115 CASE( 'no' ) 2116 2116 ztmp3(:,:,:) = 0.0 ; ztmp4(:,:,:) = 0.0 2117 2117 DO jl=1,jpl 2118 ztmp3(:,:,1) = ztmp3(:,:,1) + h t_i(:,:,jl) * a_i(:,:,jl)2119 ztmp4(:,:,1) = ztmp4(:,:,1) + h t_s(:,:,jl) * a_i(:,:,jl)2118 ztmp3(:,:,1) = ztmp3(:,:,1) + h_i(:,:,jl) * a_i(:,:,jl) 2119 ztmp4(:,:,1) = ztmp4(:,:,1) + h_s(:,:,jl) * a_i(:,:,jl) 2120 2120 ENDDO 2121 2121 CASE default ; CALL ctl_stop( 'sbc_cpl_snd: wrong definition of sn_snd_thick%clcat' ) … … 2124 2124 SELECT CASE( sn_snd_thick%clcat ) 2125 2125 CASE( 'yes' ) 2126 ztmp3(:,:,1:jpl) = h t_i(:,:,1:jpl)2127 ztmp4(:,:,1:jpl) = h t_s(:,:,1:jpl)2126 ztmp3(:,:,1:jpl) = h_i(:,:,1:jpl) 2127 ztmp4(:,:,1:jpl) = h_s(:,:,1:jpl) 2128 2128 CASE( 'no' ) 2129 2129 WHERE( SUM( a_i, dim=3 ) /= 0. ) 2130 ztmp3(:,:,1) = SUM( h t_i * a_i, dim=3 ) / SUM( a_i, dim=3 )2131 ztmp4(:,:,1) = SUM( h t_s * a_i, dim=3 ) / SUM( a_i, dim=3 )2130 ztmp3(:,:,1) = SUM( h_i * a_i, dim=3 ) / SUM( a_i, dim=3 ) 2131 ztmp4(:,:,1) = SUM( h_s * a_i, dim=3 ) / SUM( a_i, dim=3 ) 2132 2132 ELSEWHERE 2133 2133 ztmp3(:,:,1) = 0.
Note: See TracChangeset
for help on using the changeset viewer.