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.
trdmld_rst.F90 in trunk/NEMO/OPA_SRC/TRD – NEMO

source: trunk/NEMO/OPA_SRC/TRD/trdmld_rst.F90 @ 1409

Last change on this file since 1409 was 1229, checked in by smasson, 16 years ago

define restart names in the namelist, see ticket:258
  • Property svn:eol-style set to native
  • Property svn:keywords set to Id
File size: 9.3 KB
Line 
1MODULE trdmld_rst
2   !!=================================================================================
3   !!                       ***  MODULE  trdmld_rst  ***
4   !! Ocean dynamic :  Input/Output files for restart on mixed-layer diagnostics
5   !!=================================================================================
6   !! History :  9.0  ! 05-05 (C. Deltel) Original code
7   !!---------------------------------------------------------------------------------
8#if defined key_trdmld
9   !!---------------------------------------------------------------------------------
10   USE dom_oce         ! ocean space and time domain
11   USE trdmod_oce      ! ocean variables for trend diagnostics (i.e. icp/mixed-layer/vorticity)
12   USE in_out_manager  ! I/O manager
13   USE daymod          ! calendar
14   USE iom             ! I/O module
15   USE restart         ! only for lrst_oce
16
17   IMPLICIT NONE
18   PRIVATE
19 
20   PUBLIC   trd_mld_rst_read    ! routine called by trd_mld_init
21   PUBLIC   trd_mld_rst_write   ! routine called by step.F90
22 
23   INTEGER ::   nummldw         ! logical unit for mld restart
24
25   !!---------------------------------------------------------------------------------
26   !! OPA 9.0 , LOCEAN-IPSL (2006)
27   !! $Id$
28   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
29   !!---------------------------------------------------------------------------------
30 
31CONTAINS
32 
33   SUBROUTINE trd_mld_rst_write( kt )
34      !!--------------------------------------------------------------------------------
35      !!                  ***  SUBROUTINE trd_mld_rst_wri  ***
36      !!               
37      !! ** Purpose :   Write mixed-layer diagnostics restart fields.
38      !!--------------------------------------------------------------------------------
39      INTEGER, INTENT( in ) ::   kt   ! ocean time-step index
40      !
41      CHARACTER (len=35) :: charout
42      INTEGER ::   jk                 ! loop indice
43      CHARACTER(LEN=20)   ::   clkt     ! ocean time-step deine as a character
44      CHARACTER(LEN=50)   ::   clname   ! ice output restart file name
45      !!--------------------------------------------------------------------------------
46
47      ! to get better performances with NetCDF format:
48      ! we open and define the ocean restart_mld file one time step before writing the data (-> at nitrst - 1)
49      ! except if we write ocean restart_mld files every time step or if an ocean restart_mld file was writen at nitend - 1
50      IF( kt == nitrst - 1 .OR. nstock == 1 .OR. ( kt == nitend .AND. MOD( nitend - 1, nstock ) == 0 ) ) THEN
51         ! beware of the format used to write kt (default is i8.8, that should be large enough...)
52         IF( nitrst > 999999999 ) THEN   ;   WRITE(clkt, *       ) nitrst
53         ELSE                            ;   WRITE(clkt, '(i8.8)') nitrst
54         ENDIF
55         ! create the file
56         clname = TRIM(cexper)//"_"//TRIM(ADJUSTL(clkt))//"_"//TRIM(cn_trdrst_out)
57         IF(lwp) THEN
58            WRITE(numout,*)
59            SELECT CASE ( jprstlib )
60            CASE ( jprstdimg )   ;   WRITE(numout,*) '             open ocean restart_mld binary file: '//clname
61            CASE DEFAULT         ;   WRITE(numout,*) '             open ocean restart_mld NetCDF file: '//clname
62            END SELECT
63            IF( kt == nitrst - 1 ) THEN   ;   WRITE(numout,*) '             kt = nitrst - 1 = ', kt,' date= ', ndastp
64            ELSE                          ;   WRITE(numout,*) '             kt = '             , kt,' date= ', ndastp
65            ENDIF
66         ENDIF
67
68         CALL iom_open( clname, nummldw, ldwrt = .TRUE., kiolib = jprstlib )
69      ENDIF
70
71      IF( kt == nitrst .AND. lwp ) THEN
72         WRITE(numout,*)
73         WRITE(numout,*) 'trdmld_rst: output for ML diags. restart, with trd_mld_rst_write routine kt =', kt
74         WRITE(numout,*) '~~~~~~~~~~'
75         WRITE(numout,*)
76      ENDIF
77
78      IF( ln_trdmld_instant ) THEN 
79         !-- Temperature
80         CALL iom_rstput( kt, nitrst, nummldw, 'tmlbb'           , tmlbb           )
81         CALL iom_rstput( kt, nitrst, nummldw, 'tmlbn'           , tmlbn           )
82         CALL iom_rstput( kt, nitrst, nummldw, 'tmlatfb'         , tmlatfb         )
83
84         !-- Salinity
85         CALL iom_rstput( kt, nitrst, nummldw, 'smlbb'           , smlbb           )
86         CALL iom_rstput( kt, nitrst, nummldw, 'smlbn'           , smlbn           )
87         CALL iom_rstput( kt, nitrst, nummldw, 'smlatfb'         , smlatfb         )
88      ELSE
89         CALL iom_rstput( kt, nitrst, nummldw, 'rmldbn'          , rmldbn          )
90
91         !-- Temperature
92         CALL iom_rstput( kt, nitrst, nummldw, 'tmlbn'           , tmlbn           )
93         CALL iom_rstput( kt, nitrst, nummldw, 'tml_sumb'        , tml_sumb        )
94         DO jk = 1, jpltrd
95            IF( jk < 10 ) THEN   ;   WRITE(charout,FMT="('tmltrd_csum_ub_', I1)") jk
96            ELSE                 ;   WRITE(charout,FMT="('tmltrd_csum_ub_', I2)") jk
97            ENDIF
98            CALL iom_rstput( kt, nitrst, nummldw, charout,  tmltrd_csum_ub(:,:,jk) )
99         ENDDO
100         CALL iom_rstput( kt, nitrst, nummldw, 'tmltrd_atf_sumb' , tmltrd_atf_sumb )
101
102         !-- Salinity
103         CALL iom_rstput( kt, nitrst, nummldw, 'smlbn'           , smlbn           )
104         CALL iom_rstput( kt, nitrst, nummldw, 'sml_sumb'        , sml_sumb        )
105         DO jk = 1, jpltrd
106            IF( jk < 10 ) THEN   ;   WRITE(charout,FMT="('smltrd_csum_ub_', I1)") jk
107            ELSE                 ;   WRITE(charout,FMT="('smltrd_csum_ub_', I2)") jk
108            ENDIF
109            CALL iom_rstput( kt, nitrst, nummldw, charout , smltrd_csum_ub(:,:,jk) )
110         ENDDO
111         CALL iom_rstput( kt, nitrst, nummldw, 'smltrd_atf_sumb' , smltrd_atf_sumb )
112      ENDIF
113      !
114      IF( kt == nitrst ) THEN
115         CALL iom_close( nummldw )     ! close the restart file (only at last time step)
116         lrst_oce = .FALSE.
117      ENDIF
118      !
119      !   
120   END SUBROUTINE trd_mld_rst_write
121
122
123   SUBROUTINE trd_mld_rst_read
124    !!----------------------------------------------------------------------------
125    !!                   ***  SUBROUTINE trd_mld_rst_lec  ***
126    !!                   
127    !! ** Purpose :   Read file for mixed-layer diagnostics restart.
128    !!----------------------------------------------------------------------------
129    INTEGER  ::  inum       ! temporary logical unit
130    !
131    CHARACTER (len=35) :: charout
132    INTEGER ::   jk         ! loop indice
133    !!-----------------------------------------------------------------------------
134
135    IF(lwp)  THEN
136       WRITE(numout,*)
137       WRITE(numout,*) ' trd_mld_rst_read : read the NetCDF MLD restart file'
138       WRITE(numout,*) ' ~~~~~~~~~~~~~~~~'
139    ENDIF
140
141    CALL iom_open( cn_trdrst_in, inum, kiolib = jprstlib ) 
142
143    IF( ln_trdmld_instant ) THEN 
144       !-- Temperature
145       CALL iom_get( inum, jpdom_autoglo, 'tmlbb'           , tmlbb          )
146       CALL iom_get( inum, jpdom_autoglo, 'tmlbn'           , tmlbn          )
147       CALL iom_get( inum, jpdom_autoglo, 'tmlatfb'         , tmlatfb        )
148
149       !-- Salinity
150       CALL iom_get( inum, jpdom_autoglo, 'smlbb'           , smlbb          )
151       CALL iom_get( inum, jpdom_autoglo, 'smlbn'           , smlbn          )
152       CALL iom_get( inum, jpdom_autoglo, 'smlatfb'         , smlatfb        )
153    ELSE
154       CALL iom_get( inum, jpdom_autoglo, 'rmldbn'          , rmldbn         ) ! needed for rmld_sum
155
156       !-- Temperature
157       CALL iom_get( inum, jpdom_autoglo, 'tmlbn'           , tmlbn          ) ! needed for tml_sum
158       CALL iom_get( inum, jpdom_autoglo, 'tml_sumb'        , tml_sumb       )
159       DO jk = 1, jpltrd
160          IF( jk < 10 )   THEN
161             WRITE(charout,FMT="('tmltrd_csum_ub_', I1)") jk
162          ELSE
163             WRITE(charout,FMT="('tmltrd_csum_ub_', I2)") jk
164          ENDIF
165          CALL iom_get( inum, jpdom_autoglo, charout, tmltrd_csum_ub(:,:,jk) )
166       ENDDO
167       CALL iom_get( inum, jpdom_autoglo, 'tmltrd_atf_sumb' , tmltrd_atf_sumb)
168
169       !-- Salinity
170       CALL iom_get( inum, jpdom_autoglo, 'smlbn'           , smlbn          ) ! needed for sml_sum
171       CALL iom_get( inum, jpdom_autoglo, 'sml_sumb'        , sml_sumb       )
172       DO jk = 1, jpltrd
173          IF( jk < 10 )   THEN
174             WRITE(charout,FMT="('smltrd_csum_ub_', I1)") jk
175          ELSE
176             WRITE(charout,FMT="('smltrd_csum_ub_', I2)") jk
177          ENDIF
178          CALL iom_get( inum, jpdom_autoglo, charout, smltrd_csum_ub(:,:,jk) )
179       ENDDO
180       CALL iom_get( inum, jpdom_autoglo, 'smltrd_atf_sumb' , smltrd_atf_sumb)
181
182       CALL iom_close( inum )
183    ENDIF
184
185  END SUBROUTINE trd_mld_rst_read
186 
187#else
188  !!=================================================================================
189  !!                       ***  MODULE  trdmld_rst  ***
190  !! Ocean dynamic :  Input/Output files for restart on mixed-layer diagnostics
191  !!=================================================================================
192CONTAINS
193  SUBROUTINE trd_mld_rst_write( kt )           !  No ML diags ==> empty routine
194    WRITE(*,*) 'trd_mld_rst_wri: You should not have seen this print! error?', kt
195  END SUBROUTINE trd_mld_rst_write
196  SUBROUTINE trd_mld_rst_read                  !  No ML Diags ==> empty routine
197  END SUBROUTINE trd_mld_rst_read
198#endif
199
200  !!=================================================================================
201END MODULE trdmld_rst
Note: See TracBrowser for help on using the repository browser.