Changeset 7025 for branches/2016/dev_r6999_CONFIGMAN_1/NEMOGCM/TOOLS/SIREN/src/create_coord.f90

Timestamp:

2016-10-13T11:08:38+02:00 (8 years ago)

Author:

jpaul

Message:

see ticket #1781

File:

• branches/2016/dev_r6999_CONFIGMAN_1/NEMOGCM/TOOLS/SIREN/src/create_coord.f90 (modified) (8 diffs)

branches/2016/dev_r6999_CONFIGMAN_1/NEMOGCM/TOOLS/SIREN/src/create_coord.f90

@@ -76 +76 @@
 !>
 !>    * _nesting namelist (namnst)_:<br/>
+!>       you could define sub domain with coarse grid indices or with coordinates.
 !>       - in_imin0 : i-direction lower left  point indice 
 !> of coarse grid subdomain to be used
@@ -84 +85 @@
 !>       - in_jmax0 : j-direction upper right point indice
 !> of coarse grid subdomain to be used
+!>       - rn_lonmin0 : lower left  longitude of coarse grid subdomain to be used
+!>       - rn_lonmax0 : upper right longitude of coarse grid subdomain to be used
+!>       - rn_latmin0 : lower left  latitude  of coarse grid subdomain to be used
+!>       - rn_latmax0 : upper right latitude  of coarse grid subdomain to be used
 !>       - in_rhoi  : refinement factor in i-direction
 !>       - in_rhoj  : refinement factor in j-direction<br/>
 !>
 !>       \image html  grid_zoom_40.png 
-!>       \image latex grid_zoom_40.png
+!>       <center> \image latex grid_zoom_40.png 
+!>       </center>
 !>
 !>    * _output namelist (namout)_:
@@ -102 +108 @@
 !> @date September, 2015
 !> - manage useless (dummy) variable, attributes, and dimension
+!> @date September, 2016
+!> - allow to use coordinate to define subdomain
 !>
 !> @note Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt)
@@ -139 +147 @@
    INTEGER(i4)                                          :: il_nvar
    INTEGER(i4)                                          :: il_ew
+   INTEGER(i4)                                          :: il_imin0
+   INTEGER(i4)                                          :: il_imax0
+   INTEGER(i4)                                          :: il_jmin0
+   INTEGER(i4)                                          :: il_jmax0
+
    INTEGER(i4)      , DIMENSION(ip_maxdim)              :: il_rho
+   INTEGER(i4)      , DIMENSION(2)                      :: il_index
    INTEGER(i4)      , DIMENSION(2,2,ip_npoint)          :: il_offset
 
@@ -177 +191 @@
 
    !namnst
+   REAL(sp)          :: rn_lonmin0 = -360.
+   REAL(sp)          :: rn_lonmax0 = -360.
+   REAL(sp)          :: rn_latmin0 = -360.
+   REAL(sp)          :: rn_latmax0 = -360.
    INTEGER(i4)       :: in_imin0 = 0
    INTEGER(i4)       :: in_imax0 = 0
@@ -207 +225 @@
    
    NAMELIST /namnst/ &  !  nesting namelist
+   &  rn_lonmin0, &     !< lower left  coarse grid longitude
+   &  rn_lonmax0, &     !< upper right coarse grid longitude
+   &  rn_latmin0, &     !< lower left  coarse grid latitude
+   &  rn_latmax0, &     !< upper right coarse grid latitude
    &  in_imin0,   &     !< i-direction lower left  point indice 
    &  in_imax0,   &     !< i-direction upper right point indice
@@ -301 +323 @@
 
    ! check nesting parameters
-   IF( in_imin0 < 0 .OR. in_imax0 < 0 .OR. in_jmin0 < 0 .OR. in_jmax0 < 0)THEN
+   il_index(:)=0
+   IF( rn_lonmin0 >= -180. .AND. rn_lonmin0 <= 360 .AND. &
+     & rn_latmin0 >= -90.  .AND. rn_latmin0 <= 90. )THEN
+
+      il_index(:)=grid_get_closest(tl_coord0, &
+         &                         REAL(rn_lonmin0,dp), REAL(rn_latmin0,dp), &
+         &                         cd_pos='ll') 
+      il_imin0=il_index(1)
+      il_jmin0=il_index(2)
+   ELSE
+      il_imin0=in_imin0
+      il_jmin0=in_jmin0
+   ENDIF
+
+   il_index(:)=0
+   IF( rn_lonmax0 >= -180. .AND. rn_lonmax0 <= 360 .AND. &
+     & rn_latmax0 >= -90.  .AND. rn_latmax0 <= 90. )THEN
+      il_index(:)=grid_get_closest(tl_coord0, &
+         &                         REAL(rn_lonmax0,dp), REAL(rn_latmax0,dp), &
+         &                         cd_pos='ur') 
+      il_imax0=il_index(1)
+      il_jmax0=il_index(2)
+   ELSE
+      il_imax0=in_imax0
+      il_jmax0=in_jmax0
+   ENDIF
+
+   ! forced indices for east west cyclic domain
+   IF( rn_lonmin0 == rn_lonmax0 .AND. &
+     & rn_lonmin0 /= -360. )THEN
+      il_imin0=0
+      il_imax0=0
+   ENDIF
+
+   IF( il_imin0 < 0 .OR. il_imax0 < 0 .OR. il_jmin0 < 0 .OR. il_jmax0 < 0)THEN
       CALL logger_fatal("CREATE COORD: invalid points indices."//&
       &  " check namelist "//TRIM(cl_namelist))
@@ -318 +374 @@
 
    ! check domain validity
-   CALL grid_check_dom(tl_coord0, in_imin0, in_imax0, in_jmin0, in_jmax0 )
+   CALL grid_check_dom(tl_coord0, il_imin0, il_imax0, il_jmin0, il_jmax0 )
 
    ! compute domain
    tl_dom=dom_init( tl_coord0,         &
-   &                in_imin0, in_imax0,&
-   &                in_jmin0, in_jmax0 )
+   &                il_imin0, il_imax0,&
+   &                il_jmin0, il_jmax0 )
 
    ! add extra band (if need be) to compute interpolation