- Timestamp:
- 2018-12-18T16:57:34+01:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2018/dev_r9947_SI3_advection/src/ICE/icevar.F90
r9943 r10411 555 555 !!------------------------------------------------------------------- 556 556 ! 557 ! 558 DO jl = 1, jpl !== loop over the categories ==! 559 ! 560 !---------------------------------------- 561 ! zap ice energy and send it to the ocean 562 !---------------------------------------- 563 DO jk = 1, nlay_i 564 DO jj = 1 , jpj 565 DO ji = 1 , jpi 566 IF( pe_i(ji,jj,jk,jl) < 0._wp .OR. pa_i(ji,jj,jl) < 0._wp ) THEN 567 hfx_res(ji,jj) = hfx_res(ji,jj) - pe_i(ji,jj,jk,jl) * r1_rdtice ! W.m-2 <0 568 pe_i(ji,jj,jk,jl) = 0._wp 569 ENDIF 570 END DO 571 END DO 572 END DO 573 ! 574 DO jk = 1, nlay_s 575 DO jj = 1 , jpj 576 DO ji = 1 , jpi 577 IF( pe_s(ji,jj,jk,jl) < 0._wp .OR. pa_i(ji,jj,jl) < 0._wp ) THEN 578 hfx_res(ji,jj) = hfx_res(ji,jj) - pe_s(ji,jj,jk,jl) * r1_rdtice ! W.m-2 <0 579 pe_s(ji,jj,jk,jl) = 0._wp 580 ENDIF 581 END DO 582 END DO 583 END DO 584 ! 585 !----------------------------------------------------- 586 ! zap ice and snow volume, add water and salt to ocean 587 !----------------------------------------------------- 588 DO jj = 1 , jpj 589 DO ji = 1 , jpi 590 IF( pv_i(ji,jj,jl) < 0._wp .OR. pa_i(ji,jj,jl) < 0._wp ) THEN 591 wfx_res(ji,jj) = wfx_res(ji,jj) + pv_i (ji,jj,jl) * rhoi * r1_rdtice 592 pv_i (ji,jj,jl) = 0._wp 593 ENDIF 594 IF( pv_s(ji,jj,jl) < 0._wp .OR. pa_i(ji,jj,jl) < 0._wp ) THEN 595 wfx_res(ji,jj) = wfx_res(ji,jj) + pv_s (ji,jj,jl) * rhos * r1_rdtice 596 pv_s (ji,jj,jl) = 0._wp 597 ENDIF 598 IF( psv_i(ji,jj,jl) < 0._wp .OR. pa_i(ji,jj,jl) < 0._wp ) THEN 599 sfx_res(ji,jj) = sfx_res(ji,jj) + psv_i(ji,jj,jl) * rhoi * r1_rdtice 600 psv_i (ji,jj,jl) = 0._wp 601 ENDIF 602 END DO 603 END DO 604 ! 605 END DO 606 ! 557 607 WHERE( pato_i(:,:) < 0._wp ) pato_i(:,:) = 0._wp 558 608 WHERE( poa_i (:,:,:) < 0._wp ) poa_i (:,:,:) = 0._wp … … 561 611 WHERE( pv_ip (:,:,:) < 0._wp ) pv_ip (:,:,:) = 0._wp ! in theory one should change wfx_pnd(-) and wfx_sum(+) 562 612 ! but it does not change conservation, so keep it this way is ok 563 !564 DO jl = 1, jpl !== loop over the categories ==!565 !566 !----------------------------------------567 ! zap ice energy and send it to the ocean568 !----------------------------------------569 DO jk = 1, nlay_i570 DO jj = 1 , jpj571 DO ji = 1 , jpi572 IF( pe_i(ji,jj,jk,jl) < 0._wp ) THEN573 hfx_res(ji,jj) = hfx_res(ji,jj) - pe_i(ji,jj,jk,jl) * r1_rdtice ! W.m-2 <0574 pe_i(ji,jj,jk,jl) = 0._wp575 ENDIF576 END DO577 END DO578 END DO579 !580 DO jk = 1, nlay_s581 DO jj = 1 , jpj582 DO ji = 1 , jpi583 IF( pe_s(ji,jj,jk,jl) < 0._wp ) THEN584 hfx_res(ji,jj) = hfx_res(ji,jj) - pe_s(ji,jj,jk,jl) * r1_rdtice ! W.m-2 <0585 pe_s(ji,jj,jk,jl) = 0._wp586 ENDIF587 END DO588 END DO589 END DO590 !591 !-----------------------------------------------------592 ! zap ice and snow volume, add water and salt to ocean593 !-----------------------------------------------------594 DO jj = 1 , jpj595 DO ji = 1 , jpi596 IF( pv_i(ji,jj,jl) < 0._wp ) THEN597 wfx_res(ji,jj) = wfx_res(ji,jj) + pv_i (ji,jj,jl) * rhoi * r1_rdtice598 pv_i (ji,jj,jl) = 0._wp599 ENDIF600 IF( pv_s(ji,jj,jl) < 0._wp ) THEN601 wfx_res(ji,jj) = wfx_res(ji,jj) + pv_s (ji,jj,jl) * rhos * r1_rdtice602 pv_s (ji,jj,jl) = 0._wp603 ENDIF604 IF( psv_i(ji,jj,jl) < 0._wp ) THEN605 sfx_res(ji,jj) = sfx_res(ji,jj) + psv_i(ji,jj,jl) * rhoi * r1_rdtice606 psv_i (ji,jj,jl) = 0._wp607 ENDIF608 END DO609 END DO610 !611 END DO612 613 ! 613 614 END SUBROUTINE ice_var_zapneg
Note: See TracChangeset
for help on using the changeset viewer.