Changeset 8374 for branches/2017/dev_r8183_ICEMODEL/NEMOGCM
- Timestamp:
- 2017-07-25T19:48:56+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/limsbc.F90
r8341 r8374 107 107 REAL(wp) :: zqmass ! Heat flux associated with mass exchange ice->ocean (W.m-2) 108 108 REAL(wp) :: zqsr ! New solar flux received by the ocean 109 REAL(wp), POINTER, DIMENSION(:,:,:) :: zalb_cs, zalb_os ! 3D workspace110 REAL(wp), POINTER, DIMENSION(:,:):: zalb ! 2D workspace109 REAL(wp), DIMENSION(jpi,jpj,jpl) :: zalb_cs, zalb_os ! 3D workspace 110 REAL(wp), DIMENSION(jpi,jpj) :: zalb ! 2D workspace 111 111 !!--------------------------------------------------------------------- 112 112 … … 122 122 123 123 ! albedo output 124 CALL wrk_alloc( jpi,jpj, zalb )125 126 124 zalb(:,:) = 0._wp 127 125 WHERE ( at_i_b <= epsi06 ) ; zalb(:,:) = rn_alb_oce … … 133 131 IF( iom_use('albedo' ) ) CALL iom_put( "albedo" , zalb(:,:) ) ! surface albedo output 134 132 135 CALL wrk_dealloc( jpi,jpj, zalb )136 133 137 134 DO jj = 1, jpj … … 217 214 ! Snow/ice albedo (only if sent to coupler, useless in forced mode) ! 218 215 !------------------------------------------------------------------------! 219 CALL wrk_alloc( jpi,jpj,jpl, zalb_cs, zalb_os )220 216 CALL albedo_ice( t_su, ht_i, ht_s, a_ip_frac, h_ip, ln_pnd_rad, zalb_cs, zalb_os ) ! cloud-sky and overcast-sky ice albedos 221 217 222 218 alb_ice(:,:,:) = ( 1. - cldf_ice ) * zalb_cs(:,:,:) + cldf_ice * zalb_os(:,:,:) 223 CALL wrk_dealloc( jpi,jpj,jpl, zalb_cs, zalb_os )224 219 225 220 ! conservation test
Note: See TracChangeset
for help on using the changeset viewer.