Changeset 2531 for CONFIG/UNIFORM/v6_rc0/IPSLCM6/SOURCES

Timestamp:

05/29/15 11:16:25 (8 years ago)

Author:

aclsce

Message:

Modifications to have correct NEMO diaptr files with mask.

File:

• CONFIG/UNIFORM/v6_rc0/IPSLCM6/SOURCES/NEMO/iom.F90 (modified) (4 diffs)

CONFIG/UNIFORM/v6_rc0/IPSLCM6/SOURCES/NEMO/iom.F90

@@ -127 +127 @@
          CALL set_grid( "V", glamv, gphiv )
          CALL set_grid( "W", glamt, gphit )
+         CALL set_grid_znl( gphit )
       ENDIF
 
@@ -136 +137 @@
          CALL set_grid( "V", glamv_crs, gphiv_crs ) 
          CALL set_grid( "W", glamt_crs, gphit_crs ) 
+         CALL set_grid_znl( gphit_crs )
           !
          CALL dom_grid_glo   ! Return to parent grid domain
@@ -1212 +1214 @@
 
 
+   SUBROUTINE set_grid_znl( plat )
+      !!----------------------------------------------------------------------
+      !!                     ***  ROUTINE set_grid_znl  ***
+      !!
+      !! ** Purpose :   define 2D grids for zonal mean
+      !!
+      !!----------------------------------------------------------------------
+      REAL(wp), DIMENSION(jpi,jpj), INTENT(in) ::   plat
+      !
+      REAL(wp), DIMENSION(:), ALLOCATABLE  ::   zlon
+      INTEGER  :: ni,nj, ix, iy
+
+      
+      ni=nlei-nldi+1 ; nj=nlej-nldj+1
+      ALLOCATE( zlon(ni*nj) )       ;       zlon(:) = 0.
+
+      CALL iom_set_domain_attr("gznl", ni_glo=jpiglo, nj_glo=jpjglo, ibegin=nimpp+nldi-1, jbegin=njmpp+nldj-1, ni=ni, nj=nj)
+      CALL iom_set_domain_attr("gznl", data_dim=2, data_ibegin = 1-nldi, data_ni = jpi, data_jbegin = 1-nldj, data_nj = jpj)
+      CALL iom_set_domain_attr("gznl", lonvalue = zlon,   &
+         &                             latvalue = RESHAPE(plat(nldi:nlei, nldj:nlej),(/ ni*nj /)))  
+      CALL dom_ngb( 180., 90., ix, iy, 'T' ) !  i-line that passes near the North Pole : Reference latitude (used in plots)
+      CALL iom_set_domain_attr ('ptr', zoom_ibegin=ix, zoom_nj=jpjglo)
+      CALL iom_update_file_name('ptr')
+      !
+   END SUBROUTINE set_grid_znl
+
+
    SUBROUTINE set_scalar
       !!----------------------------------------------------------------------
@@ -1291 +1320 @@
       zlatpira = (/ -19.0, -14.0,  -8.0, 0.0, 4.0, 8.0, 12.0, 15.0, 20.0 /)
       CALL set_mooring( zlonpira, zlatpira )
-    
-      ! diaptr : zonal mean 
-      CALL iom_set_domain_attr ('ptr', zoom_ibegin=nxline, zoom_nj=jpjglo)
-      CALL iom_update_file_name('ptr')
       !
    END SUBROUTINE set_xmlatt