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.
trcrst.F90 in trunk/NEMO/TOP_SRC – NEMO

source: trunk/NEMO/TOP_SRC/trcrst.F90 @ 561

Last change on this file since 561 was 561, checked in by opalod, 17 years ago

nemo_v1_bugfix_073:CE:well used of iom module to read variables in CFC's restart files
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 18.9 KB
Line 
1MODULE trcrst
2   !!======================================================================
3   !!
4   !!                       *** MODULE trcrst ***
5   !!
6   !!   Read the restart files for passive tracers
7   !!
8   !!======================================================================
9   !!  TOP 1.0,  LOCEAN-IPSL (2005)
10   !! $Header$
11   !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt
12   !!----------------------------------------------------------------------
13#if defined key_passivetrc   
14   !!----------------------------------------------------------------------
15   !! * Modules used
16   !! ==============
17   USE oce_trc
18   USE trc
19   USE sms
20   USE trctrp_lec   
21   USE lib_mpp
22   
23   IMPLICIT NONE
24   PRIVATE
25   
26   !! * Accessibility
27   PUBLIC trc_rst
28   PUBLIC trc_wri
29   
30   !! * Module variables
31   CHARACTER (len=48) ::   &
32      trestart = 'initial.trc.nc'   ! restart file name
33
34   !! * Substitutions
35#  include "passivetrc_substitute.h90"
36   
37CONTAINS
38
39#if defined key_fdir
40   !!----------------------------------------------------------------------
41   !!   'key_fdir'                                       direct access file
42   !!----------------------------------------------------------------------
43#include "trcrst_fdir.h90"
44   
45#else
46
47   SUBROUTINE trc_rst 
48      !!===========================================================================================
49      !!
50      !!                       ROUTINE trc_rst
51      !!                     *******************
52      !!
53      !!  PURPOSE :
54      !!  ---------
55      !!     READ files for restart for passive tracer
56      !!
57      !!   METHOD :
58      !!   -------
59      !!      READ the previous fields on the FILE nutrst
60      !!      the first record indicates previous characterics
61      !!      after control with the present run, we READ :
62      !!      - prognostic variables on the second and more record
63      !!
64      !!   History:
65      !!   --------
66      !!  original  : 96-11
67      !!  00-05 (A. Estublier) TVD Limiter Scheme key_trc_tvd
68      !!  00-12 (O. Aumont, E. Kestenare) read restart file for sediments
69      !!  01-05 (O. Aumont, E. Kestenare) read restart file for calcite and silicate sediments
70      !!  05-03 (O. Aumont and A. El Moussaoui) F90           
71      !!------------------------------------------------------------------------
72      !! * Modules used
73      USE iom
74
75
76      !! local declarations
77      !! ==================
78
79      INTEGER :: &
80         ji, jj, jk, jn 
81      INTEGER  ::   &
82         inum, iarak0                 ! temporary logical unit
83      REAL(wp), DIMENSION(1, 1, 3)  ::   zinfo
84
85      CHARACTER (len=32) :: clname1,clname2,clname
86      REAL(wp) :: caralk,bicarb
87
88#if defined key_trc_pisces 
89#   if ! defined key_cfg_1d && ( defined key_orca_r4 || defined key_orca_r2 || defined key_orca_r05 || defined key_orca_r025 )
90      REAL(wp) ::   zareatot, ztrasum
91#   endif
92#endif
93
94      !!---------------------------------------------------------------------
95      !!  OPA.9 03-2005 
96      !!---------------------------------------------------------------------
97      !! 0. initialisations
98      !!------------------
99
100
101      IF( ln_trcadv_cen2 .OR. ln_trcadv_tvd ) THEN
102         iarak0 = 1
103      ELSE
104         iarak0 = 0
105      ENDIF
106
107
108      IF(lwp) WRITE(numout,*) ' '
109      IF(lwp) WRITE(numout,*) ' *** trc_rst beginning of restart for'
110      IF(lwp) WRITE(numout,*) ' passive tracer'
111      IF(lwp) WRITE(numout,*) ' the present run :'
112      IF(lwp) WRITE(numout,*) '   number job is  : ',no
113      IF(lwp) WRITE(numout,*) '   with the time nit000 : ',nit000
114      IF(lwp) THEN
115         IF( iarak0 == 1 ) THEN
116            WRITE(numout,*) '   and before fields for Arakawa sheme '
117         ENDIF
118         WRITE(numout,*) ' '
119      ENDIF
120
121      ! Time domain : restart
122      ! -------------------------
123
124      IF(lwp) WRITE(numout,*)
125      IF(lwp) WRITE(numout,*)
126      IF(lwp) WRITE(numout,*) ' *** passive tracer restart option'
127      SELECT CASE ( nrsttr )
128      CASE ( 0 )
129         IF(lwp) WRITE(numout,*) ' nrsttr = 0 no control of nit000'
130      CASE ( 1 )
131         IF(lwp) WRITE(numout,*) ' nrsttr = 1 we control the date of nit000'
132      CASE ( 2 )
133         IF(lwp) WRITE(numout,*) ' nrsttr = 2 the date adatrj is read in restart file'
134      CASE DEFAULT
135         IF(lwp) WRITE(numout,*) '  ===>>>> nrsttr not equal 0, 1 or 2 : no control of the date'
136         IF(lwp) WRITE(numout,*) ' =======                   ========='
137      END SELECT
138
139
140
141      CALL iom_open ( 'restart.trc', inum )
142
143      CALL iom_get ( inum, jpdom_unknown, 'info', zinfo )
144
145      IF(lwp) WRITE(numout,*)
146      IF(lwp) WRITE(numout,*) ' Info on the restart file read : '
147      IF(lwp) WRITE(numout,*) '   job number          : ', NINT( zinfo(1, 1, 1) )
148      IF(lwp) WRITE(numout,*) '   time-step           : ', NINT( zinfo(1, 1, 2) )
149      IF(lwp) WRITE(numout,*) '   arakawa option      : ', NINT( zinfo(1, 1, 3) )
150      IF(lwp) WRITE(numout,*)
151
152
153      !! control of date
154      !! -------------------
155
156      IF( nittrc000 - NINT( zinfo( 1, 1, 2 ) ) /= 1 .AND. nrsttr /= 0 )  &
157           & CALL ctl_stop( ' ===>>>> : problem with nit000 for the restart', &
158           & ' verify the restart file or rerun with nrstdt = 0 (namelist)' )
159
160      !! Control of the scheme
161      !! ------------------------
162
163      IF( iarak0 /= NINT( zinfo(1, 1, 3 ) ) ) &
164           & CALL ctl_stop( ' ===>>>> : problem with advection scheme', &
165           & ' it must be the same type for both restart and previous run', &
166           & ' centered or euler '  )
167
168
169      !! ... READ prognostic variables and computes diagnostic variable
170      !! ---------------------------------------------------------------
171
172      DO jn = 1, jptra
173         clname = 'TRN'//ctrcnm(jn)
174         CALL iom_get( inum, jpdom_local, clname, trn(:,:,:,jn)   ) 
175      END DO
176
177      DO jn = 1, jptra
178         clname = 'TRB'//ctrcnm(jn)
179         CALL iom_get( inum, jpdom_local, clname, trb(:,:,:,jn)   ) 
180      END DO
181
182#if defined key_trc_lobster1
183      clname  = 'SEDB'//ctrcnm(jpdet)
184      clname1 = 'SEDN'//ctrcnm(jpdet)
185      CALL iom_get( inum, jpdom_local, clname , sedpocb(:,:) ) 
186      CALL iom_get( inum, jpdom_local, clname1, sedpocn(:,:) ) 
187
188#elif defined key_trc_pisces
189      clname = 'Silicalim'
190      CALL iom_get( inum, jpdom_local, clname, xksi(:,:) ) 
191      xksimax = xksi
192
193      clname  = 'SED'//ctrcnm(jppoc)
194      clname1 = 'SED'//ctrcnm(jpcal)
195      clname2 = 'SED'//ctrcnm(jpsil)
196
197      CALL iom_get( inum, jpdom_local, clname , sedpoc(:,:) ) 
198      CALL iom_get( inum, jpdom_local, clname1, sedcal(:,:) ) 
199      CALL iom_get( inum, jpdom_local, clname2, sedsil(:,:) ) 
200
201#elif defined key_cfc
202      clname  = 'qint'
203      clname1 = 'qtr'
204      DO jn = 1, jptra
205         CALL iom_get( inum, jpdom_unknown, clname,  qint(:,:,jn), jn )
206         CALL iom_get( inum, jpdom_unknown, clname1, qtr (:,:,jn), jn )
207      ENDDO
208#endif
209
210#if defined key_trc_pisces 
211
212#if ! defined key_cfg_1d && ( defined key_orca_r4 || defined key_orca_r2 || defined key_orca_r05 || defined key_orca_r025 )
213
214      zareatot = 0.
215      DO jk = 1, jpk
216         DO jj = 1, jpj
217            DO ji = 1, jpi
218               zareatot = zareatot + tmask(ji,jj,jk) * tmask_i(ji,jj)  &
219#if defined key_off_degrad
220                  &                * facvol(ji,jj,jk)    &
221#endif
222
223                  &                * e1t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,jk) 
224            END DO
225         END DO
226      END DO
227
228      IF( lk_mpp ) THEN
229         CALL mpp_sum( zareatot )     ! sum over the global domain 
230      END IF
231
232      ztrasum = 0.
233      DO jk = 1, jpk
234         DO jj = 1, jpj
235            DO ji = 1, jpi
236               ztrasum = ztrasum + trn(ji,jj,jk,jptal) * tmask(ji,jj,jk) * tmask_i(ji,jj)    &
237#if defined key_off_degrad
238                  &              * facvol(ji,jj,jk)   &
239#endif
240
241                  &              * e1t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,jk)
242            END DO
243         END DO
244      END DO
245
246      IF( lk_mpp ) THEN
247         CALL mpp_sum( ztrasum )     ! sum over the global domain 
248      END IF
249
250      WRITE(0,*) 'TALK moyen ', ztrasum/zareatot*1E6
251      ztrasum = ztrasum/zareatot*1E6
252      trn(:,:,:,jptal) = trn(:,:,:,jptal)*2391./ztrasum
253
254      ztrasum = 0.
255      DO jk = 1, jpk
256         DO jj = 1, jpj
257            DO ji = 1, jpi
258               ztrasum = ztrasum + trn(ji,jj,jk,jppo4) * tmask(ji,jj,jk) * tmask_i(ji,jj)  &
259#if defined key_off_degrad
260                  &              * facvol(ji,jj,jk)   &
261#endif
262
263                  &              * e1t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,jk)
264            END DO
265         END DO
266      END DO
267
268      IF( lk_mpp ) THEN
269         CALL mpp_sum( ztrasum )     ! sum over the global domain 
270      END IF
271
272
273      WRITE(0,*) 'PO4 moyen ', ztrasum/zareatot*1E6/122.
274      ztrasum = ztrasum/zareatot*1E6/122.
275      trn(:,:,:,jppo4) = trn(:,:,:,jppo4)*2.165/ztrasum
276
277      ztrasum = 0.
278      DO jk = 1, jpk
279         DO jj = 1, jpj
280            DO ji = 1, jpi
281               ztrasum = ztrasum + trn(ji,jj,jk,jpno3) * tmask(ji,jj,jk) * tmask_i(ji,jj)  &
282#if defined key_off_degrad
283                  &              * facvol(ji,jj,jk)   &
284#endif
285
286                  &              * e1t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,jk)
287            END DO
288         END DO
289      END DO
290
291      IF( lk_mpp ) THEN
292         CALL mpp_sum( ztrasum )     ! sum over the global domain 
293      END IF
294
295
296      WRITE(0,*) 'NO3 moyen ', ztrasum/zareatot*1E6/7.6
297      ztrasum = ztrasum/zareatot*1E6/7.6
298      trn(:,:,:,jpno3) = trn(:,:,:,jpno3)*30.9/ztrasum
299
300      ztrasum = 0.
301      DO jk = 1, jpk
302         DO jj = 1, jpj
303            DO ji = 1, jpi
304               ztrasum = ztrasum + trn(ji,jj,jk,jpsil) * tmask(ji,jj,jk) * tmask_i(ji,jj)   &
305#if defined key_off_degrad
306                  &              * facvol(ji,jj,jk)   &
307#endif
308
309                  &              * e1t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,jk)
310            END DO
311         END DO
312      END DO
313
314      IF( lk_mpp ) THEN
315         CALL mpp_sum( ztrasum )     ! sum over the global domain 
316      END IF
317
318      WRITE(0,*) 'SiO3 moyen ', ztrasum/zareatot*1E6
319      ztrasum = ztrasum/zareatot*1E6
320      trn(:,:,:,jpsil) = MIN( 400E-6,trn(:,:,:,jpsil)*91.51/ztrasum) 
321
322#endif
323      !!  Initialization of chemical variables of the carbon cycle
324      !!  --------------------------------------------------------
325      DO jk = 1,jpk
326         DO jj = 1,jpj
327            DO ji = 1,jpi
328               caralk = trn(ji,jj,jk,jptal)-       &
329                  &        borat(ji,jj,jk)/(1.+1.E-8/(rtrn+akb3(ji,jj,jk)))
330               co3(ji,jj,jk)=(caralk-trn(ji,jj,jk,jpdic))*tmask(ji,jj,jk)   &
331                  &        +(1.-tmask(ji,jj,jk))*.5e-3
332               bicarb = (2.*trn(ji,jj,jk,jpdic)-caralk)
333               hi(ji,jj,jk) = (ak23(ji,jj,jk)*bicarb/co3(ji,jj,jk))     &
334                  &  *tmask(ji,jj,jk)+(1.-tmask(ji,jj,jk))*1.e-9
335               h2co3(ji,jj) = 1.e-5
336            ENDDO
337         ENDDO
338      ENDDO
339#endif
340      trb(:,:,:,:) = trn(:,:,:,:)
341
342      CALL iom_close( inum )
343
344
345   END SUBROUTINE trc_rst
346
347   SUBROUTINE trc_wri(kt)
348      !! ==================================================================================
349      !!
350      !!                       ROUTINE trc_wri
351      !!                     ******************
352      !!
353      !!  PURPOSE :
354      !!  ---------
355      !!     WRITE restart fields in nutwrs
356      !!   METHOD :
357      !!   -------
358      !!
359      !!   nutwrs FILE:
360      !!   each nstock time step , SAVE fields which are necessary for
361      !!   passive tracer restart
362      !!
363      !!
364      !!   INPUT :
365      !!   -----
366      !!      argument
367      !!              kt              : time step
368      !!      COMMON
369      !!            /cottrc/          : passive tracers fields (before,now
370      !!                                  ,after)
371      !!
372      !!   OUTPUT :
373      !!   ------
374      !!      FILE
375      !!           nutwrs          : standard restart fields OUTPUT
376      !!
377      !!   WORKSPACE :
378      !!   ---------
379      !!      ji,jj,jk,jl,ino0,it0,iarak0
380      !!
381      !!   History:
382      !!   --------
383      !!      original : 96-12
384      !!      addition : 99-12 (M.-A. Foujols) NetCDF FORMAT with ioipsl
385      !!      additions : 00-05 (A. Estublier)
386      !!                  TVD Limiter Scheme : key_trc_tvd
387      !!      additions : 01-01 (M.A Foujols, E. Kestenare) bug fix: restclo
388      !!      additions : 01-01 (O. Aumont, E. Kestenare)
389      !!                  write restart file for sediments
390      !!      additions : 01-05 (O. Aumont, E. Kestenare)
391      !!                  write restart file for calcite and silicate sediments
392      !!   05-03 (O. Aumont and A. El Moussaoui) F90
393      !!========================================================================================!
394      !! * Modules used
395      USE ioipsl
396
397      !! * Arguments
398      !! -----------
399      INTEGER, INTENT( in ) :: kt
400
401      !! * local declarations
402      !! ====================
403
404      LOGICAL :: clbon         !!!
405      CHARACTER (len=50) :: clname,clname1,clname2,cln
406
407      INTEGER :: jn,   &
408         ino0,it0,iarak0,     &
409         ic,jc,ji,jj,jk,      &
410         itime
411
412      REAL(wp) :: zdate0, zinfo(3),zdiag_var,    &
413         zdiag_varmin, zdiag_varmax, zdiag_tot, zder
414
415
416      !! 1. OUTPUT of restart fields (nutwrs)
417      !! ---------------------------
418
419      IF( kt == nittrc000 ) THEN
420         IF(lwp) WRITE(numout,*)
421         IF(lwp) WRITE(numout,*) 'trc_wri : write passive tracers restart.output NetCDF file'
422         IF(lwp) WRITE(numout,*) '~~~~~~~'
423       
424
425         areatot = 0.
426         DO jk = 1, jpk
427            DO jj = 1, jpj
428               DO ji = 1, jpi
429                  areatot = areatot + tmask(ji,jj,jk)*tmask_i(ji,jj)*e1t(ji,jj)*e2t(ji,jj)*fse3t(ji,jj,jk) 
430               END DO
431            END DO
432         END DO
433
434         IF( lk_mpp ) THEN
435             CALL mpp_sum(areatot)     ! sum over the global domain 
436         END IF
437
438         trai = 0.
439         DO jn = 1, jptra
440            DO jk = 1,jpk
441               DO jj = 1, jpj
442                  DO ji = 1, jpi
443                     trai=trai+tmask(ji,jj,jk)*trn(ji,jj,jk,jn)*     &
444                        &    tmask_i(ji,jj)* e1t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,jk) 
445                  END DO
446               END DO
447            END DO
448         END DO
449
450         IF( lk_mpp ) THEN
451             CALL mpp_sum(trai)         ! sum over the global domain 
452         END IF
453
454         IF (lwp) WRITE(numout,*) 'Integral of all tracers over the full domain at NIT000 =',trai
455
456      ENDIF
457
458
459      IF( MOD(kt,nstock) == 0 .OR. kt == nitend ) THEN
460
461         !! 0. initialisations
462         !! ------------------
463
464         IF(lwp) WRITE(numout,*) ' '
465         IF(lwp) WRITE(numout,*) 'trc_wri : write the passive tracer restart file in NetCDF format ',   &
466            'at it= ',kt,' date= ',ndastp
467         IF(lwp) WRITE(numout,*) '~~~~~~~~~'
468
469
470         ino0 =no
471         it0  =kt
472         IF( ln_trcadv_cen2 .OR. ln_trcadv_tvd ) THEN
473            iarak0 = 1
474         ELSE
475            iarak0=0
476         ENDIF
477
478         zinfo(1)=FLOAT(ino0)
479         zinfo(2)=FLOAT(it0)
480         zinfo(3)=FLOAT(iarak0)
481
482         !! 1. WRITE in nutwrs
483         !! ------------------
484         !!... first information
485
486         INQUIRE (FILE=trestart,EXIST=clbon)
487         IF(clbon) THEN
488            OPEN(UNIT=nutwrs,FILE=trestart,STATUS='old')
489            CLOSE(nutwrs,STATUS='delete')
490         ENDIF
491
492         ic=1
493         DO jc=1,16
494            IF(cexper(jc:jc) /= ' ') ic = jc
495         END DO
496         WRITE(cln,'("_",i4.4,i2.2,i2.2,"_restart.trc")') nyear, nmonth, nday
497         clname=cexper(1:ic)//cln
498         ic=1
499         DO jc=1,48
500            IF(clname(jc:jc) /= ' ') ic=jc
501         END DO
502         trestart=clname(1:ic)//".nc"
503         itime=0
504         CALL ymds2ju(nyear,nmonth,nday,0.0,zdate0)
505         CALL restini('NONE',jpi,jpj,glamt,gphit,jpk,gdept_0,clname           &
506            &        ,itime,zdate0,rdt*nstock,nutwrs,domain_id=nidom)
507
508         CALL restput(nutwrs,'info',1,1,3,0,zinfo)
509
510         ! prognostic variables
511         ! --------------------
512
513         IF (lwp) WRITE(numout,*) '----TRACER STAT----'
514         zdiag_tot=0.
515         DO jn=1,jptra
516            clname='TRN'//ctrcnm(jn)
517            CALL restput(nutwrs,clname,jpi,jpj,jpk,0,trn(:,:,:,jn))
518
519            zdiag_var=0.
520            zdiag_varmin=0.
521            zdiag_varmax=0.
522
523            DO ji=1, jpi
524               DO jj=1, jpj
525                  DO jk=1,jpk
526                    zdiag_var=zdiag_var+tmask(ji,jj,jk)*trn(ji,jj,jk,jn)*     &
527                               tmask_i(ji,jj)* e1t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,jk) 
528
529                  END DO
530               END DO
531            END DO
532
533            zdiag_varmin=MINVAL(trn(:,:,:,jn), mask= ((tmask*SPREAD(tmask_i,DIM=3,NCOPIES=jpk).NE.0.)))
534            zdiag_varmax=MAXVAL(trn(:,:,:,jn), mask= ((tmask*SPREAD(tmask_i,DIM=3,NCOPIES=jpk).NE.0.)))
535
536            IF( lk_mpp ) THEN
537               CALL mpp_min(zdiag_varmin)      ! min over the global domain 
538               CALL mpp_max(zdiag_varmax)      ! max over the global domain 
539               CALL mpp_sum(zdiag_var)         ! sum over the global domain 
540            END IF
541
542            zdiag_tot=zdiag_tot+zdiag_var
543            zdiag_var=zdiag_var/areatot
544
545            IF (lwp) WRITE(numout,*) 'MEAN NO ',jn,ctrcnm(jn),' =',zdiag_var,'MIN= '  &
546               ,zdiag_varmin,'MAX= ',zdiag_varmax
547
548         END DO
549
550         zdiag_tot=zdiag_tot
551         zder=((zdiag_tot-trai)/trai)*100._wp
552         IF (lwp) WRITE(numout,*) 'Integral of all tracers over the full domain  =',zdiag_tot 
553         IF (lwp) WRITE(numout,*) 'Drift of the sum of all tracers =',zder, '%' 
554
555         DO jn=1,jptra
556            clname='TRB'//ctrcnm(jn)
557            CALL restput(nutwrs,clname,jpi,jpj,jpk,0,trb(:,:,:,jn))
558         END DO
559
560#if defined key_trc_lobster1
561         clname='SEDB'//ctrcnm(jpdet)
562         clname1='SEDN'//ctrcnm(jpdet)
563         CALL restput(nutwrs,clname,jpi,jpj,1,0,sedpocb(:,:))
564         CALL restput(nutwrs,clname1,jpi,jpj,1,0,sedpocn(:,:))
565#elif defined key_trc_pisces
566         clname='SED'//ctrcnm(jppoc)
567         clname1='SED'//ctrcnm(jpcal)
568         clname2='SED'//ctrcnm(jpsil)
569         CALL restput(nutwrs,clname1,jpi,jpj,1,0,sedcal(:,:))
570         CALL restput(nutwrs,clname2,jpi,jpj,1,0,sedsil(:,:))
571         CALL restput(nutwrs,clname,jpi,jpj,1,0,sedpoc(:,:))
572
573         clname='Silicalim'
574         CALL restput(nutwrs,clname,jpi,jpj,1,0,xksi(:,:))
575#elif defined key_cfc
576         clname='qint'
577         CALL restput(nutwrs,clname,jpi,jpj,jptra,0,qint(:,:,:))
578         clname1='qtr'
579         CALL restput(nutwrs,clname1,jpi,jpj,jptra,0,qtr(:,:,:))
580#endif
581
582
583         CALL restclo(nutwrs)
584
585      ENDIF
586
587   END SUBROUTINE trc_wri
588
589#endif
590
591#else
592   !!======================================================================
593   !!  Empty module : No passive tracer
594   !!======================================================================
595CONTAINS
596
597   SUBROUTINE trc_rst
598      !! no passive tracers
599   END SUBROUTINE trc_rst
600
601   SUBROUTINE trc_wri(kt)
602      !! no passive tracers
603      INTEGER, INTENT ( in ) :: kt
604      WRITE(*,*) 'trc_wri: You should not have seen this print! error?', kt
605   END SUBROUTINE trc_wri
606   
607#endif
608   
609END MODULE trcrst
Note: See TracBrowser for help on using the repository browser.