Changeset 1530 for trunk/LMDZ.COMMON/libf/dynlonlat_phylonlat

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

Location:

trunk/LMDZ.COMMON/libf/dynlonlat_phylonlat

Files:

• phytitan/readstartphy.F (moved) (moved from trunk/LMDZ.TITAN/libf/phytitan/readstartphy.F)
• phytitan/writerestartphy.F (moved) (moved from trunk/LMDZ.TITAN/libf/phytitan/writerestartphy.F)
• phyvenus/readstartphy.F (moved) (moved from trunk/LMDZ.VENUS/libf/phyvenus/readstartphy.F) (17 diffs)
• phyvenus/writerestartphy.F (moved) (moved from trunk/LMDZ.VENUS/libf/phyvenus/writerestartphy.F)

trunk/LMDZ.COMMON/libf/dynlonlat_phylonlat/phyvenus/readstartphy.F

@@ -15 +15 @@
 c======================================================================
       use dimphy
+      USE mod_grid_phy_lmdz, only: klon_glo
       IMPLICIT none
-#include "dimensions.h"
 #include "netcdf.inc"
 #include "dimsoil.h"
 c======================================================================
-      integer ngridmx
-      parameter (ngridmx=(2+(jjm-1)*iim - 1/jjm))
 
       CHARACTER*(*) fichnom
-      REAL rlat(ngridmx), rlon(ngridmx)
-      REAL tsol(ngridmx)
-      REAL tsoil(ngridmx,nsoilmx)
-      REAL albe(ngridmx)
-      REAL radsol(ngridmx)
-      REAL sollw(ngridmx)
-      real solsw(ngridmx)
-      real fder(ngridmx)
-      real dlw(ngridmx)
-      real sollwdown(ngridmx)
-      REAL zmea(ngridmx), zstd(ngridmx)
-      REAL zsig(ngridmx), zgam(ngridmx), zthe(ngridmx)
-      REAL zpic(ngridmx), zval(ngridmx)
+      REAL rlat(klon_glo), rlon(klon_glo)
+      REAL tsol(klon_glo)
+      REAL tsoil(klon_glo,nsoilmx)
+      REAL albe(klon_glo)
+      REAL radsol(klon_glo)
+      REAL sollw(klon_glo)
+      real solsw(klon_glo)
+      real fder(klon_glo)
+      real dlw(klon_glo)
+      real sollwdown(klon_glo)
+      REAL zmea(klon_glo), zstd(klon_glo)
+      REAL zsig(klon_glo), zgam(klon_glo), zthe(klon_glo)
+      REAL zpic(klon_glo), zval(klon_glo)
       INTEGER length
       PARAMETER (length=100)
@@ -129 +127 @@
          xmin = 1.0E+20
          xmax = -1.0E+20
-         DO i = 1, ngridmx
+         DO i = 1, klon_glo
             xmin = MIN(tsol(i),xmin)
             xmax = MAX(tsol(i),xmax)
@@ -148 +146 @@
          PRINT*, "phyetat0: Le champ <Tsoil"//str2//"> est absent"
          PRINT*, "          Il prend donc la valeur de surface"
-         DO i=1, ngridmx
+         DO i=1, klon_glo
              tsoil(i,isoil)=tsol(i)
          ENDDO
@@ -185 +183 @@
          xmin = 1.0E+20
          xmax = -1.0E+20
-         DO i = 1, ngridmx
+         DO i = 1, klon_glo
             xmin = MIN(albe(i),xmin)
             xmax = MAX(albe(i),xmax)
@@ -213 +211 @@
       xmin = 1.0E+20
       xmax = -1.0E+20
-      DO i = 1, ngridmx
+      DO i = 1, klon_glo
          xmin = MIN(solsw(i),xmin)
          xmax = MAX(solsw(i),xmax)
@@ -239 +237 @@
       xmin = 1.0E+20
       xmax = -1.0E+20
-      DO i = 1, ngridmx
+      DO i = 1, klon_glo
          xmin = MIN(sollw(i),xmin)
          xmax = MAX(sollw(i),xmax)
@@ -266 +264 @@
       xmin = 1.0E+20
       xmax = -1.0E+20
-      DO i = 1, ngridmx
+      DO i = 1, klon_glo
          xmin = MIN(fder(i),xmin)
          xmax = MAX(fder(i),xmax)
@@ -293 +291 @@
       xmin = 1.0E+20
       xmax = -1.0E+20
-      DO i = 1, ngridmx
+      DO i = 1, klon_glo
          xmin = MIN(dlw(i),xmin)
          xmax = MAX(dlw(i),xmax)
@@ -320 +318 @@
       xmin = 1.0E+20
       xmax = -1.0E+20
-      DO i = 1, ngridmx
+      DO i = 1, klon_glo
          xmin = MIN(sollwdown(i),xmin)
          xmax = MAX(sollwdown(i),xmax)
@@ -345 +343 @@
       xmin = 1.0E+20
       xmax = -1.0E+20
-      DO i = 1, ngridmx
+      DO i = 1, klon_glo
          xmin = MIN(radsol(i),xmin)
          xmax = MAX(radsol(i),xmax)
@@ -372 +370 @@
       xmin = 1.0E+20
       xmax = -1.0E+20
-      DO i = 1, ngridmx
+      DO i = 1, klon_glo
          xmin = MIN(zmea(i),xmin)
          xmax = MAX(zmea(i),xmax)
@@ -396 +394 @@
       xmin = 1.0E+20
       xmax = -1.0E+20
-      DO i = 1, ngridmx
+      DO i = 1, klon_glo
          xmin = MIN(zstd(i),xmin)
          xmax = MAX(zstd(i),xmax)
@@ -420 +418 @@
       xmin = 1.0E+20
       xmax = -1.0E+20
-      DO i = 1, ngridmx
+      DO i = 1, klon_glo
          xmin = MIN(zsig(i),xmin)
          xmax = MAX(zsig(i),xmax)
@@ -444 +442 @@
       xmin = 1.0E+20
       xmax = -1.0E+20
-      DO i = 1, ngridmx
+      DO i = 1, klon_glo
          xmin = MIN(zgam(i),xmin)
          xmax = MAX(zgam(i),xmax)
@@ -468 +466 @@
       xmin = 1.0E+20
       xmax = -1.0E+20
-      DO i = 1, ngridmx
+      DO i = 1, klon_glo
          xmin = MIN(zthe(i),xmin)
          xmax = MAX(zthe(i),xmax)
@@ -492 +490 @@
       xmin = 1.0E+20
       xmax = -1.0E+20
-      DO i = 1, ngridmx
+      DO i = 1, klon_glo
          xmin = MIN(zpic(i),xmin)
          xmax = MAX(zpic(i),xmax)
@@ -516 +514 @@
       xmin = 1.0E+20
       xmax = -1.0E+20
-      DO i = 1, ngridmx
+      DO i = 1, klon_glo
          xmin = MIN(zval(i),xmin)
          xmax = MAX(zval(i),xmax)