Changeset 1738 for trunk/NEMO

Timestamp:

2009-11-18T15:37:06+01:00 (14 years ago)

Author:

smasson

Message:

iom_put for scalar values, see ticket:604

File:

• trunk/NEMO/OPA_SRC/IOM/iom.F90 (modified) (6 diffs)

trunk/NEMO/OPA_SRC/IOM/iom.F90

@@ -59 +59 @@
    END INTERFACE
   INTERFACE iom_put
-     MODULE PROCEDURE iom_p2d, iom_p3d
+     MODULE PROCEDURE iom_p0d, iom_p2d, iom_p3d
   END INTERFACE
 #if defined key_iomput
@@ -99 +99 @@
 
       ! horizontal grid definition
+      CALL set_scalar
       CALL set_grid( "grid_T", glamt, gphit )
       CALL set_grid( "grid_U", glamu, gphiu )
@@ -910 +911 @@
    !!                   INTERFACE iom_put
    !!----------------------------------------------------------------------
+   SUBROUTINE iom_p0d( cdname, pfield0d )
+      CHARACTER(LEN=*), INTENT(in) ::   cdname
+      REAL(wp)        , INTENT(in) ::   pfield0d
+#if defined key_iomput
+      CALL event__write_field2D( cdname, RESHAPE( (/pfield0d/), (/1,1/) ) )
+#else
+      IF( .FALSE. )   WRITE(numout,*) cdname, pfield0d   ! useless test to avoid compilation warnings
+#endif
+   END SUBROUTINE iom_p0d
+
    SUBROUTINE iom_p2d( cdname, pfield2d )
       CHARACTER(LEN=*)            , INTENT(in) ::   cdname
@@ -953 +964 @@
 
 
+   SUBROUTINE set_scalar
+      !!----------------------------------------------------------------------
+      !!                     ***  ROUTINE   ***
+      !!
+      !! ** Purpose :   define fake grids for scalar point
+      !!
+      !!----------------------------------------------------------------------
+      REAL(wp), DIMENSION(1,1) ::   zz = 1.
+      !!----------------------------------------------------------------------
+      CALL event__set_grid_dimension( 'scalarpoint', jpnij, 1)
+      CALL event__set_grid_domain   ( 'scalarpoint', 1, 1, narea, 1, zz, zz )
+      CALL event__set_grid_type_nemo( 'scalarpoint' )
+
+   END SUBROUTINE set_scalar
+
+
    SUBROUTINE set_xmlatt
       !!----------------------------------------------------------------------
@@ -960 +987 @@
       !!
       !!----------------------------------------------------------------------
-      CHARACTER(len=6),DIMENSION( 5) ::   clsuff                   ! suffix name
+      CHARACTER(len=6),DIMENSION( 6) ::   clsuff                   ! suffix name
       CHARACTER(len=1),DIMENSION( 3) ::   clgrd                    ! suffix name
       CHARACTER(len=50)              ::   clname                   ! file name
@@ -988 +1015 @@
       
       ! output file names (attribut: name)
-      clsuff(:) = (/ 'grid_T', 'grid_U', 'grid_V', 'grid_W', 'icemod' /)      
+      clsuff(:) = (/ 'grid_T', 'grid_U', 'grid_V', 'grid_W', 'icemod', 'scalar' /)      
       DO jg = 1, SIZE(clsuff)                                                                  ! grid type
          DO jh = 1, 12                                                                         ! 1, 2, 3, 4, 6, 12 hours