Changeset 834 for trunk/NEMO/LIM_SRC_3/limthd_sal.F90
- Timestamp:
- 2008-03-07T18:11:35+01:00 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/LIM_SRC_3/limthd_sal.F90
r825 r834 1 1 MODULE limthd_sal 2 2 #if defined key_lim3 3 !!---------------------------------------------------------------------- 4 !! 'key_lim3' LIM3 sea-ice model 5 !!---------------------------------------------------------------------- 3 6 !!====================================================================== 4 7 !! *** MODULE limthd_sal *** … … 20 23 USE limvar 21 24 USE par_ice 22 USE limicepoints23 25 24 26 IMPLICIT NONE … … 60 62 !! 61 63 !! ** History : 64 !! 65 !! "Je ne suis pour l'instant qu'a 80% de ma condition, mais c'est 66 !! les 30% qui restent qui seront les plus difficiles" 67 !! E. Mpenza 68 !! 62 69 !!------------------------------------------------------------------- 63 70 !! History : … … 65 72 !! 3.0 ! 05-12 Routine rewritten for the 3-D version 66 73 !!--------------------------------------------------------------------- 74 !! 67 75 !! * Local variables 68 76 INTEGER, INTENT(in) :: & … … 74 82 75 83 REAL(wp) :: & 76 ! ds_i_gd, & !: salt loss by gravity drainage77 ! ds_i_seg, & !: salt gain by segregation of seawater at the bottom78 ! ds_i_flush, & !: summer salt loss by flushing79 zargtemp, & !: used to compute Scharzacher profile80 84 zsold, & !: old salinity 81 ! ds_i_snice, & !: increase of salinity due to snow ice82 85 zeps=1.0e-06 , & !: very small 83 86 iflush , & !: flushing (1) or not (0) … … 87 90 i_ice_switch , & !: ice thickness above a certain treshold or not 88 91 ztmelts , & !: freezing point of sea ice 89 zs_snowice , & !: salinity of forming snow ice90 92 zaaa , & !: dummy factor 91 93 zbbb , & !: dummy factor … … 95 97 REAL(wp), DIMENSION(jpij) :: & 96 98 ze_init , & !initial total enthalpy 97 ze_end , & !final total enthalpy98 99 zhiold , & 99 100 zsiold … … 157 158 END DO ! ji 158 159 END DO ! jk 159 160 ! IF (jiindex_1d .GT. 0 ) THEN161 ! WRITE(numout,*) ' s_i_b : ', s_i_b(jiindex_1d,1:nlay_i)162 ! WRITE(numout,*) ' t_i_b : ', t_i_b(jiindex_1d,1:nlay_i)163 ! WRITE(numout,*) ' q_i_b : ', q_i_b(jiindex_1d,1:nlay_i)164 ! ENDIF165 160 166 161 DO ji = kideb, kiut … … 215 210 CALL lim_var_salprof1d(kideb,kiut) 216 211 217 !--------------------------------------218 ! Energy of melting for each ice layer219 !--------------------------------------220 ! q=q(S,T) J.m-3221 ! DO jk = 1, nlay_i222 ! DO ji = kideb, kiut223 ! ! New heat content after desalination224 ! ztmelts = - tmut * s_i_b(ji,jk) + rtt ! Tfi in K225 ! q_i_b(ji,jk) = rhoic * &226 ! ( cpic * ( ztmelts-t_i_b(ji,jk) ) &227 ! + lfus * ( 1.0 - (ztmelts-rtt) / &228 ! MIN( ( t_i_b(ji,jk) - rtt ) , -zeps ) ) &229 ! - rcp * ( ztmelts-rtt ) )230 ! END DO ! ji231 ! END DO232 233 !--------------------234 ! Total heat content235 !--------------------236 ! ze_end(:) = 0.0237 ! DO jk = 1, nlay_i238 ! DO ji = kideb, kiut239 ! ze_end(ji) = ze_end(ji) + q_i_b(ji,jk) * ht_i_b(ji) / nlay_i240 ! END DO241 ! END DO242 243 ! IF (jiindex_1d .GT. 0 ) THEN244 ! WRITE(numout,*) ' ze_init : ', ze_init(jiindex_1d) / rdt_ice245 ! WRITE(numout,*) ' ze_end : ', ze_end (jiindex_1d) / rdt_ice246 ! WRITE(numout,*) ' s_i_b : ', s_i_b(jiindex_1d,1:nlay_i)247 ! WRITE(numout,*) ' t_i_b : ', t_i_b(jiindex_1d,1:nlay_i)248 ! WRITE(numout,*) ' q_i_b : ', q_i_b(jiindex_1d,1:nlay_i)249 ! ENDIF250 251 212 !---------------------------- 252 213 ! Heat flux - brine drainage 253 214 !---------------------------- 254 255 ! WRITE(numout,*) ' jiindex, jjindex : ', jiindex, jjindex256 ! WRITE(numout,*) ' jiindex_1d : ', jiindex_1d257 ! WRITE(numout,*) ' kideb, kiut: ', kideb, kiut258 215 259 216 DO ji = kideb, kiut … … 265 222 iaccrbo = MAX( 0.0 , SIGN ( 1.0 , dh_i_bott(ji) ) ) 266 223 267 ! IF ( ( iflush .EQ. 1 ) .OR. ( iaccrbo .EQ. 1 ) ) THEN268 ! fhbri_1d(ji) = fhbri_1d(ji) + ( ze_end(ji) - ze_init(ji) ) * &269 ! a_i_b(ji) / at_i_b(ji) / rdt_ice270 ! ze_end GT ze_init... ?271 ! ENDIF272 224 fhbri_1d(ji) = 0.0 273 274 !++++++ 275 ! zji = MOD( npb(ji) - 1, jpi ) + 1 276 ! zjj = ( npb(ji) - 1 ) / jpi + 1 277 ! fhbricat(zji,zjj,jl) = ( ze_end(ji) - ze_init(ji) ) / & 278 ! rdt_ice 279 !++++++ 280 ! ftotal_fin(ji) = ftotal_fin(ji) + fhbricat(zji,zjj,jl) 281 ! IF ( ji .EQ. jiindex_1d) WRITE(numout,*) ' fhbri : ', & 282 ! fhbricat(zji,zjj,jl) * rdt_ice 283 ! IF ( ji .EQ. jiindex_1d ) & 284 ! WRITE(numout,*) ' ftotal_fin : ', ftotal_fin(jiindex_1d) 285 ! !+++++ 286 287 END DO ! ji 288 289 ! IF (jiindex_1d .GT. 0 ) THEN 290 ! WRITE(numout,*) ' Category : ', jl 291 ! WRITE(numout,*) ' jiindex_1d : ', jiindex_1d 292 ! WRITE(numout,*) ' fhbricat : ', fhbricat(jiindex,jjindex,jl) 293 ! WRITE(numout,*) ' fhbri_1d : ', fhbri_1d(jiindex_1d) 294 ! WRITE(numout,*) ' ze_end : ', ze_end(jiindex_1d) 295 ! WRITE(numout,*) ' ze_init : ', ze_init(jiindex_1d) 296 ! WRITE(numout,*) ' a_i_b : ', a_i_b(jiindex_1d) 297 ! WRITE(numout,*) ' at_i_b : ', at_i_b(jiindex_1d) 298 ! WRITE(numout,*) ' rdt_ice : ', rdt_ice 299 ! ENDIF 225 END DO ! ji 300 226 301 227 !---------------------------- … … 377 303 ENDIF 378 304 379 !!--multiyear sea-ice salinity380 !!-- ca va pas alors on vire bordel de djeuille381 ! if (ht.gt.360) then382 ! do jk = 1, nlay_i383 ! s_i_b(ji,jk) = 1.58 + 0.18*ht_i_b(ji)384 ! end do385 ! endif386 305 DO jk = 1, nlay_i 387 306 s_i_b(ji,jk) = sm_i_b(ji) … … 415 334 END DO ! ji 416 335 ENDIF 417 418 !+++++419 ! IF ( (zji.EQ.jiindex) .AND. (zjj.EQ.jjindex) ) THEN420 ! WRITE(numout,*) ' *** FS_EQV *** '421 ! WRITE(numout,*) ' limthd_sal '422 ! WRITE(numout,*) ' fseqv ', fseqv_1d(ji)423 ! WRITE(numout,*) ' sss_io ', sss_io(zji,zjj)424 ! WRITE(numout,*) ' s_i_new ', s_i_new(ji)425 ! WRITE(numout,*) ' dh_i_bott ', dh_i_bott(ji)426 ! WRITE(numout,*) ' a_i_b ', a_i_b(ji)427 ! WRITE(numout,*) ' *** FS_BRI *** '428 ! WRITE(numout,*) ' fsbri ', fsbri_1d(ji)429 ! WRITE(numout,*) ' ht_i ', ht_i_b(ji)430 ! WRITE(numout,*) ' dsm_i_gd ', dsm_i_gd_1d(ji)431 ! WRITE(numout,*) ' dsm_i_fl ', dsm_i_fl_1d(ji)432 ! ENDIF433 434 !+++++435 336 436 337 !-- End of salinity computations
Note: See TracChangeset
for help on using the changeset viewer.