Ignore:
Timestamp:
2018-12-18T16:57:34+01:00 (23 months ago)
Author:
clem
Message:

make all the sette tests successful

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2018/dev_r9947_SI3_advection/src/ICE/icevar.F90

    r9943 r10411  
    555555      !!------------------------------------------------------------------- 
    556556      ! 
     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      ! 
    557607      WHERE( pato_i(:,:)   < 0._wp )   pato_i(:,:)   = 0._wp 
    558608      WHERE( poa_i (:,:,:) < 0._wp )   poa_i (:,:,:) = 0._wp 
     
    561611      WHERE( pv_ip (:,:,:) < 0._wp )   pv_ip (:,:,:) = 0._wp ! in theory one should change wfx_pnd(-) and wfx_sum(+) 
    562612      !                                                        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 ocean 
    568          !---------------------------------------- 
    569          DO jk = 1, nlay_i 
    570             DO jj = 1 , jpj 
    571                DO ji = 1 , jpi 
    572                   IF( pe_i(ji,jj,jk,jl) < 0._wp ) THEN 
    573                      hfx_res(ji,jj)   = hfx_res(ji,jj) - pe_i(ji,jj,jk,jl) * r1_rdtice ! W.m-2 <0 
    574                      pe_i(ji,jj,jk,jl) = 0._wp 
    575                   ENDIF 
    576                END DO 
    577             END DO 
    578          END DO 
    579          ! 
    580          DO jk = 1, nlay_s 
    581             DO jj = 1 , jpj 
    582                DO ji = 1 , jpi 
    583                   IF( pe_s(ji,jj,jk,jl) < 0._wp ) THEN 
    584                      hfx_res(ji,jj)   = hfx_res(ji,jj) - pe_s(ji,jj,jk,jl) * r1_rdtice ! W.m-2 <0 
    585                      pe_s(ji,jj,jk,jl) = 0._wp 
    586                   ENDIF 
    587                END DO 
    588             END DO 
    589          END DO 
    590          ! 
    591          !----------------------------------------------------- 
    592          ! zap ice and snow volume, add water and salt to ocean 
    593          !----------------------------------------------------- 
    594          DO jj = 1 , jpj 
    595             DO ji = 1 , jpi 
    596               IF( pv_i(ji,jj,jl) < 0._wp ) THEN 
    597                   wfx_res(ji,jj)    = wfx_res(ji,jj) + pv_i (ji,jj,jl) * rhoi * r1_rdtice 
    598                   pv_i   (ji,jj,jl) = 0._wp 
    599                ENDIF 
    600                IF( pv_s(ji,jj,jl) < 0._wp ) THEN 
    601                   wfx_res(ji,jj)    = wfx_res(ji,jj) + pv_s (ji,jj,jl) * rhos * r1_rdtice 
    602                   pv_s   (ji,jj,jl) = 0._wp 
    603                ENDIF 
    604                IF( psv_i(ji,jj,jl) < 0._wp ) THEN 
    605                   sfx_res(ji,jj)    = sfx_res(ji,jj) + psv_i(ji,jj,jl) * rhoi * r1_rdtice 
    606                   psv_i  (ji,jj,jl) = 0._wp 
    607                ENDIF 
    608             END DO 
    609          END DO 
    610          ! 
    611       END DO  
    612613      ! 
    613614   END SUBROUTINE ice_var_zapneg 
Note: See TracChangeset for help on using the changeset viewer.