limwri_2.F90: wrong average for EVP (NEMO_3.3)

[molines]

I think that when using LIM2_EVP (standard) or LIM2_VP (option) the way the average is done should not be the same, when wriing u_ice or v_ice. Therefore I suggest to have something like

            zcmo(ji,jj,4)  = frld  (ji,jj)
            zcmo(ji,jj,5)  = sist  (ji,jj)
            zcmo(ji,jj,6)  = fbif  (ji,jj)
           IF (lk_lim2_vp) THEN
            zcmo(ji,jj,7)  = zindb * (  u_ice(ji,jj  ) * tmu(ji,jj  ) + u_ice(ji+1,jj  ) * tmu(ji+1,jj  )   &
                                      + u_ice(ji,jj+1) * tmu(ji,jj+1) + u_ice(ji+1,jj+1) * tmu(ji+1,jj+1) ) &
                                  / ztmu 

            zcmo(ji,jj,8)  = zindb * (  v_ice(ji,jj  ) * tmu(ji,jj  ) + v_ice(ji+1,jj  ) * tmu(ji+1,jj  )   &
                                      + v_ice(ji,jj+1) * tmu(ji,jj+1) + v_ice(ji+1,jj+1) * tmu(ji+1,jj+1) ) &
                                  / ztmu
           ELSE

            zcmo(ji,jj,7)  = zindb * (  u_ice(ji,jj  ) * tmu(ji,jj)                       &
             &                        + u_ice(ji-1,jj) * tmu(ji-1,jj) )                   &
             &                    / 2.0
            zcmo(ji,jj,8)  = zindb * (  v_ice(ji,jj  ) * tmv(ji,jj)                       &
             &                        + v_ice(ji,jj-1) * tmv(ji,jj-1) )                   &
             &                    / 2.0

           ENDIF
            zcmo(ji,jj,9)  = sst_m(ji,jj)
            zcmo(ji,jj,10) = sss_m(ji,jj)
            zcmo(ji,jj,11) = qns(ji,jj) + qsr(ji,jj)

instead of zcmo(ji,jj,4) = frld (ji,jj)

zcmo(ji,jj,5) = sist (ji,jj)
zcmo(ji,jj,6) = fbif (ji,jj)
zcmo(ji,jj,7) = zindb * ( u_ice(ji,jj ) * tmu(ji,jj ) + u_ice(ji+1,jj ) * tmu(ji+1,jj ) &

+ u_ice(ji,jj+1) * tmu(ji,jj+1) + u_ice(ji+1,jj+1) * tmu(ji+1,jj+1) ) &

/ ztmu

zcmo(ji,jj,8) = zindb * ( v_ice(ji,jj ) * tmu(ji,jj ) + v_ice(ji+1,jj ) * tmu(ji+1,jj ) &

+ v_ice(ji,jj+1) * tmu(ji,jj+1) + v_ice(ji+1,jj+1) * tmu(ji+1,jj+1) ) &

/ ztmu

zcmo(ji,jj,9) = sst_m(ji,jj)
zcmo(ji,jj,10) = sss_m(ji,jj)
zcmo(ji,jj,11) = qns(ji,jj) + qsr(ji,jj)

}}}

Of course the same apply to the counter part in limwri_dimg_2.F90

[smasson]

We keep this piece of of code only for compatibility (or as "souvenir" of the old time?). It doesn't exist anymore when key_iomput is activated (the default).
When rewriting outputs with iom_put, I must admit that I didn't really understood the idea/the need/the pertinence of this weighted interpolation of u/v_ice by zindb (with its incredible definition in a pure LIM2 style) and tmu...
With iom_put we decided to output the u/v_ice on its original grid (knowing that, soon, we will switch to VP on C grid). Why should we output ice velocity on T point whereas we never do it for the ocean? We decided to only weight these fields by the "ice presence" variable, i.e. for example if you have ice only half of the time, its speed is not divided by 2.

btw, I will not be able to accept you correction in the current form (with the horrible "/2.") the Master would kill me :-)

[rblod]

Considering your constance in this request, I've fixed it in 3.4 trunk for next 3.4.1, despite Sebastien opposition (but with the agreement of the Master). Note you may so be a source of tension inside the NEMO team Paris . Cheers. Rachid