Changeset 4105 for LMDZ6/trunk/libf/phylmd

Timestamp:

Mar 17, 2022, 11:24:22 PM (3 years ago)

Author:

jyg

Message:

Bug fix in dyn1d/old_lmdz1d.F90: compute the geopotential before it is used

File:

• LMDZ6/trunk/libf/phylmd/dyn1d/old_lmdz1d.F90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmd/dyn1d/old_lmdz1d.F90

@@ -1050 +1050 @@
 
 !---------------------------------------------------------------------
-! Interpolation of forcings in time and onto model levels
-!---------------------------------------------------------------------
-
-#include "old_1D_interp_cases.h"
-
-      if (forcing_GCM2SCM) then
-        write (*,*) 'forcing_GCM2SCM not yet implemented'
-        stop 'in time loop'
-      endif ! forcing_GCM2SCM
-
-!---------------------------------------------------------------------
 !  Geopotential :
 !---------------------------------------------------------------------
@@ -1069 +1058 @@
      &    (play(l)-play(l+1))/(play(l)+play(l+1))
         enddo
+
+!---------------------------------------------------------------------
+! Interpolation of forcings in time and onto model levels
+!---------------------------------------------------------------------
+
+#include "old_1D_interp_cases.h"
+
+      if (forcing_GCM2SCM) then
+        write (*,*) 'forcing_GCM2SCM not yet implemented'
+        stop 'in time loop'
+      endif ! forcing_GCM2SCM
+
+!!!!---------------------------------------------------------------------
+!!!!  Geopotential :
+!!!!---------------------------------------------------------------------
+!!!
+!!!        phi(1)=RD*temp(1)*(plev(1)-play(1))/(.5*(plev(1)+play(1)))
+!!!        do l = 1, llm-1
+!!!          phi(l+1)=phi(l)+RD*(temp(l)+temp(l+1))*                           &
+!!!     &    (play(l)-play(l+1))/(play(l)+play(l+1))
+!!!        enddo
 
 !---------------------------------------------------------------------