- Timestamp:
- 2017-10-04T20:08:33+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90
r8563 r8592 32 32 USE geo2ocean ! 33 33 USE oce , ONLY : tsn, un, vn, sshn, ub, vb, sshb, fraqsr_1lev 34 USE albedooce!34 USE ocealb ! 35 35 USE eosbn2 ! 36 36 USE sbcrnf, ONLY : l_rnfcpl … … 177 177 TYPE( DYNARR ), SAVE, DIMENSION(jprcv) :: frcv ! all fields recieved from the atmosphere 178 178 179 REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: alb edo_oce_mix ! ocean albedo sent to atmosphere (mix clear/overcast sky)179 REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: alb_oce_mix ! ocean albedo sent to atmosphere (mix clear/overcast sky) 180 180 181 181 REAL(wp) :: rpref = 101000._wp ! reference atmospheric pressure[N/m2] … … 201 201 ierr(:) = 0 202 202 ! 203 ALLOCATE( alb edo_oce_mix(jpi,jpj), nrcvinfo(jprcv), STAT=ierr(1) )203 ALLOCATE( alb_oce_mix(jpi,jpj), nrcvinfo(jprcv), STAT=ierr(1) ) 204 204 205 205 #if ! defined key_lim3 && ! defined key_cice … … 735 735 ! 2. receiving mixed oce-ice solar radiation 736 736 IF ( TRIM ( sn_snd_alb%cldes ) == 'mixed oce-ice' .OR. TRIM ( sn_rcv_qsr%cldes ) == 'mixed oce-ice' ) THEN 737 CALL albedo_oce( zaos, zacs )737 CALL oce_alb( zaos, zacs ) 738 738 ! Due to lack of information on nebulosity : mean clear/overcast sky 739 alb edo_oce_mix(:,:) = ( zacs(:,:) + zaos(:,:) ) * 0.5739 alb_oce_mix(:,:) = ( zacs(:,:) + zaos(:,:) ) * 0.5 740 740 ENDIF 741 741 … … 1884 1884 ! ( see OASIS3 user guide, 5th edition, p39 ) 1885 1885 zqsr_ice(:,:,1) = frcv(jpr_qsrmix)%z3(:,:,1) * ( 1.- palbi(:,:,1) ) & 1886 & / ( 1.- ( alb edo_oce_mix(:,: ) * ziceld(:,:) &1886 & / ( 1.- ( alb_oce_mix(:,: ) * ziceld(:,:) & 1887 1887 & + palbi (:,:,1) * picefr(:,:) ) ) 1888 1888 END SELECT … … 2051 2051 ztmp1(:,:) = SUM( alb_ice (:,:,1:jpl) * a_i(:,:,1:jpl), dim=3 ) / SUM( a_i(:,:,1:jpl), dim=3 ) 2052 2052 ELSEWHERE 2053 ztmp1(:,:) = alb edo_oce_mix(:,:)2053 ztmp1(:,:) = alb_oce_mix(:,:) 2054 2054 END WHERE 2055 2055 CASE default ; CALL ctl_stop( 'sbc_cpl_snd: wrong definition of sn_snd_alb%clcat' ) … … 2079 2079 2080 2080 IF( ssnd(jps_albmix)%laction ) THEN ! mixed ice-ocean 2081 ztmp1(:,:) = alb edo_oce_mix(:,:) * zfr_l(:,:)2081 ztmp1(:,:) = alb_oce_mix(:,:) * zfr_l(:,:) 2082 2082 DO jl=1,jpl 2083 2083 ztmp1(:,:) = ztmp1(:,:) + alb_ice(:,:,jl) * a_i(:,:,jl)
Note: See TracChangeset
for help on using the changeset viewer.