Newer
Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
MODULE trcini_cfc
!!======================================================================
!! *** MODULE trcini_cfc ***
!! TOP : initialisation of the CFC tracers
!!======================================================================
!! History : 2.0 ! 2007-12 (C. Ethe, G. Madec)
!!----------------------------------------------------------------------
!!----------------------------------------------------------------------
!! trc_ini_cfc : CFC model initialisation
!!----------------------------------------------------------------------
USE oce_trc ! Ocean variables
USE par_trc ! TOP parameters
USE trc ! TOP variables
USE trcnam_cfc ! CFC SMS namelist
USE trcsms_cfc ! CFC sms trends
IMPLICIT NONE
PRIVATE
PUBLIC trc_ini_cfc ! called by trcini.F90 module
INTEGER :: inum ! unit number
REAL(wp) :: ylats = -10. ! 10 degrees south
REAL(wp) :: ylatn = 10. ! 10 degrees north
!! * Substitutions
# include "do_loop_substitute.h90"
!!----------------------------------------------------------------------
!! NEMO/TOP 4.0 , NEMO Consortium (2018)
!! $Id: trcini_cfc.F90 15090 2021-07-06 14:25:18Z cetlod $
!! Software governed by the CeCILL license (see ./LICENSE)
!!----------------------------------------------------------------------
CONTAINS
SUBROUTINE trc_ini_cfc( Kmm )
!!----------------------------------------------------------------------
!! *** trc_ini_cfc ***
!!
!! ** Purpose : initialization for cfc model
!!
!! ** Method : - Read the namcfc namelist and check the parameter values
!!----------------------------------------------------------------------
INTEGER, INTENT(in) :: Kmm ! time level indices
INTEGER :: ji, jj, jn, jl, jm, js, io, ierr
INTEGER :: iskip = 6 ! number of 1st descriptor lines
REAL(wp) :: zyy, zyd
CHARACTER(len = 20) :: cltra
!!----------------------------------------------------------------------
!
CALL trc_nam_cfc
!
IF(lwp) WRITE(numout,*)
IF(lwp) WRITE(numout,*) ' trc_ini_cfc: initialisation of CFC chemical model'
IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~~'
!
IF(lwp) WRITE(numout,*) 'Read annual atmospheric concentratioins from formatted file : ' // TRIM(clname)
CALL ctl_opn( inum, clname, 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE. )
REWIND(inum)
! compute the number of year in the file
! file starts in 1931 do jn represent the year in the century
jn = 31
DO
READ(inum,'(1x)',END=100)
jn = jn + 1
END DO
100 jpyear = jn - 1 - iskip
IF ( lwp) WRITE(numout,*) ' ---> ', jpyear ,' years read'
! ! Allocate CFC arrays
ALLOCATE( p_cfc(jpyear,jphem,3), STAT=ierr )
IF( ierr > 0 ) THEN
CALL ctl_stop( 'trc_ini_cfc: unable to allocate p_cfc array' ) ; RETURN
ENDIF
IF( trc_sms_cfc_alloc() /= 0 ) CALL ctl_stop( 'STOP', 'trc_ini_cfc: unable to allocate CFC arrays' )
! Initialization of boundaries conditions
! ---------------------------------------
xphem (:,:) = 0._wp
p_cfc(:,:,:) = 0._wp
! Initialization of qint in case of no restart
!----------------------------------------------
qtr_cfc(:,:,:) = 0._wp
IF( .NOT. ln_rsttr ) THEN
IF(lwp) THEN
WRITE(numout,*)
WRITE(numout,*) 'Initialisation of qint ; No restart : qint equal zero '
ENDIF
qint_cfc(:,:,:) = 0._wp
DO jl = 1, jp_cfc
jn = jp_cfc0 + jl - 1
tr(:,:,:,jn,Kmm) = 0._wp
END DO
ENDIF
REWIND(inum)
DO jm = 1, iskip ! Skip over 1st six descriptor lines
READ(inum,'(1x)')
END DO
! file starts in 1931 do jn represent the year in the century.jhh
! Read file till the end
jn = 31
DO
READ(inum,*, IOSTAT=io) zyy, p_cfc(jn,1:2,1), p_cfc(jn,1:2,2), p_cfc(jn,1:2,3)
IF( io < 0 ) exit
jn = jn + 1
END DO
!p_cfc(32,1:2,1) = 5.e-4 ! modify the values of the first years
!p_cfc(33,1:2,1) = 8.e-4
!p_cfc(34,1:2,1) = 1.e-6
!p_cfc(35,1:2,1) = 2.e-3
!p_cfc(36,1:2,1) = 4.e-3
!p_cfc(37,1:2,1) = 6.e-3
!p_cfc(38,1:2,1) = 8.e-3
!p_cfc(39,1:2,1) = 1.e-2
IF(lwp) THEN ! Control print
WRITE(numout,*)
WRITE(numout,*) ' Year c11NH c11SH c12NH c12SH SF6NH SF6SH'
DO jn = 30, jpyear
WRITE(numout, '( 1I4, 6F10.4)') jn, p_cfc(jn,1:2,1), p_cfc(jn,1:2,2), p_cfc(jn,1:2,3)
END DO
ENDIF
! Interpolation factor of atmospheric partial pressure
! Linear interpolation between 2 hemispheric function of latitud between ylats and ylatn
!---------------------------------------------------------------------------------------
zyd = ylatn - ylats
DO_2D( nn_hls, nn_hls, nn_hls, nn_hls )
IF( gphit(ji,jj) >= ylatn ) THEN ; xphem(ji,jj) = 1.e0
ELSEIF( gphit(ji,jj) <= ylats ) THEN ; xphem(ji,jj) = 0.e0
ELSE ; xphem(ji,jj) = ( gphit(ji,jj) - ylats) / zyd
ENDIF
END_2D
!
IF(lwp) WRITE(numout,*) 'Initialization of CFC tracers done'
IF(lwp) WRITE(numout,*) ' '
!
END SUBROUTINE trc_ini_cfc
!!======================================================================
END MODULE trcini_cfc