Changeset 1530 for trunk/LMDZ.TITAN/libf/phytitan/ballon.F

Timestamp:

Apr 6, 2016, 8:51:34 AM (9 years ago)

Author:

emillour

Message:

Venus and Titan GCMs:
Updates in the physics to keep up with updates in LMDZ5 (up to
LMDZ5 trunk, rev 2350) concerning dynamics/physics separation:

• Adapted makelmdz and makelmdz_fcm script to stop if trying to compile 1d model or newstart or start2archive in parallel.
• got rid of references to "dimensions.h" in physics. Within physics packages, use nbp_lon (=iim), nbp_lat (=jjmp1) and nbp_lev (=llm) from module mod_grid_phy_lmdz (in phy_common) instead. Only partially done for Titan, because of many hard-coded commons; a necessary first step will be to clean these up (using modules).

EM

File:

• trunk/LMDZ.TITAN/libf/phytitan/ballon.F (modified) (13 diffs)

trunk/LMDZ.TITAN/libf/phytitan/ballon.F

@@ -3 +3 @@
 
       use dimphy
+      use mod_grid_phy_lmdz, only: nbp_lon, nbp_lat
       implicit none
 
@@ -94 +95 @@
 c======================================================================
 
-#include "dimensions.h"
 #include "YOMCST.h"
 c
@@ -110 +110 @@
       INTEGER jj,ii,ll
 
-      REAL zlon(iim+1),zlat(jjm+1)
-      save zlon,zlat
+      REAL,SAVE,ALLOCATABLE :: zlon(:),zlat(:)
 
       REAL time
@@ -141 +140 @@
       print*,"BALLOONS ACTIVATED"
 
+      allocate(zlon(nbp_lon+1))
+      allocate(zlat(nbp_lat))
+
 C Latitudes:
       zlat(1)=plat(1)*RPI/180.
-      do j = 2,jjm
-         k=(j-2)*iim+2
+      do j = 2,nbp_lat-1
+         k=(j-2)*nbp_lon+2
          zlat(j)=plat(k)*RPI/180.
       enddo
-      zlat(jjm+1)=plat(klon)*RPI/180.
+      zlat(nbp_lat)=plat(klon)*RPI/180.
 
 C Longitudes:
-      do i = 1,iim
+      do i = 1,nbp_lon
          k=i+1
          zlon(i)=plon(k)*RPI/180.
       enddo
-      zlon(iim+1)=zlon(1)+2.*RPI
+      zlon(nbp_lon+1)=zlon(1)+2.*RPI
 
 c verif init     lat de 90 à -90, lon de -180 à 180
@@ -256 +258 @@
 
       use dimphy
+      use mod_grid_phy_lmdz, only: nbp_lon, nbp_lat
       implicit none
 
@@ -283 +286 @@
 c======================================================================
 
-#include "dimensions.h"
 #include "YOMCST.h"
 c
@@ -289 +291 @@
 c
       real map_u(klon),map_v(klon),map_a(klon)
-      real latit(jjm+1),longit(iim)
+      real latit(nbp_lat),longit(nbp_lon)
       real phi,lam,ubal,vbal,abal
 c
@@ -296 +298 @@
       INTEGER i,j,k
       INTEGER jj,ii
-      REAL    ujj(iim+1),vjj(iim+1),ajj(iim+1)
+      REAL    ujj(nbp_lon+1),vjj(nbp_lon+1),ajj(nbp_lon+1)
       REAL    factlat,factlon
 
@@ -302 +304 @@
 c-------------------------------------------------
         jj=1  ! POLE NORD 
-        do j=2,jjm
+        do j=2,nbp_lat-1
           if (phi.lt.latit(j)) jj=j
         enddo
@@ -309 +311 @@
 c pole nord
         if (jj.eq.1) then
-         do i=1,iim
+         do i=1,nbp_lon
           ujj(i) = map_u(i+1)*factlat + map_u(1)*(1-factlat)
           vjj(i) = map_v(i+1)*factlat + map_v(1)*(1-factlat)
@@ -315 +317 @@
          enddo
 c pole sud
-        elseif (jj.eq.jjm) then
-         do i=1,iim
-          k = (jj-2)*iim+1+i
+        elseif (jj.eq.nbp_lat-1) then
+         do i=1,nbp_lon
+          k = (jj-2)*nbp_lon+1+i
           ujj(i) = map_u(klon)*factlat + map_u(k)*(1-factlat)
           vjj(i) = map_v(klon)*factlat + map_v(k)*(1-factlat)
@@ -324 +326 @@
 c autres latitudes
         else
-         do i=1,iim
-          k = (jj-2)*iim+1+i
-          ujj(i) = map_u(k+iim)*factlat + map_u(k)*(1-factlat)
-          vjj(i) = map_v(k+iim)*factlat + map_v(k)*(1-factlat)
-          ajj(i) = map_a(k+iim)*factlat + map_a(k)*(1-factlat)
+         do i=1,nbp_lon
+          k = (jj-2)*nbp_lon+1+i
+          ujj(i) = map_u(k+nbp_lon)*factlat + map_u(k)*(1-factlat)
+          vjj(i) = map_v(k+nbp_lon)*factlat + map_v(k)*(1-factlat)
+          ajj(i) = map_a(k+nbp_lon)*factlat + map_a(k)*(1-factlat)
          enddo
         endif
-        ujj(iim+1)=ujj(1)
-        vjj(iim+1)=vjj(1)
-        ajj(iim+1)=ajj(1)
+        ujj(nbp_lon+1)=ujj(1)
+        vjj(nbp_lon+1)=vjj(1)
+        ajj(nbp_lon+1)=ajj(1)
 
 c Interpolation in longitudes
 c-------------------------------------------------
         ii=1  ! lon=-180
-        do i=2,iim
+        do i=2,nbp_lon
           if (lam.gt.longit(i)) ii=i
         enddo
@@ -370 +372 @@
 c======================================================================
 
-#include "dimensions.h"
 #include "YOMCST.h"
 c