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.
trdmxl_rst.F90 in branches/UKMO/dev_r7750_GO6_package_oasis_timers/NEMOGCM/NEMO/OPA_SRC/TRD – NEMO

source: branches/UKMO/dev_r7750_GO6_package_oasis_timers/NEMOGCM/NEMO/OPA_SRC/TRD/trdmxl_rst.F90 @ 8895

Last change on this file since 8895 was 8895, checked in by andmirek, 6 years ago

#1978 new timers for NEMO restart write
File size: 9.9 KB
Line 
1MODULE trdmxl_rst
2   !!=================================================================================
3   !!                       ***  MODULE  trdmxl_rst  ***
4   !! Ocean dynamic :  Input/Output files for restart on mixed-layer diagnostics
5   !!=================================================================================
6   !! History :  1.0  ! 2005-05 (C. Deltel)  Original code
7   !!---------------------------------------------------------------------------------
8
9   !!---------------------------------------------------------------------------------
10   !!  trd_mxl_rst_write : write mixed layer trend restart
11   !!  trd_mxl_rst_read  : read  mixed layer trend restart
12   !!---------------------------------------------------------------------------------
13   USE dom_oce         ! ocean space and time domain
14   USE trd_oce         ! trends: ocean variables
15   USE in_out_manager  ! I/O manager
16   USE iom             ! I/O module
17   USE restart         ! only for lrst_oce
18   USE timing
19
20   IMPLICIT NONE
21   PRIVATE
22 
23   PUBLIC   trd_mxl_rst_read    ! routine called by trd_mxl_init
24   PUBLIC   trd_mxl_rst_write   ! routine called by step.F90
25 
26   INTEGER ::   nummxlw         ! logical unit for mxl restart
27
28   !!---------------------------------------------------------------------------------
29   !! NEMO/OPA 3.3 , NEMO Consortium (2010)
30   !! $Id$
31   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
32   !!---------------------------------------------------------------------------------
33CONTAINS
34 
35   SUBROUTINE trd_mxl_rst_write( kt )
36      !!--------------------------------------------------------------------------------
37      !!                  ***  SUBROUTINE trd_mxl_rst_wri  ***
38      !!               
39      !! ** Purpose :   Write mixed-layer diagnostics restart fields.
40      !!--------------------------------------------------------------------------------
41      INTEGER, INTENT( in ) ::   kt   ! ocean time-step index
42      !
43      CHARACTER (len=35) :: charout
44      INTEGER ::   jk                 ! loop indice
45      CHARACTER(LEN=20)   ::   clkt     ! ocean time-step deine as a character
46      CHARACTER(LEN=50)   ::   clname   ! output restart file name
47      CHARACTER(LEN=256)  ::   clpath   ! full path to restart file
48      !!--------------------------------------------------------------------------------
49
50      ! to get better performances with NetCDF format:
51      ! we open and define the ocean restart_mxl file one time step before writing the data (-> at nitrst - 1)
52      ! except if we write ocean restart_mxl files every time step or if an ocean restart_mxl file was writen at nitend - 1
53      IF( kt == nitrst - 1 .OR. nstock == 1 .OR. ( kt == nitend .AND. MOD( nitend - 1, nstock ) == 0 ) ) THEN
54         ! beware of the format used to write kt (default is i8.8, that should be large enough...)
55         IF( nitrst > 999999999 ) THEN   ;   WRITE(clkt, *       ) nitrst
56         ELSE                            ;   WRITE(clkt, '(i8.8)') nitrst
57         ENDIF
58         ! create the file
59         clname = TRIM(cexper)//"_"//TRIM(ADJUSTL(clkt))//"_"//TRIM(cn_trdrst_out)
60         clpath = TRIM(cn_ocerst_outdir)
61         IF( clpath(LEN_TRIM(clpath):) /= '/' ) clpath = TRIM(clpath) // '/'
62         IF(lwp) THEN
63            WRITE(numout,*)
64            SELECT CASE ( jprstlib )
65            CASE ( jprstdimg )   ;   WRITE(numout,*) '             open ocean restart_mxl binary file: '//clname
66            CASE DEFAULT         ;   WRITE(numout,*) '             open ocean restart_mxl NetCDF file: '//clname
67            END SELECT
68            IF( kt == nitrst - 1 ) THEN   ;   WRITE(numout,*) '             kt = nitrst - 1 = ', kt,' date= ', ndastp
69            ELSE                          ;   WRITE(numout,*) '             kt = '             , kt,' date= ', ndastp
70            ENDIF
71         ENDIF
72
73         CALL iom_open( TRIM(clpath)//TRIM(clname), nummxlw, ldwrt = .TRUE., kiolib = jprstlib )
74      ENDIF
75
76      IF( kt == nitrst .AND. lwp ) THEN
77         WRITE(numout,*)
78         WRITE(numout,*) 'trdmxl_rst: output for ML diags. restart, with trd_mxl_rst_write routine kt =', kt
79         WRITE(numout,*) '~~~~~~~~~~'
80         WRITE(numout,*)
81      ENDIF
82
83      IF( ln_trdmxl_instant ) THEN 
84         !-- Temperature
85         IF(nn_timing == 2)  CALL timing_start('rst_put')
86         CALL iom_rstput( kt, nitrst, nummxlw, 'tmlbb'           , tmlbb           )
87         CALL iom_rstput( kt, nitrst, nummxlw, 'tmlbn'           , tmlbn           )
88         CALL iom_rstput( kt, nitrst, nummxlw, 'tmlatfb'         , tmlatfb         )
89
90         !-- Salinity
91         CALL iom_rstput( kt, nitrst, nummxlw, 'smlbb'           , smlbb           )
92         CALL iom_rstput( kt, nitrst, nummxlw, 'smlbn'           , smlbn           )
93         CALL iom_rstput( kt, nitrst, nummxlw, 'smlatfb'         , smlatfb         )
94      ELSE
95         CALL iom_rstput( kt, nitrst, nummxlw, 'hmxlbn'          , hmxlbn          )
96
97         !-- Temperature
98         CALL iom_rstput( kt, nitrst, nummxlw, 'tmlbn'           , tmlbn           )
99         CALL iom_rstput( kt, nitrst, nummxlw, 'tml_sumb'        , tml_sumb        )
100         DO jk = 1, jpltrd
101            IF( jk < 10 ) THEN   ;   WRITE(charout,FMT="('tmltrd_csum_ub_', I1)")   jk
102            ELSE                 ;   WRITE(charout,FMT="('tmltrd_csum_ub_', I2)")   jk
103            ENDIF
104            CALL iom_rstput( kt, nitrst, nummxlw, charout,  tmltrd_csum_ub(:,:,jk) )
105         ENDDO
106         CALL iom_rstput( kt, nitrst, nummxlw, 'tmltrd_atf_sumb' , tmltrd_atf_sumb )
107
108         !-- Salinity
109         CALL iom_rstput( kt, nitrst, nummxlw, 'smlbn'           , smlbn           )
110         CALL iom_rstput( kt, nitrst, nummxlw, 'sml_sumb'        , sml_sumb        )
111         DO jk = 1, jpltrd
112            IF( jk < 10 ) THEN   ;   WRITE(charout,FMT="('smltrd_csum_ub_', I1)")   jk
113            ELSE                 ;   WRITE(charout,FMT="('smltrd_csum_ub_', I2)")   jk
114            ENDIF
115            CALL iom_rstput( kt, nitrst, nummxlw, charout , smltrd_csum_ub(:,:,jk) )
116         ENDDO
117         CALL iom_rstput( kt, nitrst, nummxlw, 'smltrd_atf_sumb' , smltrd_atf_sumb )
118         IF(nn_timing == 2)  CALL timing_stop('rst_put')
119      ENDIF
120      !
121      IF( kt == nitrst ) THEN
122         CALL iom_close( nummxlw )     ! close the restart file (only at last time step)
123         lrst_oce = .FALSE.
124      ENDIF
125      !
126   END SUBROUTINE trd_mxl_rst_write
127
128
129   SUBROUTINE trd_mxl_rst_read
130      !!----------------------------------------------------------------------------
131      !!                   ***  SUBROUTINE trd_mxl_rst_lec  ***
132      !!                   
133      !! ** Purpose :   Read file for mixed-layer diagnostics restart.
134      !!----------------------------------------------------------------------------
135      INTEGER  ::  inum       ! temporary logical unit
136      !
137      CHARACTER (len=35) :: charout
138      INTEGER ::   jk         ! loop indice
139      INTEGER ::   jlibalt = jprstlib
140      LOGICAL ::   llok
141      CHARACTER(LEN=256)  ::   clpath   ! full path to restart file
142      !!-----------------------------------------------------------------------------
143
144      IF(lwp)  THEN
145         WRITE(numout,*)
146         WRITE(numout,*) ' trd_mxl_rst_read : read the NetCDF mixed layer trend restart file'
147         WRITE(numout,*) ' ~~~~~~~~~~~~~~~~'
148      ENDIF
149
150      clpath = TRIM(cn_ocerst_indir)
151      IF( clpath(LEN_TRIM(clpath):) /= '/' ) clpath = TRIM(clpath) // '/'
152
153      IF ( jprstlib == jprstdimg ) THEN
154         ! eventually read netcdf file (monobloc)  for restarting on different number of processors
155         ! if {cn_trdrst_in}.nc exists, then set jlibalt to jpnf90
156         INQUIRE( FILE = TRIM(clpath)//TRIM(cn_trdrst_in)//'.nc', EXIST = llok )
157         IF ( llok ) THEN   ;   jlibalt = jpnf90   
158         ELSE               ;   jlibalt = jprstlib   
159         ENDIF
160      ENDIF
161
162      CALL iom_open( TRIM(clpath)//TRIM(cn_trdrst_in), inum, kiolib = jlibalt ) 
163
164      IF( ln_trdmxl_instant ) THEN 
165         !-- Temperature
166         CALL iom_get( inum, jpdom_autoglo, 'tmlbb'           , tmlbb          )
167         CALL iom_get( inum, jpdom_autoglo, 'tmlbn'           , tmlbn          )
168         CALL iom_get( inum, jpdom_autoglo, 'tmlatfb'         , tmlatfb        )
169         !
170         !-- Salinity
171         CALL iom_get( inum, jpdom_autoglo, 'smlbb'           , smlbb          )
172         CALL iom_get( inum, jpdom_autoglo, 'smlbn'           , smlbn          )
173         CALL iom_get( inum, jpdom_autoglo, 'smlatfb'         , smlatfb        )
174      ELSE
175         CALL iom_get( inum, jpdom_autoglo, 'hmxlbn'          , hmxlbn         ) ! needed for hmxl_sum
176         !
177         !-- Temperature
178         CALL iom_get( inum, jpdom_autoglo, 'tmlbn'           , tmlbn          ) ! needed for tml_sum
179         CALL iom_get( inum, jpdom_autoglo, 'tml_sumb'        , tml_sumb       )
180         DO jk = 1, jpltrd
181            IF( jk < 10 ) THEN   ;   WRITE(charout,FMT="('tmltrd_csum_ub_', I1)")   jk
182            ELSE                 ;   WRITE(charout,FMT="('tmltrd_csum_ub_', I2)")   jk
183            ENDIF
184            CALL iom_get( inum, jpdom_autoglo, charout, tmltrd_csum_ub(:,:,jk) )
185         END DO
186         CALL iom_get( inum, jpdom_autoglo, 'tmltrd_atf_sumb' , tmltrd_atf_sumb)
187         !
188         !-- Salinity
189         CALL iom_get( inum, jpdom_autoglo, 'smlbn'           , smlbn          ) ! needed for sml_sum
190         CALL iom_get( inum, jpdom_autoglo, 'sml_sumb'        , sml_sumb       )
191         DO jk = 1, jpltrd
192            IF( jk < 10 ) THEN   ;   WRITE(charout,FMT="('smltrd_csum_ub_', I1)")   jk
193            ELSE                 ;   WRITE(charout,FMT="('smltrd_csum_ub_', I2)")   jk
194            ENDIF
195            CALL iom_get( inum, jpdom_autoglo, charout, smltrd_csum_ub(:,:,jk) )
196         END DO
197         CALL iom_get( inum, jpdom_autoglo, 'smltrd_atf_sumb' , smltrd_atf_sumb)
198         !
199         CALL iom_close( inum )
200      ENDIF
201      !
202   END SUBROUTINE trd_mxl_rst_read
203 
204  !!=================================================================================
205END MODULE trdmxl_rst
Note: See TracBrowser for help on using the repository browser.