Changeset 2346 for LMDZ5/trunk/libf/phylmd/oasis.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/oasis.F90 (modified) (13 diffs)

LMDZ5/trunk/libf/phylmd/oasis.F90

@@ -100 +100 @@
 #endif
     USE print_control_mod, ONLY: lunout
-
-    INCLUDE "dimensions.h"
+    USE mod_grid_phy_lmdz, ONLY nbp_lon, nbp_lat
 
 ! Local variables
@@ -228 +227 @@
 !************************************************************************************
     ig_paral(1) = 1                            ! apple partition for //
-    ig_paral(2) = (jj_begin-1)*iim+ii_begin-1  ! offset
-    ig_paral(3) = (jj_end*iim+ii_end) - (jj_begin*iim+ii_begin) + 1
-
-    IF (mpi_rank==mpi_size-1) ig_paral(3)=ig_paral(3)+iim-1
+    ig_paral(2) = (jj_begin-1)*nbp_lon+ii_begin-1  ! offset
+    ig_paral(3) = (jj_end*nbp_lon+ii_end) - (jj_begin*nbp_lon+ii_begin) + 1
+
+    IF (mpi_rank==mpi_size-1) ig_paral(3)=ig_paral(3)+nbp_lon-1
     WRITE(lunout,*) mpi_rank,'ig_paral--->',ig_paral(2),ig_paral(3)
     
@@ -248 +247 @@
 
     il_var_actual_shape(1) = 1
-    il_var_actual_shape(2) = iim
+    il_var_actual_shape(2) = nbp_lon
     il_var_actual_shape(3) = 1
-    il_var_actual_shape(4) = jjm+1
+    il_var_actual_shape(4) = nbp_lat
    
     il_var_type = PRISM_Real
@@ -320 +319 @@
 !
     USE print_control_mod, ONLY: lunout
-    INCLUDE "dimensions.h"
+    USE mod_grid_phy_lmdz, ONLY nbp_lon, nbp_lat
 ! Input arguments
 !************************************************************************************
@@ -327 +326 @@
 ! Output arguments
 !************************************************************************************
-    REAL, DIMENSION(iim, jj_nb,maxrecv), INTENT(OUT) :: tab_get
+    REAL, DIMENSION(nbp_lon, jj_nb,maxrecv), INTENT(OUT) :: tab_get
 
 ! Local variables
@@ -335 +334 @@
     CHARACTER (len = 20)          :: modname = 'fromcpl'
     CHARACTER (len = 80)          :: abort_message 
-    REAL, DIMENSION(iim*jj_nb)    :: field
+    REAL, DIMENSION(nbp_lon*jj_nb)    :: field
 
 !************************************************************************************
@@ -344 +343 @@
     istart=ii_begin
     IF (is_south_pole) THEN
-       iend=(jj_end-jj_begin)*iim+iim
+       iend=(jj_end-jj_begin)*nbp_lon+nbp_lon
     ELSE
-       iend=(jj_end-jj_begin)*iim+ii_end
+       iend=(jj_end-jj_begin)*nbp_lon+ii_end
     ENDIF
     
@@ -353 +352 @@
           field(:) = -99999.
           CALL prism_get_proto(inforecv(i)%nid, ktime, field(istart:iend), ierror)
-          tab_get(:,:,i) = RESHAPE(field(:),(/iim,jj_nb/))
+          tab_get(:,:,i) = RESHAPE(field(:),(/nbp_lon,jj_nb/))
        
           IF (ierror .NE. PRISM_Ok .AND. ierror.NE.PRISM_Recvd .AND. &
@@ -382 +381 @@
 ! 
     USE print_control_mod, ONLY: lunout
-    INCLUDE "dimensions.h"
+    USE mod_grid_phy_lmdz, ONLY nbp_lon, nbp_lat
 ! Input arguments
 !************************************************************************************
     INTEGER, INTENT(IN)                              :: ktime
     LOGICAL, INTENT(IN)                              :: last
-    REAL, DIMENSION(iim, jj_nb, maxsend), INTENT(IN) :: tab_put
+    REAL, DIMENSION(nbp_lon, jj_nb, maxsend), INTENT(IN) :: tab_put
 
 ! Local variables
@@ -395 +394 @@
     INTEGER                          :: wstart,wend
     INTEGER                          :: ierror, i
-    REAL, DIMENSION(iim*jj_nb)       :: field
+    REAL, DIMENSION(nbp_lon*jj_nb)       :: field
     CHARACTER (len = 20),PARAMETER   :: modname = 'intocpl'
     CHARACTER (len = 80)             :: abort_message 
@@ -410 +409 @@
     istart=ii_begin
     IF (is_south_pole) THEN
-       iend=(jj_end-jj_begin)*iim+iim
+       iend=(jj_end-jj_begin)*nbp_lon+nbp_lon
     ELSE
-       iend=(jj_end-jj_begin)*iim+ii_end
+       iend=(jj_end-jj_begin)*nbp_lon+ii_end
     ENDIF
     
@@ -418 +417 @@
        wstart=istart
        wend=iend
-       IF (is_north_pole) wstart=istart+iim-1
-       IF (is_south_pole) wend=iend-iim+1
+       IF (is_north_pole) wstart=istart+nbp_lon-1
+       IF (is_south_pole) wend=iend-nbp_lon+1
        
        DO i = 1, maxsend
           IF (infosend(i)%action) THEN
-             field = RESHAPE(tab_put(:,:,i),(/iim*jj_nb/))
+             field = RESHAPE(tab_put(:,:,i),(/nbp_lon*jj_nb/))
              CALL writefield_phy(infosend(i)%name,field(wstart:wend),1)
           END IF
@@ -435 +434 @@
     DO i = 1, maxsend
       IF (infosend(i)%action .AND. infosend(i)%nid .NE. -1 ) THEN
-          field = RESHAPE(tab_put(:,:,i),(/iim*jj_nb/))
+          field = RESHAPE(tab_put(:,:,i),(/nbp_lon*jj_nb/))
           CALL prism_put_proto(infosend(i)%nid, ktime, field(istart:iend), ierror)