Changeset 2346 for LMDZ5/trunk/libf/phylmd/aaam_bud.F90

Timestamp:

Aug 21, 2015, 5:13:46 PM (9 years ago)

Author:

Ehouarn Millour

Message:

Physics/dynamics separation:

• remove all references to dimensions.h from physics. nbp_lon (==iim) , nbp_lat (==jjm+1) and nbp_lev (==llm) from mod_grid_phy_lmdz should be used instead.
• added module regular_lonlat_mod in phy_common to store information about the global (lon-lat) grid cell boundaries and centers.

EM

File:

• LMDZ5/trunk/libf/phylmd/aaam_bud.F90 (modified) (13 diffs)

LMDZ5/trunk/libf/phylmd/aaam_bud.F90

@@ -6 +6 @@
 
   USE dimphy
+  USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat
   IMPLICIT NONE
   ! ======================================================================
@@ -58 +59 @@
   ! ===================
 
-  ! iim--common-I: Number of longitude intervals
-  ! jjm--common-I: Number of latitude intervals
+  ! nbp_lon--common-I: Number of longitude intervals
+  ! (nbp_lat-1)--common-I: Number of latitude intervals
   ! klon-common-I: Number of points seen by the physics
-  ! iim*(jjm-1)+2 for instance
+  ! nbp_lon*(nbp_lat-2)+2 for instance
   ! klev-common-I: Number of vertical layers
   ! ======================================================================
@@ -84 +85 @@
   ! ======================================================================
 
-  include "dimensions.h"
-  ! cc#include "dimphy.h"
-
   ! ARGUMENTS
 
@@ -121 +119 @@
   ! PUT AAM QUANTITIES AT ZERO:
 
-  IF (iim+1>801 .OR. jjm+1>401) THEN
+  IF (nbp_lon+1>801 .OR. nbp_lat>401) THEN
     abort_message = 'Pb de dimension dans aaam_bud'
     CALL abort_physic(modname, abort_message, 1)
@@ -129 +127 @@
   hadley = 1.E18
   hadday = 1.E18*24.*3600.
-  dlat = xpi/real(jjm)
-  dlon = 2.*xpi/real(iim)
+  dlat = xpi/real(nbp_lat-1)
+  dlon = 2.*xpi/real(nbp_lon)
 
   DO iax = 1, 3
@@ -155 +153 @@
   zlat(1) = plat(l)*xpi/180.
 
-  DO i = 1, iim + 1
+  DO i = 1, nbp_lon + 1
 
     zs(i, 1) = phis(l)/rg
@@ -169 +167 @@
 
 
-  DO j = 2, jjm
+  DO j = 2, nbp_lat-1
 
     ! Values at Greenwich (Periodicity)
 
-    zs(iim+1, j) = phis(l+1)/rg
-    ps(iim+1, j) = p(l+1, 1)
-    ssou(iim+1, j) = dragu(l+1) + liftu(l+1)
-    ssov(iim+1, j) = dragv(l+1) + liftv(l+1)
-    blsu(iim+1, j) = phyu(l+1) - dragu(l+1) - liftu(l+1)
-    blsv(iim+1, j) = phyv(l+1) - dragv(l+1) - liftv(l+1)
-    zlon(iim+1) = -plon(l+1)*xpi/180.
+    zs(nbp_lon+1, j) = phis(l+1)/rg
+    ps(nbp_lon+1, j) = p(l+1, 1)
+    ssou(nbp_lon+1, j) = dragu(l+1) + liftu(l+1)
+    ssov(nbp_lon+1, j) = dragv(l+1) + liftv(l+1)
+    blsu(nbp_lon+1, j) = phyu(l+1) - dragu(l+1) - liftu(l+1)
+    blsv(nbp_lon+1, j) = phyv(l+1) - dragv(l+1) - liftv(l+1)
+    zlon(nbp_lon+1) = -plon(l+1)*xpi/180.
     zlat(j) = plat(l+1)*xpi/180.
 
-    ub(iim+1, j) = 0.
-    vb(iim+1, j) = 0.
+    ub(nbp_lon+1, j) = 0.
+    vb(nbp_lon+1, j) = 0.
     DO k = 1, nlev
-      ub(iim+1, j) = ub(iim+1, j) + u(l+1, k)*(p(l+1,k)-p(l+1,k+1))/rg
-      vb(iim+1, j) = vb(iim+1, j) + v(l+1, k)*(p(l+1,k)-p(l+1,k+1))/rg
-    END DO
-
-
-    DO i = 1, iim
+      ub(nbp_lon+1, j) = ub(nbp_lon+1, j) + u(l+1, k)*(p(l+1,k)-p(l+1,k+1))/rg
+      vb(nbp_lon+1, j) = vb(nbp_lon+1, j) + v(l+1, k)*(p(l+1,k)-p(l+1,k+1))/rg
+    END DO
+
+
+    DO i = 1, nbp_lon
 
       l = l + 1
@@ -215 +213 @@
   ! South Pole
 
-  IF (jjm>1) THEN
+  IF (nbp_lat-1>1) THEN
     l = l + 1
-    ub(1, jjm+1) = 0.
-    vb(1, jjm+1) = 0.
+    ub(1, nbp_lat) = 0.
+    vb(1, nbp_lat) = 0.
     DO k = 1, nlev
-      ub(1, jjm+1) = ub(1, jjm+1) + u(l, k)*(p(l,k)-p(l,k+1))/rg
-      vb(1, jjm+1) = vb(1, jjm+1) + v(l, k)*(p(l,k)-p(l,k+1))/rg
-    END DO
-    zlat(jjm+1) = plat(l)*xpi/180.
-
-    DO i = 1, iim + 1
-      zs(i, jjm+1) = phis(l)/rg
-      ps(i, jjm+1) = p(l, 1)
-      ssou(i, jjm+1) = dragu(l) + liftu(l)
-      ssov(i, jjm+1) = dragv(l) + liftv(l)
-      blsu(i, jjm+1) = phyu(l) - dragu(l) - liftu(l)
-      blsv(i, jjm+1) = phyv(l) - dragv(l) - liftv(l)
-      ub(i, jjm+1) = ub(1, jjm+1)
-      vb(i, jjm+1) = vb(1, jjm+1)
+      ub(1, nbp_lat) = ub(1, nbp_lat) + u(l, k)*(p(l,k)-p(l,k+1))/rg
+      vb(1, nbp_lat) = vb(1, nbp_lat) + v(l, k)*(p(l,k)-p(l,k+1))/rg
+    END DO
+    zlat(nbp_lat) = plat(l)*xpi/180.
+
+    DO i = 1, nbp_lon + 1
+      zs(i, nbp_lat) = phis(l)/rg
+      ps(i, nbp_lat) = p(l, 1)
+      ssou(i, nbp_lat) = dragu(l) + liftu(l)
+      ssov(i, nbp_lat) = dragv(l) + liftv(l)
+      blsu(i, nbp_lat) = phyu(l) - dragu(l) - liftu(l)
+      blsv(i, nbp_lat) = phyv(l) - dragv(l) - liftv(l)
+      ub(i, nbp_lat) = ub(1, nbp_lat)
+      vb(i, nbp_lat) = vb(1, nbp_lat)
     END DO
   END IF
@@ -240 +238 @@
   ! MOMENT ANGULAIRE
 
-  DO j = 1, jjm
-    DO i = 1, iim
+  DO j = 1, nbp_lat-1
+    DO i = 1, nbp_lon
 
       raam(1) = raam(1) - rea**3*dlon*dlat*0.5*(cos(zlon(i))*sin(zlat(j))*cos &
@@ -274 +272 @@
 
 
-  DO j = 1, jjm
-    DO i = 1, iim
+  DO j = 1, nbp_lat-1
+    DO i = 1, nbp_lon
       tmou(1) = tmou(1) - rea**2*dlon*0.5*sin(zlon(i))*(zs(i,j)-zs(i,j+1))*( &
         cos(zlat(j+1))*ps(i,j+1)+cos(zlat(j))*ps(i,j))
@@ -283 +281 @@
   END DO
 
-  DO j = 2, jjm
-    DO i = 1, iim
+  DO j = 2, nbp_lat-1
+    DO i = 1, nbp_lon
       tmou(1) = tmou(1) + rea**2*dlat*0.5*sin(zlat(j))*(zs(i+1,j)-zs(i,j))*( &
         cos(zlon(i+1))*ps(i+1,j)+cos(zlon(i))*ps(i,j))
@@ -298 +296 @@
 
   l = 1
-  DO j = 2, jjm
-    DO i = 1, iim
+  DO j = 2, nbp_lat-1
+    DO i = 1, nbp_lon
       l = l + 1
       tsso(1) = tsso(1) - rea**3*cos(zlat(j))*dlon*dlat*ssou(i, j)*sin(zlat(j &
@@ -341 +339 @@
 100 FORMAT (F12.5, 15(1X,F12.5))
 
-  ! write(iam+1,*)((zs(i,j),i=1,iim),j=1,jjm+1)
-  ! write(iam+1,*)((ps(i,j),i=1,iim),j=1,jjm+1)
-  ! write(iam+1,*)((ub(i,j),i=1,iim),j=1,jjm+1)
-  ! write(iam+1,*)((vb(i,j),i=1,iim),j=1,jjm+1)
-  ! write(iam+1,*)((ssou(i,j),i=1,iim),j=1,jjm+1)
-  ! write(iam+1,*)((ssov(i,j),i=1,iim),j=1,jjm+1)
-  ! write(iam+1,*)((blsu(i,j),i=1,iim),j=1,jjm+1)
-  ! write(iam+1,*)((blsv(i,j),i=1,iim),j=1,jjm+1)
+  ! write(iam+1,*)((zs(i,j),i=1,nbp_lon),j=1,nbp_lat)
+  ! write(iam+1,*)((ps(i,j),i=1,nbp_lon),j=1,nbp_lat)
+  ! write(iam+1,*)((ub(i,j),i=1,nbp_lon),j=1,nbp_lat)
+  ! write(iam+1,*)((vb(i,j),i=1,nbp_lon),j=1,nbp_lat)
+  ! write(iam+1,*)((ssou(i,j),i=1,nbp_lon),j=1,nbp_lat)
+  ! write(iam+1,*)((ssov(i,j),i=1,nbp_lon),j=1,nbp_lat)
+  ! write(iam+1,*)((blsu(i,j),i=1,nbp_lon),j=1,nbp_lat)
+  ! write(iam+1,*)((blsv(i,j),i=1,nbp_lon),j=1,nbp_lat)
 
   aam = raam(3)