Changeset 5087 for LMDZ6/branches/Amaury_dev/libf/phylmd/yamada4.F90
- Timestamp:
- Jul 20, 2024, 12:00:23 PM (4 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/yamada4.F90
r5082 r5087 393 393 tkeprov=q2(ig,k)/ydeux 394 394 tkeprov= tkeprov* & 395 &(tkeprov+dt*(shear(ig,k)+max(0.,buoy(ig,k))))/ &396 &(tkeprov+dt*((-1.)*min(0.,buoy(ig,k))+dissip(ig,k)+drgpro(ig,k)*tkeprov))395 (tkeprov+dt*(shear(ig,k)+max(0.,buoy(ig,k))))/ & 396 (tkeprov+dt*((-1.)*min(0.,buoy(ig,k))+dissip(ig,k)+drgpro(ig,k)*tkeprov)) 397 397 q2(ig,k)=tkeprov*ydeux 398 398 ENDDO … … 430 430 tkeprov= tkeprov+dt*(shear(ig,k)+max(0.,buoy(ig,k))) 431 431 tkeprov= tkeprov* & 432 &tkeprov/ &433 &(tkeprov+dt*((-1.)*min(0.,buoy(ig,k))+dissip(ig,k)))432 tkeprov/ & 433 (tkeprov+dt*((-1.)*min(0.,buoy(ig,k))+dissip(ig,k))) 434 434 q2(ig,k)=tkeprov*ydeux 435 435 ! En cas stable, on traite la flotabilite comme la … … 453 453 winds(ig,k)=sqrt(u(ig,k)**2+v(ig,k)**2) 454 454 tkeprov= (shear(ig,k)+ & 455 &drgpro(ig,k)*(winds(ig,k))**3)*tkeprov/disseff*(1.-tkeexp)+tkeprov*tkeexp455 drgpro(ig,k)*(winds(ig,k))**3)*tkeprov/disseff*(1.-tkeexp)+tkeprov*tkeexp 456 456 q2(ig,k)=tkeprov*ydeux 457 457 ! En cas stable, on traite la flotabilite comme la
Note: See TracChangeset
for help on using the changeset viewer.