Changeset 1373 for LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/physiq.F

Timestamp:

Apr 28, 2010, 4:27:15 PM (14 years ago)

Author:

musat

Message:

Some changes for new physics :

• new formulation for wakes' alp_offset (physiq)
• bug fix for lalim in thermals plume
• add iflag_pbl=9 option (yamada4)

FH/IM

File:

• LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/physiq.F (modified) (1 diff)

LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/physiq.F

@@ -2105 +2105 @@
        else
          IF(prt_level>9)WRITE(lunout,*)'ALE et ALP couples au thermique'
-          do i = 1,klon
-              ALE(i) = max(ale_wake(i),Ale_bl(i))
-c avant        ALP(i) = alp_wake(i) + Alp_bl(i)
-              ALP(i) = alp_wake(i) + Alp_bl(i) + alp_offset ! modif sb
-c         write(20,*)'ALE',ALE(i),Ale_bl(i),ale_wake(i)
-c         write(21,*)'ALP',ALP(i),Alp_bl(i),alp_wake(i)
-          enddo
+!         do i = 1,klon
+!             ALE(i) = max(ale_wake(i),Ale_bl(i))
+! avant        ALP(i) = alp_wake(i) + Alp_bl(i)
+!             ALP(i) = alp_wake(i) + Alp_bl(i) + alp_offset ! modif sb
+!         write(20,*)'ALE',ALE(i),Ale_bl(i),ale_wake(i)
+!         write(21,*)'ALP',ALP(i),Alp_bl(i),alp_wake(i)
+!         enddo
+
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+! Modif FH 2010/04/27. Sans doute temporaire.
+! Deux options pour le alp_offset : constant si >Â 0 ou proportionnel Ãa
+! w si <0
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+       do i = 1,klon
+          ALE(i) = max(ale_wake(i),Ale_bl(i))
+          if (alp_offset>=0.) then
+            ALP(i) = alp_wake(i) + Alp_bl(i) + alp_offset ! modif sb
+          else
+            ALP(i)=alp_wake(i)+Alp_bl(i)+alp_offset*min(omega(i,6),0.)
+            if (alp(i)<0.) then
+               print*,'ALP ',alp(i),alp_wake(i)
+     s         ,Alp_bl(i),alp_offset*min(omega(i,6),0.)
+            endif
+          endif
+       enddo
+!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+
        endif
        do i=1,klon