Changeset 2098 for trunk/LMDZ.TITAN/libf/phytitan/physiq_mod.F90

Timestamp:

Feb 11, 2019, 10:31:18 PM (6 years ago)

Author:

jvatant

Message:

Fix forgotten do loop in previous commit
--JVO

File:

• trunk/LMDZ.TITAN/libf/phytitan/physiq_mod.F90 (modified) (1 diff)

trunk/LMDZ.TITAN/libf/phytitan/physiq_mod.F90

@@ -765 +765 @@
       else !  if not moyzon
       
-         DO l=1,nlayer   
-           zzlay_eff(ig,l)=g*rad*rad/(g*rad-(pphi(ig,l)+phisfi(ig)))-rad ! reference = geoid
-         ENDDO
-         zzlev_eff(ig,1)=phisfi(ig)/g
-         DO l=2,nlayer
+        DO ig=1,ngrid
+          DO l=1,nlayer   
+            zzlay_eff(ig,l)=g*rad*rad/(g*rad-(pphi(ig,l)+phisfi(ig)))-rad ! reference = geoid
+          ENDDO
+          zzlev_eff(ig,1)=phisfi(ig)/g
+          DO l=2,nlayer
             z1=(pplay(ig,l-1)+pplev(ig,l))/ (pplay(ig,l-1)-pplev(ig,l))
             z2=(pplev(ig,l)  +pplay(ig,l))/(pplev(ig,l)  -pplay(ig,l))
             zzlev_eff(ig,l)=(z1*zzlay_eff(ig,l-1)+z2*zzlay_eff(ig,l))/(z1+z2)
-         ENDDO
-         zzlev_eff(ig,nlayer+1)=zzlay_eff(ig,nlayer)+(zzlay_eff(ig,nlayer)-zzlev_eff(ig,nlayer))
+          ENDDO
+          zzlev_eff(ig,nlayer+1)=zzlay_eff(ig,nlayer)+(zzlay_eff(ig,nlayer)-zzlev_eff(ig,nlayer))
+        ENDDO
 
       endif  ! moyzon