- Timestamp:
- Aug 2, 2024, 2:12:03 PM (3 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/dyn3d_common/traceurpole.f90
r5136 r5158 30 30 sommemasses = 0 31 31 sommemqs = 0 32 dol = 1, llm33 doi = 1, iip132 DO l = 1, llm 33 DO i = 1, iip1 34 34 sommemasses(l) = sommemasses(l) + masse(i, jjp1, l) 35 35 sommemqs(l) = sommemqs(l) + masse(i, jjp1, l) * q(i, jjp1, l) … … 41 41 sommemassen = 0 42 42 sommemqn = 0 43 dol = 1, llm44 doi = 1, iip143 DO l = 1, llm 44 DO i = 1, iip1 45 45 sommemassen(l) = sommemassen(l) + masse(i, 1, l) 46 46 sommemqn(l) = sommemqn(l) + masse(i, 1, l) * q(i, 1, l) … … 50 50 51 51 ! On force le traceur à prendre cette valeur aux pôles 52 dol = 1, llm53 doi = 1, iip152 DO l = 1, llm 53 DO i = 1, iip1 54 54 q(i, 1, l) = qpolen(l) 55 55 q(i, jjp1, l) = qpoles(l)
Note: See TracChangeset
for help on using the changeset viewer.