Changeset 13644


Ignore:
Timestamp:
2020-10-20T12:17:11+02:00 (6 months ago)
Author:
techene
Message:

#2385 come back to previous e3 weighting in ene and ens when qco not activated (SETTE pass)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/dev_r13327_KERNEL-06_2_techene_e3/src/OCE/DYN/dynvor.F90

    r13627 r13644  
    386386         END SELECT 
    387387         ! 
    388          !                                   !==  horizontal fluxes and potential vorticity ==! 
     388#if defined key_qco 
     389         !                                   !==  horizontal fluxes and potential vorticity  ==! 
    389390         zwx(:,:) = e2u(:,:) * e3u(:,:,jk,Kmm) * pu(:,:,jk) 
    390391         zwy(:,:) = e1v(:,:) * e3v(:,:,jk,Kmm) * pv(:,:,jk) 
    391392         zwz(:,:) = zwz(:,:) / e3f(:,:,jk) 
     393#else 
     394         !                                   !==  horizontal fluxes and potential vorticity  ==! 
     395         IF( ln_sco ) THEN                        ! weighted by e3 
     396            zwz(:,:) = zwz(:,:) / e3f(:,:,jk) 
     397            zwx(:,:) = e2u(:,:) * e3u(:,:,jk,Kmm) * pu(:,:,jk) 
     398            zwy(:,:) = e1v(:,:) * e3v(:,:,jk,Kmm) * pv(:,:,jk) 
     399         ELSE                                     ! ugly fix: unweighted by e3 
     400            zwx(:,:) = e2u(:,:) * pu(:,:,jk) 
     401            zwy(:,:) = e1v(:,:) * pv(:,:,jk) 
     402         ENDIF 
     403#endif 
    392404         ! 
    393405         !                                   !==  compute and add the vorticity term trend  =! 
     
    484496         ! 
    485497         ! 
    486          !                                   !==  horizontal fluxes and potential vorticity ==! 
     498#if defined key_qco 
     499         !                                   !==  horizontal fluxes and potential vorticity  ==! 
    487500         zwx(:,:) = e2u(:,:) * e3u(:,:,jk,Kmm) * pu(:,:,jk) 
    488501         zwy(:,:) = e1v(:,:) * e3v(:,:,jk,Kmm) * pv(:,:,jk) 
    489502         zwz(:,:) = zwz(:,:) / e3f(:,:,jk) 
     503#else 
     504         !                                   !==  horizontal fluxes and potential vorticity  ==! 
     505         IF( ln_sco ) THEN                        ! weighted by e3 
     506            zwz(:,:) = zwz(:,:) / e3f(:,:,jk) 
     507            zwx(:,:) = e2u(:,:) * e3u(:,:,jk,Kmm) * pu(:,:,jk) 
     508            zwy(:,:) = e1v(:,:) * e3v(:,:,jk,Kmm) * pv(:,:,jk) 
     509         ELSE                                     ! ugly fix: unweighted by e3 
     510            zwx(:,:) = e2u(:,:) * pu(:,:,jk) 
     511            zwy(:,:) = e1v(:,:) * pv(:,:,jk) 
     512         ENDIF 
     513#endif 
    490514         ! 
    491515         !                                   !==  compute and add the vorticity term trend  =! 
Note: See TracChangeset for help on using the changeset viewer.