Changeset 914 for trunk/NEMO
- Timestamp:
- 2008-04-30T16:30:23+02:00 (16 years ago)
- Location:
- trunk/NEMO/OPA_SRC/SBC
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/SBC/sbcblk_clio.F90
r892 r914 270 270 REAL(wp) :: zdtetar, ztvmoyr, zlxins, zchcm, zclcm ! - - 271 271 REAL(wp) :: zmt1, zmt2, zmt3, ztatm3, ztamr, ztaevbk ! - - 272 REAL(wp) :: zsst, ztatm, zcco1, zpatm , zinda! - -272 REAL(wp) :: zsst, ztatm, zcco1, zpatm ! - - 273 273 REAL(wp) :: zrhoa, zev, zes, zeso, zqatm, zevsqr ! - - 274 274 !! … … 296 296 297 297 CALL blk_clio_qsr_oce( qsr ) 298 299 ! CAUTION: ocean shortwave radiation sets to zero if more than 50% of sea-ice !!gm to be removed300 DO jj = 1, jpj301 DO ji = 1, jpi302 zinda = MAX( 0.e0, SIGN( 1.e0, -( -1.5 - freeze(ji,jj) ) ) )303 qsr(ji,jj) = zinda * qsr(ji,jj)304 END DO305 END DO306 307 298 308 299 !------------------------! -
trunk/NEMO/OPA_SRC/SBC/sbcice_lim.F90
r906 r914 52 52 USE in_out_manager ! I/O manager 53 53 USE prtctl ! Print control 54 USE ocfzpt ! ocean freezing point 54 55 55 56 IMPLICIT NONE … … 99 100 INTEGER :: alert_id ! number of the current alert 100 101 REAL(wp) :: ztmelts ! ice layer melting point 101 INTEGER , DIMENSION(20) :: numal ! number of alerts positive 102 REAL(wp) :: zinda 103 INTEGER , DIMENSION(20) :: numal ! number of alerts positive 102 104 CHARACTER (len=30), DIMENSION(20) :: alname ! name of alert 103 105 REAL(wp), DIMENSION(jpi,jpj,jpl) :: alb_ice_os ! albedo of the ice under overcast sky … … 159 161 & tprecip , sprecip , & 160 162 & fr1_i0 , fr2_i0 , cl_grid ) 163 164 ! CAUTION: ocean shortwave radiation sets to zero if more than 50% of sea-ice !!gm to be removed 165 DO jj = 1, jpj 166 DO ji = 1, jpi 167 zinda = MAX( 0.e0, SIGN( 1.e0, -( -1.5 - freeze(ji,jj) ) ) ) 168 qsr(ji,jj) = zinda * qsr(ji,jj) 169 END DO 170 END DO 171 161 172 CASE( 4 ) ! CORE bulk formulation 162 173 CALL blk_ice_core( t_su , u_ice , v_ice , alb_ice_cs, & -
trunk/NEMO/OPA_SRC/SBC/sbcice_lim_2.F90
r900 r914 46 46 USE in_out_manager ! I/O manager 47 47 USE prtctl ! Print control 48 USE ocfzpt ! ocean freezing point 48 49 49 50 IMPLICIT NONE … … 89 90 !! 90 91 INTEGER :: ji, jj ! dummy loop indices 92 REAL(wp) :: zinda 91 93 REAL(wp), DIMENSION(jpi,jpj,1) :: alb_ice_os ! albedo of the ice under overcast sky 92 94 REAL(wp), DIMENSION(jpi,jpj,1) :: alb_ice_cs ! albedo of ice under clear sky … … 152 154 & tprecip , sprecip , & 153 155 & fr1_i0 , fr2_i0 , cl_grid ) 156 157 ! CAUTION: ocean shortwave radiation sets to zero if more than 50% of sea-ice !!gm to be removed 158 DO jj = 1, jpj 159 DO ji = 1, jpi 160 zinda = MAX( 0.e0, SIGN( 1.e0, -( -1.5 - freeze(ji,jj) ) ) ) 161 qsr(ji,jj) = zinda * qsr(ji,jj) 162 END DO 163 END DO 164 154 165 CASE( 4 ) ! CORE bulk formulation 155 166 CALL blk_ice_core( zsist , ui_ice , vi_ice , alb_ice_cs , &
Note: See TracChangeset
for help on using the changeset viewer.