- Timestamp:
- 2017-06-21T16:27:21+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_GO6M_dev/NEMOGCM/NEMO/TOP_SRC/CFC/trcini_cfc.F90
r6486 r8201 5 5 !!====================================================================== 6 6 !! History : 2.0 ! 2007-12 (C. Ethe, G. Madec) 7 !! ! 2017-04 (A. Yool) Add SF6 7 8 !!---------------------------------------------------------------------- 8 9 #if defined key_cfc … … 22 23 PUBLIC trc_ini_cfc ! called by trcini.F90 module 23 24 24 CHARACTER (len=34) :: clname = 'cfc1112 .atm' ! ???25 CHARACTER (len=34) :: clname = 'cfc1112sf6.atm' ! ??? 25 26 26 27 INTEGER :: inum ! unit number … … 44 45 !!---------------------------------------------------------------------- 45 46 INTEGER :: ji, jj, jn, jl, jm, js, io, ierr 46 INTEGER :: iskip = 6! number of 1st descriptor lines47 INTEGER :: iskip = 7 ! number of 1st descriptor lines 47 48 REAL(wp) :: zyy, zyd 48 49 !!---------------------------------------------------------------------- … … 53 54 54 55 55 IF(lwp) WRITE(numout,*) 'read of formatted file cfc1112 atm'56 IF(lwp) WRITE(numout,*) 'read of formatted file cfc1112sf6.atm' 56 57 57 58 CALL ctl_opn( inum, clname, 'OLD', 'FORMATTED', 'SEQUENTIAL', -1, numout, .FALSE. ) … … 69 70 ! ! Allocate CFC arrays 70 71 71 ALLOCATE( p_cfc(jpyear,jphem, 2), STAT=ierr )72 ALLOCATE( p_cfc(jpyear,jphem,3), STAT=ierr ) 72 73 IF( ierr > 0 ) THEN 73 74 CALL ctl_stop( 'trc_ini_cfc: unable to allocate p_cfc array' ) ; RETURN … … 90 91 ENDIF 91 92 qint_cfc(:,:,:) = 0._wp 92 DO jl = 1, jp_cfc 93 jn = jp_cfc0 + jl - 1 94 trn(:,:,:,jn) = 0._wp 95 END DO 93 trn(:,:,:,jp_cfc0:jp_cfc1) = 0._wp 96 94 ENDIF 97 95 … … 105 103 jn = 31 106 104 DO 107 READ(inum,*, IOSTAT=io) zyy, p_cfc(jn,1,1), p_cfc(jn,1,2), p_cfc(jn,2,1), p_cfc(jn,2,2) 105 READ(inum,*, IOSTAT=io) zyy, p_cfc(jn,1,1), p_cfc(jn,1,2), & 106 & p_cfc(jn,1,3), p_cfc(jn,2,1), & 107 & p_cfc(jn,2,2), p_cfc(jn,2,3) 108 108 IF( io < 0 ) exit 109 109 jn = jn + 1 110 110 END DO 111 111 112 p_cfc(32,1:2,1) = 5.e-4 ! modify the values of the first years 113 p_cfc(33,1:2,1) = 8.e-4 114 p_cfc(34,1:2,1) = 1.e-6 115 p_cfc(35,1:2,1) = 2.e-3 116 p_cfc(36,1:2,1) = 4.e-3 117 p_cfc(37,1:2,1) = 6.e-3 118 p_cfc(38,1:2,1) = 8.e-3 119 p_cfc(39,1:2,1) = 1.e-2 112 ! AXY (25/04/17): do not adjust 113 ! p_cfc(32,1:2,1) = 5.e-4 ! modify the values of the first years 114 ! p_cfc(33,1:2,1) = 8.e-4 115 ! p_cfc(34,1:2,1) = 1.e-6 116 ! p_cfc(35,1:2,1) = 2.e-3 117 ! p_cfc(36,1:2,1) = 4.e-3 118 ! p_cfc(37,1:2,1) = 6.e-3 119 ! p_cfc(38,1:2,1) = 8.e-3 120 ! p_cfc(39,1:2,1) = 1.e-2 120 121 121 122 IF(lwp) THEN ! Control print 122 123 WRITE(numout,*) 123 WRITE(numout,*) ' Year p11HN p11HS p12HN p12HS '124 WRITE(numout,*) ' Year p11HN p11HS p12HN p12HS pSF6N pSF6S ' 124 125 DO jn = 30, jpyear 125 WRITE(numout, '( 1I4, 4F9.2)') jn, p_cfc(jn,1,1), p_cfc(jn,2,1), p_cfc(jn,1,2), p_cfc(jn,2,2) 126 WRITE(numout, '( 1I4, 6F9.2)') jn, p_cfc(jn,1,1), p_cfc(jn,2,1), & 127 & p_cfc(jn,1,2), p_cfc(jn,2,2), & 128 & p_cfc(jn,1,3), p_cfc(jn,2,3) 126 129 END DO 127 130 ENDIF 128 129 131 130 132 ! Interpolation factor of atmospheric partial pressure
Note: See TracChangeset
for help on using the changeset viewer.