Changeset 1293 for LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/limit_netcdf.F

Timestamp:

Dec 22, 2009, 12:07:26 PM (14 years ago)

Author:

lguez

Message:

1) Replaced calls to "float" by calls to "real" and "dble" in
"phys_cosp". "call construct_cosp_gridbox(float(itap),..." was a bug
since the corresponding dummy argument has the type "double
precision". (And "float" is not in the Fortran standard.)

2) Modifications for the program "create_etat0_limit" only, "gcm" is
not touched:

2.1) Removed generic interface "startget" in module "startvar". Replaced
calls to "startget" by calls to "startget_phys2d", "startget_dyn" and
"startget_phys1d".

2.2) Simplified the interface of "startget_dyn" and made it more secure by
removing arguments for sizes of arrays and using assumed-shape array
arguments.

2.3) Corrected bug in call to "startget_dyn" for northward velocity. See
ticket #25 in Trac.

2.4) Collected procedures "inter_barxy", "inter_barx", "inter_bary",
"ord_coord" and "ord_coordm" in a module. "inter_barxy" is the only
public procedure in the module. Rewrote those five procedures: removed
arguments for sizes of arrays; used assumed-shape array arguments;
used array expressions; translated some spaghetti code in "inter_bary"
into structured code; corrected bug in "inter_bary", described by
ticket #26 in Trac.

2.5) Corrected a bug in "inter_bary": "y0" was initialized at 0
instead of -90. This bug made values at the south pole wrong. (This bug
has been acknowledeged by P. Le Van.)

2.6) Made the variable "champ" in procedure "limit_netcdf" an array of
rank 2. Thus it can be an argument in the call to the newly-secured
"inter_barxy".

2.7) The files "start.nc", "startphy.nc" and "limit.nc" are impacted
by this revision. There is a significant change in the variable "vcov"
of "start.nc". See
​http://web.lmd.jussieu.fr/~lglmd/Rev_1293/index.html. Note also the
changed value of "vcov" at the south pole. There is very little change
in other variables, in all three files.

File:

• LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/limit_netcdf.F (modified) (11 diffs)

LMDZ4/branches/LMDZ4V5.0-dev/libf/dyn3dpar/limit_netcdf.F

@@ -9 +9 @@
       USE dimphy
       use phys_state_var_mod , ONLY : pctsrf
+      use inter_barxy_m, only: inter_barxy
+
       IMPLICIT none
 c
@@ -70 +72 @@
       REAL , ALLOCATABLE :: dlon(:), dlat(:)
       REAL , ALLOCATABLE :: dlon_ini(:), dlat_ini(:)
-      REAL , ALLOCATABLE :: champ_msk(:), champ(:)
+      REAL , ALLOCATABLE :: champ_msk(:), champ(:, :)
       REAL , ALLOCATABLE :: work(:,:)
 
@@ -335 +337 @@
       ENDIF
 c
-      ALLOCATE( champ(imdep*jmdep) )
+      ALLOCATE( champ(imdep, jmdep) )
 
       DO  200 l = 1, lmdep
@@ -364 +366 @@
 
        IF ( interbar )   THEN
-         DO j = 1, imdep * jmdep
-           champ(j) = LOG(champ(j))
-         ENDDO
-
         IF( l.EQ.1 )  THEN
          WRITE(6,*) '-------------------------------------------------',
@@ -376 +374 @@
      ,'------------------------'
         ENDIF
-        CALL inter_barxy ( imdep,jmdep -1,dlon,dlat,champ ,
-     ,                  iim,jjm,rlonu,rlatv, jjp1,champint )
+        CALL inter_barxy(dlon, dlat(:jmdep -1), log(champ), rlonu(:iim),
+     $       rlatv, champint)
          DO j=1,jjp1
           DO i=1,iim
@@ -556 +554 @@
       ENDIF 
 c
-      ALLOCATE ( champ(imdep*jmdep) )
+      ALLOCATE ( champ(imdep, jmdep) )
 
       DO l = 1, lmdep
@@ -604 +602 @@
 cIM   ENDDO
 cIM22/02/2002
-         CALL inter_barxy ( imdep,jmdep -1,dlon, dlat ,
-     ,     champ, iim, jjm, rlonu, rlatv, jjp1, champint )
+       CALL inter_barxy (dlon, dlat(:jmdep -1), champ, rlonu(:iim),
+     $      rlatv, champint)
       ELSE
          CALL sea_ice(imdep, jmdep, dlon, dlat, champ,
@@ -881 +879 @@
       ENDIF 
 
-       ALLOCATE( champ(imdep*jmdep) )
+       ALLOCATE( champ(imdep, jmdep) )
        IF( extrap )   THEN
          ALLOCATE ( work(imdep,jmdep) )
@@ -924 +922 @@
      ,'------------------------'
         ENDIF
-       CALL inter_barxy ( imdep,jmdep -1,dlon, dlat ,
-     , champ, iim, jjm, rlonu, rlatv, jjp1, champint )
+        CALL inter_barxy (dlon, dlat(:jmdep -1), champ, rlonu(:iim),
+     $       rlatv, champint)
       ELSE
        CALL grille_m(imdep, jmdep, dlon, dlat, champ,
@@ -1082 +1080 @@
       ENDIF 
 c
-      ALLOCATE ( champ(imdep*jmdep) )
+      ALLOCATE ( champ(imdep, jmdep) )
 
       DO l = 1, lmdep
@@ -1120 +1118 @@
         ENDIF
 
-       CALL inter_barxy ( imdep,jmdep -1,dlon, dlat ,
-     , champ, iim, jjm, rlonu, rlatv, jjp1, champint )
+       CALL inter_barxy(dlon, dlat(:jmdep-1), champ, rlonu(:iim),
+     $       rlatv, champint)
       ELSE
        CALL grille_m(imdep, jmdep, dlon, dlat, champ,