New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
domain.F90 in utils/tools/DOMAINcfg/src – NEMO

source: utils/tools/DOMAINcfg/src/domain.F90 @ 14449

Last change on this file since 14449 was 14199, checked in by mathiot, 4 years ago

ticket #2588: various bug in domain_cfg (undefined variable before use, isf, closed sea re-activated ...)
File size: 23.5 KB
Line 
1MODULE domain
2   !!==============================================================================
3   !!                       ***  MODULE domain   ***
4   !! Ocean initialization : domain initialization
5   !!==============================================================================
6   !! History :  OPA  !  1990-10  (C. Levy - G. Madec)  Original code
7   !!                 !  1992-01  (M. Imbard) insert time step initialization
8   !!                 !  1996-06  (G. Madec) generalized vertical coordinate
9   !!                 !  1997-02  (G. Madec) creation of domwri.F
10   !!                 !  2001-05  (E.Durand - G. Madec) insert closed sea
11   !!   NEMO     1.0  !  2002-08  (G. Madec)  F90: Free form and module
12   !!            2.0  !  2005-11  (V. Garnier) Surface pressure gradient organization
13   !!            3.3  !  2010-11  (G. Madec)  initialisation in C1D configuration
14   !!            3.6  !  2013     ( J. Simeon, C. Calone, G. Madec, C. Ethe ) Online coarsening of outputs
15   !!            3.7  !  2015-11  (G. Madec, A. Coward)  time varying zgr by default
16   !!----------------------------------------------------------------------
17   
18   !!----------------------------------------------------------------------
19   !!   dom_init       : initialize the space and time domain
20   !!   dom_nam        : read and contral domain namelists
21   !!   dom_ctl        : control print for the ocean domain
22   !!----------------------------------------------------------------------
23   USE dom_oce         ! domain: ocean
24   USE phycst          ! physical constants
25   USE domhgr          ! domain: set the horizontal mesh
26   USE domzgr          ! domain: set the vertical mesh
27   USE dommsk          ! domain: set the mask system
28   USE domclo          ! domain: set closed sea mask
29   !
30   USE lib_mpp         !
31   USE in_out_manager  ! I/O manager
32   USE iom             !
33
34   IMPLICIT NONE
35   PRIVATE
36
37   PUBLIC   dom_init   ! called by opa.F90
38   PUBLIC   dom_nam  ! called by opa.F90
39   PUBLIC   cfg_write   ! called by opa.F90
40
41   !!-------------------------------------------------------------------------
42   !! NEMO/OCE 4.0 , NEMO Consortium (2018)
43   !! $Id: domain.F90 6140 2015-12-21 11:35:23Z timgraham $
44   !! Software governed by the CeCILL licence        (./LICENSE)
45   !!-------------------------------------------------------------------------
46CONTAINS
47
48   SUBROUTINE dom_init
49      !!----------------------------------------------------------------------
50      !!                  ***  ROUTINE dom_init  ***
51      !!                   
52      !! ** Purpose :   Domain initialization. Call the routines that are
53      !!              required to create the arrays which define the space
54      !!              and time domain of the ocean model.
55      !!
56      !! ** Method  : - dom_msk: compute the masks from the bathymetry file
57      !!              - dom_hgr: compute or read the horizontal grid-point position
58      !!                         and scale factors, and the coriolis factor
59      !!              - dom_zgr: define the vertical coordinate and the bathymetry
60      !!              - dom_stp: defined the model time step
61      !!              - dom_wri: create the meshmask file if nmsh=1
62      !!              - 1D configuration, move Coriolis, u and v at T-point
63      !!----------------------------------------------------------------------
64      INTEGER ::   jk          ! dummy loop indices
65      INTEGER ::   iconf = 0   ! local integers
66      REAL(wp), POINTER, DIMENSION(:,:) ::   z1_hu_0, z1_hv_0
67      INTEGER , DIMENSION(jpi,jpj) ::   ik_top , ik_bot       ! top and bottom ocean level
68      !!----------------------------------------------------------------------
69      !
70      IF(lwp) THEN
71         WRITE(numout,*)
72         WRITE(numout,*) 'dom_init : domain initialization'
73         WRITE(numout,*) '~~~~~~~~'
74      ENDIF
75      !
76      !                       !==  Reference coordinate system  ==!
77      !
78      CALL dom_nam               ! read namelist ( namrun, namdom )
79         
80      CALL dom_hgr               ! Horizontal mesh
81      CALL dom_zgr( ik_top, ik_bot )  ! Vertical mesh and bathymetry
82      CALL dom_msk( ik_top, ik_bot )  ! Masks
83      IF ( ln_domclo ) CALL dom_clo               ! Closed seas and lake
84      !
85      CALL dom_ctl                  ! print extrema of masked scale factors
86      !
87#if ! defined key_agrif
88      CALL cfg_write                ! create the configuration file
89#endif
90      !
91   END SUBROUTINE dom_init
92
93   SUBROUTINE dom_nam
94      !!----------------------------------------------------------------------
95      !!                     ***  ROUTINE dom_nam  ***
96      !!                   
97      !! ** Purpose :   read domaine namelists and print the variables.
98      !!
99      !! ** input   : - namrun namelist
100      !!              - namdom namelist
101      !!              - namnc4 namelist   ! "key_netcdf4" only
102      !!----------------------------------------------------------------------
103      USE ioipsl
104      NAMELIST/namrun/ cn_exp   ,    &         
105         &             nn_it000, nn_itend , nn_date0    , nn_time0     , nn_leapy  ,     &
106         &             ln_mskland  , ln_clobber   , nn_chunksz,     &
107         &             ln_cfmeta, ln_iscpl
108
109      NAMELIST/namdom/ ln_read_cfg, nn_bathy, cn_domcfg, cn_topo, cn_bath, cn_lon, cn_lat, rn_scale, nn_interp, &
110         &             cn_topolvl, cn_fisfd, cn_visfd, cn_bathlvl, cn_fcoord,                       & 
111         &             rn_bathy , rn_e3zps_min, rn_e3zps_rat, nn_msh, rn_hmin,                       &
112         &             rn_atfp , rn_rdt   ,  ln_crs      , jphgr_msh ,                               &
113         &             ppglam0, ppgphi0, ppe1_deg, ppe2_deg, ppe1_m, ppe2_m,                         &
114         &             ppsur, ppa0, ppa1, ppkth, ppacr, ppdzmin, pphmax, ldbletanh,                  &
115         &             ppa2, ppkth2, ppacr2
116
117      INTEGER  ::   ios                 ! Local integer output status for namelist read
118      !!----------------------------------------------------------------------
119
120      REWIND( numnam_ref )              ! Namelist namrun in reference namelist : Parameters of the run
121      READ  ( numnam_ref, namrun, IOSTAT = ios, ERR = 901)
122901   IF( ios /= 0 )   CALL ctl_nam ( ios , 'namrun in reference namelist', lwp )
123
124      REWIND( numnam_cfg )              ! Namelist namrun in configuration namelist : Parameters of the run
125      READ  ( numnam_cfg, namrun, IOSTAT = ios, ERR = 902 )
126902   IF( ios /= 0 )   CALL ctl_nam ( ios , 'namrun in configuration namelist', lwp )
127      IF(lwm) WRITE ( numond, namrun )
128      !
129      IF(lwp) THEN                  ! control print
130         WRITE(numout,*)
131         WRITE(numout,*) 'dom_nam  : domain initialization through namelist read'
132         WRITE(numout,*) '~~~~~~~ '
133         WRITE(numout,*) '   Namelist namrun'
134         WRITE(numout,*) '      experiment name for output      cn_exp     = ', cn_exp
135         WRITE(numout,*) '      mask land points                ln_mskland = ', ln_mskland
136         WRITE(numout,*) '      additional CF standard metadata ln_cfmeta  = ', ln_cfmeta
137         WRITE(numout,*) '      overwrite an existing file      ln_clobber = ', ln_clobber
138         WRITE(numout,*) '      NetCDF chunksize (bytes)        nn_chunksz = ', nn_chunksz
139      ENDIF
140
141      cexper = cn_exp
142      nit000 = nn_it000
143      nitend = nn_itend
144      ndate0 = nn_date0
145      nleapy = nn_leapy
146
147      !
148      cn_topo =''
149      cn_bath =''
150      cn_lon  =''
151      cn_lat  =''
152      rn_scale = 1.
153
154      REWIND( numnam_ref )              ! Namelist namdom in reference namelist : space & time domain (bathymetry, mesh, timestep)
155      READ  ( numnam_ref, namdom, IOSTAT = ios, ERR = 903)
156903   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namdom in reference namelist', lwp )
157 
158      !
159      REWIND( numnam_cfg )              ! Namelist namdom in configuration namelist : space & time domain (bathymetry, mesh, timestep)
160      READ  ( numnam_cfg, namdom, IOSTAT = ios, ERR = 904 )
161904   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namdom in configuration namelist', lwp )
162      IF(lwm) WRITE ( numond, namdom )
163      !
164
165
166
167      IF(lwp) THEN
168         WRITE(numout,*)
169         WRITE(numout,*) '   Namelist namdom : space & time domain'
170         WRITE(numout,*) '      flag read/compute bathymetry      nn_bathy     = ', nn_bathy
171         IF( nn_bathy == 1 ) THEN
172            WRITE(numout,*) '   read bathymetry from file      cn_topo      = ' ,TRIM(cn_topo)
173            WRITE(numout,*) '   bathymetry name in file        cn_bath      = ' ,TRIM(cn_bath)
174            WRITE(numout,*) '   read isf draft from file       cn_fisfd     = ' ,TRIM(cn_fisfd)
175            WRITE(numout,*) '   isf draft name in file         cn_visfd     = ' ,TRIM(cn_visfd)
176         ELSE IF( nn_bathy == 2 ) THEN
177            WRITE(numout,*) '   compute bathymetry from file      cn_topo      = ' , cn_topo
178            WRITE(numout,*) '   bathymetry name in file           cn_bath      = ' , cn_bath
179            WRITE(numout,*) '   longitude name in file            cn_lon       = ' , cn_lon
180            WRITE(numout,*) '   latitude  name in file            cn_lat       = ' , cn_lat
181            WRITE(numout,*) '   bathmetry scale factor            rn_scale     = ' , rn_scale 
182         ENDIF   
183         WRITE(numout,*) '      Depth (if =0 bathy=jpkm1)         rn_bathy     = ', rn_bathy
184         WRITE(numout,*) '      min depth of the ocean    (>0) or    rn_hmin   = ', rn_hmin
185         WRITE(numout,*) '      min number of ocean level (<0)       '
186         WRITE(numout,*) '      minimum thickness of partial      rn_e3zps_min = ', rn_e3zps_min, ' (m)'
187         WRITE(numout,*) '         step level                     rn_e3zps_rat = ', rn_e3zps_rat
188         WRITE(numout,*) '      create mesh/mask file(s)          nn_msh       = ', nn_msh
189         WRITE(numout,*) '           = 0   no file created           '
190         WRITE(numout,*) '           = 1   mesh_mask                 '
191         WRITE(numout,*) '           = 2   mesh and mask             '
192         WRITE(numout,*) '           = 3   mesh_hgr, msh_zgr and mask'
193         WRITE(numout,*) '      type of horizontal mesh jphgr_msh           = ', jphgr_msh
194         WRITE(numout,*) '      longitude of first raw and column T-point ppglam0 = ', ppglam0
195         WRITE(numout,*) '      latitude  of first raw and column T-point ppgphi0 = ', ppgphi0
196         WRITE(numout,*) '      zonal      grid-spacing (degrees) ppe1_deg        = ', ppe1_deg
197         WRITE(numout,*) '      meridional grid-spacing (degrees) ppe2_deg        = ', ppe2_deg
198         WRITE(numout,*) '      zonal      grid-spacing (degrees) ppe1_m          = ', ppe1_m
199         WRITE(numout,*) '      meridional grid-spacing (degrees) ppe2_m          = ', ppe2_m
200         WRITE(numout,*) '      ORCA r4, r2 and r05 coefficients  ppsur           = ', ppsur
201         WRITE(numout,*) '                                        ppa0            = ', ppa0
202         WRITE(numout,*) '                                        ppa1            = ', ppa1
203         WRITE(numout,*) '                                        ppkth           = ', ppkth
204         WRITE(numout,*) '                                        ppacr           = ', ppacr
205         WRITE(numout,*) '      Minimum vertical spacing ppdzmin                  = ', ppdzmin
206         WRITE(numout,*) '      Maximum depth pphmax                              = ', pphmax
207         WRITE(numout,*) '      Use double tanf function for vertical coordinates ldbletanh = ', ldbletanh
208         WRITE(numout,*) '      Double tanh function parameters ppa2              = ', ppa2
209         WRITE(numout,*) '                                      ppkth2            = ', ppkth2
210         WRITE(numout,*) '                                      ppacr2            = ', ppacr2
211      ENDIF
212      !
213      ntopo     = nn_bathy          ! conversion DOCTOR names into model names (this should disappear soon)
214      e3zps_min = rn_e3zps_min
215      e3zps_rat = rn_e3zps_rat
216      nmsh      = nn_msh
217      atfp      = rn_atfp
218      rdt       = rn_rdt
219
220      snc4set%luse = .FALSE.        ! No NetCDF 4 case
221      !
222   END SUBROUTINE dom_nam
223
224
225   SUBROUTINE dom_ctl
226      !!----------------------------------------------------------------------
227      !!                     ***  ROUTINE dom_ctl  ***
228      !!
229      !! ** Purpose :   Domain control.
230      !!
231      !! ** Method  :   compute and print extrema of masked scale factors
232      !!----------------------------------------------------------------------
233      INTEGER ::   iimi1, ijmi1, iimi2, ijmi2, iima1, ijma1, iima2, ijma2
234      INTEGER, DIMENSION(2) ::   iloc   !
235      REAL(wp) ::   ze1min, ze1max, ze2min, ze2max
236      !!----------------------------------------------------------------------
237      !
238#undef CHECK_DOM
239#ifdef CHECK_DOM
240      IF(lk_mpp) THEN
241         CALL mpp_minloc( 'dom_ctl', e1t(:,:), tmask_i(:,:), ze1min, iloc )
242         iimi1 = iloc(1) ; ijmi1 = iloc(2)
243         CALL mpp_minloc( 'dom_ctl', e2t(:,:), tmask_i(:,:), ze2min, iloc )
244         iimi2 = iloc(1) ; ijmi2 = iloc(2)
245         CALL mpp_maxloc( 'dom_ctl', e1t(:,:), tmask_i(:,:), ze1max, iloc )
246         iima1 = iloc(1) ; ijma1 = iloc(2)
247         CALL mpp_maxloc( 'dom_ctl', e2t(:,:), tmask_i(:,:), ze2max, iloc )
248         iima2 = iloc(1) ; ijma2 = iloc(2)
249      ELSE
250         ze1min = MINVAL( e1t(:,:), mask = tmask_i(:,:) == 1._wp )   
251         ze2min = MINVAL( e2t(:,:), mask = tmask_i(:,:) == 1._wp )   
252         ze1max = MAXVAL( e1t(:,:), mask = tmask_i(:,:) == 1._wp )   
253         ze2max = MAXVAL( e2t(:,:), mask = tmask_i(:,:) == 1._wp )   
254
255         iloc  = MINLOC( e1t(:,:), mask = tmask_i(:,:) == 1._wp )
256         iimi1 = iloc(1) + nimpp - 1
257         ijmi1 = iloc(2) + njmpp - 1
258         iloc  = MINLOC( e2t(:,:), mask = tmask_i(:,:) == 1._wp )
259         iimi2 = iloc(1) + nimpp - 1
260         ijmi2 = iloc(2) + njmpp - 1
261         iloc  = MAXLOC( e1t(:,:), mask = tmask_i(:,:) == 1._wp )
262         iima1 = iloc(1) + nimpp - 1
263         ijma1 = iloc(2) + njmpp - 1
264         iloc  = MAXLOC( e2t(:,:), mask = tmask_i(:,:) == 1._wp )
265         iima2 = iloc(1) + nimpp - 1
266         ijma2 = iloc(2) + njmpp - 1
267      ENDIF
268      IF(lwp) THEN
269         WRITE(numout,*)
270         WRITE(numout,*) 'dom_ctl : extrema of the masked scale factors'
271         WRITE(numout,*) '~~~~~~~'
272         WRITE(numout,"(14x,'e1t maxi: ',1f10.2,' at i = ',i5,' j= ',i5)") ze1max, iima1, ijma1
273         WRITE(numout,"(14x,'e1t mini: ',1f10.2,' at i = ',i5,' j= ',i5)") ze1min, iimi1, ijmi1
274         WRITE(numout,"(14x,'e2t maxi: ',1f10.2,' at i = ',i5,' j= ',i5)") ze2max, iima2, ijma2
275         WRITE(numout,"(14x,'e2t mini: ',1f10.2,' at i = ',i5,' j= ',i5)") ze2min, iimi2, ijmi2
276      ENDIF
277#endif
278      !
279      ! check that all processes are still there... If some process have an error,
280      ! they will never enter in cfg_write
281      IF( lk_mpp )   CALL mpp_max( 'nemogcm',nstop )
282      IF (nstop /= 0) THEN
283         WRITE(numout,*) ''
284         WRITE(numout,*) '========================================================'
285         WRITE(numout,*) 'E R R O R : ',nstop, ' error have been found'
286         WRITE(numout,*) '========================================================'
287         WRITE(numout,*) ''
288         IF ( lk_mpp ) THEN
289            CALL mppstop()
290         ELSE
291            STOP 123
292         END IF
293      END IF
294      !
295   END SUBROUTINE dom_ctl
296
297
298   SUBROUTINE cfg_write
299      !!----------------------------------------------------------------------
300      !!                  ***  ROUTINE cfg_write  ***
301      !!                   
302      !! ** Purpose :   Create the "domain_cfg" file, a NetCDF file which
303      !!              contains all the ocean domain informations required to
304      !!              define an ocean configuration.
305      !!
306      !! ** Method  :   Write in a file all the arrays required to set up an
307      !!              ocean configuration.
308      !!
309      !! ** output file :   domain_cfg.nc : domain size, characteristics,horizontal mesh,
310      !!                              Coriolis parameter, and vertical scale factors
311      !!                              NB: also contains ORCA family information (if cp_cfg = "ORCA")
312      !!                              and depths (ln_e3_dep=F)
313      !!----------------------------------------------------------------------
314      INTEGER           ::   ji, jj, jk   ! dummy loop indices
315      INTEGER           ::   izco, izps, isco, icav
316      INTEGER           ::   inum     ! temprary units for 'domain_cfg.nc' file
317      CHARACTER(len=21) ::   clnam    ! filename (mesh and mask informations)
318      REAL(wp), DIMENSION(jpi,jpj) ::   z2d   ! workspace
319      !!----------------------------------------------------------------------
320      !
321      IF(lwp) WRITE(numout,*)
322      IF(lwp) WRITE(numout,*) 'cfg_write : create the "domain_cfg.nc" file containing all required configuration information'
323      IF(lwp) WRITE(numout,*) '~~~~~~~~~'
324      !
325      !                       ! ============================= !
326      !                       !  create 'domain_cfg.nc' file  !
327      !                       ! ============================= !
328      !         
329      clnam = 'domain_cfg'  ! filename (configuration information)
330      CALL iom_open( TRIM(clnam), inum, ldwrt = .TRUE.)!, kiolib = jprstlib )
331     
332      !
333      !                             !==  ORCA family specificities  ==!
334      IF( cp_cfg == "ORCA" ) THEN
335         CALL iom_rstput( 0, 0, inum, 'ORCA'      , 1._wp            , ktype = jp_i4 )
336         CALL iom_rstput( 0, 0, inum, 'ORCA_index', REAL( jp_cfg, wp), ktype = jp_i4 )         
337      ENDIF
338      !                             !==  global domain size  ==!
339      !
340      CALL iom_rstput( 0, 0, inum, 'jpiglo', REAL( jpiglo, wp), ktype = jp_i4 )
341      CALL iom_rstput( 0, 0, inum, 'jpjglo', REAL( jpjglo, wp), ktype = jp_i4 )
342      CALL iom_rstput( 0, 0, inum, 'jpkglo', REAL( jpk   , wp), ktype = jp_i4 )
343      !
344      !                             !==  domain characteristics  ==!
345      !
346      !                                   ! lateral boundary of the global
347      !                                   domain
348      CALL iom_rstput( 0, 0, inum, 'jperio', REAL( jperio, wp), ktype = jp_i4 )
349      !
350      !                                   ! type of vertical coordinate
351      IF( ln_zco    ) THEN   ;   izco = 1   ;   ELSE   ;   izco = 0   ;   ENDIF
352      IF( ln_zps    ) THEN   ;   izps = 1   ;   ELSE   ;   izps = 0   ;   ENDIF
353      IF( ln_sco    ) THEN   ;   isco = 1   ;   ELSE   ;   isco = 0   ;   ENDIF
354      CALL iom_rstput( 0, 0, inum, 'ln_zco'   , REAL( izco, wp), ktype = jp_i4 )
355      CALL iom_rstput( 0, 0, inum, 'ln_zps'   , REAL( izps, wp), ktype = jp_i4 )
356      CALL iom_rstput( 0, 0, inum, 'ln_sco'   , REAL( isco, wp), ktype = jp_i4 )
357      !
358      !                                   ! ocean cavities under iceshelves
359      IF( ln_isfcav ) THEN   ;   icav = 1   ;   ELSE   ;   icav = 0   ;   ENDIF
360      CALL iom_rstput( 0, 0, inum, 'ln_isfcav', REAL( icav, wp), ktype = jp_i4 )
361      !
362      !                             !==  horizontal mesh  !
363      !
364      CALL iom_rstput( 0, 0, inum, 'glamt', glamt, ktype = jp_r8 )   ! latitude
365      CALL iom_rstput( 0, 0, inum, 'glamu', glamu, ktype = jp_r8 )
366      CALL iom_rstput( 0, 0, inum, 'glamv', glamv, ktype = jp_r8 )
367      CALL iom_rstput( 0, 0, inum, 'glamf', glamf, ktype = jp_r8 )
368      !                               
369      CALL iom_rstput( 0, 0, inum, 'gphit', gphit, ktype = jp_r8 )   ! longitude
370      CALL iom_rstput( 0, 0, inum, 'gphiu', gphiu, ktype = jp_r8 )
371      CALL iom_rstput( 0, 0, inum, 'gphiv', gphiv, ktype = jp_r8 )
372      CALL iom_rstput( 0, 0, inum, 'gphif', gphif, ktype = jp_r8 )
373      !                               
374      CALL iom_rstput( 0, 0, inum, 'e1t'  , e1t  , ktype = jp_r8 )   ! i-scale factors (e1.)
375      CALL iom_rstput( 0, 0, inum, 'e1u'  , e1u  , ktype = jp_r8 )
376      CALL iom_rstput( 0, 0, inum, 'e1v'  , e1v  , ktype = jp_r8 )
377      CALL iom_rstput( 0, 0, inum, 'e1f'  , e1f  , ktype = jp_r8 )
378      !
379      CALL iom_rstput( 0, 0, inum, 'e2t'  , e2t  , ktype = jp_r8 )   ! j-scale factors (e2.)
380      CALL iom_rstput( 0, 0, inum, 'e2u'  , e2u  , ktype = jp_r8 )
381      CALL iom_rstput( 0, 0, inum, 'e2v'  , e2v  , ktype = jp_r8 )
382      CALL iom_rstput( 0, 0, inum, 'e2f'  , e2f  , ktype = jp_r8 )
383      !
384      CALL iom_rstput( 0, 0, inum, 'ff_f' , ff_f , ktype = jp_r8 )   ! coriolis factor
385      CALL iom_rstput( 0, 0, inum, 'ff_t' , ff_t , ktype = jp_r8 )
386      !
387      !                             !==  vertical mesh  ==!
388      !                                                     
389      CALL iom_rstput( 0, 0, inum, 'e3t_1d'  , e3t_1d  , ktype = jp_r8 )   !  reference 1D-coordinate
390      CALL iom_rstput( 0, 0, inum, 'e3w_1d'  , e3w_1d  , ktype = jp_r8 )
391      !
392      CALL iom_rstput( 0, 0, inum, 'e3t_0'   , e3t_0   , ktype = jp_r8 )   !  vertical scale factors (e
393      CALL iom_rstput( 0, 0, inum, 'e3u_0'   , e3u_0   , ktype = jp_r8 )
394      CALL iom_rstput( 0, 0, inum, 'e3v_0'   , e3v_0   , ktype = jp_r8 )
395      CALL iom_rstput( 0, 0, inum, 'e3f_0'   , e3f_0   , ktype = jp_r8 )
396      CALL iom_rstput( 0, 0, inum, 'e3w_0'   , e3w_0   , ktype = jp_r8 )
397      CALL iom_rstput( 0, 0, inum, 'e3uw_0'  , e3uw_0  , ktype = jp_r8 )
398      CALL iom_rstput( 0, 0, inum, 'e3vw_0'  , e3vw_0  , ktype = jp_r8 )
399      !
400      IF(.NOT.ln_e3_dep ) THEN                                             !  depth (t- & w-points)
401         CALL iom_rstput( 0, 0, inum, 'gdept_1d', gdept_1d, ktype = jp_r8 )   ! required only with 
402         CALL iom_rstput( 0, 0, inum, 'gdepw_1d', gdepw_1d, ktype = jp_r8 )   ! the old e3. definition
403         CALL iom_rstput( 0, 0, inum, 'gdept_0' , gdept_0 , ktype = jp_r8 )
404         CALL iom_rstput( 0, 0, inum, 'gdepw_0' , gdepw_0 , ktype = jp_r8 )
405      ENDIF
406      !                                         
407      !                             !==  ocean top and bottom level  ==!
408      !
409      CALL iom_rstput( 0, 0, inum, 'bottom_level' , REAL( mbkt, wp )*ssmask , ktype = jp_i4 )   ! nb of ocean T-points
410      CALL iom_rstput( 0, 0, inum, 'top_level'    , REAL( mikt, wp )*ssmask , ktype = jp_i4 )   ! nb of ocean T-points (ISF)
411      CALL iom_rstput( 0, 0, inum, 'isf_draft'    , risfdep , ktype = jp_r8 )
412      DO jj = 1,jpj
413         DO ji = 1,jpi
414            z2d (ji,jj) = SUM ( e3t_0(ji,jj, 1:mbkt(ji,jj) ) ) * ssmask(ji,jj) 
415         END DO
416      END DO
417      CALL iom_rstput( 0, 0, inum, 'bathy_metry'   , z2d , ktype = jp_r8 )
418      !
419      !                              !== closed sea ==!
420      IF (ln_domclo) THEN
421         ! mask for the open sea
422         CALL iom_rstput( 0, 0, inum, 'mask_opensea' , msk_opnsea  , ktype = jp_i4 )
423         ! mask for all the under closed sea
424         CALL iom_rstput( 0, 0, inum, 'mask_csundef' , msk_csundef , ktype = jp_i4 )
425         ! mask for global, local net precip, local net precip and evaporation correction
426         CALL iom_rstput( 0, 0, inum, 'mask_csglo'   , msk_csglo   , ktype = jp_i4 )
427         CALL iom_rstput( 0, 0, inum, 'mask_csemp'   , msk_csemp   , ktype = jp_i4 )
428         CALL iom_rstput( 0, 0, inum, 'mask_csrnf'   , msk_csrnf   , ktype = jp_i4 )
429         ! mask for the various river mouth (in case multiple lake in the same outlet)
430         CALL iom_rstput( 0, 0, inum, 'mask_csgrpglo', msk_csgrpglo, ktype = jp_i4 )
431         CALL iom_rstput( 0, 0, inum, 'mask_csgrpemp', msk_csgrpemp, ktype = jp_i4 )
432         CALL iom_rstput( 0, 0, inum, 'mask_csgrprnf', msk_csgrprnf, ktype = jp_i4 )
433      END IF
434      !
435      !                                ! ============================
436      !                                !        close the files
437      !                                ! ============================
438      CALL iom_close( inum )
439      !
440   END SUBROUTINE cfg_write
441
442   !!======================================================================
443END MODULE domain
Note: See TracBrowser for help on using the repository browser.